Molecular Phylogenetic Relationships and Unveiling Novel Genetic Diversity among Slow and Pygmy Lorises, including Resurrection of Xanthonycticebus intermedius.

Genetic analysis of historical museum collections presents an opportunity to clarify the evolutionary history of understudied primate groups, improve taxonomic inferences, and inform conservation efforts. Among the most understudied primate groups, slow and pygmy lorises (genera Nycticebus and Xanthonycticebus) are nocturnal strepsirrhines found in South and Southeast Asia. Previous molecular studies have supported five species, but studies using morphological data suggest the existence of at least nine species. We sequenced four mitochondrial loci, CO1, cytb, d-loop, and ND4, for a total of 3324 aligned characters per sample from 41 historical museum specimens for the most comprehensive geographic coverage to date for these genera. We then combined these sequences with a larger dataset composed of samples collected in Vietnam as well as previously published sequences (total sample size N = 62). We inferred phylogenetic relationships using Bayesian inference and maximum likelihood methods based on data from each locus and on concatenated sequences. We also inferred divergence dates for the most recent common ancestors of major lineages using a BEAST analysis. Consistent with previous studies, we found support for Xanthonycticebus pygmaeus as a basal taxon to the others in the group. We also confirmed the separation between lineages of X. pygmaeus from northern Vietnam/Laos/China and southern Vietnam/Cambodia and included a taxonomic revision recognizing a second taxon of pygmy loris, X. intermedius. Our results found support for multiple reciprocally monophyletic taxa within Borneo and possibly Java. The study will help inform conservation management of these trade-targeted animals as part of a genetic reference database for determining the taxonomic unit and provenance of slow and pygmy lorises confiscated from illegal wildlife trade activities.

Functional genomics analysis reveals the evolutionary adaptation and demographic history of pygmy lorises.

Lorises are a group of globally threatened strepsirrhine primates that exhibit many unusual physiological and behavioral features, including a low metabolic rate, slow movement, and hibernation. Here, we assembled a chromosome-level genome sequence of the pygmy loris (Xanthonycticebus pygmaeus) and resequenced whole genomes from 50 pygmy lorises and 6 Bengal slow lorises (Nycticebus bengalensis). We found that many gene families involved in detoxification have been specifically expanded in the pygmy loris, including the GSTA gene family, with many newly derived copies functioning specifically in the liver. We detected many genes displaying evolutionary convergence between pygmy loris and koala, including PITRM1. Significant decreases in PITRM1 enzymatic activity in these two species may have contributed to their characteristic low rate of metabolism. We also detected many evolutionarily convergent genes and positively selected genes in the pygmy loris that are involved in muscle development. Functional assays demonstrated the decreased ability of one positively selected gene, MYOF, to up-regulate the fast-type muscle fiber, consistent with the lower proportion of fast-twitch muscle fibers in the pygmy loris. The protein product of another positively selected gene in the pygmy loris, PER2, exhibited weaker binding to the key circadian core protein CRY, a finding that may be related to this species' unusual circadian rhythm. Finally, population genomics analysis revealed that these two extant loris species, which coexist in the same habitat, have exhibited an inverse relationship in terms of their demography over the past 1 million years, implying strong interspecies competition after speciation.

Variation in predicted COVID-19 risk among lemurs and lorises.

The novel coronavirus SARS-CoV-2, which in humans leads to the disease COVID-19, has caused global disruption and more than 2 million fatalities since it first emerged in late 2019. As we write, infection rates are at their highest point globally and are rising extremely rapidly in some areas due to more infectious variants. The primary target of SARS-CoV-2 is the cellular receptor angiotensin-converting enzyme-2 (ACE2). Recent sequence analyses of the ACE2 gene predict that many nonhuman primates are also likely to be highly susceptible to infection. However, the anticipated risk is not equal across the Order. Furthermore, some taxonomic groups show high ACE2 amino acid conservation, while others exhibit high variability at this locus. As an example of the latter, analyses of strepsirrhine primate ACE2 sequences to date indicate large variation among lemurs and lorises compared to other primate clades despite low sampling effort. Here, we report ACE2 gene and protein sequences for 71 individual strepsirrhines, spanning 51 species and 19 genera. Our study reinforces previous results while finding additional variability in other strepsirrhine species, and suggests several clades of lemurs have high potential susceptibility to SARS-CoV-2 infection. Troublingly, some species, including the rare and endangered aye-aye (Daubentonia madagascariensis), as well as those in the genera Avahi and Propithecus, may be at high risk. Given that lemurs are endemic to Madagascar and among the primates at highest risk of extinction globally, further understanding of the potential threat of COVID-19 to their health should be a conservation priority. All feasible actions should be taken to limit their exposure to SARS-CoV-2.

