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.

Lemurs are a well-known example of adaptive radiation. Since colonizing Madagascar, more than 100 extant lemur species have evolved to fill the variety of ecological niches on the island. However, recent work suggests that lemurs do not exhibit one of the hallmarks of adaptive radiations: explosive speciation rates that decline over time. We test this idea using a phylogenomic dataset with broad taxonomic sampling of lemurs and their sister group, the lorisiforms of Asia and continental Africa. We find higher rates of speciation in Madagascar’s lemurs compared to lorisiforms and we confirm that lemurs did not experience an “early burst” of speciation after colonizing Madagascar. Instead, we identify three independent bursts of speciation approximately 15 million years ago that underly much of today’s lemur diversity. We demonstrate that the lemur clades with exceptionally high diversification rates have higher rates of introgression. This suggests that hybridization in these primates is not an evolutionary dead- end, but a driving force for diversification. Considering the conservation crisis affecting strepsirrhine primates, with approximately 95% of species being threatened with extinction, this phylogenomic study offers a new perspective for explaining Madagascar’s exceptional primate diversity and reveals patterns of speciation, extinction, and gene flow that will help inform future conservation decisions.

Summary1) Research Aims — The extent of Madagascar’s grasslands prior to human colonization is unresolved. We used population genetic analyses of a broadly dominant C4fire-adapted grass,Loudetia simplex, as a proxy for estimating grassland change through time. In the absence of population genomic resources, we used target-enrichment data. We carefully examined the utility of these data for population genetics to make recommendations on cost-effective strategies for conservation genetics. We explored the potential of estimating individual ploidy levels from target-enrichment data and how assumptions about ploidy could affect analyses.2) Methods — We developed a novel bioinformatic pipeline based to estimate ploidy and genotypes from target-enrichment data. We estimated standard population genetic summary statistics in addition to species trees and population structure. Extended Bayesian skyline plots provided estimates of population size through time for empirical and simulated data.3) Key Result — All MalagasyLoudetia simplexindividuals sampled in this study formed a clade and possibly indicated an ancestral high-altitude distribution near the present-day Itremo protected area. Demographic models suggested grassland expansions occurred prior to the Last Interglacial Period and supported extensive grasslands prior to human colonization. Though there are limitations to target-enrichment data for population genetic studies, we find that analyses of population structure are reliable.4) Key Point —Genetic variation inLoudetia simplexsupports widespread grasslands in Madagascar prior to the more recent periods of notable paleoclimatic change. However, the methods explored here could not differentiate between paleoclimatic change near the Last Glacial Maximum and anthropogenic effects. Target-enrichment data can be a valuable tool for analyses of population structure in the absence of reference genomic resources.

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 gut microbiome is gaining recognition for its role in primate nutrition, but we stand to benefit from microbiome comparisons across diverse hosts and environmental conditions. We compared gut microbiome structure in four lemur species from four phylogenetic lineages, including 9 individual mouse lemurs (Microcebus danfossi), 6 brown lemurs (Eulemur fulvus), 20 sifakas (Propithecus coquereli), and a single sportive lemur (Lepilemur grewcockorum). In northwestern Madagascar, these species are sympatric, but use different feeding strategies to cope with environmental challenges, including relying on tree gums and insects (mouse lemurs), and some vs. significant leaf matter (brown lemurs vs. sifakas and sportive lemurs). From one fecal sample collected per lemur in the dry season in the Anjajavy Forest, we determined gut microbiome diversity, variability, and membership via 16S rRNA sequencing. The lemurs harbored strongly species-specific gut microbiomes. Brown lemurs showed more diverse and generalized consortia; mouse lemurs, sifakas, and the sportive lemur had less diverse consortia with more distinct memberships. Consistent with their fallback foods, mouse lemur microbiomes included taxa putatively associated with gum and insect digestion, whereas those of sifakas and the sportive lemur showed stronger and distinct signatures of leaf fiber and secondary compound metabolism. These results point to feeding strategy, intertwined with host phylogeny, as a driver of gut microbiome composition, but highlight real-time dietary specificity as a contributing driver of microbiome diversity. While illuminating how gut microbiomes facilitate host nutrition on challenging foods, these results help explain how ecologically diverse primates living in sympatry may differentially cope with seasonal or stochastic lean times.

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.

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

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.

