Chromosome-length genome assemblies and cytogenomic analyses of pangolins reveal remarkable chromosome counts and plasticity.

We report the first chromosome-length genome assemblies for three species in the mammalian order Pholidota: the white-bellied, Chinese, and Sunda pangolins. Surprisingly, we observe extraordinary karyotypic plasticity within this order and, in female white-bellied pangolins, the largest number of chromosomes reported in a Laurasiatherian mammal: 2n = 114. We perform the first karyotype analysis of an African pangolin and report a Y-autosome fusion in white-bellied pangolins, resulting in 2n = 113 for males. We employ a novel strategy to confirm the fusion and identify the autosome involved by finding the pseudoautosomal region (PAR) in the female genome assembly and analyzing the 3D contact frequency between PAR sequences and the rest of the genome in male and female white-bellied pangolins. Analyses of genetic variability show that white-bellied pangolins have intermediate levels of genome-wide heterozygosity relative to Chinese and Sunda pangolins, consistent with two moderate declines of historical effective population size. Our results reveal a remarkable feature of pangolin genome biology and highlight the need for further studies of these unique and endangered mammals.

Waking the undead: Implications of a soft explosive model for the timing of placental mammal diversification.

The explosive, long fuse, and short fuse models represent competing hypotheses for the timing of placental mammal diversification. Support for the explosive model, which posits both interordinal and intraordinal diversification after the KPg mass extinction, derives from morphological cladistic studies that place Cretaceous eutherians outside of crown Placentalia. By contrast, most molecular studies favor the long fuse model wherein interordinal cladogenesis occurred in the Cretaceous followed by intraordinal cladogenesis after the KPg boundary. Phillips (2016) proposed a soft explosive model that allows for the emergence of a few lineages (Xenarthra, Afrotheria, Euarchontoglires, Laurasiatheria) in the Cretaceous, but otherwise agrees with the explosive model in positing the majority of interordinal diversification after the KPg mass extinction. Phillips (2016) argues that rate transference errors associated with large body size and long lifespan have inflated previous estimates of interordinal divergence times, and further suggests that most interordinal divergences are positioned after the KPg boundary when rate transference errors are avoided through the elimination of calibrations in large-bodied and/or long lifespan clades. Here, we show that rate transference errors can also occur in the opposite direction and drag forward estimated divergence dates when calibrations in large-bodied/long lifespan clades are omitted. This dragging forward effect results in the occurrence of more than half a billion years of 'zombie lineages' on Phillips' preferred timetree. By contrast with ghost lineages, which are a logical byproduct of an incomplete fossil record, zombie lineages occur when estimated divergence dates are younger than the minimum age of the oldest crown fossils. We also present the results of new timetree analyses that address the rate transference problem highlighted by Phillips (2016) by deleting taxa that exceed thresholds for body size and lifespan. These analyses recover all interordinal divergence times in the Cretaceous and are consistent with the long fuse model of placental diversification. Finally, we outline potential problems with morphological cladistic analyses of higher-level relationships among placental mammals that may account for the perceived discrepancies between molecular and paleontological estimates of placental divergence times.

Resolution of a concatenation/coalescence kerfuffle: partitioned coalescence support and a robust family-level tree for Mammalia.

Recent phylogenetic analyses of a large dataset for mammalian families (169 taxa, 26 loci) portray contrasting results. Supermatrix (concatenation) methods support a generally robust tree with only a few inconsistently resolved polytomies, whereas MP-EST coalescence analysis of the same dataset yields a weakly supported tree that conflicts with many traditionally recognized clades. Here, we evaluate this discrepancy via improved coalescence analyses with reference to the rich history of phylogenetic studies on mammals. This integration clearly demonstrates that both supermatrix and coalescence analyses of just 26 loci yield a congruent, well-supported phylogenetic hypothesis for Mammalia. Discrepancies between published studies are explained by implementation of overly simple DNA substitution models, inadequate tree-search routines and limitations of the MP-EST method. We develop a simple measure, partitioned coalescence support (PCS), which summarizes the distribution of support and conflict among gene trees for a given clade. Extremely high PCS scores for outlier gene trees at two nodes in the mammalian tree indicate a troubling bias in the MP-EST method. We conclude that in this age of phylogenomics, a solid understanding of systematics fundamentals, choice of valid methodology and a broad knowledge of a clade's taxonomic history are still required to yield coherent phylogenetic inferences.