Species Identification and Unlocking Hidden Genetic Diversity of Confiscated Slow Lorises (Nycticebus spp.) Based on Mitochondrial DNA Markers.

Confiscated slow lorises (Nycticebus spp.) at Bangpra Water-Bird Breeding Center (BWBC) in Thailand provided an opportunity to demonstrate the application of noninvasive genetic approaches for species identification when morphology of the animals was ambiguous. The slow lorises at BWBC had been assigned to either N. bengalensis or N. pygmaeus, based on body size. However, the morphology of N. bengalensis is highly variable and overlaps with that of N. coucang (sensu stricto). Phylogenetic analysis of cytochrome b and d-loop mitochondrial regions placed all confiscated N. pygmaeus with the published sequences of N. pygmaeus and distinguished them from other Nycticebus. All other confiscated individuals formed a monophyletic clade, most individuals grouping with published N. bengalensis sequences from wild populations in Vietnam and distinct from Peninsular Malaysian and Sumatran N. coucang, Javan N. javanicus and Bornean N. menagensis. Six individuals within the N. bengalensis clade formed a separate subgroup that did not group with any reference material as indicated by phylogenetic and haplotype network analyses. Whether these trafficked individuals are undiscovered wild populations will require further investigation. Additional genetic studies of wild slow loris populations in different regions are therefore urgently required for reference to aid the protection and conservation of these threatened species.

Using a multi-gene approach to infer the complicated phylogeny and evolutionary history of lorises (Order Primates: Family Lorisidae).

Extensive phylogenetic studies have found robust phylogenies are modeled by using a multi-gene approach and sampling from the majority of the taxa of interest. Yet, molecular studies focused on the lorises, a cryptic primate family, have often relied on one gene, or just mitochondrial DNA, and many were unable to include all four genera in the analyses, resulting in inconclusive phylogenies. Past phylogenetic loris studies resulted in lorises being monophyletic, paraphyletic, or an unresolvable trichotomy with the closely related galagos. The purpose of our study is to improve our understanding of loris phylogeny and evolutionary history by using a multi-gene approach. We used the mitochondrial genes cytochrome b, and cytochrome c oxidase subunit 1, along with a nuclear intron (recombination activating gene 2) and nuclear exon (the melanocortin 1 receptor). Maximum Likelihood and Bayesian phylogenetic analyses were conducted based on data from each locus, as well as on the concatenated sequences. The robust, concatenated results found lorises to be a monophyletic family (Lorisidae) (PP ≥ 0.99) with two distinct subfamilies: the African Perodictinae (PP ≥ 0.99) and the Asian Lorisinae (PP ≥ 0.99). Additionally, from these analyses all four genera were all recovered as monophyletic (PP ≥ 0.99). Some of our single-gene analyses recovered monophyly, but many had discordances, with some showing paraphyly or a deep-trichotomy. Bayesian partitioned analyses inferred the most recent common ancestors of lorises emerged ∼42 ±â€¯6 million years ago (mya), the Asian Lorisinae separated ∼30 ±â€¯9 mya, and Perodictinae arose ∼26 ±â€¯10 mya. These times fit well with known historical tectonic shifts of the area, as well as with the sparse loris fossil record. Additionally, our results agree with previous multi-gene studies on Lorisidae which found lorises to be monophyletic and arising ∼40 mya (Perelman et al., 2011; Pozzi et al., 2014). By taking a multi-gene approach, we were able to recover a well-supported, monophyletic loris phylogeny and inferred the evolutionary history of this cryptic family.

