Female reproductive aging in seven primate species: Patterns and consequences.

Age-related changes in fertility have increasingly been documented in wild animal populations: In many species the youngest and oldest reproducers are disadvantaged relative to prime adults. How do these effects evolve, and what explains their diversity across species? Tackling this question requires detailed data on patterns of age-related reproductive performance in multiple animal species. Here, we compare patterns and consequences of age-related changes in female reproductive performance in seven primate populations that have been subjects of long-term continuous study for 29 to 57 y. We document evidence of age effects on fertility and on offspring performance in most, but not all, of these primate species. Specifically, females of six species showed longer interbirth intervals in the oldest age classes, youngest age classes, or both, and the oldest females also showed relatively fewer completed interbirth intervals. In addition, five species showed markedly lower survival among offspring born to the oldest mothers, and two species showed reduced survival for offspring born to both the youngest and the oldest mothers. In contrast, we found mixed evidence that maternal age affects the age at which daughters first reproduce: Only in muriquis and to some extent in chimpanzees, the only two species with female-biased dispersal, did relatively young mothers produce daughters that tended to have earlier first reproduction. Our findings demonstrate shared patterns as well as contrasts in age-related changes in female fertility across species of nonhuman primates and highlight species-specific behavior and life-history patterns as possible explanations for species-level differences.

Developing evolutionary anthropology in local ecosystems.

The traditional regional focus of evolutionary anthropology-typically defined as places where hominin fossils, nonhuman primates, and non-western populations reside-forms the basis of much evolutionary anthropological research. Using the highly biodiverse temperate region of Appalachia as an example, we suggest that evolutionary anthropologists have much to gain by stepping outside of this traditional geographic area. Being purposely provocative, we argue that evolutionary anthropologists might also benefit from conducting research in Appalachia and other temperate ecosystems. We briefly discuss multiple areas of study-including studies of seed dispersal, functional redundancy, convergent evolution, human behavioral ecology, and conservation-and how they can be considered within the purview of integrative and evolutionary anthropology. We also highlight broader impacts to higher education that evolutionary anthropologists can help promote by working in local ecosystems.

Maternal death and offspring fitness in multiple wild primates.

Primate offspring often depend on their mothers well beyond the age of weaning, and offspring that experience maternal death in early life can suffer substantial reductions in fitness across the life span. Here, we leverage data from eight wild primate populations (seven species) to examine two underappreciated pathways linking early maternal death and offspring fitness that are distinct from direct effects of orphaning on offspring survival. First, we show that, for five of the seven species, offspring face reduced survival during the years immediately preceding maternal death, while the mother is still alive. Second, we identify an intergenerational effect of early maternal loss in three species (muriquis, baboons, and blue monkeys), such that early maternal death experienced in one generation leads to reduced offspring survival in the next. Our results have important implications for the evolution of slow life histories in primates, as they suggest that maternal condition and survival are more important for offspring fitness than previously realized.

Inflammatory Cytokine Profile in Individuals with Inherited Chromosomally Integrated Human Herpesvirus 6.

Acute graft-versus-host-disease (aGVHD) is a major complication following hematopoietic cell transplantations (HCTs). We have shown that HCT recipients in whom either the donor or patient had inherited chromosomally integrated human herpesvirus 6 (iciHHV-6) have a higher incidence of developing more severe aGVHD. Previous studies established that increased proinflammatory cytokines are associated with increased risk for aGVHD and nonrelapse mortality post-HCT. We hypothesized that HCT recipients with donor or recipient iciHHV-6 (iciHHV-6pos HCT cases) will have higher cytokine levels compared with HCT recipients without iciHHV-6 (iciHHV-6neg HCT controls). We identified 64 iciHHV-6pos HCT cases with plasma from days 7, 14, and/or 21 post-HCT and before aGVHD onset in patients who developed aGVHD. We identified 64 iciHHV-6neg HCT controls matched for aGVHD risk factors. We also identified 28 donors with iciHHV-6 and 56 matched donors without iciHHV-6. We measured plasma cytokine concentrations for IL-6, suppression of tumorigenicity 2, T cell immunoglobulin and mucin-domain containing 3, TNFα, soluble TNF receptor 1 (TNFRp55), and C-reactive protein (CRP). We used Mann-Whitney tests and repeated-measures models to compare cytokine levels. iciHHV-6pos HCT cases had higher CRP levels on day 7 and day 21 and higher TNFRp55 levels on day 14 and day 21 compared with iciHHV-6neg HCT controls. These findings were recapitulated in a repeated-measures model. The differences were most evident among patients who subsequently developed aGVHD grades 2 to 4. Additionally, iciHHV-6pos HCT cases had earlier-onset aGVHD (median, 20 versus 27 days post-HCT; P = .02). There were no differences in cytokine levels among healthy donors with or without iciHHV-6. This study demonstrates that HCT recipients with iciHHV-6 have higher proinflammatory cytokines that may be associated with increased risk for aGVHD.

