Molecular phylogeny and divergence times of Malagasy tenrecs: influence of data partitioning and taxon sampling on dating analyses

Phylogenomics reveals an almost perfect polytomy among the almost ungulates (Paenungulata)

Phylogenetic studies have resolved most relationships among Eutherian Orders. However, the branching order of elephants (Proboscidea), hyraxes (Hyracoidea), and sea cows (Sirenia) (i.e., the Paenungulata) has remained uncertain since at least 1758, when Linnaeus grouped elephants and manatees into a single Order (Bruta) to the exclusion of hyraxes. Subsequent morphological, molecular, and large-scale phylogenomic datasets have reached conflicting conclusions on the branching order within Paenungulates. We use a phylogenomic dataset of alignments from 13,388 protein-coding genes across 261 Eutherian mammals to infer phylogenetic relationships within Paenungulates. We find that gene trees almost equally support the three alternative resolutions of Paenungulate relationships and that despite strong support for a Proboscidea+Hyracoidea split in the multispecies coalescent (MSC) tree, there is significant evidence for gene tree uncertainty, incomplete lineage sorting, and introgression among Proboscidea, Hyracoidea, and Sirenia. Indeed, only 8-10% of genes have statistically significant phylogenetic signal to reject the hypothesis of a Paenungulate polytomy. These data indicate little support for any resolution for the branching order Proboscidea, Hyracoidea, and Sirenia within Paenungulata and suggest that Paenungulata may be as close to a real, or at least unresolvable, polytomy as possible.

Testing for adaptive radiation: A new approach applied to Madagascar frogs

Adaptive radiation is a key topic at the intersection of ecology and evolutionary biology. Yet the definition and identification of adaptive radiation both remain contentious. Here, we introduce a new approach for identifying adaptive radiations that combines key aspects of two widely used definitions. Our approach compares evolutionary rates in morphology, performance, and diversification between the candidate radiation and other clades. We then apply this approach to a putative adaptive radiation of frogs from Madagascar (Mantellidae). We present new data on morphology and performance from mantellid frogs, then compare rates of diversification and multivariate evolution of size, shape, and performance between mantellids and other frogs. We find that mantellids potentially pass our test for accelerated rates of evolution for shape, but not for size, performance, or diversification. Our results demonstrate that clades can have accelerated phenotypic evolution without rapid diversification (dubbed "adaptive non-radiation"). We also highlight general issues in testing for adaptive radiation, including taxon sampling and the problem of including another adaptive radiation among the comparison clades. Finally, we suggest that similar tests should be conducted on other putative adaptive radiations on Madagascar, comparing their evolutionary rates to those of related clades outside Madagascar. Based on our results, we speculate that older Madagascar clades may show evolutionary patterns more similar to those on a continent than an island.

Jaw shape and mechanical advantage are indicative of diet in Mesozoic mammals

Jaw morphology is closely linked to both diet and biomechanical performance, and jaws are one of the most common Mesozoic mammal fossil elements. Knowledge of the dietary and functional diversity of early mammals informs on the ecological structure of palaeocommunities throughout the longest era of mammalian evolution: the Mesozoic. Here, we analyse how jaw shape and mechanical advantage of the masseter (MAM) and temporalis (MAT) muscles relate to diet in 70 extant and 45 extinct mammals spanning the Late Triassic-Late Cretaceous. In extant mammals, jaw shape discriminates well between dietary groups: insectivores have long jaws, carnivores intermediate to short jaws, and herbivores have short jaws. Insectivores have low MAM and MAT, carnivores have low MAM and high MAT, and herbivores have high MAM and MAT. These traits are also informative of diet among Mesozoic mammals (based on previous independent determinations of diet) and set the basis for future ecomorphological studies.

Speciation and gene flow in two sympatric small mammals from Madagascar, Microgale fotsifotsy and M. soricoides (Mammalia: Tenrecidae)

Madagascar's shrew tenrecs (Mammalia: Tenrecidae; Microgale, Nesogale) represent an excellent system for studying speciation. Most species are endemic to the island's eastern humid forests, a region renowned for high levels of biodiversity and a high rate of in situ diversification. We set out to understand the speciation dynamics in a clade of recently described taxa: Microgale fotsifotsy and M. soricoides, which have nearly identical distributions in the moist evergreen forest, and M. nasoloi, which occurs in the western dry deciduous forest. A phylogenetic analysis using mitochondrial DNA data recovered two distinct clades of M. fotsifotsy: a south clade that is sister to, and broadly sympatric with, M. soricoides, and a north clade that is sister to the dry-forest and distantly allopatric species M. nasoloi. To better understand this result, we analysed cranioskeletal measurements and performed demographic analyses using nuclear sequence data from ultraconserved elements. Nuclear data did not support a sister relationship between M. soricoides and the south clade of M. fotsifotsy but did demonstrate introgression between these clades, which probably explains the discordance between nuclear and mitochondrial phylogenies. Demographic analyses also revealed the absence of gene flow between the north and south clades of M. fotsifotsy. Morphometric data revealed several major differences between M. soricoides and M. fotsifotsy, as well as more subtle differences between the two clades of M. fotsifotsy. In light of these results, we treat the south clade of M. fotsifotsy as a new candidate species. Our findings demonstrate the utility of integrating multiple data types to understand complex speciation histories, and contribute to a growing body of evidence that species diversity on Madagascar is underestimated.

Pathogenic Leptospira and their animal reservoirs: testing host specificity through experimental infection

Leptospirosis is caused by pathogenic Leptospira transmitted through contact with contaminated environments. Most mammalian species are infectable by Leptospira but only few act as efficient reservoir being capable of establishing long term kidney colonization and shedding Leptospira in urine. In Madagascar, a large diversity of pathogenic Leptospira display a tight specificity towards their endemic volant or terrestrial mammalian hosts. The basis of this specificity is unknown: it may indicate some genetically determined compatibility between host cells and bacteria or only reflect ecological constraints preventing contacts between specific hosts. In this study, Rattus norvegicus was experimentally infected with either Leptospira interrogans, Leptospira borgpetersenii or Leptospira mayottensis isolated from rats, bats or tenrecs, respectively. Leptospira borgpetersenii and L. mayottensis do not support renal colonization as featured by no shedding of live bacteria in urine and low level and sporadic detection of Leptospira DNA in kidneys. In contrast 2 out of the 7 R. norvegicus challenged with L. interrogans developed renal colonization and intense Leptospira shedding in urine throughout the 3 months of experimental infection. These data suggest that host-Leptospira specificity in this biodiversity hotspot is driven at least in part by genetic determinants likely resulting from long-term co-diversification processes.

Sympatric lineages in the Mantidactylus ambreensis complex of Malagasy frogs originated allopatrically rather than by in-situ speciation