chromDraw: an R package for visualization of linear and circular karyotypes.

Species-specific sets of chromosomes-karyotypes-are traditionally depicted as linear ideograms with individual chromosomes represented by vertical bars. However, linear visualization has its limitations when the shared collinearity and/or chromosomal rearrangements differentiating two or more karyotypes need to be demonstrated. In these instances, circular visualization might provide easier comprehension and interpretation of inter-species chromosomal collinearity. The chromDraw graphical tool was developed as a user-friendly graphical tool for visualizing both linear and circular karyotypes based on the same input data matrix. The output graphics, saved in two different formats (EPS and SVG), can be easily imported to and modified in presentation and image-editing computer programs. The tool is freely distributed under GNU General Public License (GPL) and can be installed from Bioconductor or from the chromDraw home page.

Range-Wide Snow Leopard Phylogeography Supports Three Subspecies.

The snow leopard, Panthera uncia, is an elusive high-altitude specialist that inhabits vast, inaccessible habitat across Asia. We conducted the first range-wide genetic assessment of snow leopards based on noninvasive scat surveys. Thirty-three microsatellites were genotyped and a total of 683 bp of mitochondrial DNA sequenced in 70 individuals. Snow leopards exhibited low genetic diversity at microsatellites (AN = 5.8, HO = 0.433, HE = 0.568), virtually no mtDNA variation, and underwent a bottleneck in the Holocene (∼8000 years ago) coinciding with increased temperatures, precipitation, and upward treeline shift in the Tibetan Plateau. Multiple analyses supported 3 primary genetic clusters: (1) Northern (the Altai region), (2) Central (core Himalaya and Tibetan Plateau), and (3) Western (Tian Shan, Pamir, trans-Himalaya regions). Accordingly, we recognize 3 subspecies, Panthera uncia irbis (Northern group), Panthera uncia uncia (Western group), and Panthera uncia uncioides (Central group) based upon genetic distinctness, low levels of admixture, unambiguous population assignment, and geographic separation. The patterns of variation were consistent with desert-basin "barrier effects" of the Gobi isolating the northern subspecies (Mongolia), and the trans-Himalaya dividing the central (Qinghai, Tibet, Bhutan, and Nepal) and western subspecies (India, Pakistan, Tajikistan, and Kyrgyzstan). Hierarchical Bayesian clustering analysis revealed additional subdivision into a minimum of 6 proposed management units: western Mongolia, southern Mongolia, Tian Shan, Pamir-Himalaya, Tibet-Himalaya, and Qinghai, with spatial autocorrelation suggesting potential connectivity by dispersing individuals up to ∼400 km. We provide a foundation for global conservation of snow leopard subspecies, and set the stage for in-depth landscape genetics and genomic studies.

A cytogenetic and comparative map of camelid chromosome 36 and the minute in alpacas.

