Search

Cross-species introgression can have significant impacts on phylogenomic reconstruction of species divergence events. Here, we used simulations to show how the presence of even a small amount of introgression can bias divergence time estimates when gene flow is ignored in the analysis. Using advances in analytical methods under the multispecies coalescent (MSC) model, we demonstrate that by accounting for incomplete lineage sorting and introgression using large phylogenomic data sets this problem can be avoided. The multispecies-coalescent-with-introgression (MSci) model is capable of accurately estimating both divergence times and ancestral effective population sizes, even when only a single diploid individual per species is sampled. We characterize some general expectations for biases in divergence time estimation under three different scenarios: 1) introgression between sister species, 2) introgression between non-sister species, and 3) introgression from an unsampled (i.e., ghost) outgroup lineage. We also conducted simulations under the isolation-with-migration (IM) model and found that the MSci model assuming episodic gene flow was able to accurately estimate species divergence times despite high levels of continuous gene flow. We estimated divergence times under the MSC and MSci models from two published empirical datasets with previous evidence of introgression, one of 372 target-enrichment loci from baobabs (Adansonia), and another of 1000 transcriptome loci from 14 species of the tomato relative, Jaltomata. The empirical analyses not only confirm our findings from simulations, demonstrating that the MSci model can reliably estimate divergence times but also show that divergence time estimation under the MSC can be robust to the presence of small amounts of introgression in empirical datasets with extensive taxon sampling. [divergence time; gene flow; hybridization; introgression; MSci model; multispecies coalescent]

The occurrence of natural hybridization has been reported in a wide range of organisms, including primates. The present study focuses on the endemic lemurs of Madagascar, primates for which only a few species occur in sympatry or parapatry with congeners, thereby creating limited opportunity for natural hybridization. This study examines RADseq data from 480 individuals to investigate whether the recent expansion of Microcebus murinus towards the northwest and subsequent secondary contact with Microcebus ravelobensis has resulted in the occurrence of hybridization between the two species. Admixture analysis identified one individual with 26% of nuclear admixture, which may correspond to an F2- or F3-hybrid. A composite-likelihood approach was subsequently used to test the fit of alternative phylogeographic scenarios to the genomic data and to date introgression. The simulations yielded support for low levels of gene flow (2Nm0 = 0.063) between the two species starting before the Last Glacial Maximum (between 54 and 142 kyr). Since M. murinus most likely colonized northwestern Madagascar during the Late Pleistocene, the rather recent secondary contact with M. ravelobensis has likely created the opportunity for occasional hybridization. Although reproductive isolation between these distantly related congeners is not complete, it is effective in maintaining species boundaries.

Duke Lemur Centre publication #1510; International audience; Madagascar's Central Highlands are largely composed of grasslands, interspersed with patches of forest. The historical perspective was that Madagascar's grasslands had anthropogenic origins, but emerging evidence suggests that grasslands were a component of the pre-human Central Highlands vegetation. Consequently, there is now vigorous debate regarding the extent to which these grasslands have expanded due to anthropogenic pressures. Here, we shed light on the temporal dynamics of Madagascar's vegetative composition by conducting a population genomic investigation of Goodman's mouse lemur (Microcebus lehilahytsara; Cheirogaleidae). These small-bodied primates occur both in Madagascar's eastern rainforests and in the Central Highlands, making them a valuable indicator species. Population divergences among forest-dwelling mammals will reflect changes to their habitat, including fragmentation, whereas patterns of post-divergence gene flow can reveal formerly wooded migration corridors. To explore these patterns, we used RADseq data to infer population genetic structure, demographic models of post-divergence gene flow, and population size change through time. The results offer evidence that open habitats are an ancient component of the Central Highlands, and that widespread forest fragmentation occurred naturally during a period of decreased precipitation near the last glacial maximum. Models of gene flow suggest that migration across the Central Highlands has been possible from the Pleistocene through the recent Holocene via riparian corridors. Though our findings support the hypothesis that Central Highland grasslands predate human arrival, we also find evidence for human-mediated population declines. This highlights the extent to which species imminently threatened by human-mediated deforestation may already be vulnerable from paleoclimatic conditions.

Sensory gene families are of special interest for both what they can tell us about molecular evolution and what they imply as mediators of social communication. The vomeronasal type-1 receptors (V1Rs) have often been hypothesized as playing a fundamental role in driving or maintaining species boundaries given their likely function as mediators of intraspecific mate choice, particularly in nocturnal mammals. Here, we employ a comparative genomic approach for revealing patterns of V1R evolution within primates, with a special focus on the small-bodied nocturnal mouse and dwarf lemurs of Madagascar (genera Microcebus and Cheirogaleus, respectively). By doubling the existing genomic resources for strepsirrhine primates (i.e. the lemurs and lorises), we find that the highly speciose and morphologically cryptic mouse lemurs have experienced an elaborate proliferation of V1Rs that we argue is functionally related to their capacity for rapid lineage diversification. Contrary to a previous study that found equivalent degrees of V1R diversity in diurnal and nocturnal lemurs, our study finds a strong correlation between nocturnality and V1R elaboration, with nocturnal lemurs showing elaborate V1R repertoires and diurnal lemurs showing less diverse repertoires. Recognized subfamilies among V1Rs show unique signatures of diversifying positive selection, as might be expected if they have each evolved to respond to specific stimuli. Furthermore, a detailed syntenic comparison of mouse lemurs with mouse (genus Mus) and other mammalian outgroups shows that orthologous mammalian subfamilies, predicted to be of ancient origin, tend to cluster in a densely populated region across syntenic chromosomes that we refer to as a V1R “hotspot.”

The published article [...]

Critically endangered sifaka lemurs exhibit high genetic diversity and genomic evidence for dietary adaptation., Sifakas (genus Propithecus) are critically endangered, large-bodied diurnal lemurs that eat leaf-based diets and show corresponding anatomical and microbial adaptations to folivory. We report on the genome assembly of Coquerel’s sifaka (P. coquereli) and the resequenced genomes of Verreaux’s (P. verreauxi), the golden-crowned (P. tattersalli), and the diademed (P. diadema) sifakas. We find high heterozygosity in all sifakas compared with other primates and endangered mammals. Demographic reconstructions nevertheless suggest declines in effective population size beginning before human arrival on Madagascar. Comparative genomic analyses indicate pervasive accelerated evolution in the ancestral sifaka lineage affecting genes in several complementary pathways relevant to folivory, including nutrient absorption and xenobiotic and fatty acid metabolism. Sifakas show convergent evolution at the level of the pathway, gene family, gene, and amino acid substitution with other folivores. Although sifakas have relatively generalized diets, the physiological challenges of habitual folivory likely led to strong selection.

Germline mutations are the raw material for natural selection, driving species evolution and the generation of earth's biodiversity. Without this driver of genetic diversity, life on earth would stagnate. Yet, it is a double-edged sword. An excess of mutations can have devastating effects on fitness and population viability. It is therefore one of the great challenges of molecular ecology to determine the rate and mechanisms by which these mutations accrue across the tree of life. Advances in high-throughput sequencing technologies are providing new opportunities for characterizing the rates and mutational spectra within species and populations thus informing essential evolutionary parameters such as the timing of speciation events, the intricacies of historical demography, and the degree to which lineages are subject to the burdens of mutational load. Here, we will focus on both the challenge and promise of whole-genome comparisons among parents and their offspring from known pedigrees for the detection of germline mutations as they arise in a single generation. The potential of these studies is high, but the field is still in its infancy and much uncertainty remains. Namely, the technical challenges are daunting given that pedigree-based genome comparisons are essentially searching for needles in a haystack given the very low signal to noise ratio. Despite the challenges, we predict that rapidly developing methods for whole-genome comparisons hold great promise for integrating empirically derived estimates of de novo mutation rates and mutation spectra across many molecular ecological applications.