Epigenetic Clocks.

Recent research has revealed clock-like patterns of epigenetic change across the life span in humans. Models describing these epigenetic changes have been dubbed "epigenetic clocks," and they can not only predict chronological age but also reveal biological age, which measures physiological homeostasis and deterioration over the life span. Comparative studies of the epigenetic clocks of different primate species are likely to provide insights into the evolution of life history schedules, as well as shed light on the physiological and genetic bases of aging and aging-related diseases. Chronological age estimation using clock-based calculators may also offer biological anthropologists a useful tool for applying to forensic and demographic studies.

Predictive Value of Clinical Findings and Plasma Biomarkers after Fourteen Days of Prednisone Treatment for Acute Graft-versus-host Disease.

We examined the hypothesis that plasma biomarkers and concomitant clinical findings after initial glucocorticoid therapy can accurately predict failure of graft-versus-host-disease (GVHD) treatment and mortality. We analyzed plasma samples and clinical data in 165 patients after 14 days of glucocorticoid therapy and used logistic regression and areas under receiver-operating characteristic curves (AUC) to evaluate associations with treatment failure and nonrelapse mortality (NRM). Initial treatment of GVHD was unsuccessful in 49 patients (30%). For predicting GVHD treatment failure, the best clinical combination (total serum bilirubin and skin GVHD stage: AUC, .70) was competitive with the best biomarker combination (T cell immunoglobulin and mucin domain 3 [TIM3] and [interleukin 1 receptor family encoded by the IL1RL1 gene, ST2]: AUC, .73). The combination of clinical features and biomarker results offered only a slight improvement (AUC, .75). For predicting NRM at 1 year, the best clinical predictor (total serum bilirubin: AUC, .81) was competitive with the best biomarker combination (TIM3 and soluble tumor necrosis factor receptor-1 [sTNFR1]: AUC, .85). The combination offered no improvement (AUC, .85). Infection was the proximate cause of death in virtually all patients. We conclude that after 14 days of glucocorticoid therapy, clinical findings (serum bilirubin, skin GVHD) and plasma biomarkers (TIM3, ST2, sTNFR1) can predict failure of GVHD treatment and NRM. These biomarkers reflect counter-regulatory mechanisms and provide insight into the pathophysiology of GVHD reactions after glucocorticoid treatment. The best predictive models, however, exhibit inadequate positive predictive values for identifying high-risk GVHD cohorts for investigational trials, as only a minority of patients with high-risk GVHD would be identified and most patients would be falsely predicted to have adverse outcomes.

Plasma biomarkers of acute GVHD and nonrelapse mortality: predictive value of measurements before GVHD onset and treatment.

We identified plasma biomarkers that presaged outcomes in patients with gastrointestinal graft-versus-host disease (GVHD) by measuring 23 biomarkers in samples collected before initiation of treatment. Six analytes with the greatest accuracy in predicting grade 3-4 GVHD in the first cohort (74 patients) were then tested in a second cohort (76 patients). The same 6 analytes were also tested in samples collected at day 14 ± 3 from 167 patients free of GVHD at the time. Logistic regression and calculation of an area under a receiver-operating characteristic (ROC) curve for each analyte were used to determine associations with outcome. Best models in the GVHD onset and landmark analyses were determined by forward selection. In samples from the second cohort, collected a median of 4 days before start of treatment, levels of TIM3, IL6, and sTNFR1 had utility in predicting development of peak grade 3-4 GVHD (area under ROC curve, 0.88). Plasma ST2 and sTNFR1 predicted nonrelapse mortality within 1 year after transplantation (area under ROC curve, 0.90). In the landmark analysis, plasma TIM3 predicted subsequent grade 3-4 GVHD (area under ROC curve, 0.76). We conclude that plasma levels of TIM3, sTNFR1, ST2, and IL6 are informative in predicting more severe GVHD and nonrelapse mortality.

