Search

Hibernation, a metabolic strategy, allows individuals to reduce energetic demands in times of energetic deficits. Hibernation is pervasive in nature, occurring in all major mammalian lineages and geographical regions; however, its expression is variable across species, populations, and individuals, suggesting that trade-offs are at play. Whereas hibernation reduces energy expenditure, energetically expensive arousals may impose physiological burdens. The torpor optimization hypothesis posits that hibernation should be expressed according to energy availability. The greater the energy surplus, the lower the hibernation output. The thrifty female hypothesis, a variation of the torpor optimization hypothesis, states that females should conserve more energy because of their more substantial reproductive costs. Contrarily, if hibernation's benefits offset its costs, hibernation may be maximized rather than optimized (e.g., hibernators with greater fat reserves could afford to hibernate longer). We assessed torpor expression in captive dwarf lemurs, primates that are obligate, seasonal, and tropical hibernators. Across 4.5 mo in winter, we subjected eight individuals at the Duke Lemur Center to conditions conducive to hibernation, recorded estimates of skin temperature hourly (a proxy for torpor), and determined body mass and tail fat reserves bimonthly. Across and between consecutive weigh-ins, heavier dwarf lemurs spent less time in torpor and lost more body mass. At equivalent body mass, females spent more time torpid and better conserved energy than did males. Although preliminary, our results support the torpor optimization and thrifty female hypotheses, suggesting that individuals optimize rather than maximize torpor according to body mass. These patterns are consistent with hibernation phenology in Madagascar, where dwarf lemurs hibernate longer in more seasonal habitats.

Papillomaviruses (PVs) are host-species-specific and tissue-specific viruses that infect a diverse array of vertebrate hosts, including humans and non-human primates, with varying pathogenic outcomes. Although primate PVs have been studied extensively, no complete genome sequences of PVs from lemurs have been determined to date. Saliva samples from three critically endangered, captive black-and-white ruffed lemurs (Varecia variegata variegata) at the Duke Lemur Center (USA) were analyzed, using high-throughput sequencing, for the presence of oral papillomaviruses. We identified three PVs from two individuals, one of which had a coinfection with two different PVs. Two of the three PVs share 99.6% nucleotide sequence identity, and we have named these isolates "Varecia variegata papillomavirus 1" (VavPV1). The third PV shares ~63% nucleotide sequence identity with VavPV1, and thus, we have named it "Varecia variegata papillomavirus 2" (VavPV2). Based on their E1 + E2 + L1 protein sequence phylogeny, the VavPVs form a distinct clade. This clade likely represents a novel genus, with VavPV1 and VavPV2 belonging to two distinct species. Our findings represent the first complete genome sequences of PVs found in lemuriform primates, with their presence suggesting the potential existence of diverse PVs across the over 100 species of lemurs.

Microsatellites have been a workhorse of evolutionary genetic studies for decades and are still commonly in use for estimating signatures of genetic diversity at the population and species level across a multitude of taxa. Yet, the very high mutation rate of these loci is a double-edged sword, conferring great sensitivity at shallow levels of analysis (e.g. paternity analysis) but yielding considerable uncertainty for deeper evolutionary comparisons. For the present study, we used reduced representation genome-wide data (restriction site-associated DNA sequencing (RADseq)) to test for patterns of interspecific hybridization previously characterized using microsatellite data in a contact zone between two closely related mouse lemur species in Madagascar ( Microcebus murinus and Microcebus griseorufus ). We revisit this system by examining populations in, near, and far from the contact zone, including many of the same individuals that had previously been identified as hybrids with microsatellite data. Surprisingly, we find no evidence for admixed nuclear ancestry. Instead, re-analyses of microsatellite data and simulations suggest that previously inferred hybrids were false positives and that the program N ew H ybrids can be particularly sensitive to erroneously inferring hybrid ancestry. Combined with results from coalescent-based analyses and evidence for local syntopic co-occurrence, we conclude that the two mouse lemur species are in fact completely reproductively isolated, thus providing a new understanding of the evolutionary rate whereby reproductive isolation can be achieved in a primate.