Hypervariability of Nucleolus Organizer Regions in Bengal Slow Lorises, Nycticebus bengalensis (Primates, Lorisidae).

Slow lorises are a cryptic species complex, and thus genetic markers are needed to identify distinct evolutionary lineages or species. We examined the nucleolus organizer regions (NORs) of Bengal slow lorises (Nycticebus bengalensis) using FISH with 18S rDNA (rDNA-FISH) and silver nitrate staining (Ag-NOR stain). Ten individuals of the putatively single species N. bengalensis showed higher variability in localization than 3 other congeners, though their overall karyotypes were similar. The rDNA-FISH analysis detected a total of 18 loci, in contrast to previous studies of other slow loris species that revealed far fewer (6-10) loci. Eight of the 18 loci detected in the present analysis were found to be semi-stable localizations at 4 different chromosomes, while 10 were found to be unstable localizations at 5 other chromosomes. The semi-stable locations showed occasional presence/absence of variations for rDNA-FISH, and unstable locations were polymorphic among individuals, contributing to the higher variability of NORs in this taxon. We hypothesize that the larger numbers of rDNA loci found in N. bengalensis were introduced by genomic dispersion through ectopic recombination in association with terminal regions including rDNA. Such differences are potentially very powerful chromosomal markers to be used in species identification and conservation.

Molecular cytogenetic studies in strepsirrhine primates, Dermoptera and Scandentia.

Since the first chromosome painting study between human and strepsirrhine primates was performed in 1996, nearly 30 species in Strepsirrhini, Dermoptera and Scandentia have been analyzed by cross-species chromosome painting. Here, the contribution of chromosome painting data to our understanding of primate genome organization, chromosome evolution and the karyotype phylogenetic relationships within strepsirrhine primates, Dermoptera and Scandentia is reviewed. Twenty-six to 43 homologous chromosome segments have been revealed in different species with human chromosome-specific paint probes. Various landmark rearrangements characteristic for each different lineage have been identified, as cytogenetic signatures that potentially unite certain lineages within strepsirrhine primates, Dermoptera and Scandentia.

Phylogenetic relationships among the Lorisoidea as indicated by craniodental morphology and mitochondrial sequence data.

The phylogeny of the Afro-Asian Lorisoidea is controversial. While postcranial data attest strongly to the monophyly of the Lorisidae, most molecular analyses portray them as paraphyletic and group the Galagidae alternately with the Asian or African lorisids. One of the problems that has bedevilled phylogenetic analysis of the group in the past is the limited number of taxa sampled for both ingroup families. We present the results of a series of phylogenetic analyses based on 635 base pairs (bp) from two mitochondrial genes (12S and 16S rRNA) with and without 36 craniodental characters, for 11 galagid and five lorisid taxa. The outgroup was the gray mouse lemur (Microcebus murinus). Analyses of the molecular data included maximum parsimony (MP), neighbor joining (NJ), maximum likelihood (ML), and Bayesian methods. The model-based analyses and the combined "molecules+morphology" analyses supported monophyly of the Lorisidae and Galagidae. The lorisids form two geographically defined clades. We find no support for the taxonomy of Galagidae as proposed recently by Groves [Primate Taxonomy, Washington, DC: Smithsonian Institution Press. 350 p, 2001]. The taxonomy of Nash et al. [International Journal of Primatology 10:57-80, 1989] is supported by the combined "molecules+morphology" analysis; however, the model-based analyses suggest that Galagoides may be an assemblage of species united by plesiomorphic craniodental characters.

Multidirectional chromosome painting reveals a remarkable syntenic homology between the greater galagos and the slow loris.

