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

A chromosome-level, haplotype-resolved genome assembly and annotation for the Eurasian minnow (Leuciscidae: Phoxinus phoxinus) provide evidence of haplotype diversity

BACKGROUND

In this study, we present an in-depth analysis of the Eurasian minnow (Phoxinus phoxinus) genome, highlighting its genetic diversity, structural variations, and evolutionary adaptations. We generated an annotated haplotype-phased, chromosome-level genome assembly (2n = 50) by integrating high-fidelity (HiFi) long reads and chromosome conformation capture data (Hi-C).

RESULTS

We achieved a haploid size of 940 megabase pairs (Mbp) for haplome 1 and 929 Mbp for haplome 2 with high scaffold N50 values of 36.4 Mb and 36.6 Mb and BUSCO scores of 96.9% and 97.2%, respectively, indicating a highly complete genome assembly. We detected notable heterozygosity (1.43%) and a high repeat content (approximately 54%), primarily consisting of DNA transposons, which contribute to genome rearrangements and variations. We found substantial structural variations within the genome, including insertions, deletions, inversions, and translocations. These variations affect genes enriched in functions such as dephosphorylation, developmental pigmentation, phagocytosis, immunity, and stress response. In the annotation of protein-coding genes, 30,980 messenger RNAs and 23,497 protein-coding genes were identified with a high completeness score, which further underpins the high contiguity of our genome assemblies. We performed a gene family evolution analysis by comparing our proteome to 10 other teleost species, which identified immune system gene families that prioritize histone-based disease prevention over NB-LRR-related-based immune responses. Additionally, demographic analysis indicates historical fluctuations in the effective population size of P. phoxinus, likely correlating with past climatic changes.

CONCLUSIONS

This annotated, phased reference genome provides a crucial resource for resolving the taxonomic complexity within the genus Phoxinus and highlights the importance of haplotype-phased assemblies in understanding haplotype diversity in species characterized by high heterozygosity.

De Novo Genome Assembly for an Endangered Lemur Using Portable Nanopore Sequencing in Rural Madagascar

As one of the most threatened mammalian taxa, lemurs of Madagascar are facing unprecedented anthropogenic pressures. To address conservation imperatives such as this, researchers have increasingly relied on conservation genomics to identify populations of particular concern. However, many of these genomic approaches necessitate high-quality genomes. While the advent of next-generation sequencing technologies and the resulting reduction in associated costs have led to the proliferation of genomic data and high-quality reference genomes, global discrepancies in genomic sequencing capabilities often result in biological samples from biodiverse host countries being exported to facilities in the Global North, creating inequalities in access and training within genomic research. Here, we present the first published reference genome for the endangered red-fronted brown lemur (Eulemur rufifrons) from sequencing efforts conducted entirely within the host country using portable Oxford Nanopore sequencing. Using an archived E. rufifrons specimen, we conducted long-read, nanopore sequencing at the Centre ValBio Research Station near Ranomafana National Park, in rural Madagascar, generating over 750 Gb of sequencing data from 10 MinION flow cells. Exclusively using this long-read data, we assembled 2.157 gigabase, 2980-contig nuclear assembly with an N50 of 101.6 Mb and a 17,108 bp mitogenome. The nuclear assembly had 30× average coverage and was comparable in completeness to other primate reference genomes, with a 96.1% BUSCO completeness score for primate-specific genes. As the first published reference genome for E. rufifrons and the only annotated genome available for the speciose Eulemur genus, this resource will prove vital for conservation genomic studies while our efforts exhibit the potential of this protocol to address research inequalities and build genomic capacity.

Effects of Social Structure on Effective Population Size Change Estimates