Recent advances in camelid genomics have provided draft sequence assemblies and the first comparative and gene maps for the dromedary (CDR) and the alpaca (LPA). However, no map information is currently available for the smallest camelid autosome-chr36. The chromosome is also of clinical interest because of its involvement in the minute chromosome syndrome (MCS) in infertile alpacas. Here, we developed molecular markers for camelid chr36 by direct sequencing CDR36 and LPA minute and by bioinformatics analysis of alpaca unplaced sequence scaffolds. We constructed a cytogenetic map for chr36 in the alpaca, llama, and dromedary and showed its homology to human chromosome 7 (HSA7) at 49.8-55.5 Mb. The chr36 map comprised seven markers, including two genes-ZPBP and WVC2. Comparative status of HSA7 was further refined by cytogenetic mapping of 16 HSA7 orthologs in camelid chromosomes 7 and 18 and by the analysis of HSA7-conserved synteny blocks across 11 vertebrate species. Finally, mapping chr36 markers in infertile alpacas confirmed that the minute chromosome was a derivative of chr36, but the small size was not a result of a large deletion or a translocation. Instead, cytogenetic mapping of 5.8S, 18S, and 28S rRNA genes (nucleolus organizer region (NOR)) revealed that the size difference between chr36 homologs in infertile alpacas was due to a heterozygous presence of NOR, whereas chr36 in fertile alpacas had no NOR. We theorized that the heterozygous NOR might affect chr36 pairing, recombination, and segregation in meiosis and, thus fertility.

Genetically based low oxygen affinities of felid hemoglobins: lack of biochemical adaptation to high-altitude hypoxia in the snow leopard.

Genetically based modifications of hemoglobin (Hb) function that increase blood-O2 affinity are hallmarks of hypoxia adaptation in vertebrates. Among mammals, felid Hbs are unusual in that they have low intrinsic O2 affinities and reduced sensitivities to the allosteric cofactor 2,3-diphosphoglycerate (DPG). This combination of features compromises the acclimatization capacity of blood-O2 affinity and has led to the hypothesis that felids have a restricted physiological niche breadth relative to other mammals. In seeming defiance of this conjecture, the snow leopard (Panthera uncia) has an extraordinarily broad elevational distribution and occurs at elevations above 6000 m in the Himalayas. Here, we characterized structural and functional variation of big cat Hbs and investigated molecular mechanisms of Hb adaptation and allosteric regulation that may contribute to the extreme hypoxia tolerance of the snow leopard. Experiments revealed that purified Hbs from snow leopard and African lion exhibited equally low O2 affinities and DPG sensitivities. Both properties are primarily attributable to a single amino acid substitution, ß2His→Phe, which occurred in the common ancestor of Felidae. Given the low O2 affinity and reduced regulatory capacity of feline Hbs, the extreme hypoxia tolerance of snow leopards must be attributable to compensatory modifications of other steps in the O2-transport pathway.

Range-wide multilocus phylogeography of the red fox reveals ancient continental divergence, minimal genomic exchange and distinct demographic histories.

Widely distributed taxa provide an opportunity to compare biogeographic responses to climatic fluctuations on multiple continents and to investigate speciation. We conducted the most geographically and genomically comprehensive study to date of the red fox (Vulpes vulpes), the world's most widely distributed wild terrestrial carnivore. Analyses of 697 bp of mitochondrial sequence in ~1000 individuals suggested an ancient Middle Eastern origin for all extant red foxes and a 400 kya (SD = 139 kya) origin of the primary North American (Nearctic) clade. Demographic analyses indicated a major expansion in Eurasia during the last glaciation (~50 kya), coinciding with a previously described secondary transfer of a single matriline (Holarctic) to North America. In contrast, North American matrilines (including the transferred portion of Holarctic clade) exhibited no signatures of expansion until the end of the Pleistocene (~12 kya). Analyses of 11 autosomal loci from a subset of foxes supported the colonization time frame suggested by mtDNA (and the fossil record) but, in contrast, reflected no detectable secondary transfer, resulting in the most fundamental genomic division of red foxes at the Bering Strait. Endemic continental Y-chromosome clades further supported this pattern. Thus, intercontinental genomic exchange was overall very limited, consistent with long-term reproductive isolation since the initial colonization of North America. Based on continental divergence times in other carnivoran species pairs, our findings support a model of peripatric speciation and are consistent with the previous classification of the North American red fox as a distinct species, V. fulva.

Accelerated evolution of CES7, a gene encoding a novel major urinary protein in the cat family.