Madagascar's biota is characterized by a high degree of microendemism at different taxonomic levels, but how colonization and in-situ speciation contribute to the assembly of local species communities has rarely been studied on this island. Here we analyze the phylogenetic relationships of riparian frogs of the Mantidactylus ambreensis species complex, which is distributed in the north of Madagascar and was originally described from Montagne d'Ambre, an isolated mountain of volcanic origin, currently protected within Montagne d'Ambre National Park (MANP). Data from mitochondrial DNA, and phylogenomic data from FrogCap, a sequence capture method, independently confirm that this species complex is monophyletic within the subgenus Ochthomantis, and identify two main clades within it. These two clades are separated by 5.6-6.8% pairwise distance in the mitochondrial 16S rRNA gene and co-occur in MANP, with one distributed at high elevations (940-1375 m a.s.l.) and the other at lower elevations (535-1010 m a.s.l.), but show almost no haplotype sharing in the nuclear RAG1 gene. This occurrence in syntopy without admixture confirms them as independent evolutionary lineages that merit recognition as separate species, and we here refer to them as high-elevation (HE) and low-elevation (LE) lineage; they will warrant taxonomic assessment to confidently assign the name ambreensis to one or the other. Populations of the M. ambreensis complex from elsewhere in northern Madagascar all belong to the LE lineage, although they do occur over a larger elevational range than in Montagne d'Ambre (285-1040 m a.s.l.). Within LE there are several phylogroups (LE1-LE4) of moderately deep divergence (1.5-2.8% in 16S), but phylogroup LE4 that occurs in MANP has a deeply nested phylogenetic position, as recovered separately by mitochondrial and sequence capture datasets. This suggests that HE and LE did not diverge by a local fission of lower and upper populations, but instead arose through a more complex biogeographic scenario. The branching pattern of phylogroups LE1-LE4 shows a clear south-to-north phylogeographic pattern. We derive from these results a testable hypothesis of vicariant speciation that restricted the HE lineage to MANP and the LE candidate species to a climatic refugium further south, with subsequent northwards range expansion and secondary colonization of MANP by LE. These results provide an example for complex assembly of local microendemic amphibian faunas on Madagascar.

Molecular phylogeny and biogeography of the land snail subfamily Leptaxinae (Gastropoda: Hygromiidae)

The subfamily Leptaxinae is included within the highly diverse land snail family Hygromiidae. In the absence of clear diagnostic morphological differences, the subfamily status is currently based solely on molecular information and includes three disjunctly distributed tribes, Leptaxini, Cryptosaccini and Metafruticicolini. However, the phylogenetic relationships among these tribes are not fully resolved and the clustering of some of the genera to the tribes is not statistically supported. To resolve the relationships within Leptaxinae and their position within Hygromiidae, we reconstructed their phylogeny using a multi-locus approach with two mitochondrial genes and eight nuclear markers. The phylogeny was further calibrated and an analysis of ancestral area estimation was carried out to infer the biogeographic history of the group. We elevated Metafruticicolini to subfamily level (Metafruticicolinae) and we restricted Leptaxinae to Cryptosaccini and Leptaxini. The Lusitanian genus Portugala was moved to Leptaxini, previously containing only the Macaronesian genus Leptaxis. Within Cryptosaccini, a new genus strictly confined to the Sierra de la Cabrera (Spain) is described, Fractanella gen. nov. According to our results, Leptaxinae originated in the Early Miocene in the Iberian Peninsula, from which the Macaronesian Islands were colonized. Due to the old split recovered for the divergence between Macaronesian and Iberian lineages, we hypothesize that this colonization may have occurred via the once emerged seamounts located between the archipelagos and the European and African continents, although this could also have occurred through the oldest now emerged islands of Macaronesia. In the Iberian Peninsula, the climatic shift that began during the Middle Miocene, changing progressively from subtropical climate towards the present-day Mediterranean climate, was identified as an important factor shaping the subfamily's diversification, along with Pleistocene climatic fluctuations.

Review of the status and conservation of tenrecs (Mammalia: Afrotheria: Tenrecidae)

The mammal family Tenrecidae (Afrotheria: Afrosoricida) is endemic to Madagascar. Here we present the conservation priorities for the 31 species of tenrec that were assessed or reassessed in 2015–2016 for the IUCN Red List of Threatened Species. Six species (19.4%) were found to be threatened (4 Vulnerable, 2 Endangered) and one species was categorized as Data Deficient. The primary threat to tenrecs is habitat loss, mostly as a result of slash-and-burn agriculture, but some species are also threatened by hunting and incidental capture in fishing traps. In the longer term, climate change is expected to alter tenrec habitats and ranges. However, the lack of data for most tenrecs on population size, ecology and distribution, together with frequent changes in taxonomy (with many cryptic species being discovered based on genetic analyses) and the poorly understood impact of bushmeat hunting on spiny species (Tenrecinae), hinders conservation planning. Priority conservation actions are presented for Madagascar's tenrecs for the first time since 1990 and focus on conserving forest habitat (especially through improved management of protected areas) and filling essential knowledge gaps. Tenrec research, monitoring and conservation should be integrated into broader sustainable development objectives and programmes targeting higher profile species, such as lemurs, if we are to see an improvement in the conservation status of tenrecs in the near future.

An Evaluation of Different Partitioning Strategies for Bayesian Estimation of Species Divergence Times

The explosive growth of molecular sequence data has made it possible to estimate species divergence times under relaxed-clock models using genome-scale data sets with many gene loci. In order to improve both model realism and to best extract information about relative divergence times in the sequence data, it is important to account for the heterogeneity in the evolutionary process across genes or genomic regions. Partitioning is a commonly used approach to achieve those goals. We group sites that have similar evolutionary characteristics into the same partition and those with different characteristics into different partitions, and then use different models or different values of model parameters for different partitions to account for the among-partition heterogeneity. However, how to partition data in practical phylogenetic analysis, and in particular in relaxed-clock dating analysis, is more art than science. Here, we use computer simulation and real data analysis to study the impact of the partition scheme on divergence time estimation. The partition schemes had relatively minor effects on the accuracy of posterior time estimates when the prior assumptions were correct and the clock was not seriously violated, but showed large differences when the clock was seriously violated, when the fossil calibrations were in conflict or incorrect, or when the rate prior was mis-specified. Concatenation produced the widest posterior intervals with the least precision. Use of many partitions increased the precision, as predicted by the infinite-sites theory, but the posterior intervals might fail to include the true ages because of the conflicting fossil calibrations or mis-specified rate priors. We analyzed a data set of 78 plastid genes from 15 plant species with serious clock violation and showed that time estimates differed significantly among partition schemes, irrespective of the rate drift model used. Multiple and precise fossil calibrations reduced the differences among partition schemes and were important to improving the precision of divergence time estimates. While the use of many partitions is an important approach to reducing the uncertainty in posterior time estimates, we do not recommend its general use for the present, given the limitations of current models of rate drift for partitioned data and the challenges of interpreting the fossil evidence to construct accurate and informative calibrations.

