Genetic structure in the southernmost populations of black-and-gold howler monkeys (Alouatta caraya) and its conservation implications

Phylogeography, taxonomy, and conservation of the endangered brown howler monkey, Alouatta guariba (Primates, Atelidae), of the Atlantic Forest

The brown howler, Alouatta guariba, endemic to the Atlantic Forest of Brazil and Argentina, is threatened by habitat loss and fragmentation, hunting, and its susceptibility to yellow fever. Two subspecies have been recognized, but their names, validity, and geographic ranges have been controversial. We obtained samples covering the species' entire distribution in Brazil and Argentina to clarify these issues by investigating their genetic diversity and structure and assessing their evolutionary history. We analyzed, for the first time, a set of ten microsatellite markers (N = 153), plus mitochondrial DNA (mtDNA) segments of the control region (N = 207) and cytochrome b gene (N = 116). The microsatellite data support two to three genetic clusters with biological significance. The southern populations (Argentina, Santa Catarina, and Rio Grande do Sul) presented a homogeneous genetic component, and populations from São Paulo (SP) to the north presented another component, although most presented ∼20% of the southern component. With K = 3, SP emerged as a third component while sharing some ancestry with Rio de Janeiro and Argentina. The mtDNA phylogenies revealed three main clades that diverged almost simultaneously around 250 thousand years ago (kya). Clades A and B are from central SP to the north and east, while clade C is from SP to the south and southwest. Samples from SP presented haplotypes in all three clades, sometimes in the same population. The demographic history of the species estimated with the Bayesian skyline plot of the mtDNA showed a strong expansion ∼40-20 kya and a strong reduction over the last ∼4-2 kya. Although the genetic clusters identified here deserve appropriate management strategies as conservation units, the absence of (i) concordance between the mtDNA and microsatellite data, (ii) reciprocal monophyly in the mtDNA, and (iii) clear-cut non-genetic diagnostic characters advises against considering them as different taxonomic entities. None of the previous taxonomic proposals were corroborated by our data. Our results elucidate the taxonomy of the Atlantic Forest brown howler, indicating it should be considered a monotypic species, A. guariba. We also clarify the evolutionary history of the species regarding its intraspecific genetic diversity, which is crucial information for its conservation and population management.

Human and Non-Human Primate Coexistence in Argentina: Conflicts and Solutions

There are five different primate species inhabiting widely distinct ecoregions in Argentina. Each of them faces various threats in terms of conservation and conflicts that hamper their ability to coexist with human populations. We present here some of the drivers known to be the causes of conflicts between humans and primates in the southernmost area of distribution of Latin American primates. We focus our synthesis on two of the biggest sources of conflict: the effects of different anthropogenic disturbances, and human misconceptions concerning the role of primates in the ecosystem. In each section, we briefly characterize the conflicts worldwide and then provide specific cases and examples from Argentina. In the last part of the manuscript, we further describe some ongoing national and regional educational, research, and conservation approaches to mitigate those effects.

Influence of landscape structure on previous exposure to Leptospira spp. and Brucella abortus in free-living neotropical primates from southern Brazil