DNA sequences from three mitochondrial genes and one nuclear gene were analyzed to determine the phylogeny of the Malagasy primate family Lemuridae. Whether analyzed separately or in combination, the data consistently indicate that Eulemur species comprise a clade that is sister to a Lemur catta plus Hapalemur clade. The genus Varecia is basal to both. Resolution of cladogenic events within Eulemur was found to be extremely problematic with a total of six alternative arrangements offered by various data sets and weighting regimes. We attempt to determine the best arrangement of Eulemur taxa through a variety of character and taxon sampling strategies. Because our study includes all but one Eulemur species, increased taxon sampling is probably not an option for enhancing phylogenetic accuracy. We find, however, that the combined genetic data set is more robust to changes in taxon sample than are any of the individual data sets, suggesting that increased character sampling stabilizes phylogenetic resolution. Nonetheless, due to the difficult nature of the problem, we may have to accept certain aspects of Eulemur interrelationships as uncertain.

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 [...]

In the past decade, several studies have estimated the human per-generation germline mutation rate using large pedigrees. More recently, estimates for various non-human 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 how to appropriately account 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 two-fold 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.

Despite being one of the most fundamental biological processes, the process of speciation remains poorly understood in many groups of organisms. Mouse lemurs are a species-rich genus of small primates endemic to Madagascar, whose diversity has only recently been uncovered using genetic data and is primarily found among morphologically cryptic, allopatric populations. To assess to what extent described species represent reproductively isolated entities, studies are needed in areas where mouse lemur taxa come into contact. Hybridization has previously been reported in a contact zone between two closely related mouse lemur species (Microcebus murinus and M. griseorufus) based on microsatellite data. Here, we revisit this system using RADseq data for populations in, near, and far from the contact zone, including many of the individuals that had previously been identified as hybrids. Surprisingly, we find no evidence for admixed nuclear ancestry in any of the individuals. Re-analyses of microsatellite data and simulations suggest that previously inferred hybrids were false positives and that the program NewHybrids can be particularly sensitive to erroneously inferring hybrid ancestry. Using coalescent-bases analyses, we also show an overall lack of recent gene flow between the two species, and low levels of ancestral gene flow. Combined with evidence for local syntopic occurrence, these data indicate that M. murinus and M. griseorufus are reproductively isolated. Finally, we estimate that they diverged less than a million years ago, suggesting that completion of speciation is relatively rapid in mouse lemurs. Future work should focus on the underpinnings of reproductive isolation in this cryptic primate radiation, which are mostly unknown. Our study also provides a cautionary tale for the inference of hybridization with microsatellite data.

Folivory evolved independently at least three times over the last 40 million years among Madagascar’s lemurs. Many extant lemuriform folivores exist in sympatry in Madagascar’s remaining forests. These species avoid feeding competition by adopting different dietary strategies within folivory, reflected in behavioral, morphological, and microbiota diversity across species. These conditions make lemurs an ideal study system for understanding adaptation to leaf-eating. Most folivorous lemurs are also highly endangered. The significance of folivory for conservation outlook is complex. Though generalist folivores may be relatively well equipped to survive habitat disturbance, specialist folivores occupying narrow dietary niches may be less resilient. Characterizing the genetic bases of adaptation to folivory across species and lineages can provide insights into their differential physiology and potential to resist habitat change. We recently reported accelerated genetic change in RNASE1, a gene encoding an enzyme (RNase 1) involved in molecular adaptation in mammalian folivores, including various monkeys and sifakas (genus Propithecus; family Indriidae). Here, we sought to assess whether other lemurs, including phylogenetically and ecologically diverse folivores, might show parallel adaptive change in RNASE1 that could underlie a capacity for efficient folivory. We characterized RNASE1 in 21 lemur species representing all five families and members of the three extant folivorous lineages: 1) bamboo lemurs (family Lemuridae), 2) sportive lemurs (family Lepilemuridae), and 3) indriids (family Indriidae). We found pervasive sequence change in RNASE1 across all indriids, a dN/dS value > 3 in this clade, and evidence for shared change in isoelectric point, indicating altered enzymatic function. Sportive and bamboo lemurs, in contrast, showed more modest sequence change. The greater change in indriids may reflect a shared strategy emphasizing complex gut morphology and microbiota to facilitate folivory. This case study illustrates how genetic analysis may reveal differences in functional traits that could influence species’ ecology and, in turn, their resilience to habitat change. Moreover, our results support the contention that not all primate folivores are built the same and highlight the need to avoid generalizations about dietary guild in considering conservation outlook, particularly in lemurs where such diversity in folivory has probably led to extensive specialization via niche partitioning.

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.