Extreme physiological plasticity in a hibernating basoendothermic mammal, Tenrec ecaudatus

Physiological plasticity allows organisms to respond to diverse conditions. However, can being too plastic actually be detrimental? Malagasy common tenrecs, Tenrec ecaudatus, have many plesiomorphic traits and may represent a basal placental mammal. We established a laboratory population of T. ecaudatus and found extreme plasticity in thermoregulation and metabolism, a novel hibernation form, variable annual timing, and remarkable growth and reproductive biology. For instance, tenrec body temperature (Tb) may approximate ambient temperature to as low as 12°C even when tenrecs are fully active. Conversely, tenrecs can hibernate with Tbs of 28°C. During the active season, oxygen consumption may vary 25-fold with little or no changes in Tb. During the Austral winter, tenrecs are consistently torpid but the depth of torpor may be variable. A righting assay revealed that Tb contributes to but does not dictate activity status. Homeostatic processes are not always linked e.g. a hibernating tenrec experienced a ∼34% decrease in heart rate while maintaining constant body temperature and oxygen consumption rates. Tenrec growth rates vary but young may grow ∼40-fold in the 5 weeks until weaning and may possess indeterminate growth as adults. Despite all of this profound plasticity, tenrecs are surprisingly intolerant to extremes in ambient temperature (<8 or >34°C). We contend that while plasticity may confer numerous energetic advantages in consistently moderate environments, environmental extremes may have limited the success and distribution of plastic basal mammals.

Muroid rodent phylogenetics: 900-species tree reveals increasing diversification rates

We combined new sequence data for more than 300 muroid rodent species with our previously published sequences for up to five nuclear and one mitochondrial genes to generate the most widely and densely sampled hypothesis of evolutionary relationships across Muroidea. An exhaustive screening procedure for publically available sequences was implemented to avoid the propagation of taxonomic errors that are common to supermatrix studies. The combined data set of carefully screened sequences derived from all available sequences on GenBank with our new data resulted in a robust maximum likelihood phylogeny for 900 of the approximately 1,620 muroids. Several regions that were equivocally resolved in previous studies are now more decisively resolved, and we estimated a chronogram using 28 fossil calibrations for the most integrated age and topological estimates to date. The results were used to update muroid classification and highlight questions needing additional data. We also compared the results of multigene supermatrix studies like this one with the principal published supertrees and concluded that the latter are unreliable for any comparative study in muroids. In addition, we explored diversification patterns as an explanation for why muroid rodents represent one of the most species-rich groups of mammals by detecting evidence for increasing net diversification rates through time across the muroid tree. We suggest the observation of increasing rates may be due to a combination of parallel increases in rate across clades and high average extinction rates. Five increased diversification-rate-shifts were inferred, suggesting that multiple, but perhaps not independent, events have led to the remarkable species diversity in the superfamily. Our results provide a phylogenetic framework for comparative studies that is not highly dependent upon the signal from any one gene.

A temporal banding approach for consistent taxonomic ranking above the species level

Comparable taxonomic ranks within clades can facilitate more consistent classifications and objective comparisons among taxa. Here we use a temporal approach to identify taxonomic ranks. This is an extension of the temporal banding approach including a Temporal Error Score that finds an objective cut-off for each taxonomic rank using information for the current classification. We illustrate this method using a data set of the lichenized fungal family Parmeliaceae. To assess its performance, we simulated the effect of taxon sampling and compared our method with the other temporal banding method. For our sampled phylogeny, 11 of the 12 included families remained intact and 55 genera were confirmed, whereas 32 genera were lumped and 15 genera were split. Taxon sampling impacted the method at the genus level, whereas yielded only insignificant changes at the family level. The other available temporal approach also gives a similar cutoff point to our method. Our approach to identify taxonomic ranks enables taxonomists to revise and propose classifications on an objective basis, changing ranks of clades only when inconsistent with most taxa in a phylogenetic tree. An R script to find the time point with the minimal temporal error is provided.

Transcriptomic and macroevolutionary evidence for phenotypic uncoupling between frog life history phases

Anuran amphibians undergo major morphological transitions during development, but the contribution of their markedly different life-history phases to macroevolution has rarely been analysed. Here we generate testable predictions for coupling versus uncoupling of phenotypic evolution of tadpole and adult life-history phases, and for the underlying expression of genes related to morphological feature formation. We test these predictions by combining evidence from gene expression in two distantly related frogs, Xenopus laevis and Mantidactylus betsileanus, with patterns of morphological evolution in the entire radiation of Madagascan mantellid frogs. Genes linked to morphological structure formation are expressed in a highly phase-specific pattern, suggesting uncoupling of phenotypic evolution across life-history phases. This gene expression pattern agrees with uncoupled rates of trait evolution among life-history phases in the mantellids, which we show to have undergone an adaptive radiation. Our results validate a prevalence of uncoupling in the evolution of tadpole and adult phenotypes of frogs.

Fossils matter: improved estimates of divergence times in Pinus reveal older diversification

BACKGROUND

The taxonomy of pines (genus Pinus) is widely accepted and a robust gene tree based on entire plastome sequences exists. However, there is a large discrepancy in estimated divergence times of major pine clades among existing studies, mainly due to differences in fossil placement and dating methods used. We currently lack a dated molecular phylogeny that makes use of the rich pine fossil record, and this study is the first to estimate the divergence dates of pines based on a large number of fossils (21) evenly distributed across all major clades, in combination with applying both node and tip dating methods.

RESULTS

We present a range of molecular phylogenetic trees of Pinus generated within a Bayesian framework. We find the origin of crown Pinus is likely up to 30 Myr older (Early Cretaceous) than inferred in most previous studies (Late Cretaceous) and propose generally older divergence times for major clades within Pinus than previously thought. Our age estimates vary significantly between the different dating approaches, but the results generally agree on older divergence times. We present a revised list of 21 fossils that are suitable to use in dating or comparative analyses of pines.

CONCLUSIONS

Reliable estimates of divergence times in pines are essential if we are to link diversification processes and functional adaptation of this genus to geological events or to changing climates. In addition to older divergence times in Pinus, our results also indicate that node age estimates in pines depend on dating approaches and the specific fossil sets used, reflecting inherent differences in various dating approaches. The sets of dated phylogenetic trees of pines presented here provide a way to account for uncertainties in age estimations when applying comparative phylogenetic methods.

Mammal survival at the Cretaceous-Palaeogene boundary: metabolic homeostasis in prolonged tropical hibernation in tenrecs