The environments in which neotropical primates live have been undergoing an intense fragmentation process, constituting a major threat to the species' survival and causing resource scarcity, social isolation, and difficulty in dispersal, leaving populations increasingly vulnerable. Moreover, the proximity of wild environments to anthropized landscapes can change the dynamics of pathogens and the parasite-host-environment relationship, creating conditions that favor exposure to different pathogens. To investigate the previous exposure of free-living primates in Rio Grande do Sul State (RS), southern Brazil, to the bacterial agents Leptospira spp. and Brucella abortus, we investigated agglutinating antibodies against 23 serovars of Leptospira spp. using the microscopic agglutination test and B. abortus acidified antigen test in primate serum samples; 101 samples from primates captured between 2002 and 2016 in different forest fragments were used: 63 Alouatta caraya, 36 Alouatta guariba clamitans, and 02 Sapajus nigritus cucullatus. In addition, the forest remnants where the primates were sampled were characterized in a multiscale approach in radii ranging from 200 to 1400 m to investigate the potential relationship of previous exposure to the agent with the elements that make up the landscape structure. The serological investigation indicated the presence of antibodies for at least one of the 23 serovars of Leptospira spp. in 36.6% (37/101) of the samples analyzed, with titers ranging from 100 to 1600. The most observed serovars were Panama (17.8%), Ballum (5.9%), Butembo (5.9%), Canicola (5.9%), Hardjo (4.9%), and Tarassovi (3.9%); no samples were seropositive for Brucella abortus. Decreased forest cover and edge density were the landscape factors that had a significant relationship with Leptospira spp. exposure, indicating that habitat fragmentation may influence contact with the pathogen. The data generated in this study demonstrate the importance of understanding how changes in landscape structure affect exposure to pathogenic microorganisms of zoonotic relevance. Hence, improving epidemiological research and understanding primates' ecological role in these settings can help improve environmental surveillance and conservation strategies for primate populations in different landscapes.

Terrain Ruggedness and Canopy Height Predict Short-Range Dispersal in the Critically Endangered Black-and-White Ruffed Lemur

Dispersal is a fundamental aspect of primates' lives and influences both population and community structuring, as well as species evolution. Primates disperse within an environmental context, where both local and intervening environmental factors affect all phases of dispersal. To date, research has primarily focused on how the intervening landscape influences primate dispersal, with few assessing the effects of local habitat characteristics. Here, we use a landscape genetics approach to examine between- and within-site environmental drivers of short-range black-and-white ruffed lemur (Varecia variegata) dispersal in the Ranomafana region of southeastern Madagascar. We identified the most influential drivers of short-range ruffed lemur dispersal as being between-site terrain ruggedness and canopy height, more so than any within-site habitat characteristic evaluated. Our results suggest that ruffed lemurs disperse through the least rugged terrain that enables them to remain within their preferred tall-canopied forest habitat. Furthermore, we noted a scale-dependent environmental effect when comparing our results to earlier landscape characteristics identified as driving long-range ruffed lemur dispersal. We found that forest structure drives short-range dispersal events, whereas forest presence facilitates long-range dispersal and multigenerational gene flow. Together, our findings highlight the importance of retaining high-quality forests and forest continuity to facilitate dispersal and maintain functional connectivity in ruffed lemurs.

Future simulated landscape predicts habitat loss for the Golden Langur (Trachypithecus geei): A range level analysis for an endangered primate

Trachypithecus geei Khajuria, 1956 or Golden langur are an endangered colobine primate species that are distributed in the transboundary region of Bhutan and India. The species is severely threatened because of increasing habitat fragmentation and isolation across its entire range, especially the populations in Assam, India. The distribution range of the species has not been updated for the last two decades, nor is there any proper evaluation of the habitat requirements for the species. Therefore, we mapped the habitat suitability for the species across its entire distribution and projected its habitat suitability on the simulated landscape for the future (2031). The results indicate that out of the total range extent (66,320 km²), only 12,265 km² (18.49%) is suitable for the species at present, which will further be reduced to 8884 km² by the year 2031, indicating major range contraction. These suitable habitats are largely scattered and fragmented in southern range of the species. Among the predictors used, the distance to evergreen and deciduous broadleaf forest was the strongest predictor out of the 35 used for model building. Moreover, land use and land cover were found to be more informative than the climatic variables. Much of the suitable habitats of the species are located outside the protected area network in the landscape. Therefore, we identified landscape configurations and suitable habitat areas for the future conservation and monitoring of Golden Langur in the protected areas of its range.

Development of 16 novel EST-SSR markers for species identification and cross-genus amplification in sambar, sika, and red deer