We report on the first reciprocal chromosome painting of lorisoids and humans. The chromosome painting showed a remarkable syntenic homology between Otolemur and Nycticebus. Eight derived syntenic associations of human segments are common to both Otolemur and Nycticebus, indicative of a considerable period of common evolution between the greater galago and the slow loris. Five additional Robertsonian translocations form the slow loris karyotype, while the remaining chromosomes are syntenically equivalent, although some differ in terms of centromere position and heterochromatin additions. Strikingly, the breakpoints of the human chromosomes found fragmented in these two species are apparently identical. Only fissions of homologs to human chromosomes 1 and 15 provide significant evidence of a cytogenetic link between Lemuriformes and Lorisiformes. The association of human chromosomes 7/16 in both lorisoids strongly suggests that this chromosome was present in the ancestral primate genome.

Chromosome painting between human and lorisiform prosimians: evidence for the HSA 7/16 synteny in the primate ancestral karyotype.

Multidirectional chromosome painting with probes derived from flow-sorted chromosomes of humans (Homo sapiens, HSA, 2n = 46) and galagos (Galago moholi, GMO, 2n = 38) allowed us to map evolutionarily conserved chromosomal segments among humans, galagos, and slow lorises (Nycticebus coucang, NCO, 2n = 50). In total, the 22 human autosomal painting probes detected 40 homologous chromosomal segments in the slow loris genome. The genome of the slow loris contains 16 sytenic associations of human homologues. The ancient syntenic associations of human chromosomes such as HSA 3/21, 7/16, 12/22 (twice), and 14/15, reported in most mammalian species, were also present in the slow loris genome. Six associations (HSA 1a/19a, 2a/12a, 6a/14b, 7a/12c, 9/15b, and 10a/19b) were shared by the slow loris and galago. Five associations (HSA 1b/6b, 4a/5a, 11b/15a, 12b/19b, and 15b/16b) were unique to the slow loris. In contrast, 30 homologous chromosome segments were identified in the slow loris genome when using galago chromosome painting probes. The data showed that the karyotypic differences between these two species were mainly due to Robertsonian translocations. Reverse painting, using galago painting probes onto human chromosomes, confirmed most of the chromosome homologies between humans and galagos established previously, and documented the HSA 7/16 association in galagos, which was not reported previously. The presence of the HSA 7/16 association in the slow loris and galago suggests that the 7/16 association is an ancestral synteny for primates. Based on our results and the published homology maps between humans and other primate species, we propose an ancestral karyotype (2n = 60) for lorisiform primates.

Reconstructing the evolutionary history of the Lorisidae using morphological, molecular, and geological data.

Major aspects of lorisid phylogeny and systematics remain unresolved, despite several studies (involving morphology, histology, karyology, immunology, and DNA sequencing) aimed at elucidating them. Our study is the first to investigate the evolution of this enigmatic group using molecular and morphological data for all four well-established genera: Arctocebus, Loris, Nycticebus, and Perodicticus. Data sets consisting of 386 bp of 12S rRNA, 535 bp of 16S rRNA, and 36 craniodental characters were analyzed separately and in combination, using maximum parsimony and maximum likelihood. Outgroups, consisting of two galagid taxa (Otolemur and Galagoides) and a lemuroid (Microcebus), were also varied. The morphological data set yielded a paraphyletic lorisid clade with the robust Nycticebus and Perodicticus grouped as sister taxa, and the galagids allied with Arctocebus. All molecular analyses maximum parsimony (MP) or maximum likelihood (ML) which included Microcebus as an outgroup rendered a paraphyletic lorisid clade, with one exception: the 12S + 16S data set analyzed with ML. The position of the galagids in these paraphyletic topologies was inconsistent, however, and bootstrap values were low. Exclusion of Microcebus generated a monophyletic Lorisidae with Asian and African subclades; bootstrap values for all three clades in the total evidence tree were over 90%. We estimated mean genetic distances for lemuroids vs. lorisoids, lorisids vs. galagids, and Asian vs. African lorisids as a guide to relative divergence times. We present information regarding a temporary land bridge that linked the two now widely separated regions inhabited by lorisids that may explain their distribution. Finally, we make taxonomic recommendations based on our results.

Fluorescence in situ hybridization to chromosomes as a tool to understand human and primate genome evolution.