Most methods currently used to infer the "demographic history of species" interpret this expression as a history of population size changes. The detection, quantification, and dating of demographic changes often rely on the assumption that population structure can be neglected. However, most vertebrates are typically organized in populations subdivided into social groups that are usually ignored in the interpretation of genetic data. This could be problematic since an increasing number of studies have shown that population structure can generate spurious signatures of population size change. Here, we simulate microsatellite data from a species subdivided into social groups where reproduction occurs according to different mating systems (monogamy, polygynandry, and polygyny). We estimate the effective population size (N e) and quantify the effect of social structure on estimates of changes in N e. We analyze the simulated data with two widely used methods for demographic inference. The first approach, BOTTLENECK, tests whether the samples are at mutation-drift equilibrium and thus whether a single N e can be estimated. The second approach, msvar, aims at quantifying and dating changes in N e. We find that social structure may lead to signals of departure from mutation-drift equilibrium including signals of expansion and bottlenecks. We also find that expansion signals may be observed under simple stationary Wright-Fisher models with low diversity. Since small populations tend to characterize many endangered species, we stress that methods trying to infer N e should be interpreted with care and validated with simulated data incorporating information about structure. Spurious expansion signals due to social structure can mask critical population size changes. These can obscure true bottleneck events and be particularly problematic in endangered species.

Ecological and anthropogenic effects on the genomic diversity of lemurs in Madagascar

Ecological variation and anthropogenic landscape modification have had key roles in the diversification and extinction of mammals in Madagascar. Lemurs represent a radiation with more than 100 species, constituting roughly one-fifth of the primate order. Almost all species of lemurs are threatened with extinction, but little is known about their genetic diversity and demographic history. Here, we analyse high-coverage genome-wide resequencing data from 162 unique individuals comprising 50 species of Lemuriformes, including multiple individuals from most species. Genomic diversity varies widely across the infraorder and yet is broadly consistent among individuals within species. We show widespread introgression in multiple genera and generally high levels of genomic diversity likely resulting from allele sharing that occurred during periods of connectivity and fragmentation during climatic shifts. We find distinct patterns of demographic history in lemurs across the ecogeographic regions of Madagascar within the last million years. Within the past 2,000 years, lemurs underwent major declines in effective population size that corresponded to the timing of human population expansion in Madagascar. In multiple regions of the island, we identified chronological trajectories of inbreeding that are consistent across genera and species, suggesting localized effects of human activity. Our results show how the extraordinary diversity of these long-neglected, endangered primates has been influenced by ecological and anthropogenic factors.

Reduction of bitter taste receptor gene family in folivorous colobine primates relative to omnivorous cercopithecine primates

Bitter taste perception is important in preventing animals from ingesting potentially toxic compounds. Whole-genome assembly (WGA) data have revealed that bitter taste receptor genes (TAS2Rs) comprise a multigene family with dozens of intact and disrupted genes in primates. However, publicly available WGA data are often incomplete, especially for multigene families. In this study, we employed a targeted capture (TC) approach specifically probing TAS2Rs for ten species of cercopithecid primates with diverse diets, including eight omnivorous cercopithecine species and two folivorous colobine species. We designed RNA probes for all TAS2Rs that we modeled to be intact in the common ancestor of cercopithecids ("ancestral-cercopithecid TAS2R gene set"). The TC was followed by short-read and high-depth massive-parallel sequencing. TC retrieved more intact TAS2R genes than found in WGA databases. We confirmed a large number of gene "births" at the common ancestor of cercopithecids and found that the colobine common ancestor and the cercopithecine common ancestor had contrasting trajectories: four gene "deaths" and three gene births, respectively. The number of intact TAS2R genes was markedly reduced in colobines (25-28 detected via TC and 20-26 detected via WGA analysis) as compared with cercopithecines (27-36 via TC and 19-30 via WGA). Birth or death events occurred at almost every phylogenetic-tree branch, making the composition of intact genes variable among species. These results show that evolutionary change in intact TAS2R genes is a complex process, refute a simple general prediction that herbivory favors more TAS2R genes, and have implications for understanding dietary adaptations and the evolution of detoxification abilities.

Consistency in Verreaux's sifaka home range and core area size despite seasonal variation in resource availability as assessed by Enhanced Vegetation Index (EVI)