Deer genera around the globe are threatened by anthropogenic interference. The translocation of alien species and their subsequent genetic introgression into indigenous deer populations is particularly harmful to the species of greatest conservation concern. Products derived from deer, including venison and antler velvet, are also at risk of fraudulent labeling. The current molecular markers used to genetically identify deer species were developed from genome sequences and have limited applicability for cross-species amplification. The absence of efficacious diagnostic techniques for identifying deer species has hampered conservation and wildlife crime investigation efforts. Expressed sequence tag-simple sequence repeat (EST-SSR) markers are reliable tools for individual and species identification, especially in terms of cross-species genotyping. We conducted transcriptome sequencing of sambar (Rusa unicolor) antler velvet and acquired 11,190 EST-SSRs from 65,074 newly assembled unigenes. We identified a total of 55 unambiguous amplicons in sambar (n = 45), which were selected as markers to evaluate cross-species genotyping in sika deer (Cervus nippon, n = 30) and red deer (Cervus elaphus, n = 46), resulting in cross-species amplification rates of 94.5% and 89.1%, respectively. Based on polymorphic information content (>0.25) and genotyping fidelity, we selected 16 of these EST-SSRs for species identification. This marker set revealed significant genetic differentiation based on the fixation index and genetic distance values. Principal coordinate analysis and STRUCTURE analysis revealed distinct clusters of species and clearly identified red-sika hybrids. These markers showed applicability across different genera and proved suitable for identification and phylogenetic analyses across deer species.

Gone With the Water: The Loss of Genetic Variability in Black and Gold Howler Monkeys (Alouatta caraya) Due to Dam Construction

Black and gold howler monkeys (Alouatta caraya) inhabit several eco-regions in South America with the highest population densities in riverine forests. Dam construction for electricity production represents a severe human alteration of ecosystems with consequences for primate conservation. To evaluate the possible loss of genetic diversity in A. caraya, we analysed and compared the genetic structure of the species across 22 study sites in Argentina (14), Paraguay (1), and Brazil (7). Four of these study sites (referred to as flooded) were sampled before dam-linked flooding which most likely caused a drastic decline or functional extinction of these populations. The genetic variability of 256 individuals was evaluated using 10 autosomal microsatellites (STRs) and 112 individuals by sequencing a fragment of 507 bp of mtDNA. DNA was extracted from tissue, blood, and faecal samples. Significantly higher values of genetic variability were observed for the flooded populations both in mtDNA and STRs. Population genetic structure showed a K = 1, 2, or 5 depending on the method, separating Argentinian and Paraguayan sites from Brazilian sites and, in the case of K = 5, two clusters were mostly represented by flooded populations. Isolation-by-distance analyses showed that geographic distances influence gene flow. Analytical methods, such as Pairwise Fst’s and Nei’s and regression model of Harpending and Ward, were concordant in detecting significant genetic structuring between flooded and remaining sites examined. Although some sites have very low sample sizes, these samples are of great importance since these sampling sites are currently flooded. Our results show that the study sites where dams were built had the greatest genetic diversity. As A. caraya is currently severely threatened by yellow fever outbreaks, the remaining populations may be more vulnerable to disease outbreaks due to impoverished genetic variability. Accordingly, it is essential to implement management actions to conserve the remaining populations. Our results underline the importance for Environmental Impact Assessments (EIA) to include data on the genetic structure of species in the affected sites prior to their alteration or destruction. These genetic data are also remarkably important for determining where to relocate specific individuals to help avoid biodiversity loss.

Potential drivers of samango monkey (Cercopithecus albogularis) population subdivision in a highly fragmented mountain landscape in northern South Africa

