EMMLi: A maximum likelihood approach to the analysis of modularity

Identification of phenotypic modules, semi-autonomous sets of highly-correlated traits, can be accomplished through exploratory (e.g., cluster analysis) or confirmatory approaches (e.g., RV coefficient analysis). While statistically more robust, confirmatory approaches are generally unable to compare across different model structures. For example, RV coefficient analysis finds support for both two- and six-module models for the therian mammalian skull. Here, we present a maximum likelihood approach that takes into account model parameterization. We compare model log-likelihoods of trait correlation matrices using the finite-sample corrected Akaike Information Criterion, allowing for comparison of hypotheses across different model structures. Simulations varying model complexity and within- and between-module contrast demonstrate that this method correctly identifies model structure and parameters across a wide range of conditions. We further analyzed a dataset of 3-D data, consisting of 61 landmarks from 181 macaque (Macaca fuscata) skulls, distributed among five age categories, testing 31 models, including no modularity among the landmarks, and various partitions of 2, 3, 6, and 8 modules. Our results clearly support a complex six-module model, with separate within- and inter-module correlations. Furthermore, this model was selected for all five age categories, demonstrating that this complex pattern of integration in the macaque skull appears early and is highly conserved throughout postnatal ontogeny. Subsampling analyses demonstrate that this method is robust to relatively low sample sizes, as is commonly encountered in rare or extinct taxa. This new approach allows for the direct comparison of models with different parameterizations, providing an important tool for the analysis of modularity across diverse systems. This article is protected by copyright. All rights reserved

Complex body size trends in the evolution of sloths (Xenarthra: Pilosa)

Background Extant sloths present an evolutionary conundrum in that the two living genera are superficially similar (small-bodied, folivorous, arboreal) but diverged from one another approximately 30 million years ago and are phylogenetically separated by a radiation of medium to massive, mainly ground-dwelling, taxa. Indeed, the species in the two living genera are among the smallest, and perhaps most unusual, of the 50+ known sloth species, and must have independently and convergently evolved small size and arboreality. In order to accurately reconstruct sloth evolution, it is critical to incorporate their extinct diversity in analyses. Here, we used a dataset of 57 species of living and fossil sloths to examine changes in body mass mean and variance through their evolution, employing a general time-variable model that allows for analysis of evolutionary trends in continuous characters within clades lacking fully-resolved phylogenies, such as sloths. Results Our analyses supported eight models, all of which partition sloths into multiple subgroups, suggesting distinct modes of body size evolution among the major sloth lineages. Model-averaged parameter values supported trended walks in most clades, with estimated rates of body mass change ranging as high as 126 kg/million years for the giant ground sloth clades Megatheriidae and Nothrotheriidae. Inclusion of living sloth species in the analyses weakened reconstructed rates for their respective groups, with estimated rates for Megalonychidae (large to giant ground sloths and the extant two-toed sloth) were four times higher when the extant genus Choloepus was excluded. Conclusions Analyses based on extant taxa alone have the potential to oversimplify or misidentify macroevolutionary patterns. This study demonstrates the impact that integration of data from the fossil record can have on reconstructions of character evolution and establishes that body size evolution in sloths was complex, but dominated by trended walks towards the enormous sizes exhibited in some recently extinct forms

A symmoriiform chondrichthyan braincase and the origin of chimaeroid fishes

Chimaeroid fishes (Holocephali) are one of the four principal divisions of modern gnathostomes (jawed vertebrates). Despite only 47 described living species1, chimaeroids are the focus of resurgent interest as potential archives of genomic data2 and for the unique perspective they provide on chondrichthyan and gnathostome ancestral conditions. Chimaeroids are also noteworthy for their highly derived body plan1,3,4. However, like other living groups with distinctive anatomies5, fossils have been of limited use in unravelling their evolutionary origin, as the earliest recognized examples already exhibit many of the specializations present in modern forms6,7. Here we report the results of a computed tomography analysis of Dwykaselachus, an enigmatic chondrichthyan braincase from the ~280 million year old Karoo sediments of South Africa8. Externally, the braincase is that of a symmoriid shark9,10,11,12,13and is by far the most complete uncrushed example yet discovered. Internally, the morphology exhibits otherwise characteristically chimaeroid specializations, including the otic labyrinth arrangement and the brain space configuration relative to exceptionally large orbits. These results have important implications for our view of modern chondrichthyan origins, add robust structure to the phylogeny of early crown group gnathostomes, reveal preconditions that suggest an initial morpho-functional basis for the derived chimaeroid cranium, and shed new light on the chondrichthyan response to the extinction at the end of the Devonian period