Primates are adept at dealing with fluctuating availability of resources and display a range of responses to minimize the effects of food scarcity. An important component of primate conservation is to understand how primates adapt their foraging and ranging patterns in response to fluctuating food resources. Animals optimize resource acquisition within the home range through the selection of resource-bearing patches and choose between contrasting foraging strategies (resource-maximizing vs. area-minimizing). Our study aimed to characterize the foraging strategy of a folivorous primate, Verreaux's sifaka (Propithecus verreauxi), by evaluating whether group home range size varied between peak and lean leaf seasons within a seasonally dry tropical forest in Madagascar. We hypothesized that Verreaux's sifaka used the resource maximization strategy to select high-value resource patches so that during periods of resource depression, the home range area did not significantly change in size. We characterized resource availability (i.e., primary productivity) by season at Kirindy Mitea National Park using remotely-sensed Enhanced Vegetation Index data. We calculated group home ranges using 10 years of focal animal sampling data collected on eight groups using both 95% and 50% kernel density estimation. We used area accumulation curves to ensure each group had an adequate number of locations to reach seasonal home range asymptotes. Neither 95% home ranges nor 50% core areas differed across peak and lean leaf resource seasons, supporting the hypothesis that Verreaux's sifaka use a resource maximization strategy. With a better understanding of animal space use strategies, managers can model anticipated changes under environmental and/or anthropogenic resource depression scenarios. These findings demonstrate the value of long-term data for characterizing and understanding foraging and ranging patterns. We also illustrate the benefits of using satellite data for characterizing food resources for folivorous primates.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

Forest access restores foraging and ranging behavior in captive sifakas

Captive wildlife benefit from ecologically informed management strategies that promote natural behaviors. The Duke Lemur Center has pioneered husbandry programs rooted in species' ecology for a diversity of lemurs, including housing social groups in multiacre forest enclosures. We systematically document the foraging and ranging patterns of Coquerel's sifakas (Propithecus coquereli) living in these forest enclosures. Coquerel's sifakas are seasonal frugo-folivores that exhibit striking feeding flexibility in the wild. They are also one of the few members of the Indriidae family to persist in captivity. During all-day follows in the spring and summer of 2 consecutive years, we tracked the behavior of 14 sifakas in six forest enclosures. The sifakas' ranging and foraging patterns reflected those of wild sifakas in western Madagascar: On average, DLC sifakas occupied 3-day home ranges of 1.2 ha, traveled 473 m/day, and spent 26% of their time foraging for wild foodstuffs. The sifakas foraged most for young and mature leaves, fruits, nuts, and flowers from 39 plant species, especially red maple (Acer rubrum), tulip poplar (Liriodendron tulipifera), black locust (Robinia pseudoacacia), grapevine (Vitis rotundifolia), hickory (Carya spp.), and white oak (Quercus alba). Foraging patterns varied across seasons, enclosure areas, and groups, potentially reflecting differences in phenology, microhabitats, and individual preferences. While demonstrating that captive-bred primates express wild-like behaviors under ecologically relevant conditions, our results underscore the feeding flexibility of the Coquerel's sifaka. Captive wildlife exhibiting the range of species-specific behaviors are key resources for ecological research and might be best suited for future reintroductions.

Nuclear and Mitochondrial Phylogenomics of the Sifakas Reveal Cryptic Variation in the Diademed Sifaka

The most comprehensive phylogenomic reconstruction to date was generated on all nominal taxa within the lemur genus Propithecus. Over 200 wild-caught individuals were included in this study to evaluate the intra and interspecific relationships across this genus. Ultraconserved Elements (UCEs) resulted in well-supported phylogenomic trees. Complete mitochondrial genomes (CMGs) largely agreed with the UCEs, except where a mitochondrial introgression was detected between one clade of the diademed sifaka (Propithecus diadema) and the Milne-Edwards sifaka (P. edwardsi). Additionally, the crowned (P. coronatus) and Von der Decken’s (P. deckeni) sifakas belonged to a single admixed lineage from UCEs. Further sampling across these two species is warranted to determine if our sampling represents a hybrid zone. P. diadema recovered two well-supported clades, which were dated and estimated as being ancient as the split between the Perrier’s (P. perrierii) and silky (P. candidus) sifakas. The reconstructed demographic history of the two clades also varied over time. We then modeled the modern ecological niches of the two cryptic P. diadema clades and found that they were significantly diverged (p < 0.01). These ecological differences result in a very limited zone of geographic overlap for the P. diadema clades (<60 km2). Niche models also revealed that the Onive River acts as a potential barrier to dispersal between P. diadema and P. edwardsi. Further taxonomic work is required on P. diadema to determine if its taxonomic status should be revised. This first genomic evaluation of the genus resolved the relationships between the taxa and the recovered cryptic diversity within one species.