Forests affected by fragmentation are at risk of losing their primate populations over the long term. The impact of fragmentation on primate populations has been studied in several places in Africa, Asia and South America; however, there has been no discernible pattern of how primates react to forest disturbance and fragmentation. In fragmented habitats, the local extinction probability of a species increases due to a decrease in patch area and an increase in genetic isolation. Here we used microsatellite markers and mitochondrial DNA sequences to investigate how habitat fragmentation impacts on the genetic diversity and structure of a samango monkey population inhabiting forest patches in the Soutpansberg mountain range of northern South Africa. We sampled four local populations across the length of the mountain range and an additional outlying population from the Great Escarpment to the south. Our results indicate that local populations along the mountain range were historically more connected and less distinct than at present. In more recent times, a lack of contemporary gene flow is leading to a more pronounced genetic structure, causing population subdivision across the mountain and likely isolating the Soutpansberg population from the escarpment population to the south. Based on our results, we suggest that natural and anthropogenic fragmentation are driving population genetic differentiation, and that the matrix surrounding forests and their suitability for samango monkey utilisation play a role at the local scale. The degree of genetic isolation found for samango monkey populations in our study raises concerns about the long-term viability of populations across the mountain range.

Molecular characterization of Giardia duodenalis and evidence for cross-species transmission in Northern Argentina

Anthropogenic activities, such as human population expansion and land-use change, create ecological overlap between humans, domesticated animals, and wildlife and can exacerbate the zoonotic transmission of parasites. To improve our understanding of this dynamic, we employed multi-locus genotyping to conduct a cross-sectional study of the potential for zoonotic transmission of the protozoan parasite Giardia duodenalis among humans, household associated livestock and dogs, and black and gold howler monkeys (Alouatta caraya) in the Corrientes Province of Argentina. We found Giardia prevalence to be highest in howler monkeys (90.3% (47/52)), followed by humans (61.1% (22/36)), dogs (44.4% (16/36)), and cattle (41.9% (18/43)). We further established that howler monkeys exclusively harbored strains of assemblage B (100%) while humans were infected with either assemblage A (13.3%) or B (80%) or A and B (6.7%), and cattle and dogs were infected with either assemblage A (cattle, 94.1%; dogs, 80%)), A and C (10%), or their host-adapted assemblage (cattle, 5.9%; dogs, 10%). Our finding of G. duodenalis in both humans and domesticated animals (assemblage A) and humans and wild primates (assemblage B) suggests that cross-species transmission of multiple assemblages of G. duodenalis may occur in rural complexes such as northern Argentina where people, domesticated animals, and wildlife overlap. We further highlight the need to investigate the implications of these results for human health, the economics of livestock production, and wildlife conservation in this and similar systems.

Primate landscape genetics: A review and practical guide

Landscape genetics is an emerging field that integrates population genetics, landscape ecology, and spatial statistics to investigate how geographical and environmental features and evolutionary processes such as gene flow, genetic drift, and selection structure genetic variation at both the population and individual levels, with implications for ecology, evolution, and conservation biology. Despite being particularly well suited for primatologists, this method is currently underutilized. Here, we synthesize the current state of research on landscape genetics in primates. We begin by outlining how landscape genetics has been used to disentangle the drivers of diversity, followed by a review of how landscape genetic methods have been applied to primates. This is followed by a section highlighting special considerations when applying the methods to primates, and a practical guide to facilitate further landscape genetics studies using both existing and de novo datasets. We conclude by exploring future avenues of inquiry that could be facilitated by recent developments as well as underdeveloped applications of landscape genetics to primates.

Use of genetic tools to assess predation on reintroduced howler monkeys (Alouatta caraya) in Northeastern Argentina

Despite strong support from the media, the reintroduction of animals into natural environments does not always achieve its goal. Alouatta caraya is the primate species facing the greatest hunting pressure due to the illegal pet trade in Argentina. Confiscations of this species are common, as is the voluntary surrender of animals by owners no longer able or willing to care for them. These animals ultimately arrive at rehabilitation centers and, in many cases, are released into natural environments that may differ from the original sites where they were captured. Until recently, the lack of genetic analysis of the individuals involved led to biased relocation decisions. We followed the reintroduction of 12 A. caraya individuals in a protected area (Isla Palacio, Misiones, Argentina). The presence of potential predators such as pumas (Puma concolor) and jaguars (Panthera onca) in this area was confirmed by camera traps, footprints and feces. After the disappearance of four A. caraya at the reintroduction site, we investigated the applicability of genetic assignment tests based on genotypic data to accurately identify predated individuals. Genetic analyses allowed us to determine the predator species (P. onca) and to identify the predated individuals as two of the reintroduced animals. This procedure is promising for identifying the remains of predated individuals, and can contribute to the design of reintroduction policies based on scientific evidence.

