Mitogenomic phylogeny of Callithrix with special focus on human transferred taxa.

BACKGROUND: Callithrix marmosets are a relatively young primate radiation, whose phylogeny is not yet fully resolved. These primates are naturally para- and allopatric, but three species with highly invasive potential have been introduced into the southeastern Brazilian Atlantic Forest by the pet trade. There, these species hybridize with each other and endangered, native congeners. We aimed here to reconstruct a robust Callithrix phylogeny and divergence time estimates, and identify the biogeographic origins of autochthonous and allochthonous Callithrix mitogenome lineages. We sequenced 49 mitogenomes from four species (C. aurita, C. geoffroyi, C. jacchus, C. penicillata) and anthropogenic hybrids (C. aurita x Callithrix sp., C. penicillata x C. jacchus, Callithrix sp. x Callithrix sp., C. penicillata x C. geoffroyi) via Sanger and whole genome sequencing. We combined these data with previously published Callithrix mitogenomes to analyze five Callithrix species in total. RESULTS: We report the complete sequence and organization of the C. aurita mitogenome. Phylogenetic analyses showed that C. aurita was the first to diverge within Callithrix 3.54 million years ago (Ma), while C. jacchus and C. penicillata lineages diverged most recently 0.5 Ma as sister clades. MtDNA clades of C. aurita, C. geoffroyi, and C. penicillata show intraspecific geographic structure, but C. penicillata clades appear polyphyletic. Hybrids, which were identified by phenotype, possessed mainly C. penicillata or C. jacchus mtDNA haplotypes. The biogeographic origins of mtDNA haplotypes from hybrid and allochthonous Callithrix were broadly distributed across natural Callithrix ranges. Our phylogenetic results also evidence introgression of C. jacchus mtDNA into C. aurita. CONCLUSION: Our robust Callithrix mitogenome phylogeny shows C. aurita lineages as basal and C. jacchus lineages among the most recent within Callithrix. We provide the first evidence that parental mtDNA lineages of anthropogenic hybrid and allochthonous marmosets are broadly distributed inside and outside of the Atlantic Forest. We also show evidence of cryptic hybridization between allochthonous Callithrix and autochthonous C. aurita. Our results encouragingly show that further development of genomic resources will allow to more clearly elucidate Callithrix evolutionary relationships and understand the dynamics of Callithrix anthropogenic introductions into the Brazilian Atlantic Forest.

Unraveling the dispersal patterns and the social drivers of natal emigration of a cooperative breeding mammal, the golden lion tamarin.

The study of the social drivers of animal dispersal is key to understanding the evolution of social systems. Among the social drivers of natal emigration, the conspecific attraction, aggressive eviction, and reduced social integration hypotheses predict that sexually mature individuals who receive more aggressive behavior and are engaged in less affiliative interactions are more likely to disperse. Few reports have explored these proximate factors affecting emigration in cooperatively breeding species, particularly of Neotropical primates. In this study, we investigated the dispersal patterns and tested the social drivers of natal emigration in the golden lion tamarin (Leontopithecus rosalia) - an endangered species inhabiting Atlantic rainforests fragments in Brazil. We used behavioral and demographic data collected during 7 years from 68 groups of tamarins inhabiting 20 forest fragments. Our analyses from the 160 dispersing individuals showed that dispersal success is higher for males and for those engaged in parallel dispersal, but that males and females use different strategies to enhance their dispersal success, males immigrate into established groups while females form new groups. We did not find high levels of agonistic behavior among group members before natal emigration. Instead we found that conspecific attraction drives natal emigration in both sexes, while additionally the low level of affiliative interactions within the natal group triggers male emigration. We discuss natal emigration in the broader perspective of the cooperative breeding system and the implications of these findings for the conservation of the species.

Population Genetic Structure of an Endangered Endemic Primate (Leontopithecus chrysomelas) in a Highly Fragmented Atlantic Coastal Rain Forest.