Cauxin is a novel urinary protein recently identified in the domestic cat that regulates the excretion of felinine, a pheromone precursor involved in sociochemical communication and territorial marking of domestic and wild felids. Understanding the evolutionary history of cauxin may therefore illuminate molecular adaptations involved in the evolution of pheromone-based communication, recognition, and mate selection in wild animals. We sequenced the gene encoding cauxin, CES7, in 22 species representing all major felid lineages, and multiple outgroups and showed that it has undergone rapid evolutionary change preceding and during the diversification of the cat family. A comparison between feline cauxin and orthologous carboxylesterases from other mammalian lineages revealed evidence of strong positive Darwinian selection within and between several cat lineages, enriched at functionally important sites of the protein. The higher rate of radical amino acid replacements in small felids, coupled with the lack of felinine and extremely low levels of cauxin in the urine of the great cats (Panthera), correlates with functional divergence of this gene in Panthera, and its putative loss in the snow leopard. Expression studies found evidence for several alternatively spliced transcripts in testis and brain, suggesting additional roles in male reproductive fitness and behavior. Our work presents the first report of strong positive natural selection acting on a major urinary protein of nonrodent mammals, providing evidence for parallel selection pressure on the regulation of pheromones in different mammalian lineages, despite the use of different metabolic pathways. Our results imply that natural selection may drive rapid changes in the regulation of pheromones in urine among the different cat species, which in turn may influence social behavior, such as territorial marking and conspecific recognition, therefore serving as an important mechanism for the radiation of this group of mammals.

Pyrosequencing of 16S rRNA genes in fecal samples reveals high diversity of hindgut microflora in horses and potential links to chronic laminitis.

BACKGROUND: The nutrition and health of horses is closely tied to their gastrointestinal microflora. Gut bacteria break down plant structural carbohydrates and produce volatile fatty acids, which are a major source of energy for horses. Bacterial communities are also essential for maintaining gut homeostasis and have been hypothesized to contribute to various diseases including laminitis. We performed pyrosequencing of 16S rRNA bacterial genes isolated from fecal material to characterize hindgut bacterial communities in healthy horses and those with chronic laminitis. RESULTS: Fecal samples were collected from 10 normal horses and 8 horses with chronic laminitis. Genomic DNA was extracted and the V4-V5 segment of the 16S rRNA gene was PCR amplified and sequenced on the 454 platform generating a mean of 2,425 reads per sample after quality trimming. The bacterial communities were dominated by Firmicutes (69.21% control, 56.72% laminitis) and Verrucomicrobia (18.13% control, 27.63% laminitis), followed by Bacteroidetes, Proteobacteria, and Spirochaetes. We observed more OTUs per individual in the laminitis group than the control group (419.6 and 355.2, respectively, P = 0.019) along with a difference in the abundance of two unassigned Clostridiales genera (P = 0.03 and P = 0.01). The most abundant bacteria were Streptococcus spp., Clostridium spp., and Treponema spp.; along with unassigned genera from Subdivision 5 of Verrucomicrobia, Ruminococcaceae, and Clostridiaceae, which together constituted ~ 80% of all OTUs. There was a high level of individual variation across all taxonomic ranks. CONCLUSIONS: Our exploration of the equine fecal microflora revealed higher bacterial diversity in horses with chronic laminitis and identification of two Clostridiales genera that differed in abundance from control horses. There was large individual variation in bacterial communities that was not explained in our study. The core hindgut microflora was dominated by Streptococcus spp., several cellulytic genera, and a large proportion of uncharacterized OTUs that warrant further investigation regarding their function. Our data provide a foundation for future investigations of hindgut bacterial factors that may influence the development and progression of chronic laminitis.

Status and distribution of jaguarundi in Texas and Northeastern México: Making the case for extirpation and initiation of recovery in the United States.