Effects of anthropogenic habitat disturbance and Giardia duodenalis infection on a sentinel species' gut bacteria

Habitat disturbance, a common consequence of anthropogenic land use practices, creates human-animal interfaces where humans, wildlife, and domestic species can interact. These altered habitats can influence host-microbe dynamics, leading to potential downstream effects on host physiology and health. Here, we explored the effect of ecological overlap with humans and domestic species and infection with the protozoan parasite Giardia duodenalis on the bacteria of black and gold howler monkeys (Alouatta caraya), a key sentinel species, in northeastern Argentina. Fecal samples were screened for Giardia duodenalis infection using a nested PCR reaction, and the gut bacterial community was characterized using 16S rRNA gene amplicon sequencing. Habitat type was correlated with variation in A. caraya gut bacterial community composition but did not affect gut bacterial diversity. Giardia presence did not have a universal effect on A. caraya gut bacteria across habitats, perhaps due to the high infection prevalence across all habitats. However, some bacterial taxa were found to vary with Giardia infection. While A. caraya's behavioral plasticity and dietary flexibility allow them to exploit a range of habitat conditions, habitats are generally becoming more anthropogenically disturbed and, thus, less hospitable. Alterations in gut bacterial community dynamics are one possible indicator of negative health outcomes for A. caraya in these environments, since changes in host-microbe relationships due to stressors from habitat disturbance may lead to negative repercussions for host health. These dynamics are likely relevant for understanding organism responses to environmental change in other mammals.

Low genetic diversity and limited genetic structure across the range of the critically endangered Mexican howler monkey (Alouatta palliata mexicana)

Genetic diversity provides populations with the possibility to persist in ever-changing environments, where selective regimes change over time. Therefore, the long-term survival of a population may be affected by its level of genetic diversity. The Mexican howler monkey (Alouatta palliata mexicana) is a critically endangered primate restricted to southeast Mexico. Here, we evaluate the genetic diversity and population structure of this subspecies based on 83 individuals from 31 groups sampled across the distribution range of the subspecies, using 29 microsatellite loci. Our results revealed extremely low genetic diversity (H_O  = 0.21, H_E  = 0.29) compared to studies of other A. palliata populations and to other Alouatta species. Principal component analysis, a Bayesian clustering method, and analyses of molecular variance did not detect strong signatures of genetic differentiation among geographic populations of this subspecies. Although we detect small but significant F_ST values between populations, they can be explained by a pattern of isolation by distance. These results and the presence of unique alleles in different populations highlight the importance of implementing conservation efforts in multiple populations across the distribution range of A. p. mexicana to preserve its already low genetic diversity. This is especially important given current levels of population isolation due to the extreme habitat fragmentation across the distribution range of this primate.

Genetic assignment of illegally trafficked neotropical primates and implications for reintroduction programs

The black and gold howler monkey (Alouatta caraya) is a neotropical primate threatened by habitat loss and capture for illegal trade in Argentina. Using multilocus microsatellite genotypes from 178 A. caraya individuals sampled from 15 localities in Argentina, we built a genotype reference database (GRDB). Bayesian assignment methods applied to the GRDB allowed us to correctly re-assign 73% of individuals to their true location of origin and 93.3% to their cluster of origin. We used the GRDB to assign 22 confiscated individuals (17 of which were reintroduced), and 3 corpses to both localities and clusters of origin. We assigned with a probability >70% the locality of origin of 14 individuals and the cluster of origin of 21. We found that most of the confiscated individuals were assigned to one cluster (F-Ch-C) and two localities included in the GRDB, suggesting that trafficked A. caraya primarily originated in this area. Our results reveal that only 4 of 17 reintroduced individuals were released in sites corresponding to their cluster of origin. Our findings illustrate the applicability of genotype databases for inferring hotspots of illegal capture and for guiding future reintroduction efforts, both of which are essential elements of species protection and recovery programs.