This study evaluated the genetic structure of wild populations of the endangered primate, Leontopithecus chrysomelas. We tested the assumption that populations of L. chrysomelas, given their larger population size and a higher degree of habitat continuity, would have higher genetic diversity and less genetic structuring than other lion tamarins. We used 11 microsatellites and 122 hair samples from different locations to assess their genetic diversity and genetic structure, and to make inferences about the isolation by distance. The overall expected heterozygosity (0.51 ± 0.03) and the average number of alleles (3.6 ± 0.2) were relatively low, as is the case in other endangered lion tamarins. Genetic clustering analyses indicated two main clusters, whereas the statistical analyses based on genotype similarities and Fst suggested further substructure. A Mantel test showed that only 34% of this genetic differentiation was explained by the linear distance. In addition to linear distance, structural differences in the landscape, physical barriers and behavioural factors may be causing significant genetic structuring. Overall, this study suggests that these populations have a relatively low genetic diversity and a relatively high population genetic structure, putting in question whether the presence of agroforest systems (known locally as cabruca) is enough to fully re-establish functional landscape connectivity.

Application of PE-RADSeq to the study of genomic diversity and divergence of two Brazilian marmoset species (Callithrix jacchus and C. penicillata).

Callithrix jacchus and C. penicillata are among the smallest anthropoid primates, are highly specialized tree gougers, and largely occupy Brazil's most extreme, semi-arid biomes. However, the underlying genomic factors that underpin the evolution of these species and their unique traits are under-investigated. Additionally, exotic populations of these two species are widely established throughout Brazil and hybridize with threatened native congers. Thus, both genomic and conservation factors call for a better understanding of C. jacchus and C. penicillata evolution. Here, we applied PE-RADseq to characterize genomic variation in these two species, using six C. jacchus and seven C. penicillata individuals. We identified an average of 7,463 and 5,180 SNPs/individual in C. penicillata and C. jacchus, respectively, and also found 1,395 variable sites that were represented in both species. C. penicillata showed overall higher levels of genetic diversity than C. jacchus at the variable sites present in both species. Additionally, among these variable sites, 106 showed relative interspecific divergence levels that were significantly higher than the genome-wide average. We further compared relative and absolute divergence for C. penicillata and C. jacchus between RAD loci associated with the 106 significantly diverged variable sites and all other RAD loci present in both species. The former RAD loci set showed significantly elevated relative and absolute divergence measures in comparison to the latter set. This convergence suggests that C. jacchus and C. penicillata may have diverged under a scenario of gene flow under secondary contact. Here, we demonstrate that RADseq is an efficient method to simultaneously discover and genotype a large number of markers and validate the utility of RADseq for examining Callithrix evolution.

Hybridization effects and genetic diversity of the common and black-tufted marmoset (Callithrix jacchus and Callithrix penicillata) mitochondrial control region.

Hybridization is continually documented in primates, but effects of natural and anthropogenic hybridization on biodiversity are still unclear and differentiating between these contexts remains challenging in regards to primate evolution and conservation. Here, we examine hybridization effects on the mitochondrial DNA (mtDNA) control region of Callithrix marmosets, which provide a unique glimpse into interspecific mating under distinct anthropogenic and natural conditions. DNA was sampled from 40 marmosets along a 50-km transect from a previously uncharacterized hybrid zone in NE Brazil between the ranges of Callithrix jacchus and Callithrix penicillata. DNA was also collected from 46 marmosets along a 30-km transect in a hybrid zone in Rio de Janeiro state, Brazil, where exotic marmosets appeared in the 1980s. Combining Callithrix DNA sampled inside and outside of these hybrid zones, phylogenetic and network analyses show C. jacchus and C. penicillata being parental species to sampled hybrids. We expand limited Callithrix population genetics work by describing mtDNA diversity and demographic history of these parental species. We show ancient population expansion in C. jacchus and historically constant population size in C. penicillata, with the latter being more genetically diverse than the former. The natural hybrid zone contained higher genetic diversity relative to the anthropogenic zone. While our data suggest hybrid swarm formation within the anthropogenic zone due to removed physical reproductive barriers, this pattern is not seen in the natural hybrid zone. These results suggest different genetic dynamics within natural and anthropogenic hybridization contexts that carry important implications for primate evolution and conservation.

Bifidobacteria define gut microbiome profiles of golden lion tamarin (Leontopithecus rosalia) and marmoset (Callithrix sp.) metagenomic shotgun pools.