The gut microbiome can mediate host metabolism, including facilitating energy-saving strategies like hibernation. The dwarf lemurs of Madagascar (Cheirogaleus spp.) are the only obligate hibernators among primates. They also hibernate in the subtropics, and unlike temperate hibernators, fatten by converting fruit sugars to lipid deposits, torpor at relatively warm temperatures, and forage for a generalized diet after emergence. Despite these ecological differences, we might expect hibernation to shape the gut microbiome in similar ways across mammals. We, therefore, compare gut microbiome profiles, determined by amplicon sequencing of rectal swabs, in wild furry-eared dwarf lemurs (C. crossleyi) during fattening, hibernation, and after emergence. The dwarf lemurs exhibited reduced gut microbial diversity during fattening, intermediate diversity and increased community homogenization during hibernation, and greatest diversity after emergence. The Mycoplasma genus was enriched during fattening, whereas the Aerococcaceae and Actinomycetaceae families, and not Akkermansia, bloomed during hibernation. As expected, the dwarf lemurs showed seasonal reconfigurations of the gut microbiome; however, the patterns of microbial diversity diverged from temperate hibernators, and better resembled the shifts associated with dietary fruits and sugars in primates and model organisms. Our results thus highlight the potential for dwarf lemurs to probe microbiome-mediated metabolism in primates under contrasting conditions.

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

In the past decade, several studies have estimated the human per-generation germline mutation rate using large pedigrees. More recently, estimates for various nonhuman species have been published. However, methodological differences among studies in detecting germline mutations and estimating mutation rates make direct comparisons difficult. Here, we describe the many different steps involved in estimating pedigree-based mutation rates, including sampling, sequencing, mapping, variant calling, filtering, and appropriately accounting for false-positive and false-negative rates. For each step, we review the different methods and parameter choices that have been used in the recent literature. Additionally, we present the results from a ‘Mutationathon,’ a competition organized among five research labs to compare germline mutation rate estimates for a single pedigree of rhesus macaques. We report almost a twofold variation in the final estimated rate among groups using different post-alignment processing, calling, and filtering criteria, and provide details into the sources of variation across studies. Though the difference among estimates is not statistically significant, this discrepancy emphasizes the need for standardized methods in mutation rate estimations and the difficulty in comparing rates from different studies. Finally, this work aims to provide guidelines for computational and statistical benchmarks for future studies interested in identifying germline mutations from pedigrees.

November 2020 marked 2 y since the launch of the Earth BioGenome Project (EBP), which aims to sequence all known eukaryotic species in a 10-y timeframe. Since then, significant progress has been made across all aspects of the EBP roadmap, as outlined in the 2018 article describing the project's goals, strategies, and challenges (1). The launch phase has ended and the clock has started on reaching the EBP's major milestones. This Special Feature explores the many facets of the EBP, including a review of progress, a description of major scientific goals, exemplar projects, ethical legal and social issues, and applications of biodiversity genomics. In this Introduction, we summarize the current status of the EBP, held virtually October 5 to 9, 2020, including recent updates through February 2021. References to the nine Perspective articles included in this Special Feature are cited to guide the reader toward deeper understanding of the goals and challenges facing the EBP.

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

Habitat fragmentation is one of the major types of anthropogenic change, though fragmented landscapes predate human intervention. At present, the Central Highlands of Madagascar are covered by extensive grasslands interspersed with small discrete forest patches of unknown antiquity. Ankafobe, an actively protected site, comprises two such fragments of 12 and 30 ha, respectively, known to harbor three lemur species and other endemic wildlife. At this location, we conducted a survey of resident Goodman's mouse lemurs, Microcebus lehilahytsara, to determine baseline behavioral and ecological conditions for this isolated population. By studying primates in forest fragments, investigators can characterize the effects of shrinking habitats and decreasing connectivity on species diversity and survival, thus providing a glimpse into the potential resilience of species in the face of anthropogenic disturbance. Investigating the behavioral ecology of Goodman's mouse lemurs across their geographic range could help us understand their metabolic and ecological flexibility and predict species long-term survival prospects. We conducted night transect walks, using capture techniques and telemetry, to track eight radio-collared individuals. Preliminary density estimates based on a limited number of sightings (n = 18) were 2.19 ind/ha, and home range assessments ranged between 0.22 and 3.67 ha. Mouse lemurs traveled an average of 425 m nightly during the 5-h tracking periods and primarily fed on fruits of the mistletoe Bakerella clavata. The finding that Goodman's mouse lemurs apparently thrive in the seasonally cold and arid forest fragments in the Central Highlands indicates that they may be among the most tolerant and adaptable lemur species in Madagascar. These results point towards an exciting research program that focuses on ecological tolerance as a mechanism for long-term species survival.

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.