Free-ranging common tenrecs, Tenrec ecaudatus, from sub-tropical Madagascar, displayed long-term (nine months) hibernation which lacked any evidence of periodic interbout arousals (IBAs). IBAs are the dominant feature of the mammalian hibernation phenotype and are thought to periodically restore long-term ischaemia damage and/or metabolic imbalances (depletions and accumulations). However, the lack of IBAs in tenrecs suggests no such pathology at hibernation Tbs > 22°C. The long period of tropical hibernation that we report might explain how the ancestral placental mammal survived the global devastation that drove the dinosaurs and many other vertebrates to extinction at the Cretaceous-Palaeogene boundary following a meteorite impact. The genetics and biochemistry of IBAs are of immense interest to biomedical researchers and space exploration scientists, in the latter case, those envisioning a hibernating state in astronauts for deep space travel. Unravelling the physiological thresholds and temperature dependence of IBAs will provide new impetus to these research quests.

Morphological diversity in tenrecs (Afrosoricida, Tenrecidae): comparing tenrec skull diversity to their closest relatives

It is important to quantify patterns of morphological diversity to enhance our understanding of variation in ecological and evolutionary traits. Here, we present a quantitative analysis of morphological diversity in a family of small mammals, the tenrecs (Afrosoricida, Tenrecidae). Tenrecs are often cited as an example of an exceptionally morphologically diverse group. However, this assumption has not been tested quantitatively. We use geometric morphometric analyses of skull shape to test whether tenrecs are more morphologically diverse than their closest relatives, the golden moles (Afrosoricida, Chrysochloridae). Tenrecs occupy a wider range of ecological niches than golden moles so we predict that they will be more morphologically diverse. Contrary to our expectations, we find that tenrec skulls are only more morphologically diverse than golden moles when measured in lateral view. Furthermore, similarities among the species-rich Microgale tenrec genus appear to mask higher morphological diversity in the rest of the family. These results reveal new insights into the morphological diversity of tenrecs and highlight the importance of using quantitative methods to test qualitative assumptions about patterns of morphological diversity.

Multilocus phylogenetic and geospatial analyses illuminate diversification patterns and the biogeographic history of Malagasy endemic plated lizards (Gerrhosauridae: Zonosaurinae)

Although numerous studies have attempted to find single unifying mechanisms for generating Madagascar's unique flora and fauna, little consensus has been reached regarding the relative importance of climatic, geologic and ecological processes as catalysts of diversification of the region's unique biota. Rather, recent work has shown that both biological and physical drivers of diversification are best analysed in a case-by-case setting with attention focused on the ecological and life-history requirements of the specific phylogenetic lineage under investigation. Here, we utilize a comprehensive analytical approach to examine evolutionary drivers and elucidate the biogeographic history of Malagasy plated lizards (Zonosaurinae). Data from three genes are combined with fossil information to construct time-calibrated species trees for zonosaurines and their African relatives, which are used to test alternative diversification hypotheses. Methods are utilized for explicitly incorporating phylogenetic uncertainty into downstream analyses. Species distribution models are created for 14 of 19 currently recognized species, which are then used to estimate spatial patterns of species richness and endemicity. Spatially explicit analyses are employed to correlate patterns of diversity with both topographic heterogeneity and climatic stability through geologic time. We then use inferred geographic ranges to estimate the biogeographic history of zonosaurines within each of Madagascar's major biomes. Results suggest constant Neogene and Quaternary speciation with divergence from the African most recent common ancestor ~30 million years ago when oceanic currents and African rivers facilitated dispersal. Spatial patterns of diversity appear concentrated along coastal regions of northern and southern Madagascar. We find no relationship between either topographic heterogeneity or climatic stability and patterns of diversity. Ancestral state reconstructions suggest that western dry forests were important centres of origin with recent invasion into spiny and rain forest. These data highlight the power of combining multilocus phylogenetic and spatially explicit analyses for testing alternative diversification hypotheses within Madagascar's unique biota and more generally, particularly as applied to phylogenetically and biologically constrained systems.

The effects of partitioning on phylogenetic inference

Partitioning is a commonly used method in phylogenetics that aims to accommodate variation in substitution patterns among sites. Despite its popularity, there have been few systematic studies of its effects on phylogenetic inference, and there have been no studies that compare the effects of different approaches to partitioning across many empirical data sets. In this study, we applied four commonly used approaches to partitioning to each of 34 empirical data sets, and then compared the resulting tree topologies, branch-lengths, and bootstrap support estimated using each approach. We find that the choice of partitioning scheme often affects tree topology, particularly when partitioning is omitted. Most notably, we find occasional instances where the use of a suboptimal partitioning scheme produces highly supported but incorrect nodes in the tree. Branch-lengths and bootstrap support are also affected by the choice of partitioning scheme, sometimes dramatically so. We discuss the reasons for these effects and make some suggestions for best practice.

Modeling olfactory bulb evolution through primate phylogeny

Adaptive characterizations of primates have usually included a reduction in olfactory sensitivity. However, this inference of derivation and directionality assumes an ancestral state of olfaction, usually by comparison to a group of extant non-primate mammals. Thus, the accuracy of the inference depends on the assumed ancestral state. Here I present a phylogenetic model of continuous trait evolution that reconstructs olfactory bulb volumes for ancestral nodes of primates and mammal outgroups. Parent-daughter comparisons suggest that, relative to the ancestral euarchontan, the crown-primate node is plesiomorphic and that derived reduction in olfactory sensitivity is an attribute of the haplorhine lineage. The model also suggests a derived increase in olfactory sensitivity at the strepsirrhine node. This oppositional diversification of the strepsirrhine and haplorhine lineages from an intermediate and non-derived ancestor is inconsistent with a characterization of graded reduction through primate evolution.

Retroviral envelope syncytin capture in an ancestrally diverged mammalian clade for placentation in the primitive Afrotherian tenrecs

Syncytins are fusogenic envelope (env) genes of retroviral origin that have been captured for a function in placentation. Syncytins have been identified in Euarchontoglires (primates, rodents, Leporidae) and Laurasiatheria (Carnivora, ruminants) placental mammals. Here, we searched for similar genes in species that retained characteristic features of primitive mammals, namely the Malagasy and mainland African Tenrecidae. They belong to the superorder Afrotheria, an early lineage that diverged from Euarchotonglires and Laurasiatheria 100 Mya, during the Cretaceous terrestrial revolution. An in silico search for env genes with full coding capacity within a Tenrecidae genome identified several candidates, with one displaying placenta-specific expression as revealed by RT-PCR analysis of a large panel of Setifer setosus tissues. Cloning of this endogenous retroviral env gene demonstrated fusogenicity in an ex vivo cell-cell fusion assay on a panel of mammalian cells. Refined analysis of placental architecture and ultrastructure combined with in situ hybridization demonstrated specific expression of the gene in multinucleate cellular masses and layers at the materno-fetal interface, consistent with a role in syncytium formation. This gene, which we named "syncytin-Ten1," is conserved among Tenrecidae, with evidence of purifying selection and conservation of fusogenic activity. To our knowledge, it is the first syncytin identified to date within the ancestrally diverged Afrotheria superorder.