Novel opsin gene variation in large-bodied, diurnal lemurs.

Some primate populations include both trichromatic and dichromatic (red-green colour blind) individuals due to allelic variation at the X-linked opsin locus. This polymorphic trichromacy is well described in day-active New World monkeys. Less is known about colour vision in Malagasy lemurs, but, unlike New World monkeys, only some day-active lemurs are polymorphic, while others are dichromatic. The evolutionary pressures underlying these differences in lemurs are unknown, but aspects of species ecology, including variation in activity pattern, are hypothesized to play a role. Limited data on X-linked opsin variation in lemurs make such hypotheses difficult to evaluate. We provide the first detailed examination of X-linked opsin variation across a lemur clade (Indriidae). We sequenced the X-linked opsin in the most strictly diurnal and largest extant lemur, Indri indri, and nine species of smaller, generally diurnal indriids (Propithecus). Although nocturnal Avahi (sister taxon to Propithecus) lacks a polymorphism, at least eight species of diurnal indriids have two or more X-linked opsin alleles. Four rainforest-living taxa-I. indri and the three largest Propithecus species-have alleles not previously documented in lemurs. Moreover, we identified at least three opsin alleles in Indri with peak spectral sensitivities similar to some New World monkeys.

Non-human primates avoid the detrimental effects of prenatal androgen exposure in mixed-sex litters: combined demographic, behavioral, and genetic analyses.

Producing single versus multiple births has important life history trade-offs, including the potential benefits and risks of sharing a common in utero environment. Sex hormones can diffuse through amniotic fluid and fetal membranes, and females with male littermates risk exposure to high levels of fetal testosterone, which are shown to have masculinizing effects and negative fitness consequences in many mammals. Whereas most primates give birth to single offspring, several New World monkey and strepsirrhine species regularly give birth to small litters. We examined whether neonatal testosterone exposure might be detrimental to females in mixed-sex litters by compiling data from long-term breeding records for seven primate species (Saguinus oedipus; Varecia variegata, Varecia rubra, Microcebus murinis, Mirza coquereli, Cheirogaleus medius, Galago moholi). Litter sex ratios did not differ from the expected 1:2:1 (MM:MF:FF for twins) and 1:2:2:1 (MMM:MMF:MFF:FFF for triplets). Measures of reproductive success, including female survivorship, offspring-survivorship, and inter-birth interval, did not differ between females born in mixed-sex versus all-female litters, indicating that litter-producing non-human primates, unlike humans and rodents, show no signs of detrimental effects from androgen exposure in mixed sex litters. Although we found no evidence for CYP19A1 gene duplications-a hypothesized mechanism for coping with androgen exposure-aromatase protein evolution shows patterns of convergence among litter-producing taxa. That some primates have effectively found a way to circumvent a major cost of multiple births has implications for understanding variation in litter size and life history strategies across mammals.

Paternity in wild ring-tailed lemurs (Lemur catta): Implications for male mating strategies.

In group-living species with male dominance hierarchies where receptive periods of females do not overlap, high male reproductive skew would be predicted. However, the existence of female multiple mating and alternative male mating strategies can call into question single-male monopolization of paternity in groups. Ring-tailed lemurs (Lemur catta) are seasonally breeding primates that live in multi-male, multi-female groups. Although established groups show male dominance hierarchies, male dominance relationships can break down during mating periods. In addition, females are the dominant sex and mate with multiple males during estrus, including group residents, and extra-group males-posing the question of whether there is high or low male paternity skew in groups. In this study, we analyzed paternity in a population of wild L. catta from the Bezà Mahafaly Special Reserve in southwestern Madagascar. Paternity was determined with 80-95% confidence for 39 offspring born to nine different groups. We calculated male reproductive skew indices for six groups, and our results showed a range of values corresponding to both high and low reproductive skew. Between 21% and 33% of offspring (3 of 14 or three of nine, counting paternity assignments at the 80% or 95% confidence levels, respectively) were sired by extra-troop males. Males siring offspring within the same group during the same year appear to be unrelated. Our study provides evidence of varying male reproductive skew in different L. catta groups. A single male may monopolize paternity across one or more years, while in other groups, >1 male can sire offspring within the same group, even within a single year. Extra-group mating is a viable strategy that can result in extra-group paternity for L. catta males.