The jaguarundi (Puma yagouaroundi) is a small felid with a historical range from central Argentina through southern Texas. Information on the current distribution of this reclusive species is needed to inform recovery strategies in the United States where its last record was in 1986 in Texas. From 2003 to 2021, we conducted camera-trap surveys across southern Texas and northern Tamaulipas, México to survey for medium-sized wild cats (i.e., ocelots [Leopardus pardalis], bobcats [Lynx rufus], and jaguarundi). After 350,366 trap nights at 685 camera sites, we did not detect jaguarundis at 16 properties or along 2 highways (1050 km2) in Texas. However, we recorded 126 jaguarundi photographic detections in 15,784 trap nights on 2 properties (125.3 km2) in the northern Sierra of Tamaulipas, Tamaulipas, México. On these properties, latency to detection was 72 trap nights, with a 0.05 probability of detection per day and 0.73 photographic event rate every 100 trap nights. Due to a lack of confirmed class I sightings (e.g., specimen, photograph) in the 18 years of this study, and no other class I observations since 1986 in the United States, we conclude that the jaguarundi is likely extirpated from the United States. Based on survey effort and results from México, we would have expected to detect jaguarundis over the course of the study if still extant in Texas. We recommend that state and federal agencies consider jaguarundis as extirpated from the United States and initiate recovery actions as mandated in the federal jaguarundi recovery plan. These recovery actions include identification of suitable habitat in Texas, identification of robust populations in México, and re-introduction of the jaguarundi to Texas.

Multiscale patterns of isolation by ecology and fine-scale population structure in Texas bobcats.

Patterns of spatial genetic variation can be generated by a variety of ecological processes, including individual preferences based on habitat. These ecological processes act at multiple spatial and temporal scales, generating scale-dependent effects on gene flow. In this study, we focused on bobcats (Lynx rufus), a highly mobile, generalist felid that exhibits ecological and behavioral plasticity, high abundance, and broad connectivity across much of their range. However, bobcats also show genetic differentiation along habitat breaks, a pattern typically observed in cases of isolation-by-ecology (IBE). The IBE observed in bobcats is hypothesized to occur due to habitat-biased dispersal, but it is unknown if this occurs at other habitat breaks across their range or at what spatial scale IBE becomes most apparent. Thus, we used a multiscale approach to examine isolation by ecology (IBE) patterns in bobcats (Lynx rufus) at both fine and broad spatial scales in western Texas. We genotyped 102 individuals at nine microsatellite loci and used partial redundancy analysis (pRDA) to test if a suite of landscape variables influenced genetic variation in bobcats. Bobcats exhibited a latitudinal cline in population structure with a spatial signature of male-biased dispersal, and no clear barriers to gene flow. Our pRDA tests revealed high genetic similarity in similar habitats, and results differed by spatial scale. At the fine spatial scale, herbaceous rangeland was an important influence on gene flow whereas mixed rangeland and agriculture were significant at the broad spatial scale. Taken together, our results suggests that complex interactions between spatial-use behavior and landscape heterogeneity can create non-random gene flow in highly mobile species like bobcats. Furthermore, our results add to the growing body of data highlighting the importance of multiscale study designs when assessing spatial genetic structure.

Methane, arsenic, selenium and the origins of the DMSO reductase family.

Mononuclear molybdoenzymes of the dimethyl sulfoxide reductase (DMSOR) family catalyze a number of reactions essential to the carbon, nitrogen, sulfur, arsenic, and selenium biogeochemical cycles. These enzymes are also ancient, with many lineages likely predating the divergence of the last universal common ancestor into the Bacteria and Archaea domains. We have constructed rooted phylogenies for over 1,550 representatives of the DMSOR family using maximum likelihood methods to investigate the evolution of the arsenic biogeochemical cycle. The phylogenetic analysis provides compelling evidence that formylmethanofuran dehydrogenase B subunits, which catalyze the reduction of CO2 to formate during hydrogenotrophic methanogenesis, constitutes the most ancient lineage. Our analysis also provides robust support for selenocysteine as the ancestral ligand for the Mo/W atom. Finally, we demonstrate that anaerobic arsenite oxidase and respiratory arsenate reductase catalytic subunits represent a more ancient lineage of DMSORs compared to aerobic arsenite oxidase catalytic subunits, which evolved from the assimilatory nitrate reductase lineage. This provides substantial support for an active arsenic biogeochemical cycle on the anoxic Archean Earth. Our work emphasizes that the use of chalcophilic elements as substrates as well as the Mo/W ligand in DMSORs has indelibly shaped the diversification of these enzymes through deep time.