Gut microbiome disruptions may lead to adverse effects on wildlife fitness and viability, thus maintaining host microbiota biodiversity needs to become an integral part of wildlife conservation. The highly-endangered callitrichid golden lion tamarin (GLT-Leontopithecus rosalia) is a rare conservation success, but allochthonous callitrichid marmosets (Callithrix) serve as principle ecological GLT threats. However, incorporation of microbiome approaches to GLT conservation is impeded by limited gut microbiome studies of Brazilian primates. Here, we carried out analysis of gut metagenomic pools from 114 individuals of wild and captive GLTs and marmosets. More specifically, we analyzed the bacterial component of ultra filtered samples originally collected as part of a virome profiling study. The major findings of this study are consistent with previous studies in showing that Bifidobacterium, a bacterial species important for the metabolism of tree gums consumed by callitrichids, is an important component of the callitrichid gut microbiome - although GTLs and marmosets were enriched for different species of Bifidobacterium. Additionally, the composition of GLT and marmoset gut microbiota is sensitive to host environmental factors. Overall, our data expand baseline gut microbiome data for callitrichids to allow for the development of new tools to improve their management and conservation.

The gut microbiome of exudivorous marmosets in the wild and captivity.

Mammalian captive dietary specialists like folivores are prone to gastrointestinal distress and primate dietary specialists suffer the greatest gut microbiome diversity losses in captivity compared to the wild. Marmosets represent another group of dietary specialists, exudivores that eat plant exudates, but whose microbiome remains relatively less studied. The common occurrence of gastrointestinal distress in captive marmosets prompted us to study the Callithrix gut microbiome composition and predictive function through bacterial 16S ribosomal RNA V4 region sequencing. We sampled 59 wild and captive Callithrix across four species and their hybrids. Host environment had a stronger effect on the gut microbiome than host taxon. Wild Callithrix gut microbiomes were enriched for Bifidobacterium, which process host-indigestible carbohydrates. Captive marmoset guts were enriched for Enterobacteriaceae, a family containing pathogenic bacteria. While gut microbiome function was similar across marmosets, Enterobacteriaceae seem to carry out most functional activities in captive host guts. More diverse bacterial taxa seem to perform gut functions in wild marmosets, with Bifidobacterium being important for carbohydrate metabolism. Captive marmosets showed gut microbiome composition aspects seen in human gastrointestinal diseases. Thus, captivity may perturb the exudivore gut microbiome, which raises implications for captive exudivore welfare and calls for husbandry modifications.

The role of food transfers in wild golden lion tamarins (Leontopithecus rosalia): Support for the informational and nutritional hypothesis.

Callitrichidae is a unique primate family not only in terms of the large number of food transfers to infants but also for the prevalence of transfers that are initiated by the adults. It has been hypothesized that, as well as provisioning infants, callitrichid food transfers might function to teach the receiver what food types to eat. If food provisioning has a teaching function, we would expect successful food transfers to be more likely with food types that are novel to the juveniles. We would also expect juveniles to learn about foods from those transfers. We introduced different types of food (some familiar, some novel) to wild groups of golden lion tamarins (Leontopithecus rosalia). While novel foods were not more successfully transferred than familiar food in the experiment, transfers were more successful (i.e., the receiver obtained food) when the donor had previous experience with that food. Moreover, we found evidence suggesting that food transfers influenced the future foraging choices of juveniles. Our findings are consistent with the first and third criteria of the functional definition of teaching, which requires that tutors (the adults) modify their behavior in the presence of a naïve individual (a juvenile), and that the naïve individual learns from the modified behavior of the demonstrator. Our findings are also consistent with the provisioning function of food transfer. Social learning seems to play an important role in the development of young tamarins' foraging preferences.

Genomic skimming and nanopore sequencing uncover cryptic hybridization in one of world's most threatened primates.

The Brazilian buffy-tufted-ear marmoset (Callithrix aurita), one of the world's most endangered primates, is threatened by anthropogenic hybridization with exotic, invasive marmoset species. As there are few genetic data available for C. aurita, we developed a PCR-free protocol with minimal technical requirements to rapidly generate genomic data with genomic skimming and portable nanopore sequencing. With this direct DNA sequencing approach, we successfully determined the complete mitogenome of a marmoset that we initially identified as C. aurita. The obtained nanopore-assembled sequence was highly concordant with a Sanger sequenced version of the same mitogenome. Phylogenetic analyses unexpectedly revealed that our specimen was a cryptic hybrid, with a C. aurita phenotype and C. penicillata mitogenome lineage. We also used publicly available mitogenome data to determine diversity estimates for C. aurita and three other marmoset species. Mitogenomics holds great potential to address deficiencies in genomic data for endangered, non-model species such as C. aurita. However, we discuss why mitogenomic approaches should be used in conjunction with other data for marmoset species identification. Finally, we discuss the utility and implications of our results and genomic skimming/nanopore approach for conservation and evolutionary studies of C. aurita and other marmosets.