Genetic Evidence for Male and Female Dispersal in Wild Lemur catta.

Lemur catta has traditionally been considered a species with male-biased dispersal; however, occasional female dispersal occurs. Using molecular data, we evaluated dispersal patterns in 2 L. catta populations in southwestern Madagascar: Tsimanampesotse National Park (TNP) and Bezà Mahafaly Special Reserve (BMSR). We also investigated the genetic differentiation between the populations and dispersal partner relatedness. Results showed minor genetic differentiation between the populations (Ï´(ST) = 0.039), which may indicate gene flow historically occurring in this region, made possible by the presence of L. catta groups between the sites. Different patterns of sex-biased dispersal were found between the sites using corrected assignment indices: male-biased dispersal in TNP, and a lack of sex-biased dispersal in BMSR. Observational evidence of female dispersal in BMSR supports these results and may imply intense female resource competition in and around BMSR, because small groups of 2-3 females have been observed dispersing within BMSR and entering the reserve from outside. These dispersing groups largely consisted of mothers transferring with daughters, although we have an aunt-niece pair transferring together. Genetic data suggest that males also transfer with relatives. Our data demonstrate that dispersal partners consist of same-sexed kin for L. catta males and females, highlighting the importance of kin selection.

Converging cytokine and metabolite networks shape asymmetric T cell fate at the term human maternal-fetal interface.

Placentation presents immune conflict between mother and fetus, yet in normal pregnancy maternal immunity against infection is maintained without expense to fetal tolerance. This is believed to result from adaptations at the maternal-fetal interface (MFI) which affect T cell programming, but the identities (i.e., memory subsets and antigenic specificities) of T cells and the signals that mediate T cell fates and functions at the MFI remain poorly understood. We found intact recruitment programs as well as pro-inflammatory cytokine networks that can act on maternal T cells in an antigen-independent manner. These inflammatory signals elicit T cell expression of co-stimulatory receptors necessary for tissue retention, which can be engaged by local macrophages. Although pro-inflammatory molecules elicit T cell effector functions, we show that additional cytokine (TGF-ß1) and metabolite (kynurenine) networks may converge to tune T cell function to those of sentinels. Together, we demonstrate an additional facet of fetal tolerance, wherein T cells are broadly recruited and restrained in an antigen-independent, cytokine/metabolite-dependent manner. These mechanisms provide insight into antigen-nonspecific T cell regulation, especially in tissue microenvironments where they are enriched.

A novel soluble form of Tim-3 associated with severe graft-versus-host disease.

The T cell Ig and mucin domain 3 (Tim-3) receptor has been implicated as a negative regulator of adaptive immune responses. We have utilized a proteomic strategy to identify novel proteins associated with graft versus host disease (GVHD) after allogeneic hematopoietic cell transplantation (HCT). Mass spectrometry analysis of plasma from subjects with mid-gut and upper-gut GVHD compared with those without GVHD identified increased levels of a protein identified with high confidence as Tim-3. A follow-up validation study using an immunoassay to measure Tim-3 levels in individual plasma samples from 127 patients demonstrated significantly higher plasma Tim-3 concentrations in patients with the more severe mid-gut GVHD, compared with those with upper-gut GVHD (P = .005), patients without GVHD (P = .002), and normal controls (P < .0001). Surface expression of Tim-3 was increased on CD8(+) T cells from patients with grade 2 to 4 acute GVHD (P = .01). Mass spectrometry-based profiling of plasma from multiple subjects diagnosed with common diseases provided evidence for restricted release of soluble Tim-3 in the context of GVHD. These findings have mechanistic implications for the development of novel strategies for targeting the Tim-3 immune regulatory pathway as an approach to improving control of GVHD.