The horse Y chromosome as an informative marker for tracing sire lines.

Analysis of the Y chromosome is the best-established way to reconstruct paternal family history in humans. Here, we applied fine-scaled Y-chromosomal haplotyping in horses with biallelic markers and demonstrate the potential of our approach to address the ancestry of sire lines. We de novo assembled a draft reference of the male-specific region of the Y chromosome from Illumina short reads and then screened 5.8 million basepairs for variants in 130 specimens from intensively selected and rural breeds and nine Przewalski's horses. Among domestic horses we confirmed the predominance of a young'crown haplogroup' in Central European and North American breeds. Within the crown, we distinguished 58 haplotypes based on 211 variants, forming three major haplogroups. In addition to two previously characterised haplogroups, one observed in Arabian/Coldblooded and the other in Turkoman/Thoroughbred horses, we uncovered a third haplogroup containing Iberian lines and a North African Barb Horse. In a genealogical showcase, we distinguished the patrilines of the three English Thoroughbred founder stallions and resolved a historic controversy over the parentage of the horse 'Galopin', born in 1872. We observed two nearly instantaneous radiations in the history of Central and Northern European Y-chromosomal lineages that both occurred after domestication 5,500 years ago.

Horse Y chromosome assembly displays unique evolutionary features and putative stallion fertility genes.

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility.

Large Deletions at the SHOX Locus in the Pseudoautosomal Region Are Associated with Skeletal Atavism in Shetland Ponies.

Skeletal atavism in Shetland ponies is a heritable disorder characterized by abnormal growth of the ulna and fibula that extend the carpal and tarsal joints, respectively. This causes abnormal skeletal structure and impaired movements, and affected foals are usually killed. In order to identify the causal mutation we subjected six confirmed Swedish cases and a DNA pool consisting of 21 control individuals to whole genome resequencing. We screened for polymorphisms where the cases and the control pool were fixed for opposite alleles and observed this signature for only 25 SNPs, most of which were scattered on genome assembly unassigned scaffolds. Read depth analysis at these loci revealed homozygosity or compound heterozygosity for two partially overlapping large deletions in the pseudoautosomal region (PAR) of chromosome X/Y in cases but not in the control pool. One of these deletions removes the entire coding region of the SHOX gene and both deletions remove parts of the CRLF2 gene located downstream of SHOX. The horse reference assembly of the PAR is highly fragmented, and in order to characterize this region we sequenced bacterial artificial chromosome (BAC) clones by single-molecule real-time (SMRT) sequencing technology. This considerably improved the assembly and enabled size estimations of the two deletions to 160-180 kb and 60-80 kb, respectively. Complete association between the presence of these deletions and disease status was verified in eight other affected horses. The result of the present study is consistent with previous studies in humans showing crucial importance of SHOX for normal skeletal development.

Loss of genetic diversity among ocelots in the United States during the 20th century linked to human induced population reductions.