Background Nucleotide excision repair is the primary DNA repair mechanism that removes bulky DNA adducts such as UV-induced pyrimidine dimers. Correspondingly, genome-wide mapping of nucleotide excision repair with eXcision Repair sequencing (XR-seq), provides comprehensive profiling of DNA damage repair. A number of XR-seq experiments at a variety of conditions for different damage types revealed heterogenous repair in the human genome. Although human repair profiles were extensively studied, how repair maps vary between primates is yet to be investigated. Here, we characterized the genome-wide UV-induced damage repair in gray mouse lemur, Microcebus murinus, in comparison to human. Results We derived fibroblast cell lines from mouse lemur, exposed them to UV irradiation, and analyzed the repair events genome-wide using the XR-seq protocol. Mouse lemur repair profiles were analyzed in comparison to the equivalent human fibroblast datasets. We found that overall UV sensitivity, repair efficiency, and transcription-coupled repair levels differ between the two primates. Despite this, comparative analysis of human and mouse lemur fibroblasts revealed that genome-wide repair profiles of the homologous regions are highly correlated, and this correlation is stronger for highly expressed genes. With the inclusion of an additional XR-seq sample derived from another human cell line in the analysis, we found that fibroblasts of the two primates repair UV-induced DNA lesions in a more similar pattern than two distinct human cell lines do. Conclusion Our results suggest that mouse lemurs and humans, and possibly primates in general, share a homologous repair mechanism as well as genomic variance distribution, albeit with their variable repair efficiency. This result also emphasizes the deep homologies of individual tissue types across the eukaryotic phylogeny.

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic explosion has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we take an integrative approach to investigate species diversity in two pairs of sister lineages that occur in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among them. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species,M. mittermeieriandM. lehilahytsarais not supported – a result that is particularly striking when using the genealogical discordance index (gdi). Non-sister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages, and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC and calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated concatenated likelihood estimates, however. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case.

Mutations are the raw material on which evolution acts, and knowledge of their frequency and genomic distribution is crucial for understanding how evolution operates at both long and short timescales. At present, the rate and spectrum of de novo mutations have been directly characterized in relatively few lineages. Our study provides the first direct mutation-rate estimate for a strepsirrhine (i.e., the lemurs and lorises), which comprises nearly half of the primate clade. Using high-coverage linked-read sequencing for a focal quartet of gray mouse lemurs (Microcebus murinus), we estimated the mutation rate to be among the highest calculated for a mammal at 1.52 × 10(–8) (95% credible interval: 1.28 × 10(−8)–1.78 × 10(−8)) mutations/site/generation. Further, we found an unexpectedly low count of paternal mutations, and only a modest overrepresentation of mutations at CpG sites. Despite the surprising nature of these results, we found both the rate and spectrum to be robust to the manipulation of a wide range of computational filtering criteria. We also sequenced a technical replicate to estimate a false-negative and false-positive rate for our data and show that any point estimate of a de novo mutation rate should be considered with a large degree of uncertainty. For validation, we conducted an independent analysis of context-dependent substitution types for gray mouse lemur and five additional primate species for which de novo mutation rates have also been estimated. These comparisons revealed general consistency of the mutation spectrum between the pedigree-based and the substitution-rate analyses for all species compared.

Molecular data have been used to date species divergences ever since they were described as documents of evolutionary history in the 1960s. Yet, an inadequate fossil record and discordance between gene trees and species trees are persistently problematic. We examine how, by accommodating gene tree discordance and by scaling branch lengths to absolute time using mutation rate and generation time, multispecies coalescent (MSC) methods can potentially overcome these challenges. We find that time estimates can differ – in some cases, substantially – depending on whether MSC methods or traditional phylogenetic methods that apply concatenation are used, and whether the tree is calibrated with pedigree-based mutation rates or with fossils. We discuss the advantages and shortcomings of both approaches and provide practical guidance for data analysis when using these methods.

Hibernation is an adaptive strategy some mammals use to survive highly seasonal or unpredictable environments. We present the first investigation on the transcriptomics of hibernation in a natural population of primate hibernators: Crossley's dwarf lemurs (Cheirogaleus crossleyi). Using capture-mark-recapture techniques to track the same animals over a period of 7 months in Madagascar, we used RNA-seq to compare gene expression profiles in white adipose tissue (WAT) during three distinct physiological states. We focus on pathway analysis to assess the biological significance of transcriptional changes in dwarf lemur WAT and, by comparing and contrasting what is known in other model hibernating species, contribute to a broader understanding of genomic contributions of hibernation across Mammalia. The hibernation signature is characterized by a suppression of lipid biosynthesis, pyruvate metabolism and mitochondrial-associated functions, and an accumulation of transcripts encoding ribosomal components and iron-storage proteins. The data support a key role of pyruvate dehydrogenase kinase isoenzyme 4 (PDK4) in regulating the shift in fuel economy during periods of severe food deprivation. This pattern of PDK4 holds true across representative hibernating species from disparate mammalian groups, suggesting that the genetic underpinnings of hibernation may be ancestral to mammals.

Madagascar is a biodiversity hotspot that is facing rapid rates of deforestation, habitat destruction and poverty. Urgent action is required to document the status of biodiversity to facilitate efficacious conservation plans. Within country, new generations of Malagasy scientists and conservationists are taking on leadership roles, although many lack access to modern genetic sequencing and are underrepresented in academic publications, when compared to international counterparts.With the recent advent of portable and affordable genetic technologies, it is now possible to tackle logistical considerations. Mobile genetics labs, with the capacity for in situ DNA extraction, amplification and sequencing, can produce scientifically reproducible data under field conditions, minimizing the time between sample collection and data analysis. Additionally, mobile labs offer powerful training opportunities for in-country scientists that enable local students and researchers to actively participate and contribute fully to the research enterprise, and that further empower these communities to contribute to the conservation dialog.Here, we show “proof of concept” by deploying a miniaturized thermal cycler alongside the Oxford Nanopore MinION DNA sequencer in Madagascar, including in the newly established Anjajavy Protected Area in northwestern Madagascar. We successfully extracted DNA from tissue samples collected using minimally-invasive techniques, amplified and sequenced a phylogenetically informative mitochondrial gene (cytochrome-b; cytb), and thereby confirmed the presence of Danfoss’ mouse lemur (M. danfossi) within the Anjajavy Reserve.To demonstrate the reproducibility of our methods, we successfully performed our established molecular and analytical pipeline at two additional locations in Madagascar, where we also conducted two-day workshops at local higher-education Institutions to demonstrate the process from tissue samples to DNA sequencing. Ultimately, we show that a mobile genetics lab can provide reliable and expeditious results, become a powerful educational tool, and allow scientists to conduct genetic analyses, potentially allowing for rapid interventions under emergency conditions in situ.