Getting it right: Teaching undergraduate biology to undermine racial essentialism.

How we teach human genetics matters for social equity. The biology curriculum appears to be a crucial locus of intervention for either reinforcing or undermining students' racial essentialist views. The Mendelian genetic models dominating textbooks, particularly in combination with racially inflected language sometimes used when teaching about monogenic disorders, can increase middle and high school students' racial essentialism and opposition to policies to increase equity. These findings are of particular concern given the increasing spread of racist misinformation online and the misappropriation of human genomics research by white supremacists, who take advantage of low levels of genetics literacy in the general public. Encouragingly, however, teaching updated information about the geographical distribution of human genetic variation and the complex, multifactorial basis of most human traits, reduces students' endorsement of racial essentialism. The genetics curriculum is therefore a key tool in combating misinformation and scientific racism. Here, we describe a framework and example teaching materials for teaching students key concepts in genetics, human evolutionary history, and human phenotypic variation at the undergraduate level. This framework can be flexibly applied in biology and anthropology classes and adjusted based on time availability. Our goal is to provide undergraduate-level instructors with varying levels of expertise with a set of evidence-informed tools for teaching human genetics to combat scientific racism, including an evolving set of instructional resources, as well as learning goals and pedagogical approaches. Resources can be found at https://noto.li/YIlhZ5. Additionally, we hope to generate conversation about integrating modern genetics into the undergraduate curriculum, in light of recent findings about the risks and opportunities associated with teaching genetics.

Model selection, zero-inflated models, and predictors of primate abundance in Korup National Park, Cameroon.