Significance Based on “subfossil” skeletal remains it is known that multiple now-extinct giant lemur (primate) species with estimated body masses of up to ∼160 kg survived on Madagascar into the past millennium. In this study, we used ancient DNA methods to sequence the nuclear genome of one of these megafaunal lemurs, Megaladapis edwardsi (∼85 kg). With the power of the nuclear genome, we robustly resolved the phylogenetic relationship between Megaladapis and other lemurs, which had been a lingering uncertainty. We also identified multiple signatures of past positive natural selection across the M. edwardsi genome that support reconstructions of this taxon as a large-bodied, specialized folivore., No endemic Madagascar animal with body mass >10 kg survived a relatively recent wave of extinction on the island. From morphological and isotopic analyses of skeletal “subfossil” remains we can reconstruct some of the biology and behavioral ecology of giant lemurs (primates; up to ∼160 kg) and other extraordinary Malagasy megafauna that survived into the past millennium. Yet, much about the evolutionary biology of these now-extinct species remains unknown, along with persistent phylogenetic uncertainty in some cases. Thankfully, despite the challenges of DNA preservation in tropical and subtropical environments, technical advances have enabled the recovery of ancient DNA from some Malagasy subfossil specimens. Here, we present a nuclear genome sequence (∼2× coverage) for one of the largest extinct lemurs, the koala lemur Megaladapis edwardsi (∼85 kg). To support the testing of key phylogenetic and evolutionary hypotheses, we also generated high-coverage nuclear genomes for two extant lemurs, Eulemur rufifrons and Lepilemur mustelinus, and we aligned these sequences with previously published genomes for three other extant lemurs and 47 nonlemur vertebrates. Our phylogenetic results confirm that Megaladapis is most closely related to the extant Lemuridae (typified in our analysis by E. rufifrons) to the exclusion of L. mustelinus, which contradicts morphology-based phylogenies. Our evolutionary analyses identified significant convergent evolution between M. edwardsi and an extant folivore (a colobine monkey) and an herbivore (horse) in genes encoding proteins that function in plant toxin biodegradation and nutrient absorption. These results suggest that koala lemurs were highly adapted to a leaf-based diet, which may also explain their convergent craniodental morphology with the small-bodied folivore Lepilemur.