For the last 15 years molecular cytogenetic techniques have been extensively used to study primate evolution. Molecular probes were helpful to distinguish mammalian chromosomes and chromosome segments on the basis of their DNA content rather than solely on morphological features such as banding patterns. Various landmark rearrangements have been identified for most of the nodes in primate phylogeny while chromosome banding still provides helpful reference maps. Fluorescence in situ hybridization (FISH) techniques were used with probes of different complexity including chromosome painting probes, probes derived from chromosome sub-regions and in the size of a single gene. Since more recently, in silico techniques have been applied to trace down evolutionarily derived chromosome rearrangements by searching the human and mouse genome sequence databases. More detailed breakpoint analyses of chromosome rearrangements that occurred during higher primate evolution also gave some insights into the molecular changes in chromosome rearrangements that occurred in evolution. Hardly any "fusion genes" as known from chromosome rearrangements in cancer cells or dramatic "position effects" of genes transferred to new sites in primate genomes have been reported yet. Most breakpoint regions have been identified within gene poor areas rich in repetitive elements and/or low copy repeats (segmental duplications). The progress in various molecular and molecular-cytogenetic approaches including the recently launched chimpanzee genome project suggests that these new tools will have a significant impact on the further understanding of human genome evolution.

Odor familiarity and female preferences for males in a threatened primate, the pygmy loris Nycticebus pygmaeus: applications for genetic management of small populations.

Here we use sexual selection theory to develop a logistically simple, yet effective, method for the manipulation of female reproductive behavior for conservation goals. Mate choice leading to nonrandom mating patterns can exacerbate the loss of genetic diversity in small populations. On theoretical grounds, females should choose high-quality mates. A prediction stemming from chemical communication theory is that competitive males will be better able to saturate an area with scent marks. If this is true, females should mate preferentially with males whose odors they encounter most frequently. We tested this hypothesis with the pygmy loris, Nycticebus pygmaeus, a threatened and poorly studied nocturnal prosimian. For several weeks females were exposed repeatedly to the urine from a particular male, and were then allowed to choose between a male whose odors were familiar and one whose odors were novel. Females showed an unusually strong preference for the familiar-odor male, as indicated by several behavioral measures of mate preference. Conservation managers can use this method as a tool to obtain reproductive pairings that will maximize genetic compatibility and diversity. For example, unsuccessful males may be given the opportunity to reproduce. In captive populations, studbook managers often select pairs in order to optimize outbreeding, but these selected pairings may not coincide with the preferences of the individual animals involved. Although several authors have made theoretical arguments for manipulating mate choice for conservation, this is a novel test of a proximate mechanism that can be manipulated, cultivating applications rather than mere implications.

Molecular evolution of GH in primates: characterisation of the GH genes from slow loris and marmoset defines an episode of rapid evolutionary change.

Pituitary growth hormone (GH), like several other protein hormones, shows an unusual episodic pattern of molecular evolution in which sustained bursts of rapid change are imposed on long periods of very slow evolution (near-stasis). A marked period of rapid change occurred in the evolution of GH in primates or a primate ancestor, and gave rise to the species specificity that is characteristic of human GH. We have defined more precisely the position of this burst by cloning and sequencing the GH genes for a prosimian, the slow loris (Nycticebus pygmaeus) and a New World monkey, marmoset (Callithrix jacchus). Slow loris GH is very similar in sequence to pig GH, demonstrating that the period of rapid change occurred during primate evolution, after the separation of lines leading to prosimians and higher primates. The putative marmoset GH is similar in sequence to human GH, demonstrating that the accelerated evolution occurred before divergence of New World monkeys and Old World monkeys/apes. The burst of change was confined largely to coding sequence for mature GH, and is not marked in other components of the gene sequence including signal peptide, 5' upstream region and introns. A number of factors support the idea that this episode of rapid change was due to positive adaptive selection. Thus (1) there is no apparent loss of function of GH in man compared with non-primates, (2) after the episode of rapid change the rate of evolution fell towards the slow basal level that is seen for most mammalian GHs, (3) the accelerated rate of substitution for the exons of the GH gene significantly exceeds that for introns, and (4) the amino acids contributing to the hydrophobic core of GH are strongly conserved when higher primate and other GH sequences are compared, and for coding sequences other than that coding for hydrophobic core residues the rate of substitution for non-synonymous sites (K(A)) is significantly greater than that for synonymous sites (K(S)). In slow loris, as in most non-primate mammals, there is no evidence for duplication of the GH gene, but in marmoset, as in rhesus monkey and man, the putative GH gene is one of a cluster of closely related genes.