Accurately reconstructing the reticulate histories of polyploids remains a central challenge for understanding plant evolution. Although phylogenetic networks can provide insights into relationships among polyploid lineages, inferring networks may be hampered by the complexities of homology determination in polyploid taxa. We use simulations to show that phasing alleles from allopolyploid individuals can improve inference of phylogenetic networks under the multispecies coalescent. Phased allelic data can also improve divergence time estimates for networks, which is helpful for evaluating allopolyploid speciation hypotheses and proposing mechanisms of speciation. To achieve these outcomes, we present a novel pipeline that leverages a recently developed phasing algorithm to reliably phase alleles from polyploids. This pipeline is especially appropriate for target enrichment data, where depth of coverage is typically high enough to phase entire loci. We provide an empirical example in the North American Dryopteris fern complex that demonstrates how phasing can help reveal the mode of polyploidization and improve network inference. We establish that our pipeline (PATÉ: Phased Alleles from Target Enrichment data) is capable of recovering a high proportion of phased loci from both diploids and polyploids, and that these data improve network estimates compared to using haplotype consensus assemblies. This approach is shown to be especially effective in reticulate complexes where there are multiple hybridization events. The pipeline is available at: https://github.com/gtiley/Phasing.

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.

Germline mutations are the raw material for natural selection, driving species evolution and the creation of earth’s biodiversity. Life on earth would stagnate without this driver of genetic diversity. 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 spectrum by which these mutations accrue across the tree of life. Advances in high-throughput sequencing are providing new opportunities for characterizing these rates and patterns 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 the applications and limitations of whole-genome comparisons among closely related individuals in what are typically described as “trio” analyses for the detection of germline mutations as they arise in real time. By sequencing and comparing whole-genomes generated for individuals of known relatedness – typically, parent to offspring – investigators can ideally count and characterize mutations as they appear per generation. The promise for gaining insight into classic hypotheses of molecular evolution is high, though so too is the cost. Namely, the technical challenges are daunting given that pedigree-based studies are essentially searching for needles in a haystack. Even so, the opportunities are so enticing, and the field so young, we can say with confidence that fundamental insights have only just begun to emerge.

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.

The Central Highland Plateau of Madagascar is largely composed of grassland savanna, interspersed with patches of closed-canopy forest. Conventional wisdom has it that these grasslands are anthropogenic in nature, having been created very recently via human agricultural practices. Yet, the ancient origins of the endemic grasses suggest that the extensive savannas are natural biomes, similar to others found around the globe. We use a phylogeographic approach to compare these two competing scenarios. By sampling multiple populations of Goodman’s mouse lemur (Microcebus lehilahytsara), a small-bodied nocturnal primate, we reconstruct the phylogeographic and demographic history of these “environmental metronomes” to estimate the time at which their populations diverged, and thus proximally, when their habitats would have become fragmented. We applied coalescent methods to RADseq data to infer phylogenetic relationships, population structure, and migration corridors among sampling sites. These analyses indicate that forest fragmentation occurred rapidly during a period of decreased precipitation near the last glacial maximum and would have affected both the Central Highlands and eastern forests. Though there is clear genomic structure separating the populations of the Central Highland from those of the eastern rainforests, there is also evidence of historical migration between them. Findings support the hypothesis that the Central Highland savanna predates human arrival, indicating that it is a natural landscape that has long impacted the population dynamics of Goodman’s mouse lemur, and by extension, other forest-dwelling organisms in Madagascar.

Spontaneous germline 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 ofde novomutations have been directly characterized in only a few lineages. It is therefore critical to expand the phylogenetic scope of these studies to gain a more general understanding of observed mutation rate patterns. Our study provides the first direct mutation rate estimate for a strepsirrhine (i.e., the lemurs and lorises), which comprise 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 1.64 × 10−8(95% credible interval: 1.41 × 10−8to 1.98 × 10−8) mutations/site/generation. This estimate is higher than those measured for most previously characterized mammals. Further, we found an unexpectedly low count of paternal mutations, and only a modest overrepresentation of mutations at CpG-sites. Given the surprising nature of these observations, we conducted an independent analysis of context-dependent substitution types for gray mouse lemur and five additional primate species. This analysis yielded patterns consistent with the mutation spectrum from the pedigree mutation-rate analysis, which provides confidence in our ability to accurately identifyde novomutations with our data and bioinformatic filters.

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.

Sensory gene families are of special interest, both for what they can tell us about molecular evolution, and for 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 (generaMicrocebusandCheirogaleus, 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. Further, a detailed syntenic comparison of mouse lemurs with mouse (genusMus) 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 V1R “hotspots.”

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.