International audience; Delimitation of cryptic species is increasingly based on genetic analyses but the integration of distributional, morphological, behavioral, and ecological data offers unique complementary insights into species diversification. We surveyed communities of nocturnal mouse lemurs (Microcebus spp.) in five different sites of northeastern Madagascar, measuring a variety of morphological parameters and assessing reproductive states for 123 individuals belonging to five different lineages. We documented two different non‐sister lineages occurring in sympatry in two areas. In both cases, sympatric species pairs consisted of a locally restricted (M. macarthurii or M. sp. #3) and a more widespread lineage (M. mittermeieri or M. lehilahytsara). Estimated Extents of Occurrence (EOO) of these lineages differed remarkably with 560 and 1,500 km2 versus 9,250 and 50,700 km2, respectively. Morphometric analyses distinguished unambiguously between sympatric species and detected more subtle but significant differences among sister lineages. Tail length and body size were most informative in this regard. Reproductive schedules were highly variable among lineages, most likely impacted by phylogenetic relatedness and environmental variables. While sympatric species pairs differed in their reproductive timing (M. sp. #3/M. lehilahytsara and M. macarthurii/M. mittermeieri), warmer lowland rainforests were associated with a less seasonal reproductive schedule for M. mittermeieri and M. lehilahytsara compared with populations occurring in montane forests. Distributional, morphological, and ecological data gathered in this study support the results of genomic species delimitation analyses conducted in a companion study, which identified one lineage, M. sp. #3, as meriting formal description as a new species. Consequently, a formal species description is included. Worryingly, our data also show that geographically restricted populations of M. sp. #3 and its sister species (M. macarthurii) are at high risk of local and perhaps permanent extinction from both deforestation and habitat fragmentation.

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.

Bamboo specialization is one of the most extreme examples of convergent herbivory, yet it is unclear how this specific high-fiber diet might selectively shape the composition of the gut microbiome compared to host phylogeny. To address these questions, we used deep sequencing to investigate the nature and comparative impact of phylogenetic and dietary selection for specific gut microbial membership in three bamboo specialists-the bamboo lemur (Hapalemur griseus, Primates: Lemuridae), giant panda (Ailuropoda melanoleuca, Carnivora: Ursidae), and red panda (Ailurus fulgens, Carnivora: Musteloideadae), as well as two phylogenetic controls-the ringtail lemur (Lemur catta) and the Asian black bear (Ursus thibetanus). We detected significantly higher Shannon diversity in the bamboo lemur (10.029) compared to both the giant panda (8.256; p = 0.0001936) and the red panda (6.484; p = 0.0000029). We also detected significantly enriched bacterial taxa that distinguished each species. Our results complement previous work in finding that phylogeny predominantly governs high-level microbiome community structure. However, we also find that 48 low-abundance OTUs are shared among bamboo specialists, compared to only 8 OTUs shared by the bamboo lemur and its sister species, the ringtail lemur (Lemur catta, a generalist). Our results suggest that deep sequencing is necessary to detect low-abundance bacterial OTUs, which may be specifically adapted to a high-fiber diet. These findings provide a more comprehensive framework for understanding the evolution and ecology of the microbiome as well as the host.

Research on animal microbiomes is increasingly aimed at determining the evolutionary and ecological factors that govern host-microbiome dynamics, which are invariably intertwined and potentially synergistic. We present three empirical studies related to this topic, each of which relies on the diversity of Malagasy lemurs (representing a total of 19 species) and the comparative approach applied across scales of analysis. In Study 1, we compare gut microbial membership across 14 species in the wild to test the relative importance of host phylogeny and feeding strategy in mediating microbiome structure. Whereas host phylogeny strongly predicted community composition, the same feeding strategies shared by distant relatives did not produce convergent microbial consortia, but rather shaped microbiomes in host lineage-specific ways, particularly in folivores. In Study 2, we compare 14 species of wild and captive folivores, frugivores, and omnivores, to highlight the importance of captive populations for advancing gut microbiome research. We show that the perturbational effect of captivity is mediated by host feeding strategy and can be mitigated, in part, by modified animal management. In Study 3, we examine various scent-gland microbiomes across three species in the wild or captivity and show them to vary by host species, sex, body site, and a proxy of social status. These rare data provide support for the bacterial fermentation hypothesis in olfactory signal production and implicate steroid hormones as mediators of microbial community structure. We conclude by discussing the role of scale in comparative microbial studies, the links between feeding strategy and host-microbiome coadaptation, the underappreciated benefits of captive populations for advancing conservation research, and the need to consider the entirety of an animal's microbiota. Ultimately, these studies will help move the field from exploratory to hypothesis-driven research.

Phylogeographic analysis can be described as the study of the geological and climatological processes that have produced contemporary geographic distributions of populations and species. Here, we attempt to understand how the dynamic process of landscape change on Madagascar has shaped the distribution of a targeted clade of mouse lemurs (genus Microcebus ) and, conversely, how phylogenetic and population genetic patterns in these small primates can reciprocally advance our understanding of Madagascar's prehuman environment. The degree to which human activity has impacted the natural plant communities of Madagascar is of critical and enduring interest. Today, the eastern rainforests are separated from the dry deciduous forests of the west by a large expanse of presumed anthropogenic grassland savanna, dominated by the Family Poaceae, that blankets most of the Central Highlands. Although there is firm consensus that anthropogenic activities have transformed the original vegetation through agricultural and pastoral practices, the degree to which closed-canopy forest extended from the east to the west remains debated. Phylogenetic and population genetic patterns in a five-species clade of mouse lemurs suggest that longitudinal dispersal across the island was readily achieved throughout the Pleistocene, apparently ending at ∼55 ka. By examining patterns of both inter- and intraspecific genetic diversity in mouse lemur species found in the eastern, western, and Central Highland zones, we conclude that the natural environment of the Central Highlands would have been mosaic, consisting of a matrix of wooded savanna that formed a transitional zone between the extremes of humid eastern and dry western forest types.

Madagascar's lemurs display a diverse array of feeding strategies with complex relationships to seed dispersal mechanisms in Malagasy plants. Although these relationships have been explored previously on a case-by-case basis, we present here the first comprehensive analysis of lemuriform feeding, to our knowledge, and its hypothesized effects on seed dispersal and the long-term survival of Malagasy plant lineages. We used a molecular phylogenetic framework to examine the mode and tempo of diet evolution, and to quantify the associated morphological space occupied by Madagascar's lemurs, both extinct and extant. Using statistical models and morphometric analyses, we demonstrate that the extinction of large-bodied lemurs resulted in a significant reduction in functional morphological space associated with seed dispersal ability. These reductions carry potentially far-reaching consequences for Malagasy ecosystems, and we highlight large-seeded Malagasy plants that appear to be without extant animal dispersers. We also identify living lemurs that are endangered yet occupy unique and essential dispersal niches defined by our morphometric analyses.

Many studies have demonstrated the effects of host diet on gut microbial membership, metagenomics, and fermentation individually; but few have attempted to interpret the relationship among these biological phenomena with respect to host features (e.g. gut morphology). We quantitatively compare the fecal microbial communities, metabolic pathways, and fermentation products associated with the nutritional intake of frugivorous (fruit-eating) and folivorous (leaf-eating) lemurs. Our results provide a uniquely multidimensional and comparative perspective on the adaptive dynamics between host and microbiome. Shotgun metagenomic sequencing revealed significant differential taxonomic and metabolic pathway enrichment, tailored to digest and detoxify different diets. Frugivorous metagenomes feature pathways to degrade simple carbohydrates and host-derived glycosaminoglycans, while folivorous metagenomes are equipped to break down phytic acid and other phytochemical compounds in an anaerobic environment. We used nuclear magnetic resonance based metabolic profiling of fecal samples to link metabolic pathways to fermentation products, confirming that the dissimilar substrates provided in each diet select for specific microbial functions. Fecal samples from frugivorous lemurs contained significantly different profiles of short chain fatty acids, alcohol fermentation products, amino acids, glucose, and glycerol compared to folivorous lemurs. We present the relationships between these datasets as an integrated visual framework, which we refer to as microbial geometry. We use microbial geometry to compare empirical gut microbial profiles across different feeding strategies, and suggest additional utility as a tool for hypothesis-generation.

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.

The microbiome is now widely recognized as being important in health and disease, and makes up a substantial subset of the biome within the ecosystem of the vertebrate body. At the same time, multicellular, eukaryotic organisms such as helminths are being recognized as an important component of the biome that shaped the evolution of our genes. The absence of these macroscopic organisms during the early development and life of humans in Western culture probably leads to a wide range of human immunological diseases. However, the interaction between the microbiome and macroscopic components of the biome remains poorly characterized. In this study, the microbiome of the cecum in rats colonized for 2 generations with the small intestinal helminth Hymenolepis diminuta was evaluated. The introduction of this benign helminth, which is of considerable therapeutic interest, led to several changes in the cecal microbiome. Most of the changes were within the Firmicutes phylum, involved about 20% of the total bacteria, and generally entailed a shift from Bacilli to Clostridia species in the presence of the helminth. The results point toward ecological relationships between various components of the biome, with the observed shifts in the microbiome suggesting potential mechanisms by which this helminth might exert therapeutic effects.

Living nonhuman primates generally inhabit tropical forests, and torpor is regarded as a strategy employed by cold-adapted organisms. Yet, some primates employ daily torpor or hibernation (heterothermy) under obligatory, temporary, or emergency circumstances. Though heterothermy is present in most mammalian lineages, there are only three extant heterothermic primate lineages: bushbabies from Africa, lorises from Asia, and dwarf and mouse lemurs from Madagascar. Here, we analyze their phenotypes in the general context of tropical mammalian heterothermy. We focus on Malagasy lemurs as they have been the most intensively studied and also show an unmatched range of flexibility in their heterothermic responses. We discuss the evidence for whether heterothermy should be considered an ancestral or derived condition in primates. This consideration is particularly intriguing given that an understanding of the underlying mechanisms for hibernation in lemurs opens the possibility for insight into genotype-phenotype interactions, including those with biomedical relevance for humans.

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.

Humans first arrived on Madagascar only a few thousand years ago. Subsequent habitat destruction and hunting activities have had significant impacts on the island's biodiversity, including the extinction of megafauna. For example, we know of 17 recently extinct ‘subfossil’ lemur species, all of which were substantially larger (body mass ∼11–160 kg) than any living population of the ∼100 extant lemur species (largest body mass ∼6.8 kg). We used ancient DNA and genomic methods to study subfossil lemur extinction biology and update our understanding of extant lemur conservation risk factors by i) reconstructing a comprehensive phylogeny of extinct and extant lemurs, and ii) testing whether low genetic diversity is associated with body size and extinction risk. We recovered complete or near-complete mitochondrial genomes from five subfossil lemur taxa, and generated sequence data from population samples of two extinct and eight extant lemur species. Phylogenetic comparisons resolved prior taxonomic uncertainties and confirmed that the extinct subfossil species did not comprise a single clade. Genetic diversity estimates for the two sampled extinct species were relatively low, suggesting small historical population sizes. Low genetic diversity and small population sizes are both risk factors that would have rendered giant lemurs especially susceptible to extinction. Surprisingly, among the extant lemurs, we did not observe a relationship between body size and genetic diversity. The decoupling of these variables suggests that risk factors other than body size may have as much or more meaning for establishing future lemur conservation priorities.

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.

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.

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.

One of the most obvious physiological changes accompanying seasonal heterothermy in mammals is a fattening stage preceding periods of resource scarcity. This phenomenon reflects the interplay of both diet and physiology. Though the accrual of fat stores is known to be essential for overwintering in some species, the influence of diet on the physiology of torpor is not fully understood. Results from captive studies in heterothermic rodents and marsupials have indicated that when autumn diets are enriched with polyunsaturated fatty acids (PUFAs), animals receiving these diets experience deeper and more frequent torpor bouts than their counterparts receiving a control diet. Our study investigates this potential effect of dietary composition in animals that use daily torpor rather than prolonged torpor (i.e., hibernation). In so doing, we investigate the degree to which dietary effects on torpor are restricted to cold-adapted rodents and marsupials, or are a more general feature of mammalian heterothermy. We examined the effects of a PUFA diet and a control diet on the thermoregulation of one of the few species of primates known to use daily torpor: the grey mouse lemur (Microcebus murinus). Though the results of this study are largely inconclusive regarding the impact of dietary manipulations on torpor frequency and duration, we nonetheless find that the propensity of animals to enter torpor is directly influenced by age and seasonal changes in body mass, and thus reflect important physiological aspects of flexible thermoregulatory responses.

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.

The sagebrush lizards (Sceloporus graciosus group) consist of four taxa (S. graciosus graciosus, S. graciosus gracilis, S. graciosus vandenburgianus, and S. arenicolus) distributed in western North America. Of these, S. arenicolus is morphologically, behaviorally, and ecologically distinct as well as geographically disjunct from the other taxa, occurring only in the Mescalero-Monahans Sandhills of southeastern New Mexico and adjacent Texas. Sceloporus arenicolus is a taxon of concern because of its small range and habitat alteration due to land use practices. Understanding evolutionary relationships among members of the S. graciosus group, and especially S. arenicolus, has important implications for conservation. We examine the phylogenetic relationship of S. arenicolus relative to the three recognized subspecies of S. graciosus at mitochondrial and nuclear loci for populations sampled throughout the ranges of these taxa. Additionally, we estimate the divergence time and clade age of S. arenicolus. We find that the S. graciosus group is in need of major taxonomic revision, and also confirm that S. arenicolus is a genetically distinct and divergent lineage. These results bear important consequences for conservation and management.

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.

It is firmly believed that ancestral primates were nocturnal, with nocturnality having been maintained in most prosimian lineages. Under this traditional view, the opsin genes in all nocturnal prosimians should have undergone similar degrees of functional relaxation and accumulated similar extents of deleterious mutations. This expectation is rejected by the short-wavelength (S) opsin gene sequences from 14 representative prosimians. We found severe defects of the S opsin gene only in lorisiforms, but no defect in five nocturnal and two diurnal lemur species and only minor defects in two tarsiers and two nocturnal lemurs. Further, the nonsynonymous-to-synonymous rate ratio of the S opsin gene is highest in the lorisiforms and varies among the other prosimian branches, indicating different time periods of functional relaxation among lineages. These observations suggest that the ancestral primates were diurnal or cathemeral and that nocturnality has evolved several times in the prosimians, first in the lorisiforms but much later in other lineages. This view is further supported by the distribution pattern of the middle-wavelength (M) and long-wavelength (L) opsin genes among prosimians.

The living Malagasy lemurs constitute a spectacular radiation of >50 species that are believed to have evolved from a common ancestor that colonized Madagascar in the early Tertiary period. Yet, at least 15 additional Malagasy primate species, some of which were relative giants, succumbed to extinction within the past 2,000 years. Their existence in Madagascar is recorded predominantly in its Holocene subfossil record. To rigorously test the hypothesis that all endemic Malagasy primates constitute a monophyletic group and to determine the evolutionary relationships among living and extinct taxa, we have conducted an ancient DNA analysis of subfossil species. A total of nine subfossil individuals from the extinct genera Palaeopropithecus and Megaladapis yielded amplifiable DNA. Phylogenetic analysis of cytochrome b sequences derived from these subfossils corroborates the monophyly of endemic Malagasy primates. Our results support the close relationship of sloth lemurs to living indriids, as has been hypothesized on morphological grounds. In contrast, Megaladapis does not show a sister-group relationship with the living genus Lepilemur . Thus, the classification of the latter in the family Megaladapidae is misleading. By correlating the geographic location of subfossil specimens with relative amplification success, we reconfirm the global trend of increased success rates of ancient DNA recovery from nontropical localities.

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.

Host fitness is impacted by trillions of bacteria in the gastrointestinal tract that facilitate development and are inextricably tied to life history. During development, microbial colonization primes the gut metabolism and physiology, thereby setting the stage for adult nutrition and health. However, the ecological rules governing microbial succession are poorly understood. In this study, we examined the relationship between host lineage, captive diet, and life stage and gut microbiota characteristics in three primate species (infraorder, Lemuriformes). Fecal samples were collected from captive lemur mothers and their infants, from birth to weaning. Microbial DNA was extracted and the v4 region of 16S rDNA was sequenced on the Illumina platform using protocols from the Earth Microbiome Project. Here, we show that colonization proceeds along different successional trajectories in developing infants from species with differing dietary regimes and ecological profiles: frugivorous (fruit-eating) Varecia variegata, generalist Lemur catta, and folivorous (leaf-eating) Propithecus coquereli. Our analyses reveal community membership and succession patterns consistent with previous studies of human infants, suggesting that lemurs may serve as a useful model of microbial ecology in the primate gut. Each lemur species exhibits distinct species-specific bacterial diversity signatures correlating to life stages and life history traits, implying that gut microbial community assembly primes developing infants at species-specific rates for their respective adult feeding strategies.

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

Since its establishment in 1966, the Duke Lemur Center (DLC) has accumulated detailed records for nearly 4,200 individuals from over 40 strepsirrhine primate taxa—the lemurs, lorises, and galagos. Here we present verified data for 3,627 individuals of 27 taxa in the form of a life history table containing summarized species values for variables relating to ancestry, reproduction, longevity, and body mass, as well as the two raw data files containing direct and calculated variables from which this summary table is built. Large sample sizes, longitudinal data that in many cases span an animal’s entire life, exact dates of events, and large numbers of individuals from closely related yet biologically diverse primate taxa make these datasets unique. This single source for verified raw data and systematically compiled species values, particularly in combination with the availability of associated biological samples and the current live colony for research, will support future studies from an enormous spectrum of disciplines.

This collection of articles is focused on the evolutionary dynamics of heterothermy in mammals, specifically torpor and hibernation. Topics cover a wide range from evolutionary genetics, physiology, ecology, and applications to human health.

Debate over what is a species was already considered old hat when Darwin wrote his seminal abstract (as he called it) more than 150 years ago.(1) Endless papers, workshops, and symposia have been presented in an effort to "solve" the species problem. Yet, here we are, at it again. Has there been any progress? I believe that there has been, and that among the many advances enabled by the genomics revolution, progress on species concepts and species recognition is among them. To quote Feder and colleagues,(2) we are on the brink of a "unified theory of speciation genomics."

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.

Pattern and process are inextricably linked in biogeographic analyses, though we can observe pattern, we must infer process. Inferences of process are often based on ad hoc comparisons using a single spatial predictor. Here, we present an alternative approach that uses mixed-spatial models to measure the predictive potential of combinations of hypotheses. Biodiversity patterns are estimated from 8,362 occurrence records from 745 species of Malagasy amphibians and reptiles. By incorporating 18 spatially explicit predictions of 12 major biogeographic hypotheses, we show that mixed models greatly improve our ability to explain the observed biodiversity patterns. We conclude that patterns are influenced by a combination of diversification processes rather than by a single predominant mechanism. A 'one-size-fits-all' model does not exist. By developing a novel method for examining and synthesizing spatial parameters such as species richness, endemism and community similarity, we demonstrate the potential of these analyses for understanding the diversification history of Madagascar's biota.