Ocelots (Leopardus pardalis) in the United States currently exhibit low levels of genetic diversity. One hypothesis for this observation is that habitat fragmentation, resulting from human induced changes in the landscape during the 20(th) century, created island populations with highly reduced gene flow and increased genetic drift and inbreeding. In an effort to investigate this, we used a portion of the mitochondrial control region and 11 autosomal microsatellite loci to examine historical levels of genetic diversity and infer temporal changes in ocelot populations between 1853 and 2005. Levels of genetic diversity were higher in historical ocelot populations than in extant populations from Texas. The earliest documented loss of mitochondrial haplotype diversity occurred at Laguna Atascosa National Wildlife Refuge. The second extant population inhabiting private lands in Willacy County retained higher levels of genetic diversity through the 1990s, but subsequently lost diversity over the next decade. A similar pattern was observed for autosomal microsatellite loci. This supports the argument that low levels of genetic diversity in Texas are related to human induced population reductions and fragmentation, both of which threaten the remaining ocelots in the United States. At this time, the best means of mitigating the continued erosion of genetic variation are translocation of individuals either from larger populations in Mexico to Texas, or between the Texas populations.

Exploring the correlations between sequence evolution rate and phenotypic divergence across the Mammalian tree provides insights into adaptive evolution.

Sequence evolution behaves in a relatively consistent manner, leading to one of the fundamental paradigms in biology, the existence of a 'molecular clock'. The molecular clock can be distilled to the concept of accumulation of substitutions, through time yielding a stable rate from which we can estimate lineage divergence. Over the last 50 years, evolutionary biologists have obtained an in-depth understanding of this clock's nuances. It has been fine-tuned by taking into account the vast heterogeneity in rates across lineages and genes, leading to 'relaxed' molecular clock methods for timetree reconstruction. Sequence rate varies with life history traits including body size, generation time and metabolic rate, and we review recent studies on this topic. However, few studies have explicitly examined correlates between molecular evolution and morphological evolution. The patterns observed across diverse lineages suggest that rates of molecular and morphological evolution are largely decoupled. We discuss how identifying the molecular mechanisms behind rapid functional radiations are central to understanding evolution. The vast functional divergence within mammalian lineages that have relatively 'slow' sequence evolution refutes the hypotheses that pulses in diversification yielding major phenotypic change are the result of steady accumulation of substitutions. Patterns rather suggest phenotypic divergence is likely caused by regulatory alterations mediated through mechanisms such as insertions/deletions in functional regions. These can rapidly arise and sweep to fixation faster than predicted from a lineage's sequence neutral substitution rate, enabling species to leapfrog between phenotypic 'islands'. We suggest research directions that could illuminate mechanisms behind the functional diversity we see today.

Macroevolutionary dynamics and historical biogeography of primate diversification inferred from a species supermatrix.

Phylogenetic relationships, divergence times, and patterns of biogeographic descent among primate species are both complex and contentious. Here, we generate a robust molecular phylogeny for 70 primate genera and 367 primate species based on a concatenation of 69 nuclear gene segments and ten mitochondrial gene sequences, most of which were extracted from GenBank. Relaxed clock analyses of divergence times with 14 fossil-calibrated nodes suggest that living Primates last shared a common ancestor 71-63 Ma, and that divergences within both Strepsirrhini and Haplorhini are entirely post-Cretaceous. These results are consistent with the hypothesis that the Cretaceous-Paleogene mass extinction of non-avian dinosaurs played an important role in the diversification of placental mammals. Previous queries into primate historical biogeography have suggested Africa, Asia, Europe, or North America as the ancestral area of crown primates, but were based on methods that were coopted from phylogeny reconstruction. By contrast, we analyzed our molecular phylogeny with two methods that were developed explicitly for ancestral area reconstruction, and find support for the hypothesis that the most recent common ancestor of living Primates resided in Asia. Analyses of primate macroevolutionary dynamics provide support for a diversification rate increase in the late Miocene, possibly in response to elevated global mean temperatures, and are consistent with the fossil record. By contrast, diversification analyses failed to detect evidence for rate-shift changes near the Eocene-Oligocene boundary even though the fossil record provides clear evidence for a major turnover event ("Grande Coupure") at this time. Our results highlight the power and limitations of inferring diversification dynamics from molecular phylogenies, as well as the sensitivity of diversification analyses to different species concepts.