Madagascar’s biodiversity is notoriously threatened by deforestation and climate change. Many of these organisms are rare, cryptic, and severely threatened, making population-level sampling unrealistic. Such is the case with Madagascar’s dwarf lemurs (genusCheirogaleus), the only obligate hibernating primate. We here apply comparative genomic approaches to generate the first genome-wide estimates of genetic diversity within dwarf lemurs. We generate a reference genome for the fat-tailed dwarf lemur,Cheirogaleus medius, and use this resource to facilitate analyses of high-coverage (~30x) genome sequences for wild-caught individuals representing species:C. sp. cf. medius, C. major, C. crossleyiandC. sibreei. This study represents the largest contribution to date of novel genomic resources for Madagascar’s lemurs. We find concordant phylogenetic relationships among the four lineages ofCheirogaleusacross most of the genome, and yet detect a number of discordant genomic regions consistent with ancient admixture. We hypothesized that these regions could have resulted from adaptive introgression related to hibernation, indeed finding that genes associated with hibernation are present, though most significantly, that gene ontology categories relating to transcription are over-represented. We estimate levels of heterozygosity and find particularly low levels in an individual sampled from an isolated population ofC. mediusthat we refer to asC. sp. cf. medius. Results are consistent with a recent decline in effective population size, which is evident across species. Our study highlights the power of comparative genomic analysis for identifying species and populations of conservation concern, as well as for illuminating possible mechanisms of adaptive phenotypic evolution.

Kimura's neutral theory of molecular evolution has been essential to virtually every advance in evolutionary genetics, and by extension, is foundational to the field of conservation genetics. Conservation genetics utilizes the key concepts of neutral theory to identify species and populations at risk of losing evolutionary potential by detecting patterns of inbreeding depression and low effective population size. In turn, this information can inform the management of organisms and their habitat providing hope for the long-term preservation of both. We expand upon Avise's "inventorial" and "functional" categories of conservation genetics by proposing a third category that is linked to the coalescent and that we refer to as "process-driven." It is here that connections between Kimura's theory and conservation genetics are strongest. Process-driven conservation genetics can be especially applied to large genomic data sets to identify patterns of historical risk, such as population bottlenecks, and accordingly, yield informed intuitions for future outcomes. By examining inventorial, functional, and process-driven conservation genetics in sequence, we assess the progression from theory, to data collection and analysis, and ultimately, to the production of hypotheses that can inform conservation policies.

The rainforest biome of eastern Madagascar is renowned for its extraordinary biodiversity and restricted distribution ranges of many species, whereas the arid western region of the island is relatively species poor. We provide insight into the biogeography of western Madagascar by analyzing a multilocus phylogeographic dataset assembled for an amphibian, the widespread Malagasy bullfrog, Laliostoma labrosum. We find no cryptic species in L. labrosum (maximum 1.1% pairwise genetic distance between individuals in the 16S rRNA gene) attributable to considerable gene flow at the regional level as shown by genetic admixture in both mtDNA and three nuclear loci, especially in central Madagascar. Low breeding site fidelity, viewed as an adaptation to the unreliability of standing pools of freshwater in dry and seasonal environments, and a ubiquitous distribution within its range may underlie overall low genetic differentiation. Moreover, reductions in population size associated with periods of high aridity in western Madagascar may have purged DNA variation in this species. The mtDNA gene tree revealed seven major phylogroups within this species, five of which show mostly non-overlapping distributions. The nested positions of the northern and central mtDNA phylogroups imply a southwestern origin for all extant mtDNA lineages in L. labrosum. The current phylogeography of this species and paleo-distributions of major mtDNA lineages suggest five potential refugia in northern, western and southwestern Madagascar, likely the result of Pleistocene range fragmentation during drier and cooler climates. Lineage sorting in mtDNA and nuclear loci highlighted a main phylogeographic break between populations north and south of the Sambirano region, suggesting a role of the coastal Sambirano rainforest as a barrier to gene flow. Paleo-species distribution models and dispersal networks suggest that the persistence of some refugial populations was mainly determined by high population connectivity through space and time.

DNA quantity can be a hindrance in ecological and evolutionary research programmes due to a range of factors including endangered status of target organisms, available tissue type, and the impact of field conditions on preservation methods. A potential solution to low-quantity DNA lies in whole genome amplification (WGA) techniques that can substantially increase DNA yield. To date, few studies have rigorously examined sequence bias that might result from WGA and next-generation sequencing of nonmodel taxa. To address this knowledge deficit, we use multiple displacement amplification (MDA) and double-digest RAD sequencing on the grey mouse lemur (Microcebus murinus) to quantify bias in genome coverage and SNP calls when compared to raw genomic DNA (gDNA). We focus our efforts in providing baseline estimates of potential bias by following manufacturer's recommendations for starting DNA quantities (>100 ng). Our results are strongly suggestive that MDA enrichment does not introduce systematic bias to genome characterization. SNP calling between samples when genotyping both de-novo and with a reference genome are highly congruent (>98%) when specifying a minimum threshold of 20X stack depth to call genotypes. Relative genome coverage is also similar between MDA and gDNA, and allelic dropout is not observed. SNP concordance varies based on coverage threshold, with 95% concordance reached at ~12X coverage genotyping de-novo and ~7X coverage genotyping with the reference genome. These results suggest that MDA may be a suitable solution for next-generation molecular ecological studies when DNA quantity would otherwise be a limiting factor.

Aluelements are a highly successful family of primate-specific retrotransposons that have fundamentally shaped primate evolution, including the evolution of our own species.Alus play critical roles in the formation of neurological networks and the epigenetic regulation of biochemical processes throughout the central nervous system (CNS), and thus are hypothesized to have contributed to the origin of human cognition. Despite the benefits thatAlusprovide, deleteriousAluactivity is associated with a number of neurological and neurodegenerative disorders. In particular, neurological networks are potentially vulnerable to the epigenetic dysregulation ofAluelements operating across the suite of nuclear-encoded mitochondrial genes that are critical for both mitochondrial and CNS function. Here, we highlight the beneficial neurological aspects ofAluelements as well as their potential to cause disease by disrupting key cellular processes across the CNS. We identify at least 37 neurological and neurodegenerative disorders wherein deleteriousAluactivity has been implicated as a contributing factor for the manifestation of disease and, for many of these disorders, this activity is operating on genes that are essential for proper mitochondrial function. We conclude that the epigenetic dysregulation ofAluelements can ultimately disrupt mitochondrial homeostasis within the CNS. This mechanism is a plausible source for the incipient neuronal stress that is consistently observed across a spectrum of sporadic neurological and neurodegenerative disorders.List of AbbreviationsA-to-Iadenosine-to-inosineADAlzheimer’s DiseaseADARadenosine deaminase acting on RNAALSAmyotrophic Lateral SclerosisAMPAα-amino-3-hydroxy-5methyl-4-isoxazole propionateAPPamyloid precursor proteincircRNAscircular RNAsCNScentral nervous systemFLAMfree leftAlumonomerLINElong interspersed elementL1long interspersed element-1LTRlong-terminal repeatmRNAmessenger RNAPDParkinson’s Diseasepre-mRNAprecursor messenger RNASEDssuper-enhancer domainsSINEshort-interspersed elementTADstopologically associating domainsTOMMtranslocase of outer mitochondrial membrane