A comprehensive multilocus phylogeny of the Neotropical cotingas (Cotingidae, Aves) with a comparative evolutionary analysis of breeding system and plumage dimorphism and a revised phylogenetic classification

The Neotropical cotingas (Cotingidae: Aves) are a group of passerine birds that are characterized by extreme diversity in morphology, ecology, breeding system, and behavior. Here, we present a comprehensive phylogeny of the Neotropical cotingas based on six nuclear and mitochondrial loci (∼7500 bp) for a sample of 61 cotinga species in all 25 genera, and 22 species of suboscine outgroups. Our taxon sample more than doubles the number of cotinga species studied in previous analyses, and allows us to test the monophyly of the cotingas as well as their intrageneric relationships with high resolution. We analyze our genetic data using a Bayesian species tree method, and concatenated Bayesian and maximum likelihood methods, and present a highly supported phylogenetic hypothesis. We confirm the monophyly of the cotingas, and present the first phylogenetic evidence for the relationships of Phibalura flavirostris as the sister group to Ampelion and Doliornis, and the paraphyly of Lipaugus with respect to Tijuca. In addition, we resolve the diverse radiations within the Cotinga, Lipaugus, Pipreola, and Procnias genera. We find no support for Darwin's (1871) hypothesis that the increase in sexual selection associated with polygynous breeding systems drives the evolution of color dimorphism in the cotingas, at least when analyzed at a broad categorical scale. Finally, we present a new comprehensive phylogenetic classification of all cotinga species.

Diversification of an emerging pathogen in a biodiversity hotspot: Leptospira in endemic small mammals of Madagascar

Biodiversity hotspots and associated endemism are ideal systems for the study of parasite diversity within host communities. Here, we investigated the ecological and evolutionary forces acting on the diversification of an emerging bacterial pathogen, Leptospira spp., in communities of endemic Malagasy small mammals. We determined the infection rate with pathogenic Leptospira in 20 species of sympatric rodents (subfamily Nesomyinae) and tenrecids (family Tenrecidae) at two eastern humid forest localities. A multilocus genotyping analysis allowed the characterization of bacterial diversity within small mammals and gave insights into their genetic relationships with Leptospira infecting endemic Malagasy bats (family Miniopteridae and Vespertilionidae). We report for the first time the presence of pathogenic Leptospira in Malagasy endemic small mammals, with an overall prevalence of 13%. In addition, these hosts harbour species of Leptospira (L. kirschneri, L. borgpetersenii and L. borgpetersenii group B) which are different from those reported in introduced rats (L. interrogans) on Madagascar. The diversification of Leptospira on Madagascar can be traced millions of years into evolutionary history, resulting in the divergence of endemic lineages and strong host specificity. These observations are discussed in relation to the relative roles of endemic vs. introduced mammal species in the evolution and epidemiology of Leptospira on Madagascar, specifically how biodiversity and biogeographical processes can shape community ecology of an emerging pathogen and lead to its diversification within native animal communities.

Selecting optimal partitioning schemes for phylogenomic datasets

BACKGROUND

Partitioning involves estimating independent models of molecular evolution for different subsets of sites in a sequence alignment, and has been shown to improve phylogenetic inference. Current methods for estimating best-fit partitioning schemes, however, are only computationally feasible with datasets of fewer than 100 loci. This is a problem because datasets with thousands of loci are increasingly common in phylogenetics.

METHODS

We develop two novel methods for estimating best-fit partitioning schemes on large phylogenomic datasets: strict and relaxed hierarchical clustering. These methods use information from the underlying data to cluster together similar subsets of sites in an alignment, and build on clustering approaches that have been proposed elsewhere.

RESULTS

We compare the performance of our methods to each other, and to existing methods for selecting partitioning schemes. We demonstrate that while strict hierarchical clustering has the best computational efficiency on very large datasets, relaxed hierarchical clustering provides scalable efficiency and returns dramatically better partitioning schemes as assessed by common criteria such as AICc and BIC scores.

CONCLUSIONS

These two methods provide the best current approaches to inferring partitioning schemes for very large datasets. We provide free open-source implementations of the methods in the PartitionFinder software. We hope that the use of these methods will help to improve the inferences made from large phylogenomic datasets.

Selecting optimal partitioning schemes for phylogenomic datasets

BACKGROUND

Partitioning involves estimating independent models of molecular evolution for different subsets of sites in a sequence alignment, and has been shown to improve phylogenetic inference. Current methods for estimating best-fit partitioning schemes, however, are only computationally feasible with datasets of fewer than 100 loci. This is a problem because datasets with thousands of loci are increasingly common in phylogenetics.

METHODS

We develop two novel methods for estimating best-fit partitioning schemes on large phylogenomic datasets: strict and relaxed hierarchical clustering. These methods use information from the underlying data to cluster together similar subsets of sites in an alignment, and build on clustering approaches that have been proposed elsewhere.

RESULTS

We compare the performance of our methods to each other, and to existing methods for selecting partitioning schemes. We demonstrate that while strict hierarchical clustering has the best computational efficiency on very large datasets, relaxed hierarchical clustering provides scalable efficiency and returns dramatically better partitioning schemes as assessed by common criteria such as AICc and BIC scores.

CONCLUSIONS

These two methods provide the best current approaches to inferring partitioning schemes for very large datasets. We provide free open-source implementations of the methods in the PartitionFinder software. We hope that the use of these methods will help to improve the inferences made from large phylogenomic datasets.

Selecting optimal partitioning schemes for phylogenomic datasets

BACKGROUND

Partitioning involves estimating independent models of molecular evolution for different subsets of sites in a sequence alignment, and has been shown to improve phylogenetic inference. Current methods for estimating best-fit partitioning schemes, however, are only computationally feasible with datasets of fewer than 100 loci. This is a problem because datasets with thousands of loci are increasingly common in phylogenetics.

METHODS

We develop two novel methods for estimating best-fit partitioning schemes on large phylogenomic datasets: strict and relaxed hierarchical clustering. These methods use information from the underlying data to cluster together similar subsets of sites in an alignment, and build on clustering approaches that have been proposed elsewhere.

RESULTS

We compare the performance of our methods to each other, and to existing methods for selecting partitioning schemes. We demonstrate that while strict hierarchical clustering has the best computational efficiency on very large datasets, relaxed hierarchical clustering provides scalable efficiency and returns dramatically better partitioning schemes as assessed by common criteria such as AICc and BIC scores.

CONCLUSIONS