Howler monkey tolerance to habitat shrinking: Lifetime warranty or death sentence?

Habitat loss and fragmentation are major threats to the conservation of nonhuman primates. Given that species differ in their responses to fragmented landscapes, identifying the factors that enable them to cope with altered environments or that cause their extirpation is critical to design conservation management strategies. Howler monkeys (Alouatta spp.) are good models for studying the strategies of tolerant arboreal taxa and how they cope with spatial restriction, because they live in habitats ranging from vast pristine forests to small disturbed fragments and orchards. While some aspects of their ecology and behavior are conserved, others vary in predictable ways in response to habitat shrinking and decreasing resource availability. We argue that the ability of individual howler monkeys to inhabit low-quality environments does not guarantee the long-term persistence of the small populations that live under these conditions. Their local extirpation explains why few forest fragments below a given area threshold are frequently inhabited in landscapes where recolonization and gene flow are compromised by long isolation distances or less permeable matrices. In sum, howlers' ability to cope with habitat restriction at the individual level in the short-term may mask the inevitable fate of isolated populations, thereby compromising the persistence of the species at a regional scale in the long-term if howlers' need for protection in large forests is undervalued.

Southern extension of the geographic range of black-and-gold howler monkeys (Alouatta caraya)

The black-and-gold howler monkey (Alouatta caraya) is widely distributed in Brazil, Bolivia, Paraguay, and northeastern Argentina. Despite this wide distribution, it is locally threatened in some parts of its southern range by forest loss and fragmentation, and yellow fever outbreaks. We present 14 new localities of A. caraya occurrence in the Pampa biome of southern Brazil, extending its range southwards by approximately 100 km.

Red Fox Ancestry and Connectivity Assessments Reveal Minimal Fur Farm Introgression in Greater Yellowstone Ecosystem

Rocky Mountain red foxes Vulpes vulpes macroura potentially encounter other red fox Vulpes vulpes lineages at lower elevations, which may include nonindigenous red foxes derived from fur farms. Introgression from nonindigenous red foxes could have negative evolutionary consequences for the rare Rocky Mountain red fox subspecies. Red foxes at high elevations in the Greater Yellowstone Ecosystem exhibit lighter coat colors than those at lower elevations, potentially indicating that they represent the indigenous subspecies and that gene flow across the elevational gradient is restricted. We collected tissue samples across a 1,750-m elevation range and examined mitochondrial DNA sequences and nuclear DNA microsatellite genotypes to assess the ancestry and genetic population structure of red foxes in the northern Greater Yellowstone Ecosystem. We also used reference samples from fur farm red foxes and indigenous red foxes of the western United States to assess the extent of nonindigenous introgression across the ecosystem. We found little overlap in the elevational distribution of maternally inherited mitochondrial DNA haplotypes: above 1,600 m, we only found indigenous Rocky Mountain haplotypes (n = 4), whereas below 1,600 m, we found haplotypes not indigenous to the Rocky Mountains (n = 5) that were associated with fur farms or indigenous to the Great Plains. In contrast, biparentally inherited microsatellite variation showed little population structure across the elevational gradient. Despite this evidence of nuclear gene flow across the elevational gradient, we found little fur farm introgression in the microsatellite genotypes. It is possible that long-standing nuclear (but apparently not mitochondrial) gene flow between Rocky Mountain red foxes and indigenous red foxes on the Great Plains explained the low nuclear differentiation of these populations. Importantly, our results suggested that high elevations of the northern Greater Yellowstone Ecosystem remained free of significant fur farm introgression. Mitonuclear discordance could reflect sex-biased dispersal, which we hypothesize could be the effect of elevational differences in reproductive phenology.