Background The de novo assembly of repeat-rich mammalian genomes using only high-throughput short read sequencing data typically results in highly fragmented genome assemblies that limit downstream applications. Here, we present an iterative approach to hybrid de novo genome assembly that incorporates datasets stemming from multiple genomic technologies and methods. We used this approach to improve the gray mouse lemur (Microcebus murinus) genome from early draft status to a near chromosome-scale assembly. Methods We used a combination of advanced genomic technologies to iteratively resolve conflicts and super-scaffold the M. murinus genome. Results We improved the M. murinus genome assembly to a scaffold N50 of 93.32 Mb. Whole genome alignments between our primary super-scaffolds and 23 human chromosomes revealed patterns that are congruent with historical comparative cytogenetic data, thus demonstrating the accuracy of our de novo scaffolding approach and allowing assignment of scaffolds to M. murinus chromosomes. Moreover, we utilized our independent datasets to discover and characterize sequences associated with centromeres across the mouse lemur genome. Quality assessment of the final assembly found 96% of mouse lemur canonical transcripts nearly complete, comparable to other published high-quality reference genome assemblies. Conclusions We describe a new assembly of the gray mouse lemur (Microcebus murinus) genome with chromosome-scale scaffolds produced using a hybrid bioinformatic and sequencing approach. The approach is cost effective and produces superior results based on metrics of contiguity and completeness. Our results show that emerging genomic technologies can be used in combination to characterize centromeres of non-model species and to produce accurate de novo chromosome-scale genome assemblies of complex mammalian genomes. Electronic supplementary material The online version of this article (doi:10.1186/s12915-017-0439-6) contains supplementary material, which is available to authorized users.

Microbiologists often evaluate microbial community dynamics by formulating functional hypotheses based on ecological processes. Indeed, many of the methods and terms currently used to describe animal microbiomes derive from ecology and evolutionary biology. As our understanding of the composition and functional dynamics of "the microbiome" grows, we increasingly refer to the host as an ecosystem within which microbial processes play out. Even so, an ecosystem service framework that extends to the context of the host has thus far been lacking. Here, we argue that ecosystem services are a useful framework with which to consider the value of microbes to their hosts. We discuss those "microbiome services" in the specific context of the mammalian gut, providing a context from which to develop new hypotheses and to evaluate microbial functions in future studies and novel systems.

Primates have long been a test case for the development of phylogenetic methods for divergence time estimation. Despite a large number of studies, however, the timing of origination of crown Primates relative to the K-Pg boundary and the timing of diversification of the main crown groups remain controversial. Here we analysed a dataset of 372 taxa (367 Primates and 5 outgroups, 61 thousand base pairs) that includes nine complete primate genomes (3.4 million base pairs). We systematically explore the effect of different interpretations of fossil calibrations and molecular clock models on primate divergence time estimates. We find that even small differences in the construction of fossil calibrations can have a noticeable impact on estimated divergence times, especially for the oldest nodes in the tree. Notably, choice of molecular rate model (auto-correlated or independently distributed rates) has an especially strong effect on estimated times, with the independent rates model producing considerably more ancient estimates for the deeper nodes in the phylogeny. We implement thermodynamic integration, combined with Gaussian quadrature, in the program MCMCTree, and use it to calculate Bayes factors for clock models. Bayesian model selection indicates that the auto-correlated rates model fits the primate data substantially better, and we conclude that time estimates under this model should be preferred. We show that for eight core nodes in the phylogeny, uncertainty in time estimates is close to the theoretical limit imposed by fossil uncertainties. Thus, these estimates are unlikely to be improved by collecting additional molecular sequence data. All analyses place the origin of Primates close to the K-Pg boundary, either in the Cretaceous or straddling the boundary into the Palaeogene.

In recent years, the study of the molecular processes involved in mammalian hibernation has shifted from investigating a few carefully selected candidate genes to large-scale analysis of differential gene expression. The availability of high-throughput data provides an unprecedented opportunity to ask whether phylogenetically distant species show similar mechanisms of genetic control, and how these relate to particular genes and pathways involved in the hibernation phenotype. In order to address these questions, we compare 11 datasets of differentially expressed (DE) genes from two ground squirrel species, one bat species, and the American black bear, as well as a list of genes extracted from the literature that previously have been correlated with the drastic physiological changes associated with hibernation. We identify several genes that are DE in different species, indicating either ancestral adaptations or evolutionary convergence. When we use a network approach to expand the original datasets of DE genes to large gene networks using available interactome data, a higher agreement between datasets is achieved. This indicates that the same key pathways are important for activating and maintaining the hibernation phenotype. Functional-term-enrichment analysis identifies several important metabolic and mitochondrial processes that are critical for hibernation, such as fatty acid beta-oxidation and mitochondrial transport. We do not detect any enrichment of positive selection signatures in the coding sequences of genes from the networks of hibernation-associated genes, supporting the hypothesis that the genetic processes shaping the hibernation phenotype are driven primarily by changes in gene regulation.

Vomeronasal receptor genes have frequently been invoked as integral to the establishment and maintenance of species boundaries among mammals due to the elaborate one-to-one correspondence between semiochemical signals and neuronal sensory inputs. Here, we report the most extensive sample of vomeronasal receptor class 1 (V1R) sequences ever generated for a diverse yet phylogenetically coherent group of mammals, the tooth-combed primates (suborder Strepsirrhini). Phylogenetic analysis confirms our intensive sampling from a single V1R subfamily, apparently unique to the strepsirrhine primates. We designate this subfamily as V1Rstrep. The subfamily retains extensive repertoires of gene copies that descend from an ancestral gene duplication that appears to have occurred prior to the diversification of all lemuriform primates excluding the basal genus Daubentonia (the aye-aye). We refer to the descendent clades as V1Rstrep-α and V1Rstrep-β. Comparison of the two clades reveals different amino acid compositions corresponding to the predicted ligand-binding site and thus potentially to altered functional profiles between the two. In agreement with previous studies of the mouse lemur (genus, Microcebus), the majority of V1Rstrep gene copies appear to be intact and under strong positive selection, particularly within transmembrane regions. Finally, despite the surprisingly high number of gene copies identified in this study, it is nonetheless probable that V1R diversity remains underestimated in these nonmodel primates and that complete characterization will be limited until high-coverage assembled genomes are available. peerReviewed

It is hypothesized that retrotransposons have played a fundamental role in primate evolution and that enhanced neurologic retrotransposon activity in humans may underlie the origin of higher cognitive function. As a potential consequence of this enhanced activity, it is likely that neurons are susceptible to deleterious retrotransposon pathways that can disrupt mitochondrial function. An example is observed in the TOMM40 gene, encoding a β-barrel protein critical for mitochondrial preprotein transport. Primate-specific Alu retrotransposons have repeatedly inserted into TOMM40 introns, and at least one variant associated with late-onset Alzheimer's disease originated from an Alu insertion event. We provide evidence of enriched Alu content in mitochondrial genes and postulate that Alus can disrupt mitochondrial populations in neurons, thereby setting the stage for progressive neurologic dysfunction. This Alu neurodegeneration hypothesis is compatible with decades of research and offers a plausible mechanism for the disruption of neuronal mitochondrial homeostasis, ultimately cascading into neurodegenerative disease.

— Resolving phylogenetic relationships within the economically important genus Coffea (Rubiaceae) has proven to be difficult due to low levels of plastid and nuclear ITS sequence divergence. The recent identification of a large number of conserved, single-copy, putatively orthologous genes (COSII) provides a unique opportunity to apply regions of the nuclear genome to phylogenetic studies of Euasterid species. We expand a previously published plastid data set of Coffea with the addition of three new plastid regions and a single COSII nuclear locus. Significant phylogenetic incongruence is identified between the plastid and nuclear COSII data sets, and a combined analysis is performed after removal of incongruent taxa. Phylogenetic results from plastid, nuclear, and combined plastid/nuclear data sets allow the novel identification of early diverging lineages in Coffea. Specifically, the data show that a group of predominately Lower-Guinea/Congolian Coffea species form a clade that is sister to the ...