Failure of the ILD to determine data combinability for slow loris phylogeny.

Tests for incongruence as an indicator of among-data partition conflict have played an important role in conditional data combination. When such tests reveal significant incongruence, this has been interpreted as a rationale for not combining data into a single phylogenetic analysis. In this study of lorisiform phylogeny, we use the incongruence length difference (ILD) test to assess conflict among three independent data sets. A large morphological data set and two unlinked molecular data sets--the mitochondrial cytochrome b gene and the nuclear interphotoreceptor retinoid binding protein (exon 1)--are analyzed with various optimality criteria and weighting mechanisms to determine the phylogenetic relationships among slow lorises (Primates, Loridae). When analyzed separately, the morphological data show impressive statistical support for a monophyletic Loridae. Both molecular data sets resolve the Loridae as paraphyletic, though with different branching orders depending on the optimality criterion or character weighting used. When the three data partitions are analyzed in various combinations, an inverse relationship between congruence and phylogenetic accuracy is observed. Nearly all combined analyses that recover monophyly indicate strong data partition incongruence (P = 0.00005 in the most extreme case), whereas all analyses that recover paraphyly indicate lack of significant incongruence. Numerous lines of evidence verify that monophyly is the accurate phylogenetic result. Therefore, this study contributes to a growing body of information affirming that measures of incongruence should not be used as indicators of data set combinability.

Karyotype polymorphism in the slender loris (Loris tardigradus).

Three wild-caught Sri Lankan and two captive-bred slender lorises (Loris tardigradus) in Adelaide Zoo, South Australia, were karyotyped from leucocyte cultures. The founder male and one of the founder females were of the grey subspecies L. tardigradus grandis, of which L. tardigradus nordicus is probably a synonym. The other founder female was a small red-brown L. tardigradus tardigradus. All three wild-caught individuals were karyotypically distinct from each other and their offspring showed the expected heterozygosity. The L. tardigradus tardigradus karyotype was no more distinct from the grey subspecies karyotypes than they were from each other. On present evidence, Indian and Sri Lanka lorises appear to be chromosomally distinct.

Multiple L1 progenitors in prosimian primates: phylogenetic evidence from ORF1 sequences.

One of the uncertainties regarding the evolution of L1 elements is whether there are numerous progenitor genes. We present phylogenetic evidence from ORF1 sequences of slow loris (Nycticebus coucang) and galago (Galago crassicaudatus) that there were at least two distinct progenitors, active at the same time, in the ancestor of this family of prosimian primates. A maximum parsimony analysis that included representative L1s from human, rabbit, and rodents, along with the prosimian sequences, revealed that one of the galago L1s (Gc11) grouped very strongly with the slow loris sequences. The remaining galago elements formed their own unique and strongly supported clade. An analysis of replacement and silent site changes for each link of the most parsimonious tree indicated that during the descent of the Gc11 sequence approximately two times more synonymous than nonsynonymous substitutions had occurred, implying that the Gc11 founder was functional for some time after the split of galago and slow loris. Strong purifying selection was also evident on the galago branch of the tree. These data indicate that there were two distinct and contemporaneous L1 progenitors in the lorisoid ancestor, evolving under purifying selection, that were retained as functional L1s in the galago lineage (and presumably also in the slow loris). The prosimian ORF1 sequences could be further subdivided into subfamilies. ORF1 sequences from both the galago and slow loris have a premature termination codon near the 3' end, not shared by the other mammalian sequences, that shortens the open reading frame by 288 bp. An analysis of synonymous and nonsynonymous substitutions for the 5' and 3' portions, that included intra- and inter-subfamily comparisons, as well as comparisons among the other mammalian sequences, suggested that this premature stop codon is a prosimian acquisition that has rendered the 3' portion of ORF1 in these primates noncoding.