An Introduction to the Callithrix Genus and Overview of Recent Advances in Marmoset Research.

We provide here a current overview of marmoset (Callithrix) evolution, hybridization, species biology, basic/biomedical research, and conservation initiatives. Composed of 2 subgroups, the aurita group (C aurita and C flaviceps) and the jacchus group (C geoffroyi, C jacchus, C kuhlii, and C penicillata), this relatively young primate radiation is endemic to the Brazilian Cerrado, Caatinga, and Atlantic Forest biomes. Significant impacts on Callithrix within these biomes resulting from anthropogenic activity include (1) population declines, particularly for the aurita group; (2) widespread geographic displacement, biological invasions, and range expansions of C jacchus and C penicillata; (3) anthropogenic hybridization; and (4) epizootic Yellow Fever and Zika viral outbreaks. A number of Brazilian legal and conservation initiatives are now in place to protect the threatened aurita group and increase research about them. Due to their small size and rapid life history, marmosets are prized biomedical models. As a result, there are increasingly sophisticated genomic Callithrix resources available and burgeoning marmoset functional, immuno-, and epigenomic research. In both the laboratory and the wild, marmosets have given us insight into cognition, social group dynamics, human disease, and pregnancy. Callithrix jacchus and C penicillata are emerging neotropical primate models for arbovirus disease, including Dengue and Zika. Wild marmoset populations are helping us understand sylvatic transmission and human spillover of Zika and Yellow Fever viruses. All of these factors are positioning marmosets as preeminent models to facilitate understanding of facets of evolution, hybridization, conservation, human disease, and emerging infectious diseases.

Yellow fever in Brazil threatens successful recovery of endangered golden lion tamarins.

The golden lion tamarin is an endangered primate endemic to Brazil's Atlantic Forest. Centuries of deforestation reduced numbers to a few hundred individuals in isolated forest fragments 80 km from Rio de Janeiro city. Intensive conservation action including reintroduction of zoo-born tamarins into forest fragments 1984-2000, increased numbers to about 3,700 in 2014. Beginning in November 2016, southeastern Brazil experienced the most severe yellow fever epidemic/epizootic in the country in 80 years. In May 2018, we documented the first death of a golden lion tamarin due to yellow fever. We re-evaluated population sizes and compared them to results of a census completed in 2014. Tamarin numbers declined 32%, with ca. 2,516 individuals remaining in situ. Tamarin losses were significantly greater in forest fragments that were larger, had less forest edge and had better forest connectivity, factors that may favor the mosquito vectors of yellow fever. The future of golden lion tamarins depends on the extent of additional mortality, whether some tamarins survive the disease and acquire immunity, and the potential development of a vaccine to protect the species against yellow fever.

Estimating population sizes to evaluate progress in conservation of endangered golden lion tamarins (Leontopithecus rosalia).

Efforts to reverse the decline of endangered golden lion tamarin monkeys have been relatively successful because the Brazilian organization dedicated to the species' conservation (Associação Mico-Leão-Dourado, AMLD) relies on science-based computer modeling to determine the number of tamarins necessary to achieve demographic and genetic goals, and a process of strategic planning to achieve those goals. Accurate estimates of the numbers of tamarins in forest fragments are essential to evaluate progress in achieving goals and adapt strategies as necessary. In this report we present the results of a new method to survey the number of tamarins in the wild, a modification of the point transect with lures procedure. Using this method, we estimate that in 2014 there were approximately 3,700 golden lion tamarins in 41,400 hectares of Atlantic Forest. Of these, 59% are from remnant wild populations, 34% are descendants of captive-born reintroduced animals and 7% are descendants of wild translocated groups. The number of tamarins and amount of forest estimated in this survey exceeded values necessary to meet AMLD's definition of a viable population, determined to be 2,000 tamarins in 25,000 hectares of connected and protected forest. However, the seven forest blocks and their tamarin populations are not yet adequately connected and protected. AMLD's strategic plan to achieve a viable population of golden lion tamarins includes 12 strategies that mitigate these and other threats or contribute directly to the conservation goal. The point transect with lures survey method provides a way to evaluate progress in achieving that goal and adapt strategies as appropriate.