Understanding the remarkably high species diversity and levels of endemism found among Madagascar’s flora and fauna has been the focus of many studies. One hypothesis that has received much attention proposes that Quaternary climate fluctuations spurred diversification. While spatial patterns of distribution and phylogenetic relationships can provide support for biogeographic predictions, temporal estimates of divergence are required to determine the fit of these geospatial patterns to climatic or biogeographic mechanisms. We use multilocus DNA sequence data to test whether divergence times among Malagasy iguanid lizards of the subfamily Oplurinae are compatible with a hypotheses of Pliocene‐Pleistocene diversification. We estimate the oplurine species tree and associated divergence times under a relaxed-clock model. In addition, we examine the phylogeographic structure and population divergence times within two sister species of Oplurus primarily distributed in the north-west and south-west of Madagascar (Oplurus cuvieriandOpluruscyclurus,respectively).Wefindthatdivergenceeventsamongoplurine lineages occurred in the Oligocene and Miocene and are thus far older and incompatible with the hypothesis that recent climate fluctuations are related to current species diversity. However,thetimingofintraspecificdivergencesandspatialpatternsofpopulationgenetic structure within O. cuvieri and O. cyclurus suggest a role for both intrinsic barriers and recentclimate fluctuationsat population-leveldivergences.Integratinginformationacross spatial and temporal scales allows us to identify and better understand the mechanisms generating patterns diversity.

The systematics and speciation literature is rich with discussion relating to the potential for gene tree/species tree discordance. Numerous mechanisms have been proposed to generate discordance, including differential selection, longbranch attraction, gene duplication, genetic introgression, and/or incomplete lineage sorting. For speciose clades in which divergence has occurred recently and rapidly, recovering the true species tree can be particularly problematic due to incomplete lineage sorting. Unfortunately, the availability of multilocus or ‘‘phylogenomic’’ data sets does not simply solve the problem, particularly when the data are analyzed with standard concatenation techniques. In our study, we conduct a phylogenetic study for a nearly complete species sample of the dwarf and mouse lemur clade, Cheirogaleidae. Mouse lemurs (genus, Microcebus) have been intensively studied over the past decade for reasons relating to their high level of cryptic species diversity, and although there has been emerging consensus regarding the evolutionary diversity contained within the genus, there is no agreement as to the inter-specific relationships within the group. We attempt to resolve cheirogaleid phylogeny, focusing especially on the mouse lemurs, by employing a large multilocus data set. We compare the results of Bayesian concordance methods with those of standard gene concatenation, finding that though concatenation yields the strongest results as measured by statistical support, these results are found to be highly misleading. By employing an approach where individual alleles are treated as operational taxonomic units, we show that phylogenetic results are substantially influenced by the selection of alleles in the concatenation process.

We investigate the genetic relationships between purported island species of Pteropus fruit bat (Megachiroptera) from the western Indian Ocean islands using mitochondrial DNA sequencing in order to infer the pattern of colonisation of this biogeographic region. Most significantly, our genetic data questions the current taxonomic assignment based on morphology of many of the island species and subspecies, suggesting instead that many of the western Indian Ocean islands harbour 'races' of P. giganteus from mainland India. Our results strongly argue against a single colonisation event from mainland Asia. Evidence is presented for three colonisation events; the first to the western-most extremity of their range (Comoros and Pemba Island), the second to Rodrigues Island; and a third giving rise to the remaining extant island taxa, the latter two events occurring relatively recently and rapidly.

The chiropteran family Myzopodidae is endemic to Madagascar and is characterized by several unique morphologies, such as sessile adhesive discs on the thumb and sole. A new species, Myzopoda schliemanni, was recently described from western Madagascar that is morphologically distinct and geographically disjunct from the eastern species, M. aurita, the only other member of this family. Geographic variation within Myzopoda has only recently been studied at the morphological level and has never been addressed at the genetic level. We used a combination of phylogenetic, coalescent, and population genetic analyses to characterize the speciation history of Myzopoda and to clarify current and former patterns of gene flow within and between Myzopoda. Mitochondrial DNA sequences were used to determine whether genetic data support the morphologically distinct species M. schliemanni, to infer the distribution of the common ancestor of extant Myzopoda, to estimate effective population sizes (Ne) and levels of migration between species, and to determine patterns of population structure within species. Phylogenetic and network analyses revealed the existence of 4 well-supported clades in Myzopoda, but could not resolve relationships among those clades. Divergent haplotypes within species may result from either recent gene flow between the 2 species or more likely from incomplete lineage sorting. Multiple coalescent-based methodologies produced concordant estimates of Ne for Myzopoda, but conflicting signals for migration between the species, probably reflecting differences in the underlying models used by the methods. We found significant genetic structure within M. aurita, but no correlation with geography. This pattern may result from recent gene flow facilitated by expansion of Ravenala stands, an important day-roost tree for Myzopoda, associated with anthropogenic deforestation and the opening up of new habitat for members of this genus.

Hibernation is a complex physiological response that some mammalian species employ to evade energetic demands. Previous work in mammalian hibernators suggests that hibernation is activated not by a set of genes unique to hibernators, but by differential expression of genes that are present in all mammals. This question of universal genetic mechanisms requires further investigation and can only be tested through additional investigations of phylogenetically dispersed species. To explore this question, we use RNA-Seq to investigate gene expression dynamics as they relate to the varying physiological states experienced throughout the year in a group of primate hibernators—Madagascar’s dwarf lemurs (genus Cheirogaleus). In a novel experimental approach, we use longitudinal sampling of biological tissues as a method for capturing gene expression profiles from the same individuals throughout their annual hibernation cycle. We identify 90 candidate genes that have variable expression patterns when comparing two active states (Active 1 and Active 2) with a torpor state. These include genes that are involved in metabolic pathways, feeding behavior, and circadian rhythms, as might be expected to correlate with seasonal physiological state changes. The identified genes appear to be critical for maintaining the health of an animal that undergoes prolonged periods of metabolic depression concurrent with the hibernation phenotype. By focusing on these differentially expressed genes in dwarf lemurs, we compare gene expression patterns in previously studied mammalian hibernators. Additionally, by employing evolutionary rate analysis, we find that hibernation-related genes do not evolve under positive selection in hibernating species relative to nonhibernators., This work was supported by grants from the Society of Integrative and Comparative Biology and Sigma Xi to S.L.F., a National Science Foundation Doctoral Research Improvement Grant (NSF-BCS 1455809 to A.D.Y. and S.L.F.), a grant from the Duke University Institute for Genome Science and Policy to A.D.Y., and funding from the Ministerio de Economía y Competitividad of the Spanish Government (BIO2009-08160 to J.L.V.C. and BFU2012-36820 to J.L.V.C. and M.M.A.); and Fundació ICREA to M.M.A.