These two methods provide the best current approaches to inferring partitioning schemes for very large datasets. We provide free open-source implementations of the methods in the PartitionFinder software. We hope that the use of these methods will help to improve the inferences made from large phylogenomic datasets.

Insights into the evolution of a cryptic radiation of bats: dispersal and ecological radiation of Malagasy Miniopterus (Chiroptera: Miniopteridae)

The past decade has seen a proliferation of new species of Miniopterus bats (family Miniopteridae) recognized from Madagascar and the neighboring Comoros archipelago. The interspecific relationships of these taxa, their colonization history, and the evolution of this presumed adaptive radiation have not been sufficiently explored. Using the mitochondrial cytochrome-b gene, we present a phylogeny of the Malagasy members of this widespread Old World genus, based on 218 sequences, of which 82 are new and 136 derived from previous studies. Phylogenetic analyses recovered 18 clades, which divide into five primary lineages: (1) M. griveaudi; (2) M. mahafaliensis, M. sororculus and X3; (3) M. majori, M. gleni and M. griffithsi; (4) M. brachytragos; M. aelleniA, and M. aelleniB; and (5) M. manavi and M. petersoni recovered as sister species, which were in turn linked to a group comprising M. egeri and five genetically distinct populations referred to herein as P3, P4, P5, P6 and P7. Beast analysis indicated that the initial divergence within the Malagasy Miniopterus radiation took place 4.5 Myr; most species diverged between 4 and 2.5 Myr, and a secondary period was between 1.25 and 1 Myr. DNA K2P-distances between recognized taxa ranged from 12.9% to 2.5% and intraspecific variation was less than 1.8%. Of the 18 identified clades, Latin binomials are only associated with 11, which indicates much greater differentiation than currently recognized for Malagasy Miniopterus. These data are placed in a context of the dispersal history of this genus on the island and patterns of ecological diversity.

Imperfect isolation: factors and filters shaping Madagascar's extant vertebrate fauna

Analyses of phylogenetic topology and estimates of divergence timing have facilitated a reconstruction of Madagascar's colonization events by vertebrate animals, but that information alone does not reveal the major factors shaping the island's biogeographic history. Here, we examine profiles of Malagasy vertebrate clades through time within the context of the island's paleogeographical evolution to determine how particular events influenced the arrival of the island's extant groups. First we compare vertebrate profiles on Madagascar before and after selected events; then we compare tetrapod profiles on Madagascar to contemporary tetrapod compositions globally. We show that changes from the Mesozoic to the Cenozoic in the proportions of Madagascar's tetrapod clades (particularly its increase in the representation of birds and mammals) are tied to changes in their relative proportions elsewhere on the globe. Differences in the representation of vertebrate classes from the Mesozoic to the Cenozoic reflect the effects of extinction (i.e., the non-random susceptibility of the different vertebrate clades to purported catastrophic global events 65 million years ago), and new evolutionary opportunities for a subset of vertebrates with the relatively high potential for transoceanic dispersal potential. In comparison, changes in vertebrate class representation during the Cenozoic are minor. Despite the fact that the island's isolation has resulted in high vertebrate endemism and a unique and taxonomically imbalanced extant vertebrate assemblage (both hailed as testimony to its long isolation), that isolation was never complete. Indeed, Madagascar's extant tetrapod fauna owes more to colonization during the Cenozoic than to earlier arrivals. Madagascar's unusual vertebrate assemblage needs to be understood with reference to the basal character of clades originating prior to the K-T extinction, as well as to the differential transoceanic dispersal advantage of other, more recently arriving clades. Thus, the composition of Madagascar's endemic vertebrate assemblage itself provides evidence of the island's paleogeographic history.

A molecular phylogeny of Plesiorycteropus reassigns the extinct mammalian order 'Bibymalagasia'

Madagascar is well known for its diverse fauna and flora, being home to many species not found anywhere else in the world. However, its biodiversity in the recent past included a range of extinct enigmatic fauna, such as elephant birds, giant lemurs and dwarfed hippopotami. The 'Malagasy aardvark' (Plesiorycteropus) has remained one of Madagascar's least well-understood extinct species since its discovery in the 19(th) century. Initially considered a close relative of the aardvark (Orycteropus) within the order Tubulidentata, more recent morphological analyses challenged this placement on the grounds that the identifiably derived traits supporting this allocation were adaptations to digging rather than shared ancestry. Because the skeletal evidence showed many morphological traits diagnostic of different eutherian mammal orders, they could not be used to resolve its closest relatives. As a result, the genus was tentatively assigned its own taxonomic order 'Bibymalagasia', yet how this order relates to other eutherian mammal orders remains unclear despite numerous morphological investigations. This research presents the first known molecular sequence data for Plesiorycteropus, obtained from the bone protein collagen (I), which places the 'Malagasy aardvark' as more closely related to tenrecs than aardvarks. More specifically, Plesiorycteropus was recovered within the order Tenrecoidea (golden moles and tenrecs) within Afrotheria, suggesting that the taxonomic order 'Bibymalagasia' is obsolete. This research highlights the potential for collagen sequencing in investigating the phylogeny of extinct species as a viable alternative to ancient DNA (aDNA) sequencing, particularly in cases where aDNA cannot be recovered.

The chromosomes of Afrotheria and their bearing on mammalian genome evolution

Afrotheria is the clade of placental mammals that, together with Xenarthra, Euarchontoglires and Laurasiatheria, represents 1 of the 4 main recognized supraordinal eutherian clades. It reunites 6 orders of African origin: Proboscidea, Sirenia, Hyracoidea, Macroscelidea, Afrosoricida and Tubulidentata. The apparently unlikely relationship among such disparate morphological taxa and their possible basal position at the base of the eutherian phylogenetic tree led to a great deal of attention and research on the group. The use of biomolecular data was pivotal in Afrotheria studies, as they were the basis for the recognition of this clade. Although morphological evidence is still scarce, a plethora of molecular data firmly attests to the phylogenetic relationship among these mammals of African origin. Modern cytogenetic techniques also gave a significant contribution to the study of Afrotheria, revealing chromosome signatures for the group as a whole, as well as for some of its internal relationships. The associations of human chromosomes HSA1/19 and 5/21 were found to be chromosome signatures for the group and provided further support for Afrotheria. Additional chromosome synapomorphies were also identified linking elephants and manatees in Tethytheria (the associations HSA2/3, 3/13, 8/22, 18/19 and the lack of HSA4/8) and elephant shrews with the aardvark (HSA2/8, 3/20 and 10/17). Herein, we review the current knowledge on Afrotheria chromosomes and genome evolution. The already available data on the group suggests that further work on this apparently bizarre assemblage of mammals will provide important data to a better understanding on mammalian genome evolution.

RUNX2 tandem repeats and the evolution of facial length in placental mammals