Genomic resources for rhesus macaques (Macaca mulatta)

Rhesus macaques (Macaca mulatta) are among the most extensively studied of nonhuman primates. This species has been the subject of many investigations concerning basic primate biology and behavior, including studies of social organization, developmental psychology, physiology, endocrinology, and neurodevelopment. Rhesus macaques are also critically important as a nonhuman primate model of human health and disease, including use in studies of infectious diseases, metabolic diseases, aging, and drug or alcohol abuse. Current research addressing fundamental biological and/or applied biomedical questions benefits from various genetic and genomic analyses. As a result, the genome of rhesus macaques has been the subject of more study than most nonhuman primates. This paper briefly discusses a number of information resources that can provide interested researchers with access to genetic and genomic data describing the content of the rhesus macaque genome, available information regarding genetic variation within the species, results from studies of gene expression, and other aspects of genomic analysis. Specific online databases are discussed, including the US National Center for Biotechnology Information, the University of California Santa Cruz genome browser, Ensembl genome browser, the Macaque Genotype and Phenotype database (mGAP), Rhesusbase, and others.

Molecular Adaptation to Folivory and the Conservation Implications for Madagascar’s Lemurs

The lemurs of Madagascar include numerous species characterized by folivory across several families. 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 &gt; 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 body of work demonstrating that not all primate folivores are built the same and reiterate 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.

Maternal stress effects on infant development in wild Verreaux's sifaka (Propithecus verreauxi)

Abstract

Maternal effects mediated by nutrients or specific endocrine states of the mother can affect infant development. Specifically, pre- and postnatal maternal stress associated with elevated glucocorticoid (GC) output is known to influence the phenotype of the offspring, including their physical and behavioral development. These developmental processes, however, remain relatively poorly studied in wild vertebrates, including primates with their relatively slow life histories. Here, we investigated the effects of maternal stress, assessed by fecal glucocorticoid output, on infant development in wild Verreaux’s sifakas (Propithecus verreauxi), a group-living Malagasy primate. In a first step, we investigated factors predicting maternal fecal glucocorticoid metabolite (fGCM) concentrations, how they impact infants’ physical and behavioral development during the first 6 months of postnatal life as well as early survival during the first 1.5 years of postnatal life. We collected fecal samples of mothers for hormone assays and behavioral data of 12 infants from two birth cohorts, for which we also assessed growth rates. Maternal fGCM concentrations were higher during the late prenatal but lower during the postnatal period compared to the early/mid prenatal period and were higher during periods of low rainfall. Infants of mothers with higher prenatal fGCM concentrations exhibited faster growth rates and were more explorative in terms of independent foraging and play. Infants of mothers with high pre- and postnatal fGCM concentrations were carried less and spent more time in nipple contact. Time mothers spent carrying infants predicted infant survival: infants that were more carried had lower survival, suggesting that they were likely in poorer condition and had to be cared for longer. Thus, the physical and behavioral development of these young primates were impacted by variation in maternal fGCM concentrations during the first 6 months of their lives, presumably as an adaptive response to living in a highly seasonal, but unpredictable environment.

Significance statement

The early development of infants can be impacted by variation in maternal condition. These maternal effects can be mediated by maternal stress (glucocorticoid hormones) and are known to have downstream consequences for behavior, physiology, survival, and reproductive success well into adulthood. However, the direction of the effects of maternal physiological GC output on offspring development is highly variable, even within the same species. We contribute comparative data on maternal stress effects on infant development in a Critically Endangered primate from Madagascar. We describe variation in maternal glucocorticoid output as a function of ecological and reproductive factors and show that patterns of infant growth, behavioral development, and early survival are predicted by maternal glucocorticoids. Our study demonstrates how mothers can influence offspring fitness in response to challenging environmental conditions.

Evolutionary and phylogenetic insights from a nuclear genome sequence of the extinct, giant, "subfossil" koala lemur Megaladapis edwardsi

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.