GPI membrane anchors of cell surface glycoproteins have been shown to confer functional properties that are different from their transmembrane (TM)-anchored counterparts. For the human carcinoembryonic antigen (CEA) family, a subfamily of the immunoglobulin superfamily, conversion of the mode of membrane linkage from TM to GPI confers radical changes in function: from tumor suppression or neutrality toward inhibition of differentiation and anoikis and distortion of tissue architecture, thereby contributing to tumorigenesis. We show here that GPI anchorage in the CEA family evolved twice independently in primates, very likely from more primitive TM anchors, by different packages of mutations. Both mutational packages, one package found in many primates, including humans, and a second, novel package found only in the Cebidae radiation of New World monkeys, give rise to efficiently processed GPI-linked proteins. Both types of GPI anchors mediate inhibition of cell differentiation. The estimated rate of nonsynonymous mutations (Ka) in the anchor-determining domain for conversion from TM to GPI anchorage in the CEA family that were fixed during evolution in these primates is 7 times higher than the average Kain primates, indicating positive selection. These results suggest therefore that the functional changes mediated by CEA GPI anchors, including the inhibition of differentiation and anoikis, could be adaptive and advantageous.

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.

Background RNA-Seq has enabled high-throughput gene expression profiling to provide insight into the functional link between genotype and phenotype. Low quantities of starting RNA can be a severe hindrance for studies that aim to utilize RNA-Seq. To mitigate this bottleneck, whole transcriptome amplification (WTA) technologies have been developed to generate sufficient sequencing targets from minute amounts of RNA. Successful WTA requires accurate replication of transcript abundance without the loss or distortion of specific mRNAs. Here, we test the efficacy of NuGEN’s Ovation RNA-Seq V2 system, which uses linear isothermal amplification with a unique chimeric primer for amplification, using white adipose tissue from standard laboratory rats (Rattus norvegicus). Our goal was to investigate potential biological artifacts introduced through WTA approaches by establishing comparisons between matched raw and amplified RNA libraries derived from biological replicates. Results We found that 93% of expressed genes were identical between all unamplified versus matched amplified comparisons, also finding that gene density is similar across all comparisons. Our sequencing experiment and downstream bioinformatic analyses using the Tuxedo analysis pipeline resulted in the assembly of 25,543 high-quality transcripts. Libraries constructed from raw RNA and WTA samples averaged 15,298 and 15,253 expressed genes, respectively. Although significant differentially expressed genes (P

Implementation of the coalescent model in a Bayesian framework is an emerging strength in genetically based species delimitation studies. By providing an objective measure of species diagnosis, these methods represent a quantitative enhancement to the analysis of multilocus data, and complement more traditional methods based on phenotypic and ecological characteristics. Recognized as two species 20 years ago, mouse lemurs (genus Microcebus) now comprise more than 20 species, largely diagnosed from mtDNA sequence data. With each new species description, enthusiasm has been tempered with scientific scepticism. Here, we present a statistically justified and unbiased Bayesian approach towards mouse lemur species delimitation. We perform validation tests using multilocus sequence data and two methodologies: (i) reverse-jump Markov chain Monte Carlo sampling to assess the likelihood of different models defined a priori by a guide tree, and (ii) a Bayes factor delimitation test that compares different species-tree models without a guide tree. We assess the sensitivity of these methods using randomized individual assignments, which has been used in bpp studies, but not with Bayes factor delimitation tests. Our results validate previously diagnosed taxa, as well as new species hypotheses, resulting in support for three new mouse lemur species. As the challenge of multiple researchers using differing criteria to describe diversity is not unique to Microcebus, the methods used here have significant potential for clarifying diversity in other taxonomic groups. We echo previous studies in advocating that multiple lines of evidence, including use of the coalescent model, should be trusted to delimit new species.

Divergence time and substitution rate are seriously confounded in phylogenetic analysis, making it difficult to estimate divergence times when the molecular clock (rate constancy among lineages) is violated. This problem can be alleviated to some extent by analyzing multiple gene loci simultaneously and by using multiple calibration points. While different genes may have different patterns of evolutionary rate change, they share the same divergence times. Indeed, the fact that each gene may violate the molecular clock differently leads to the advantage of simultaneous analysis of multiple loci. Multiple calibration points provide the means for characterizing the local evolutionary rates on the phylogeny. In this paper, we extend previous likelihood models of local molecular clock for estimating species divergence times to accommodate multiple calibration points and multiple genes. Heterogeneity among different genes in evolutionary rate and in substitution process is accounted for by the models. We apply the likelihood models to analyze two mitochondrial protein-coding genes, cytochrome oxidase II and cytochrome b, to estimate divergence times of Malagasy mouse lemurs and related outgroups. The likelihood method is compared with the Bayes method of Thorne et al. (1998, Mol. Biol. Evol. 15:1647-1657), which uses a probabilistic model to describe the change in evolutionary rate over time and uses the Markov chain Monte Carlo procedure to derive the posterior distribution of rates and times. Our likelihood implementation has the drawbacks of failing to accommodate uncertainties in fossil calibrations and of requiring the researcher to classify branches on the tree into different rate groups. Both problems are avoided in the Bayes method. Despite the differences in the two methods, however, data partitions and model assumptions had the greatest impact on date estimation. The three codon positions have very different substitution rates and evolutionary dynamics, and assumptions in the substitution model affect date estimation in both likelihood and Bayes analyses. The results demonstrate that the separate analysis is unreliable, with dates variable among codon positions and between methods, and that the combined analysis is much more reliable. When the three codon positions were analyzed simultaneously under the most realistic models using all available calibration information, the two methods produced similar results. The divergence of the mouse lemurs is dated to be around 7-10 million years ago, indicating a surprisingly early species radiation for such a morphologically uniform group of primates.

The phylogeny of Rosoideae was investigated using 44 species. Here we report new sequence data from the chloroplast trnL/F region as well as an increased sample of species. The analysis of these new data, along with previously used data from the nuclear ribosomal internal transcribed spacers (ITS), significantly increased resolution as well as confidence for Rosoideae phylogeny. Using both Bayesian inference and parsimony methods, we conducted analyses on the data sets separately and in combination. The resulting phylogenies are congruent with all well‐supported clades of Rosoideae found in previous analyses of ITS or rbcL data. The support for these and other clades is improved, and we consider several clades to be supported well enough to be named. The following clades are given phylogenetic definitions: Sanguisorbeae and its subclades Agrimoniinae and Sanguisorbinae, Potentilleae and its subclades Fragariinae and Potentilla, Roperculina (Rosa + Sanguisorbeae + Potentilleae), and Sanpotina (Sanguisorbea...

Olfaction plays a critical role in both survival of the individual and in the propagation of species. Studies from across the mammalian clade have found a remarkable correlation between organismal lifestyle and molecular evolutionary properties of receptor genes in both the main olfactory system (MOS) and the vomeronasal system (VNS). When a large proportion of intact (and putatively functional) copies is observed, the inference is made that a particular mode of chemoreception is critical for an organism's fit to its environment and is thus under strong positive selection. Conversely, when the receptors in question show a disproportionately large number of pseudogene copies, this contraction is interpreted as evidence of relaxed selection potentially leading to gene family extinction. Notably, it appears that a risk factor for gene family extinction is a high rate of nonsynonymous substitution. A survey of intact versus pseudogene copies among primate vomeronasal receptor Class one genes (V1Rs) appears to substantiate this hypothesis. Molecular evolutionary complexities in the V1R gene family combine rapid rates of gene duplication, gene conversion, lineage-specific expansions, deletions, and/or pseudogenization. An intricate mix of phylogenetic footprints and current adaptive landscapes have left their mark on primate V1Rs suggesting that the primate clade offers an ideal model system for exploring the molecular evolutionary and functional properties of the VNS of mammals. Primate V1Rs tell a story of ancestral function and divergent selection as species have moved into ever diversifying adaptive regimes. The sensitivity to functional collapse in these genes, consequent to their precariously high rates of nonsynonymous substitution, confer a remarkable capacity to reveal the lifestyles of the genomes that they presently occupy as well as those of their ancestors.

High-throughput sequencing platforms are generating massive amounts of genomic data from nonmodel species, and these data sets are valuable resources that can be mined to advance a number of research areas. An example is the growing amount of transcriptome data that allow for examination of gene expression in nonmodel species. Here, we show how publicly available transcriptome data from nonmodel primates can be used to design novel research focused on immunogenomics. We mined transcriptome data from the world's most endangered group of primates, the lemurs of Madagascar, for sequences corresponding to immunoglobulins. Our results confirmed homology between strepsirrhine and haplorrhine primate immunoglobulins and allowed for high-throughput sequencing of expressed antibodies (Ig-seq) in Coquerel's sifaka (Propithecus coquereli). Using both Pacific Biosciences RS and Ion Torrent PGM sequencing, we performed Ig-seq on two individuals of Coquerel's sifaka. We generated over 150 000 sequences of expressed antibodies, allowing for molecular characterization of the antigen-binding region. Our analyses suggest that similar VDJ expression patterns exist across all primates, with sequences closely related to the human VH 3 immunoglobulin family being heavily represented in sifaka antibodies. Moreover, the antigen-binding region of sifaka antibodies exhibited similar amino acid variation with respect to haplorrhine primates. Our study represents the first attempt to characterize sequence diversity of the expressed antibody repertoire in a species of lemur. We anticipate that methods similar to ours will provide the framework for investigating the adaptive immune response in wild populations of other nonmodel organisms and can be used to advance the burgeoning field of eco-immunology.

Background Foamy viruses (FVs) are a unique subfamily of retroviruses that are widely distributed in mammals. Owing to the availability of sequences from diverse mammals coupled with their pattern of codivergence with their hosts, FVs have one of the best-understood viral evolutionary histories ever documented, estimated to have an ancient origin. Nonetheless, our knowledge of some parts of FV evolution, notably that of prosimian and afrotherian FVs, is far from complete due to the lack of sequence data. Results Here, we report the complete genome of the first extant prosimian FV (PSFV) isolated from a lorisiforme galago (PSFVgal), and a novel partial endogenous viral element with high sequence similarity to FVs, present in the afrotherian Cape golden mole genome (ChrEFV). We also further characterize a previously discovered endogenous PSFV present in the aye-aye genome (PSFVaye). Using phylogenetic methods and available FV sequence data, we show a deep divergence and stable co-evolution of FVs in eutherian mammals over 100 million years. Nonetheless, we found that the evolutionary histories of bat, aye-aye, and New World monkey FVs conflict with the evolutionary histories of their hosts. By combining sequence analysis and biogeographical knowledge, we propose explanations for these mismatches in FV-host evolutionary history. Conclusion Our discovery of ChrEFV has expanded the FV host range to cover the whole eutherian clade, and our evolutionary analyses suggest a stable mammalian FV-host co-speciation pattern which extends as deep as the exafroplacentalian basal diversification. Nonetheless, two possible cases of host switching were observed. One was among New World monkey FVs, and the other involves PSFVaye and a bat FV which may involve cross-species transmission at the level of mammalian orders. Our results highlight the value of integrating multiple sources of information to elucidate the evolutionary history of viruses, including continental and geographical histories, ancestral host locations, in addition to the natural history of host and virus. Electronic supplementary material The online version of this article (doi:10.1186/1742-4690-11-61) contains supplementary material, which is available to authorized users.

DNA sequences of the complete cytochrome b gene are shown to contain robust phylogenetic signal for the strepsirrhine primates (i.e., lemurs and lorises). The phylogeny derived from these data conforms to other molecular studies of strepsirrhine relationships despite the fact that uncorrected nucleotide distances are high for nearly all intrastrepsirrhine comparisons, with most in the 15%-20% range. Cytochrome b sequences support the hypothesis that Malagasy lemuriforms and Afro-Asian lorisiforms each comprise clades that share a sister-group relationship. A study (Adkins and Honeycutt 1994) of the cytochrome c oxidase subunit II (COII) gene placed one Malagasy primate (Daubentonia) at the base of the strepsirrhine clade, thereby suggesting a diphyletic Lemuriformes. The reanalysis of COII third-position transversions, either alone or in combination with cytochrome b third-position transversions, however, yields a tree that is congruent with phylogenetic hypotheses derived from cytochrome b and other genetic data sets.

We report the development of 18, single-copy, anonymous nuclear loci from the Malagasy plated lizard Zonosaurus madagascariensis. More than 140 clones from a genomic library were examined and 38 potential loci tested across both closely and distantly related lizards. Of the 18 loci reported here, more than half (10) work in closely related zonosaurines although only one successfully amplified a homologous fragment in the distantly related iguanid (Oplurus). Sequences of these loci revealed a high frequency of single nucleotide polymorphisms, supporting previous reports of high levels of intraspecific variation in lizards.

The field of phylogeography continues to grow in terms of power and accessibility. Initially uniting population genetics and phylogenetics, it now spans disciplines as diverse as geology, statistics, climatology, ecology, physiology, and bioinformatics to name a few. One major and recent integration driving the field forward is between “statistical phylogeography” and Geographic Information Systems (GIS) (Knowles, 2009). Merging genetic and geospatial data, and their associated methodological toolkits, is helping to bring explicit hypothesis testing to the field of phylogeography. Hypotheses derived from one approach can be reciprocally tested with data derived from the other field and the synthesis of these data can help place demographic events in an historical and spatial context, guide genetic sampling, and point to areas for further investigation. Here, we present three practical examples of empirical analysis that integrate statistical genetic and GIS tools to construct and test phylogeographic hypotheses. Insights into the evolutionary mechanisms underlying recent divergences can benefit from simultaneously considering diverse types of information to iteratively test and reformulate hypotheses. Our goal is to provide the reader with an introduction to the variety of available tools and their potential application to typical questions in phylogeography with the hope that integrative methods will be more broadly and commonly applied to other biological systems and data sets.

a b s t r a c t Approximately 20 years ago, Avise and colleagues proposed the integration of phylogenetics and popu- lation genetics for investigating the connection between micro- and macroevolutionary phenomena. The new field was termed phylogeography. Since the naming of the field, the statistical rigor of phyloge- ography has increased, in large part due to concurrent advances in coalescent theory which enabled model-based parameter estimation and hypothesis testing. The next phase will involve phylogeography increasingly becoming the integrative and comparative multi-taxon endeavor that it was originally con- ceived to be. This exciting convergence will likely involve combining spatially-explicit multiple taxon coalescent models, genomic studies of natural selection, ecological niche modeling, studies of ecological speciation, community assembly and functional trait evolution. This ambitious synthesis will allow us to determine the causal links between geography, climate change, ecological interactions and the evolution and composition of taxa across whole communities and assemblages. Although such integration presents analytical and computational challenges that will only be intensified by the growth of genomic data in non-model taxa, the rapid development of ''likelihood-free" approximate Bayesian methods should per- mit parameter estimation and hypotheses testing using complex evolutionary demographic models and genomic phylogeographic data. We first review the conceptual beginnings of phylogeography and its accomplishments and then illustrate how it evolved into a statistically rigorous enterprise with the con- current rise of coalescent theory. Subsequently, we discuss ways in which model-based phylogeography can interface with various subfields to become one of the most integrative fields in all of ecology and evo- lutionary biology.

The endemic plated lizards (Gerrhosauridae) of Madagascar are one of the most diverse groups of lizards on the island (19 species) and are found in all ecoregions. On an island that presents so many interesting biological questions, plated lizards are an ideal group for examining patterns of diversification due to their high (but tractable) diversity and wide distribution. To resolve the taxonomy and evolutionary relationships of Malagasy plated lizards, and to explore their biogeographic history, we have generated a molecular phylogeny based on >2 kb of both mitochondrial (cob, 1142 bp) and nuclear (NT-3, 419 bp; c-mos, 542 bp) DNA sequences. Phylogenetic relationships and divergence-time estimates indicate Zonosaurinae are monophyletic, arising from a single colonization event, likely from Africa to Madagascar in the Paleogene (Paleocene/early Eocene), with subsequent radiation. Furthermore, our results reveal widespread mischaracterization of species delimitation and assignment based on phenotypic characteristics. Paraphyly of a number of zonosaurine species was strongly supported by our dataset, revealing not only instances of likely cryptic species (phenotypic conservatism) but also phenotypic/genotypic discordance in the delimitation of species (phenotypic variability). These results clearly demonstrate the complex history of biotic diversification on Madagascar and provide novel insight into biogeographic patterns on the island.

Lemurs and the other strepsirrhine primates are of great interest to the primate genomics community due to their phylogenetic placement as the sister lineage to all other primates. Previous attempts to resolve the phylogeny of lemurs employed limited mitochondrial or small nuclear data sets, with many relationships poorly supported or entirely unresolved. We used genomic resources to develop 11 novel markers from nine chromosomes, representing ∼9 kb of nuclear sequence data. In combination with previously published nuclear and mitochondrial loci, this yields a data set of more than 16 kb and adds ∼275 kb of DNA sequence to current databases. Our phylogenetic analyses confirm hypotheses of lemuriform monophyly and provide robust resolution of the phylogenetic relationships among the five lemuriform families. We verify that the genus Daubentonia is the sister lineage to all other lemurs. The Cheirogaleidae and Lepilemuridae are sister taxa and together form the sister lineage to the Indriidae; this clade is the sister lineage to the Lemuridae. Divergence time estimates indicate that lemurs are an ancient group, with their initial diversification occurring around the Cretaceous-Tertiary boundary. Given the power of this data set to resolve branches in a notoriously problematic area of primate phylogeny, we anticipate that our phylogenomic toolkit will be of value to other studies of primate phylogeny and diversification. Moreover, the methods applied will be broadly applicable to other taxonomic groups where phylogenetic relationships have been notoriously difficult to resolve. [Supplemental material is available online at www.genome.org. The sequence data from this study have been submitted to GenBank under accession nos. EU057196–EU057514 and EU342218–EU342345.]

Mouse lemurs (genus Microcebus) are nocturnal primates endemic to the island of Madagascar. Until recently, they were classified as two species, one from eastern and one from western Madagascar. Previously published analyses of morphometric and mitochondrial DNA data show strong support for the recognition of more than eight species, however. Here, we test the eight-species hypothesis with DNA data derived from four independently segregating nuclear loci. We find many areas of congruence between the mitochondrial and nuclear data, but incomplete lineage sorting and low mutation rates limit the phylogenetic resolution of the nuclear data. Even so, the nuclear loci unanimously find evidence for three deeply diverged lineages within the mouse lemur radiation: one that is congruent with the mtDNA "southern clade", another that is congruent with the mtDNA "northern clade", and one monospecific branch comprised of the species Microcebus ravelobensis. The latter result in particular emphasizes the need for careful biological study of this species.

The lemurs of Madagascar are a unique radiation of primates that show an extraordinary diversity of lifestyles, morphologies and behaviours. However, very little is known about the relative antiquity of lemuriform clades due to the lack of terrestrial fossils for the Tertiary of Madagascar. Here, we employ a Bayesian method to estimate divergence dates within the lemuriform radiation using several unlinked gene loci and multiple fossil calibrations outside the lemuriform clade. Two mitochondrial genes (cytochrome oxidase II and cytochrome b ), two nuclear introns (transthyretin intron 1 and von Willebrand factor gene intron 11) and one nuclear exon (interphotoreceptor retinoid binding protein, exon 1) are used in separate and combined analyses. The genes differ in taxon sampling and evolutionary characteristics but produce congruent date estimates. Credibility intervals narrow considerably in combined analyses relative to separate analyses due to the increased amount of data. We also test the relative effects of multiple vs. single calibration points, finding that, when only single calibration points are employed, divergence dates are systematically underestimated. For the mitochondrial DNA data set, we investigate the effects of sampling density within the mouse lemur radiation (genus Microcebus ). When only two representative species are included, estimated dates throughout the phylogeny are more recent than with the complete-species sample, with basal nodes less affected than recent nodes. The difference appears to be due to the manner in which priors on node ages are constructed in the two analyses. In nearly all analyses, the age of the lemuriform clade is estimated to be approximately 62–65 Ma, with initial radiation of mouse lemurs and true lemurs (genus Eulemur ) occurring approximately 8–12 Ma. The antiquity of the mouse lemur radiation is surprising given the near uniform morphology among species. Moreover, the observation that mouse lemurs and true lemurs are of similar ages suggests discrepancies in rates of morphological, behavioural and physiological evolution in the two clades, particularly with regard to characteristics of sexual signalling. These differences appear to correlate with the nocturnal vs. diurnal lifestyles, respectively, of these two primate groups.

A chemokine receptor from the seven-transmembrane-domain G-protein-coupled receptor superfamily is an essential coreceptor for the cellular entry of human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus (SIV) strains. To investigate nonhuman primate CC-chemokine receptor 5 (CCR5) homologue structure and function, we amplified CCR5 DNA sequences from peripheral blood cells obtained from 24 representative species and subspecies of the primate suborders Prosimii (family Lemuridae) and Anthropoidea (families Cebidae, Callitrichidae, Cercopithecidae, Hylobatidae, and Pongidae) by PCR with primers flanking the coding region of the gene. Full-length CCR5 was inserted into pCDNA3.1, and multiple clones were sequenced to permit discrimination of both alleles. Compared to the human CCR5 sequence, the CCR5 sequences of the Lemuridae, Cebidae, and Cercopithecidae shared 87, 91 to 92, and 96 to 99% amino acid sequence homology, respectively. Amino acid substitutions tended to cluster in the amino and carboxy termini, the first transmembrane domain, and the second extracellular loop, with a pattern of species-specific changes that characterized CCR5 homologues from primates within a given family. At variance with humans, all primate species examined from the suborder Anthropoidea had amino acid substitutions at positions 13 (N to D) and 129 (V to I); the former change is critical for CD4-independent binding of SIV to CCR5. Within the Cebidae, Cercopithecidae, and Pongidae (including humans), CCR5 nucleotide similarities were 95.2 to 97.4, 98.0 to 99.5, and 98.3 to 99.3%, respectively. Despite this low genetic diversity, the phylogeny of the selected primate CCR5 homologue sequences agrees with present primate systematics, apart from some intermingling of species of the Cebidae and Cercopithecidae. Constructed HOS.CD4 cell lines expressing the entire CCR5 homologue protein from each of the Anthropoidea species and subspecies were tested for their ability to support HIV-1 and SIV entry and membrane fusion. Other than that of Cercopithecus pygerythrus , all CCR5 homologues tested were able to support both SIV and HIV-1 entry. Our results suggest that the shared structure and function of primate CCR5 homologue proteins would not impede the movement of primate immunodeficiency viruses between species.