Determining the ecological and anthropogenic factors that shape the abundance and distribution of wild primates is a critical component of primate conservation research. Such research is complicated, however, whenever the species under study are encountered infrequently, a characteristic of many taxa that are threatened with extinction. Typically, the resulting data sets based on surveys of such species will have a high frequency of zero counts which makes it difficult to determine the predictor variables that are associated with species abundance. In this study, we test various statistical models using survey data that was gathered on seven species of primate in Korup National Park, Cameroon. Predictor variables include hunting signs and aspects of habitat structure and floristic composition. Our statistical models include zero-inflated models that are tailored to deal with a high frequency of zero counts. First, using exploratory data analysis we found the most informative set of models as ranked by Δ-AIC (Akaike's information criterion). On the basis of this analysis, we used five predictor variables to construct several regression models including Poisson, zero-inflated Poisson, negative binomial, and zero-inflated negative binomial. Total basal area of all trees, density of secondary tree species, hunting signs, and mean basal area of all trees were significant predictors of abundance in the zero-inflated models. We discuss the statistical logic behind zero-inflated models and provide an interpretation of parameter estimates. We recommend that researchers explore a variety of models when determining the factors that correlate with primate abundance.

Evaluating ring-tailed lemurs (Lemur catta) from southwestern Madagascar for a genetic population bottleneck.

In light of historical and recent anthropogenic influences on Malagasy primate populations, in this study ring-tailed lemur (Lemur catta) samples from two sites in southwestern Madagascar, Beza Mahafaly Special Reserve (BMSR) and Tsimanampetsotsa National Park (TNP), were evaluated for the genetic signature of a population bottleneck. A total of 45 individuals (20 from BMSR and 25 from TNP) were genotyped at seven microsatellite loci. Three methods were used to evaluate these populations for evidence of a historical bottleneck: M-ratio, mode-shift, and heterozygosity excess tests. Three mutation models were used for heterozygosity excess tests: the stepwise mutation model (SMM), two-phase model (TPM), and infinite allele model (IAM). M-ratio estimations indicated a potential bottleneck in both populations under some conditions. Although mode-shift tests did not strongly indicate a population bottleneck in the recent historical past when samples from all individuals were included, a female-only analysis indicated a potential bottleneck in TNP. Heterozygosity excess was indicated under two of the three mutation models (IAM and TPM), with TNP showing stronger evidence of heterozygosity excess than BMSR. Taken together, these results suggest that a bottleneck may have occurred among L. catta in southwestern Madagascar in the recent past. Given knowledge of how current major stochastic climatic events and human-induced change can negatively impact extant lemur populations, it is reasonable that comparable events in the historical past could have caused a population bottleneck. This evaluation additionally functions to highlight the continuing environmental and anthropogenic challenges faced by lemurs in southwestern Madagascar.

Linking genotypes, phenotypes, and fitness in wild primate populations.

In the decade since the first draft of the human genome was announced, genome sequencing projects have been initiated for an additional twenty-some primate species. Within the next several years, genome sequence data will likely become available for all primate genera and for most individuals within some primate populations. At the same time, gene mapping and association studies of humans and other organisms are rapidly advancing our understanding of the genetic bases of behavioral and morphological traits. Primatologists are especially well-placed to take advantage of this coming flood of genetic data. Here we discuss what this new era of primate genomics means for field primatology and highlight some of the unprecedented opportunities it will afford, particularly with regard to examining the genetic basis of primate adaptation and diversity.

Demographic concepts and research pertaining to the study of wild primate populations.

Demography is the study of individuals as members of a population. The dynamics of a population are determined by collectively analyzing individual schedules of survival, growth, and reproduction. Together, these schedules are known as the vital rates of the population. The vital rates, along with dispersal, contribute to population structure, which refers to how the population is organized by age, sex, density, and social groups. I briefly review the history of anthropological demography as it pertains to wild primates and then I discuss basic demographic concepts and approaches for studying wild primate populations. I then turn to demographic studies of wild primate demography. Primates are generally characterized by high adult survival probabilities relative to survival at other age/stage classes and most primate populations have population growth rates near equilibrium. Changes in adult survival have the greatest impact on population growth rate (i.e., fitness) relative to other demographic traits such as juvenile/yearling survival or age at first reproduction. I discuss how these demographic patterns, and others, connect to topics and issues in behavioral ecology, life history theory, population genetics, and conservation biology. These connections help reaffirm the fact that the vital rates are both targets and agents of evolutionary change. In this regard, demographic studies of wild primates provide a critical link between the proximate socioecological processes that operate in a species and the long-term phylogenetic patterns that characterize a species.

Relevance of Plasma Matrix Metalloproteinase-9 for Bronchiolitis Obliterans Syndrome after Allogeneic Hematopoietic Cell Transplantation.

Bronchiolitis obliterans syndrome (BOS) is a highly morbid form of chronic graft-versus-host disease (GVHD) after allogeneic hematopoietic cell transplantation (HCT). Several plasma proteins have been identified as biomarkers for BOS after lung transplantation. The relevance of these biomarkers in BOS patients after allogeneic HCT has not been examined. We hypothesized that biomarkers associated with BOS after lung transplantation are also associated with BOS after allogeneic HCT. We tested plasma samples from 33 adult HCT patients who participated in a phase II multicenter study of fluticasone, azithromycin, and montelukast (FAM) treatment for new-onset BOS (NCT01307462), and matched control samples of HCT patients who had non-BOS chronic GVHD (n = 31) and those who never experienced chronic GVHD (n = 29) (NCT00637689 and NCT01902576). Candidate biomarkers included matrix metalloproteinase-9 (MMP-9), MMP-3, and chitinase-3-like-1 glycoprotein (YKL-40). MMP-9 concentrations were higher in the patients with BOS compared with those with non-BOS chronic GVHD (P = .04) or no chronic GVHD (P < .001). MMP-3 concentrations were higher in patients with BOS (P < .001) or non-BOS chronic GVHD (P < .001) compared with those with no chronic GVHD. YKL-40 concentrations did not differ statistically among the 3 groups. MMP-9 concentrations before starting FAM therapy were higher in patients who experienced treatment failure within 6 months compared with those with treatment success (P = .006), whereas MMP-3 or YKL-40 concentrations did not differ statistically between these 2 groups. Patients with an MMP-9 concentration ≥200,000 pg/mL before starting FAM therapy had worse overall survival compared with those with lower MMP-9 concentrations. Our data suggest that plasma MMP-9 concentration could serve as a relevant biomarker at diagnosis of BOS after allogeneic HCT for prognostication of survival and for prediction of treatment response. Further validation is needed to confirm our findings.

Comparative genomic analysis of sifakas (Propithecus) reveals selection for folivory and high heterozygosity despite endangered status.

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.