Mega-evolutionary dynamics of the adaptive radiation of birds

The origin and expansion of biological diversity is regulated by both developmental trajectories and limits on available ecological niches. As lineages diversify, an early and often rapid phase of species and trait proliferation gives way to evolutionary slow- downs as new species pack into ever more densely occupied regions of ecological niche space. Small clades such as Darwin’s finches demonstrate that natural selection is the driving force of adaptive radiations, but how microevolutionary processes scale up to shape the expansion of phenotypic diversity over much longer evolutionary timescales is unclear. Here we address this problem on a global scale by analysing a crowd-sourced dataset of three-dimensional scanned bill morphology from more than 2,000 species. We find that bill diversity expanded early in extant avian evolutionary history, before transitioning to a phase dominated by packing of morphological space. However, this early phenotypic diversification is decoupled from temporal variation in evolutionary rate: rates of bill evolution vary among lineages but are comparatively stable through time. We find that rare, but major, discontinuities in phenotype emerge from rapid increases in rate along single branches, sometimes leading to depauperate clades with unusual bill morphologies. Despite these jumps between groups, the major axes of within-group bill-shape evolution are remarkably consistent across birds. We reveal that macroevolutionary processes underlying global-scale adaptive radiations support Darwinian and Simpsonian ideas of microevolution within adaptive zones and accelerated evolution between distinct adaptive peaks

Multiple Determinants of Whole and Regional Brain Volume among Terrestrial Carnivorans

Mammalian brain volumes vary considerably, even after controlling for body size. Although several hypotheses have been proposed to explain this variation, most research in mammals on the evolution of encephalization has focused on primates, leaving the generality of these explanations uncertain. Furthermore, much research still addresses only one hypothesis at a time, despite the demonstrated importance of considering multiple factors simultaneously. We used phylogenetic comparative methods to investigate simultaneously the importance of several factors previously hypothesized to be important in neural evolution among mammalian carnivores, including social complexity, forelimb use, home range size, diet, life history, phylogeny, and recent evolutionary changes in body size. We also tested hypotheses suggesting roles for these variables in determining the relative volume of four brain regions measured using computed tomography. Our data suggest that, in contrast to brain size in primates, carnivoran brain size may lag behind body size over evolutionary time. Moreover, carnivore species that primarily consume vertebrates have the largest brains. Although we found no support for a role of social complexity in overall encephalization, relative cerebrum volume correlated positively with sociality. Finally, our results support negative relationships among different brain regions after accounting for overall endocranial volume, suggesting that increased size of one brain regions is often accompanied by reduced size in other regions rather than overall brain expansion

Evolution of Stress-Regulated Gene Expression in Duplicate Genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >∼0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments

Ancestral State Reconstruction Reveals Rampant Homoplasy of Diagnostic Morphological Characters in Urticaceae, Conflicting with Current Classification Schemes

Urticaceae is a family with more than 2000 species, which contains remarkable morphological diversity. It has undergone many taxonomic reorganizations, and is currently the subject of further systematic studies. To gain more resolution in systematic studies and to better understand the general patterns of character evolution in Urticaceae, based on our previous phylogeny including 169 accessions comprising 122 species across 47 Urticaceae genera, we examined 19 diagnostic characters, and analysed these employing both maximum-parsimony and maximum-likelihood approaches. Our results revealed that 16 characters exhibited multiple state changes within the family, with ten exhibiting >eight changes and three exhibiting between 28 and 40. Morphological synapomorphies were identified for many clades, but the diagnostic value of these was often limited due to reversals within the clade and/or homoplasies elsewhere. Recognition of the four clades comprising the family at subfamily level can be supported by a small number carefully chosen defining traits for each. Several non-monophyletic genera appear to be defined only by characters that are plesiomorphic within their clades, and more detailed work would be valuable to find defining traits for monophyletic clades within these. Some character evolution may be attributed to adaptive evolution in Urticaceae due to shifts in habitat or vegetation type. This study demonstrated the value of using phylogeny to trace character evolution, and determine the relative importance of morphological traits for classification