Background

When simple sequence repeats are integrated into functional genes, they can potentially act as evolutionary ‘tuning knobs’, supplying abundant genetic variation with minimal risk of pleiotropic deleterious effects. The genetic basis of variation in facial shape and length represents a possible example of this phenomenon. Runt-related transcription factor 2 (RUNX2), which is involved in osteoblast differentiation, contains a functionally-important tandem repeat of glutamine and alanine amino acids. The ratio of glutamines to alanines (the QA ratio) in this protein seemingly influences the regulation of bone development. Notably, in domestic breeds of dog, and in carnivorans in general, the ratio of glutamines to alanines is strongly correlated with facial length.

Results

In this study we examine whether this correlation holds true across placental mammals, particularly those mammals for which facial length is highly variable and related to adaptive behavior and lifestyle (e.g., primates, afrotherians, xenarthrans). We obtained relative facial length measurements and RUNX2 sequences for 41 mammalian species representing 12 orders. Using both a phylogenetic generalized least squares model and a recently-developed Bayesian comparative method, we tested for a correlation between genetic and morphometric data while controlling for phylogeny, evolutionary rates, and divergence times. Non-carnivoran taxa generally had substantially lower glutamine-alanine ratios than carnivorans (primates and xenarthrans with means of 1.34 and 1.25, respectively, compared to a mean of 3.1 for carnivorans), and we found no correlation between RUNX2 sequence and face length across placental mammals.

Conclusions

Results of our diverse comparative phylogenetic analyses indicate that QA ratio does not consistently correlate with face length across the 41 mammalian taxa considered. Thus, although RUNX2 might function as a ‘tuning knob’ modifying face length in carnivorans, this relationship is not conserved across mammals in general.

Impact of the partitioning scheme on divergence times inferred from Mammalian genomic data sets

Data partitioning has long been regarded as an important parameter for phylogenetic inference. The division of heterogeneous multigene data sets into partitions with similar substitution patterns is known to increase the performance of probabilistic phylogenetic methods. However, the effect of the partitioning scheme on divergence time estimates has generally been ignored. To investigate the impact of data partitioning on the estimation of divergence times, we have constructed two genomic data sets. The first one with 15 nuclear genes comprising 50,928 bp were selected from the OrthoMam database; the second set was composed of complete mitochondrial genomes. We studied two partitioning schemes: concatenated supermatrices and partitioned gene analysis. We have also measured the impact of taxonomic sampling on the estimates. After drawing divergence time inferences using the uncorrelated relaxed clock in BEAST, we have compared the age estimates between the partitioning schemes. Our results show that, in general, both schemes resulted in similar chronological estimates, however the concatenated data sets were more efficient than the partitioned ones in attaining suitable effective sample sizes.

The historical biogeography of Mammalia

Palaeobiogeographic reconstructions are underpinned by phylogenies, divergence times and ancestral area reconstructions, which together yield ancestral area chronograms that provide a basis for proposing and testing hypotheses of dispersal and vicariance. Methods for area coding include multi-state coding with a single character, binary coding with multiple characters and string coding. Ancestral reconstruction methods are divided into parsimony versus Bayesian/likelihood approaches. We compared nine methods for reconstructing ancestral areas for placental mammals. Ambiguous reconstructions were a problem for all methods. Important differences resulted from coding areas based on the geographical ranges of extant species versus the geographical provenance of the oldest fossil for each lineage. Africa and South America were reconstructed as the ancestral areas for Afrotheria and Xenarthra, respectively. Most methods reconstructed Eurasia as the ancestral area for Boreoeutheria, Euarchontoglires and Laurasiatheria. The coincidence of molecular dates for the separation of Afrotheria and Xenarthra at approximately 100 Ma with the plate tectonic sundering of Africa and South America hints at the importance of vicariance in the early history of Placentalia. Dispersal has also been important including the origins of Madagascar's endemic mammal fauna. Further studies will benefit from increased taxon sampling and the application of new ancestral area reconstruction methods.

A new basal caniform (Mammalia: Carnivora) from the middle Eocene of North America and remarks on the phylogeny of early carnivorans

BACKGROUND

Despite a long history of research, the phylogenetic origin and initial diversification of the mammalian crown-group Carnivora remain elusive. Well-preserved fossil materials of basal carnivorans are essential for resolving these issues, and for constraining the timing of the carnivoran origin, which constitutes an important time-calibration point in mammalian phylogenetics.

METHODOLOGY/PRINCIPAL FINDINGS

A new carnivoramorphan from the middle Eocene of southern California, Lycophocyon hutchisoni, is described. The new taxon exhibits stages of dental and basicranial evolution that are intermediate between earlier carnivoramorphans and the earliest representatives of canoid carnivorans. The evolutionary affinity of the new taxon was determined by a cladistic analysis of previously-published and newly-acquired morphological data for 30 Paleogene carnivoramorphans. The most-parsimonious trees identified L. hutchisoni as a basal caniform carnivoran, and placed (1) Tapocyon robustus, Quercygale angustidens, "Miacis" sylvestris, "M." uintensis, and "M." gracilis inside or outside the Carnivora, (2) nimravids within the Feliformia, and (3) the amphicyonid Daphoenus outside the crown-group Canoidea. Parsimony reconstructions of ancestral character states suggest that loss of the upper third molars and development of well-ossified entotympanics that are firmly fused to the basicranium (neither condition is observed in L. hutchisoni) are not associated with the origin of the Carnivora as traditionally thought, but instead occurred independently in the Caniformia and the Feliformia. A discriminant analysis of the estimated body weight and dental ecomorphology predicted a mesocarnivorous diet for L. hutchisoni, and the postcranial morphology suggests a scansorial habit.

CONCLUSIONS/SIGNIFICANCE

Lycophocyon hutchisoni illuminates the morphological evolution of early caniforms leading to the origin of crown-group canoids. Considerable uncertainty remains with respect to the phylogenetic origin of the Carnivora. The minimum date of caniform-feliform divergence is provisionally suggested to be either 47 million years ago or 38 million years ago, depending on the position of "Miacis" sylvestris within or outside the Carnivora, respectively.

Signals of recent spatial expansions in the grey mouse lemur (Microcebus murinus)

Background

Pleistocene events have shaped the phylogeography of many taxa worldwide. Their genetic signatures in tropical species have been much less explored than in those living in temperate regions. We analysed the genetic structure of a Malagasy primate species, a mouse lemur with a wide distribution (M. murinus), in order to investigate such phylogeographic processes on a large tropical island. We also evaluated the effects of anthropogenic pressures (fragmentation/deforestation) and natural features (geographic distance, rivers) on genetic structure in order to complement our understanding of past and present processes of genetic differentiation.

Results