Food-Offering Calls in Wild Golden Lion Tamarins (Leontopithecus rosalia): Evidence for Teaching Behavior?

Many animals emit calls in the presence of food, but researchers do not always know the function of these calls. Evidence suggests that adult golden lion tamarins (Leontopithecus rosalia) use food-offering calls to teach juveniles which substrate (i.e., microhabitat) to forage on, or in, for food. However, we do not yet know whether juveniles learn from this aspect of the adults' behavior. Here we examine whether juveniles learn to associate food-offering calls with a foraging substrate, as a step toward assessing whether these calls qualify as teaching behavior. We compared the performance of four wild juvenile golden lion tamarins that were introduced to a novel substrate while exposed to playbacks of food-offering calls (experimental condition) to the performance of three juveniles that were exposed to the novel substrate without the presence of food-offering playbacks (control condition). We varied the location of the novel substrate between trials. We found that food-offering calls had an immediate effect on juveniles' interactions with the novel substrate, whether they inserted their hands into the substrate and their eating behavior, and a long-term effect on eating behavior at the substrate. The findings imply that juvenile golden lion tamarins can learn through food-offering calls about the availability of food at a substrate, which is consistent with (but does not prove) teaching in golden lion tamarins through stimulus enhancement. Our findings support the hypothesis that teaching might be more likely to evolve in cooperatively breeding species with complex ecological niches.

Correction: Validation of qPCR Methods for the Detection of Mycobacterium in New World Animal Reservoirs.

[This corrects the article DOI: 10.1371/journal.pntd.0004198.].

Validation of qPCR Methods for the Detection of Mycobacterium in New World Animal Reservoirs.

Zoonotic pathogens that cause leprosy (Mycobacterium leprae) and tuberculosis (Mycobacterium tuberculosis complex, MTBC) continue to impact modern human populations. Therefore, methods able to survey mycobacterial infection in potential animal hosts are necessary for proper evaluation of human exposure threats. Here we tested for mycobacterial-specific single- and multi-copy loci using qPCR. In a trial study in which armadillos were artificially infected with M. leprae, these techniques were specific and sensitive to pathogen detection, while more traditional ELISAs were only specific. These assays were then employed in a case study to detect M. leprae as well as MTBC in wild marmosets. All marmosets were negative for M. leprae DNA, but 14 were positive for the mycobacterial rpoB gene assay. Targeted capture and sequencing of rpoB and other MTBC genes validated the presence of mycobacterial DNA in these samples and revealed that qPCR is useful for identifying mycobacterial-infected animal hosts.

Natural and Anthropogenic Hybridization in Two Species of Eastern Brazilian Marmosets (Callithrix jacchus and C. penicillata).

Animal hybridization is well documented, but evolutionary outcomes and conservation priorities often differ for natural and anthropogenic hybrids. Among primates, an order with many endangered species, the two contexts can be hard to disentangle from one another, which carries important conservation implications. Callithrix marmosets give us a unique glimpse of genetic hybridization effects under distinct natural and human-induced contexts. Here, we use a 44 autosomal microsatellite marker panel to examine genome-wide admixture levels and introgression at a natural C. jacchus and C. penicillata species border along the São Francisco River in NE Brazil and in an area of Rio de Janeiro state where humans introduced these species exotically. Additionally, we describe for the first time autosomal genetic diversity in wild C. penicillata and expand previous C. jacchus genetic data. We characterize admixture within the natural zone as bimodal where hybrid ancestry is biased toward one parental species or the other. We also show evidence that São Francisco River islands are gateways for bidirectional gene flow across the species border. In the anthropogenic zone, marmosets essentially form a hybrid swarm with intermediate levels of admixture, likely from the absence of strong physical barriers to interspecific breeding. Our data show that while hybridization can occur naturally, the presence of physical, even if leaky, barriers to hybridization is important for maintaining species genetic integrity. Thus, we suggest further study of hybridization under different contexts to set well informed conservation guidelines for hybrid populations that often fit somewhere between "natural" and "man-made."