The analysis of the mitochondrial D-loop sequences of 195 samples from 15 study sites (10 from a continuous forest and five from isolated forest fragments) from two adjacent Inter-River-Systems (IRSs) revealed that forest fragmentation and the river restrict gene flow, thereby leading to an increased genetic differentiation between populations beyond the effect of isolation-by-distance. Demographic simulations detected signals of two successive spatial expansions that could be preliminarily dated to the late Pleistocene and early Holocene. The haplotype network revealed geographic structure and showed deep molecular divergences within and between the IRSs that would be congruent with a two-step colonization scenario.

Conclusions

This study supports the hypothesis of a relatively recent spatial expansion of the grey mouse lemur in northwestern Madagascar, which may also explain why this taxon, in contrast to its congeners, has not yet undergone allopatric speciation in the studied area and possibly across its presently wide range.

Nomenclature and placental mammal phylogeny

An issue arising from recent progress in establishing the placental mammal Tree of Life concerns the nomenclature of high-level clades. Fortunately, there are now several well-supported clades among extant mammals that require unambiguous, stable names. Although the International Code of Zoological Nomenclature does not apply above the Linnean rank of family, and while consensus on the adoption of competing systems of nomenclature does not yet exist, there is a clear, historical basis upon which to arbitrate among competing names for high-level mammalian clades. Here, we recommend application of the principles of priority and stability, as laid down by G.G. Simpson in 1945, to discriminate among proposed names for high-level taxa. We apply these principles to specific cases among placental mammals with broad relevance for taxonomy, and close with particular emphasis on the Afrotherian family Tenrecidae. We conclude that no matter how reconstructions of the Tree of Life change in years to come, systematists should apply new names reluctantly, deferring to those already published and maximizing consistency with existing nomenclature.

Slow age-dependent decline of doublecortin expression and BrdU labeling in the forebrain from lesser hedgehog tenrecs

In addition to synaptic remodeling, formation of new neurons is increasingly acknowledged as an important cue for plastic changes in the central nervous system. Whereas all vertebrates retain a moderate neuroproliferative capacity, phylogenetically younger mammals become dramatically impaired in this potential during aging. The present study shows that the lesser hedgehog tenrec, an insectivore with a low encephalization index, preserves its neurogenic potential surprisingly well during aging. This was shown by quantitative analysis of 5-bromo-2'-deoxyuridine (BrdU) immunolabeling in the olfactory bulb, paleo-, archi-, and neocortices from 2- to 7-year-old animals. In addition to these newly born cells, a large number of previously formed immature neurons are present throughout adulthood as shown by doublecortin (DCX) immunostaining in various forebrain regions including archicortex, paleocortex, nucleus accumbens, and amygdala. Several ventricle-associated cells in olfactory bulb and hippocampus were double-labeled by BrdU and DCX immunoreactivity. However, most DCX cells in the paleocortex can be considered as persisting immature neurons that obviously do not enter a differentiation program since double fluorescence labeling does not reveal their co-occurrence with numerous neuronal markers, whereas only a small portion coexpresses the pan-neuronal marker HuC/D. Finally, the present study reveals tenrecs as suitable laboratory animals to study age-dependent brain alterations (e.g., of neurogenesis) or slow degenerative processes, particularly due to the at least doubled longevity of tenrecs in comparison to mice and rats.

A phylogenetic estimate for golden moles (Mammalia, Afrotheria, Chrysochloridae)

BACKGROUND

Golden moles (Chrysochloridae) are small, subterranean, afrotherian mammals from South Africa and neighboring regions. Of the 21 species now recognized, some (e.g., Chrysochloris asiatica, Amblysomus hottentotus) are relatively common, whereas others (e.g., species of Chrysospalax, Cryptochloris, Neamblysomus) are rare and endangered. Here, we use a combined analysis of partial sequences of the nuclear GHR gene and morphological characters to derive a phylogeny of species in the family Chrysochloridae.

RESULTS

Although not all nodes of the combined analysis have high support values, the overall pattern of relationships obtained from different methods of phylogeny reconstruction allow us to make several recommendations regarding the current taxonomy of golden moles. We elevate Huetia to generic status to include the species leucorhinus and confirm the use of the Linnean binomial Carpitalpa arendsi, which belongs within Amblysominae along with Amblysomus and Neamblysomus. A second group, Chrysochlorinae, includes Chrysochloris, Cryptochloris, Huetia, Eremitalpa, Chrysospalax, and Calcochloris. Bayesian methods make chrysochlorines paraphyletic by placing the root within them, coinciding with root positions favored by a majority of randomly-generated outgroup taxa. Maximum Parsimony (MP) places the root either between chrysochlorines and amblysomines (with Chlorotalpa as sister taxon to amblysomines), or at Chlorotalpa, with the former two groups reconstructed as monophyletic in all optimal MP trees.

CONCLUSIONS

The inclusion of additional genetic loci for this clade is important to confirm our taxonomic results and resolve the chrysochlorid root. Nevertheless, our optimal topologies support a division of chrysochlorids into amblysomines and chrysochlorines, with Chlorotalpa intermediate between the two. Furthermore, evolution of the chrysochlorid malleus exhibits homoplasy. The elongate malleus has evolved just once in the Cryptochloris-Chrysochloris group; other changes in shape have occurred at multiple nodes, regardless of how the root is resolved.

Performance of relaxed-clock methods in estimating evolutionary divergence times and their credibility intervals

The rapid expansion of sequence data and the development of statistical approaches that embrace varying evolutionary rates among lineages have encouraged many more investigators to use DNA and protein data to time species divergences. Here, we report results from a systematic evaluation, by means of computer simulation, of the performance of two frequently used relaxed-clock methods for estimating these times and their credibility intervals (CrIs). These relaxed-clock methods allow rates to vary in a phylogeny randomly over lineages (e.g., BEAST software) and in autocorrelated fashion (e.g., MultiDivTime software). We applied these methods for analyzing sequence data sets simulated using naturally derived parameters (evolutionary rates, sequence lengths, and base substitution patterns) and assuming that clock calibrations are known without error. We find that the estimated times are, on average, close to the true times as long as the assumed model of lineage rate changes matches the actual model. The 95% CrIs also contain the true time for >or=95% of the simulated data sets. However, the use of incorrect lineage rate model reduces this frequency to 83%, indicating that the relaxed-clock methods are not robust to the violation of underlying lineage rate model. Because these rate models are rarely known a priori and are difficult to detect empirically, we suggest building composite CrIs using CrIs produced from MultiDivTime and BEAST analysis. These composite CrIs are found to contain the true time for >or=97% data sets. Our analyses also verify the usefulness of the common practice of interpreting the congruence of times inferred from different methods as a reflection of the accuracy of time estimates. Overall, our results show that simple strategies can be used to enhance our ability to estimate times and their CrIs when using the relaxed-clock methods.