SciELO - Scientific Electronic Library Online

 
vol.34Differences between caves with and without bats in a Brazilian karst habitatTaxonomic identification using geometric morphometric approach and limited data: an example using the upper molars of two sympatric species of Calomys (Cricetidae: Rodentia) author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

Share


Zoologia (Curitiba)

Print version ISSN 1984-4670On-line version ISSN 1984-4689

Zoologia (Curitiba) vol.34  Curitiba  2017  Epub Oct 23, 2017

http://dx.doi.org/10.3897/zoologia.34.e14881 

RESEARCH ARTICLE

Morphological and genetic diversity in Callithrix hybrids in an anthropogenic area in southeastern Brazil (Primates: Cebidae: Callitrichinae)

Adrielle M. Cezar1 

Leila M. Pessôa1  * 

Cibele R. Bonvicino2  3 

1Laboratório de Mastozoologia, Departamento de Zoologia, Instituto de Biologia, Universidade Federal do Rio de Janeiro. Avenida Brigadeiro Trompowski, Ilha do Fundão, 21941-590 Rio de Janeiro, RJ, Brazil.

2Laboratório de Biologia e Parasitologia de Mamíferos Silvestres Reservatórios, Instituto Oswaldo Cruz. Avenida Brasil 4365, Manguinhos, 21040-360 Rio de Janeiro, Brazil.

3Divisão de Genética, Instituto Nacional do Câncer. Rua André Cavalcanti 37, Centro, 20231-050 Rio de Janeiro, RJ, Brazil.

ABSTRACT

Two species of Callithrix, C. jacchus (Linnaeus, 1758) and C. penicillata (É. Geoffroy, 1812), are considered invasive in Rio de Janeiro. This study determined the genetic and morphological diversity and verified the species involved in the hybridization of 10 individuals from the municipalities of Silva Jardim (N = 9) and Rio das Ostras (N = 1). We compared the external morphology and skull of C. jacchus (N = 15) and C. penicillata (N = 14) specimens deposited in the collection of the National Museum of Rio de Janeiro (MN- UFRJ). Phylogenetic (maximum likelihood and Bayesian inference) and phylogeographical analyses (network analysis) were performed based on cytochrome b sequences. These analyses included hybrids from the metropolitan region of Rio de Janeiro (N = 3), C. penicillata (N = 2), C. jacchus (N = 2), C. geoffroyi (N = 2), C. kuhlii (N = 2), C. aurita (N = 1), and as outgroups, Mico emiliae (N = 1) and Saguinus mystax (N = 1). The pelage and skull characters of most hybrids were more closely related to C. jacchus. Skull morphometric analysis revealed an intermediate state for the hybrids. Phylogenetic analyses revealed a high similarity between the hybrids and C. penicillata. Six haplotypes of hybrids were identified. Network analysis including them and C. penicillata recovered the topology generated by phylogenetic analysis. The results corroborate that C. jacchus and C. penicillata participate in the hybridization process. There was no geographic structure between hybrids from the coastal lowlands and from the metropolitan region of Rio de Janeiro.

KEY WORDS: Atlantic forest; introduced species; marmosets; morphometry; phylogeny.

INTRODUCTION

Callithrix Erxleben, 1777 has six species, all endemic to Brazil. The distribution of Callithrix species is closely associated with the Atlantic Forest. Callithrix jacchus (Linnaeus, 1758) and Callithrix penicillata (É. Geoffroy, 1812) have the largest natural geographical distribution within the genus. They are found in the Atlantic Forest and Caatinga of northeastern Brazil and in the Cerrado of central and northeastern Brazil (De Vivo 1991, Rylands et al. 1996). They are phylogenetically very close, and it has been hypothesized that their most recent common ancestor lived about 700 thousand years ago, in the Atlantic Forest, Cerrado, and Caatinga. A subsequent vicariant speciation event isolated the ancestor of C. penicillata in the Cerrado or Caatinga (Buckner et al, 2014 Malukiewicz et al. 2014.

Species of Callithrix are commonly called marmosets. The range of Callithrix species is allopatric, with some species contacting at the limits of their distribution. However, the ranges of natural species are being altered due to habitat destruction and to anthropogenic introduction of marmoset species outside their natural geographical bounds. As a result of such anthropogenic alterations, C. jacchus and C. penicillata are often found in sympatry with several other Callithrix species and in the natural ranges of other primates (Rylands et al. 1993, 2009, Ruiz-Miranda et al. 2000).

Callithrix jacchus and C. penicillata are found in the state of Rio de Janeiro, Brazil, both in forested and disturbed areas. Their introduction and settlement in the state are the result of illicit domestic and international trafficking of primates. Although the history of their introduction into the coastal lowlands is uncertain, the distribution of marmosets is increasing towards the north of the state of Rio de Janeiro in the lowlands at an estimated rate of 1.2 km per year (Ruiz-Miranda et al. 2000). There, they are found in forest fragments of costal lowlands where the golden lion tamarin, Leontopithecus rosalia (Linnaeus, 1766), is naturally distributed. The interaction between C. jacchus and C. penicillata and the native populations of golden lion tamarins is problematic because the ecology of these primates is similar, which may lead to competition for food and territory and disease transmission (Ruiz-Miranda et al. 2000.

Callithrix jacchus and C. penicillata are differentiated by the colors of the body pelage and of the auricular tufts, and by the insertion of tufts in the ear. C. penicillata has black auricular tufts arranged in front of the ear (pre-auricular) whereas C. jacchus has white auricular tufts arranged around the ear (circum-auricular) (Hershkovitz 1977). Most studies on the identification of hybrids consider only the pelage and mitochondrial DNA, excluding cranial morphological characters (e.g. Alonso et al. 1987, Fuzessy et al. 2014).

This study compares genetic hybrids from the coastal lowland with five species of the genus Callithrix and other hybrids from Rio de Janeiro’s metropolitan region, in order to evaluate morphological and genetic diversity, identify the species involved in hybridization, and verify the geographic structure. To this end, we performed cranial and pelage morphological analyses and molecular phylogenetic estimation using the cytochrome b mitochondrial gene (MT-CYB).

MATERIAL AND METHODS

The sample studied herein comprised 39 individuals: ten Callithrix hybrids from two municipalities in the coastal lowlands of the state of Rio de Janeiro, Silva Jardim (N = 9) and Rio das Ostras (N = 1); and C. penicillata (N = 14) and C. jacchus (N = 15) from localities near the type localities. This sampling strategy ensured that we handled samples of each species separately, thereby avoiding the presence of hybrids. This strategy also limited the sample size. The pure specimens analyzed were identified by the pelage description provided by Hershkovitz (1977). Hybrids were identified by the presence of intermediate characters. Voucher numbers for the specimens analyzed are available in the Appendix 1.

Auricular tuft color and disposition, and general pelage color of all individuals were analyzed to estimate the variation in pelage color and specific patterns for each species and for the hybrids.

A stereoscopic microscope was used to analyze the qualitative characters of C. jacchus and C. penicillata skulls. Differences between species and between each species and the hybrids were identified.

Ten linear cranial measurements ’were taken from the marmosets, with a digital caliper (mm). The first six were defined by Natori (1994) and De Vivo (1991) and the last four in this study: (PL) prosthion to lambda, (EE) euryon to euryon, (iFO) inside frontomalare orbitale to frontomalare orbitale, (BB) bicondylar breadth, (CM) mesial surface of the left upper canine to distal surface of the left second upper molar, (ZA) zygomatic arch breadth, (oFO) outside frontomalare orbitale to frontomalare orbitale, (MS) mandibular symphysis height, (lFM) foramen magnum length, and (bFM) foramen magnum breadth.

To analyze the morphological differences between the studied Callithrix species and the hybrids, we calculated the mean, standard deviation, maximum and minimum values of morphological measurements described above.

Student’s t-test and one-way ANOVA were used to identify differences between species and between each species and the hybrids. Principal component analysis (PCA) was performed to reveal patterns of variation between species and hybrids, and to visualize differences among them. Discriminant function analysis (DFA) was used to verify if the a priori classification of each individual as C. jacchus, C. penicillata, or hybrid using qualitative characters was correct. Analyses were performed in Statistica 8 (Statistica Software Inc.) and R 3.2.4.

DNA samples were obtained from tissue samples extracted from hybrid specimens and from the species Callithrix penicillata, C. aurita, C. kuhlii, Mico rondoni, and Saguinus mystax (Table 1). The latter two were used as outgroups. The DNA was extracted following a phenol-chloroform protocol (Sambrook and Russell 2001). Primers for L14724 (Irwin et al. 1991) and Cytb rev (Casado et al. 2010) were used to amplify MT-CYB. The PCR product of the MT-CYB gene was purified and sequenced using the same PCR primers and the internal primers Sot in1 and Sot in2 (Cassens et al. 2000); Alo aot F and Alo aot R (Menezes et al. 2010); and Citb alo (Bonvicino et al. 2001). The product was labeled with XL and BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems). Sequencing was carried out in an ABI 3130 xl platform. To improve our data, sequences of C. jacchus (accession numbers: AF295586 and AY434079), C. penicillata (accession number: KR817256.1) and C. geoffroyi (accession number: HM368005) were obtained from GenBank online database (www.ncbi.nlm.nih.gov/genbank).

Table 1 Samples used in the phylogenetic analyses. 

ID Species Locality
PRG1415 Callithrix hybrid RJ, Silva Jardim
PRG1416 Callithrix hybrid RJ, Silva Jardim
PRG1417 Callithrix hybrid RJ, Silva Jardim
PRG1454 Callithrix hybrid RJ, Silva Jardim
PRG1456 Callithrix hybrid RJ, Silva Jardim
PRG1702 Callithrix hybrid RJ, Silva Jardim
PRG1703 Callithrix hybrid RJ, Silva Jardim
PRG1706 Callithrix hybrid RJ, Silva Jardim
PRG1708 Callithrix hybrid RJ, Silva Jardim
TDX005 Callithrix hybrid RJ, Rio das Ostras
ZOOSP01031991 C. penicillata Zoológico de São Paulo
KR817256.1 C. penicillata Unavaliable
CRB2587 C. aurita SP, Cunha
CPRJ1016 C. kuhlii Centro de Primatologia do Rio de Janeiro
CRB561 Mico rondoni RO: Ariquemes
CPRJ1621 Saguinus mystax Centro de Primatologia do Rio de Janeiro
CPRJ452 C. kuhlii Centro de Primatologia do Rio de Janeiro
CRB3094 Callithrix hybrid RJ, Rio de Janeiro
CRB3095 Callithrix hybrid RJ, Rio de Janeiro
LBCE18252 Callithrix hybrid RJ, Rio de Janeiro
AF295586 C. jacchus Unavaliable
AY434079 C. jacchus Unavaliable
HM368005 C. geoffroyi Germany, Dresden Zoo

The sequences were analyzed and edited in the software ChromasPro (Mccarthy 1998) and manually aligned in MEGA 5.0 (Tamura et al. 2011).

Genetic distances were estimated with complete deletion using the Kimura 2-parameter model. The MEGA 5.0. Model Generator 0.85 (Keane et al. 2006) identified the best-fitting model for nucleotide substitution using second-order Akaike Information Criteria (AIC) (Akaike 1973).

DNAsp 5 (Librado and Rozas 2009) was used to estimate haplotype and nucleotide diversity. NETWORK was used to reconstruct a median-joining (MJ) network (Bandelt et al. 1999).

Maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees were built. The ML analysis was inferred using a TN93 + I nucleotide substitution model (Tamura and Nei 1993) and the bootstrap analysis was based on 1000 replicates with PhyML (Guindon et al. 2010). MRBAYES 3.2 (Huelsenbeck and Ronquist 2001) was used to build the Bayesian tree using a TN93 + I model.

RESULTS

Morphological analyses

Callithrix penicillata had blackish pre-auricular tufts, a dark brown neck and throat area, and grey striated hairs on the back with an orange medial band and a basal black band. C. jacchus had white circum-auricular tufts and, as C. penicillata, grey striated hairs on the back with an orange medial band and a basal black band. The general color of the body was grayish for C. jacchus, and ranged from shades of gray to brown for C. penicillata. The general color of the pelage varied widely (in light brown tones) for the hybrids. Most hybrids had white auricular tufts, similar in color but not in ear disposition to those of C. jacchus. The tufts were arranged anterior and lateral to the ear and were broken in some parts. One of the individuals had pre-auricular tufts like those of C. penicillata, but with a grayish color. Two individuals lacked tufts because they were young.

Only one cranial qualitative character could be identified as showing distinct patterns between C. penicillata and C. jacchus. Namely, the presence/absence of a space in the upper jaw after the second molar (Figs 1-4). The space was absent in C. penicillata (i.e. the maxilla ends abruptly after the last molar) and present in C. jacchus. Hybrids exhibited a pattern equivalent to that seen in C. jacchus. Table 2 summarizes the mean and standard deviation of each linear cranial measurement for all Callithrix species and hybrids. Three characters were significantly different between species (as shown by Student’s t-test and one-way ANOVA): iFO, oFO, and MS. In the hybrids, two of these characters (iFO and oFO) were more similar to those of C. penicillata and one (CM) was more similar to that of C. jacchus. The comparison between the hybrids and each species separately (Student’s t-test and one-way ANOVA) showed that EE and CM were significantly different between the hybrids and C. penicillata, whereas EE, iFO, oFO, and MS were significantly different between the hybrids and C. jacchus (Table 3).

Figures 1-4 Qualitative character differentiating Callithrix jacchus (1-2) from C. penicillata (3-4): presence/absence of a space in the upper jaw after the second molar, indicated by arrows. 

Table 2 Approximate mean and standard deviation (mm) for the linear cranial measurements obtained from samples of Callithrix species and hybrids. Measurements are identified in the left column, and species/hybrids are identified in the header. (PL) Prosthion to lambda, (EE) euryon to euryon, (iFO) inside frontomalare orbitale to frontomalare orbitale, (BB) bicondylar breadth, (CM) mesial surface of the left upper canine to distal surface of the left second upper molar, (oFO) outside frontomalare orbitale to frontomalare orbitale, (MS) mandibular symphysis height, (lFM) foramen magnum length, (bFM) foramen magnum breadth, (ZA) zygomatic arch breadth. 

Callithrix penicillata Callithrix jacchus Hybrids
PL 44.6 ± 0.86 44.32 ± 1.01 45.26 ± 1.27
EE 25.24 ± 0.65 25.30 ± 0.86 27.43 ± 0.78
iFO 22.77 ± 0.65 23.43 ± 0.90 23.89 ± 0.51
oFO 24.32 ± 0.57 25.12 ± 0.95 25.14 ± 0.59
BB 24.38 ± 0.85 24.54 ± 1.36 25.39 ± 0.67
MS 8.79 ± 0.95 8.08 ± 0.49 8.61 ± 0.61
CM 10.98 ± 0.48 10.72 ± 0.42 11.03 ± 0.29
iFM 10.98 ± 0.48 6.05 ± 0.67 6.73 ± 0.62
bFM 6.66 ± 0.17 6.20 ± 0.49 7.13 ± 0.27
ZA 10.57 ± 0.54 11.05 ± 0.43 11.13 ± 0.56

Table 3 One-way ANOVA and Student’s t-test statistical analyses of linear cranial measurements obtained from Callithrix species and hybrids. 

PL EE iFO oFO BB MS CM iFM bFM ZA
C. penicillata x C. jacchus ns ns * * ns * ns ns ns ns
C. penicillata x Hybrids ns * ns ns ns * ns ns ns ns
C. jacchus x Hybrids ns * * * ns * ns ns ns ns

ns = not significant and * = significant at p < 0.05. (PL) prosthion to lambda; (EE) euryon to euryon; (iFO) inside frontomalare orbitale to frontomalare orbitale; (BB) bicondylar breadth; (CM) mesial surface of the left upper canine to distal surface of the left second upper molar; (oFO) outside frontomalare orbitale to frontomalare orbitale; (MS) mandibular symphysis height; (lFM) foramen magnum length; (bFM) foramen magnum breadth; (ZA) zygomatic arch breadth.

The first three components of the PCA (MS, iFM, and bFM) accounted for most of the observed skull variation (PC1 = 39.9%, PC2 = 26.4%, and PC3 = 13.6%) (Figs 5-6). MS contributed positively and lFM and bFM contributed negatively to PC1; MS and bFM contributed negatively and lFM contributed slightly positively to PC2.

Figures 5-6 PCA analysis results. (5) Scatter plot of scores for principal component 1 x 2. Black circles represent Callithrix penicillata, white circles represent C. jacchus and grey circles represent hybrids. (6) Contribution of morphometric variables to the principal components. Vectors indicate the loadings of the scores for each variable on the first two principal components. 

The DFA analysis confirmed the a priori classification, revealing highly significant inter-sample variation (Wilk’s lambda = 0.065998, approximate F = 6.0744, p < 0.0001). The scatter plot showed three distinct groups of points (Figs 7-8), each group representing either one of the species or the hybrids. Measurements oFO and ZA contributed most to the first discriminant function while lFM, iFO, and EE contributed most to the second discriminant function (Figs 7-8).

Figures 7-8 DFA analysis results. (7) Scatter plot of scores for discriminant function 1 x 2. Three distinguishable groups characterize the species and hybrids analyzed: (+) Callithrix penicillata, (Δ) C. jacchus, and (○) hybrids. (8) Contribution of morphometric variables to the discriminant functions. Vectors indicate the loadings of the scores for each variable on the first two discriminant functions. 

Molecular analyses

The cytochrome b gene, comprising 1140 bp, was sequenced for all specimens. Only the hybrids shared haplotypes. The 13 hybrid sequences had six haplotypes, two of which were shared by more than one specimen. Analysis of sequences from all hybrids revealed 21 variable sites (18 transitions and three transversions), with estimates of genetic distance ranging from 0.001 to 0.01%.

The phylogenetic analyses resulted in a monophyletic Callithrix genus (Fig. 9), with C. aurita as the sister taxon to the remaining species. The hybrids were grouped into three distinct clades within the C. penicillata lineage.

Figure 9 The Bayesian and Maximum Likelihood analyses for MT-CYB of Callithrix, rooted by Saguinus. Numbers close to branches are boostrap values and posterior probability, respectively. 

The median-joining network analysis was focused on the relationship between C. penicillata and hybrids from the coastal lowlands (CL) and from the metropolitan region (MR) of Rio de Janeiro. It revealed relationships between the hybrid haplotypes of the two regions (Fig. 10). Each of the seven haplotypes in the network is separated by at least one variable site (from a total of 21 variable sites). Haplotype diversity (Hd) was 0.7582 and the nucleotide diversity (Pi) was 0.00537. Seven specimens, including CL and MR hybrids, shared haplotype 1 (H1). Each of haplotypes 2, 3, and 4 (H2, H3 and H4) included one CL individual. Haplotype 5 (H5) was shared by two CL specimens. Haplotype, H6, included two pure C. penicillata specimens and a single MR individual. Haplotype 7 (H7) included one pure C. penicillata specimen. Results from the network analysis were similar to those described for the phylogenetic analysis.

Figure 10 Haplotype network of MT-CYB sequences for Callithrix penicillata and Callithrix hybrids. Circles represent distinct haplotypes. White circles represent C. penicillata, black circles represent hybrids from the MR, and gray circles represent hybrids from the CL. The size of each circle is proportional to the number of individuals per haplotype, with the smallest circle corresponding to n = 1. The lozenge represents the medium-vector. Numbers near lines between haplotypes represent the number of mutations. 

DISCUSSION

Hybridization has been consistently documented in primates. However, the effects of natural and anthropogenic hybridization on biodiversity remain unclear. Differentiating between these two types of hybridization is a challenge in evolution and conservation studies (Malukiewicz et al. 2014). Previous reports showed that hybrids have intermediate characters between the parental species (Hershkovitz 1977, Alonso et al. 1987, Fuzessy et al .2014, a finding that is partially supported by results of the current study. One possible explanation for such finding is the fact that the parental species are phylogenetically very close. (Tagliaro et al. 1997, Buckner et al. 2014. Mallet (2005) argued that hybridization occurs in approximately 10 percent of the mammalian species, usually between groups that diverged more recently (1 to 2 million years). C. jacchus and C. penicillata are species that diverge very recently (Malukiewicz et al. 2015).

We observed that the auricular tufts (one of the main characters used to identify marmosets) of the hybrids do not match the description of any species of the genus. In the hybrids the tufts were white or gray, arranged anterior and laterally to the ear and were broken in some portions, a mosaic that may result from several generations of hybridization. In the study of Alonso et al. (1987), involving hybrids between C. jacchus and C. penicillata in natural hybrid zones, the recorded color pattern and shape of the auricular tufts suggest that there is a reproductive isolation mechanism between these species, since there was little penetration of C. jacchus characters in the C. penicillata population and vice-versa. Individuals with pure parental phenotypes were absent in the hybrid groups of that study, as was observed in the present one.

Hybrids of mixed ancestry between two marmosets species (C. penicillata and C. geoffroyi) had greater morphological variation compared with individuals of pure ancestry (Fuzessy et al. 2014), which could possibly explain the results found for the hybrid individuals in our study. Malukiewicz et al. (2014) agreed that the current situation of coastal lowland marmosets is the result of multiple introductions and that there some new genetic variations caused by new introductions. These multiple introductions may not only be from pure parental individuals but also from hybrid individuals, resulting in crosses between pure and mixed ancestries that give rise to highly variable phenotypes. The hybrids analyzed herein showed a large variation in body color, with a predominance of light brown regions throughout the body (detected even in the tail). This has not been previously observed in either of the two studied species and may be related a transgressive segregation. Studies on hybrid populations have occasionally reported the presence of phenotypes that are extreme relative to those of either parental line (Rieseberg and Ellstrand 1993, Cosse et al. 1995). The generation of these extreme phenotypes is referred to as transgressive segregation, a phenomenon specific to segregating hybrid generations and refers to the fraction of individuals that exceed parental phenotypic values in either a negative or positive direction (Rieseberg et al. 1999).

Only one qualitative character differed between C. jacchus and C. penicillata: the presence/absence of a prolongation of the maxilla, after the last molar. Hybrids had the space, being therefore similar to C. jacchus. This is an unprecedented observation in the taxonomic literature, since authors report the skulls of C. jacchus and C. penicillata as being very similar (Garbino 2015, Natori 1994).

Three measurements explained most of the variation observed in the PCA. When specimens are plotted along PC1 and PC2, the distribution is somewhat scattered. However, there are some trends according to the species and hybrids. An overlap exists between C. penicillata and C. jacchus, which may be explained by the great similarity between their skulls. Similarly, there is an overlap between hybrids and C. jacchus, indicating that hybrids are more related with C. jacchus than with C. penicillata. This result differs from that obtained by the univariate analysis.

The results of the DFA show that Wilk’s lambda was relatively small (0.065998) and the approximate F value was high (6.0744, p < 0.0001), corroborating the a priori classification.

Results of the phylogenetic analysis supported the monophyly of Callithrix, as expected, based on previous phylogenies (Tagliaro et al. 1997, Perelman et al. 2011, Garbino 2015). In the current study, C. jacchus is the sister group of C. penicillata, and these two are the most recent split in the genus. Hybrids were grouped in a clade with C. penicillata, suggesting that this species is the one involved in the maternal lineage of the individuals, being a direct or indirect parental species. Where C. penicillata was not a direct parent, crosses between lineages of hybrids with C. penicillata maternal origin may have occurred. Both cases suggest that males preferentially mate with females from the C. penicillata lineage.

Mate choice by males is primarily associated with mate availability and with variation in female quality. If males engage in paternal care, as is the case in Callithrix, the average number of available females is likely to be high compared to the capacity of males to mate with them. In this situation, males are less likely to be able to mate with all available females, rejecting some of them. If the benefit of mating with specific females exceeds the cost of assessing them, mate choice can evolve (Edward and Chapman 2011). Although fertile hybrids will be generated between the two species, some matings may result in less viable hybrids than others, influencing the process of mate choice (Coimbra-Filho et al. 1993). The relationship found between the hybrids and C. penicillata in this study suggests that hybrid males prefer C. penicillata females. This may be because hybrids that result from mating with C. penicillata females are more viable than those that result from mating with C. jacchus females.

Primates recognize potential mates (members of the same species or not) based on visual (mainly on the face), acoustic, olfactory, and other sensorial cues. Different patterns of facial color and auricular hairs have diagnostic value for each species in the taxonomy of Callithrix. One conspicuous morphologic character that distinguishes these four species is the coloration of the auricular tufts (Cavalcanti and Langguth 2008).

Cavalcanti and Langguth (2008) suggested that there is an isolation mechanism based on head color, since when two different species are together, one responds to the other with significantly less frequency than to its own species. However, this is not always the case, suggesting that, in spite of the conspicuous differences in facial coloration patterns, the evaluator species recognizes individuals of the other, “cue-bearing” species as potential sexual competitors. During speciation events in Callithrix, reproductive isolation mechanisms did not necessarily appear simultaneously in all the presently recognized species. Populations of each species have a different history. C. jacchus and C. penicillata are very closely-related and have split recently from a common ancestor. Thus, it is possible that reproductive isolation is not complete in this case. Although the primate literature is relatively rich in studies of sexual selection and mate preference in hybrid zones (Shurtliff 2013), little attention has been given to these topics in recent decades. A study on howler monkeys (Alouatta palliata and A. pigra) showed that hybridization and subsequent backcrossing are directionally biased, probably producing only fertile hybrid females and inviable or infertile males. This suggests that a process of mate choice may occur for Neotropical primates in hybrid zones (Cortés-Ortiz et al. 2007).

The low genetic distance between hybrids revealed that they are genetically very close, corroborating the hypothesis of multiple introductions suggested by the morphological results and the short passage of new genetic variations.

Haplotype H1 was found in most of the samples, containing both CL and MR individuals and precluding the existence of a geographical structure. A possible explanation is that introductions of marmosets into the coastal lowlands occurred recently with animals from the metropolitan region. De Morais et al. (2008) argued that the history of the introduction of marmosets in the coastal lowlands of Rio de Janeiro is uncertain and emphasized the occurrence of two major releases of marmosets, seized by regulatory agencies (> 60 animals), between 1983-1987, in regions close to the study area. Each of the haplotypes, H2, H3, and H4 was found in only one CL hybrid, and haplotype H5 was found in two CL individuals. This diversity in a single region may be the result of multiple introductions, with the arrival of new individuals from different localities. Haplotype H6 included one C. penicillata and one MR individual, but the fact that haplotype H1 also included an MR individual indicates that H6 is not characteristic only of the MR.

The network analysis showed close a relationship between the hybrids and C. penicillata, confirming the result from the ML and BI analyses. At least one variable site separates each haplotype, which may reflect differences due to multiple introductions or a polymorphism in the population. Twenty-one variable sites were observed, with eighteen transitions and three transversions. Haplotype diversity among hybrids and C. penicillata was high (0.7582), whereas nucleotide diversity was low (0.00537). Low nucleotide diversity may be explained by the founder effect, whereby an introduction is followed by population growth. The several haplotypes also had low diversity, possibly due to multiple anthropogenic introductions that continue to occur. The low level of genetic diversity both in the CL and between the CL and the MR suggests a recent history of population expansion, most likely associated with introductions.

These issues reflect the need to perform more detailed studies concerning hybridization and the development of morphological characters in order to obtain a better understanding of the evolution of Callithrix.

ACKNOWLEDGMENTS

We would like to thank the financing agency Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for providing a scientific initiation scholarship to Adrielle M. Cezar and research fellowships to Leila M. Pessôa (process 308505/2016-6). We also thank João A. de Oliveira and Pablo R. Gonçalves for their help with analyzing the specimens at the Museu Nacional da Universidade Federal do Rio de Janeiro (MN-UFRJ) and the Núcleo em Ecologia e Desenvolvimento Sócio-Ambiental de Macaé (NUPEM/UFRJ). We are grateful to the members of the monograph examination commission, Alcides Pissinatti, Daniel F. da Silva, Maria Lucia Lorini, and Héctor Seuánez, for their availability and contributions; to the Centro de Primatologia do Rio de Janeiro; to the São Paulo Zoo and to the Dresden Zoo for granting permission to access the materials used in this study. We are grateful to all our laboratory colleagues who helped us during this project. The English spelling and grammar of a previous version of this manuscript were edited by Publicase, PPGBBE, UFRJ financed the English editing work.

LITERATURE CITED

Akaike H (1973) Information theory and an extension of the maximum likelihood principle. In: Petrov BN, Csáki F (Eds) 2nd International Symposium on Information Theory, Tsahkadsor (Armenia), September 2-8, 1971. Akadémiai Kiadó, Budapest, 267-281. [ Links ]

Alonso C, de Faria DS, Langguth A, Santee DF (1987) Variação da pelagem na área de intergradação entre Callithrix jacchus e Callithrix penicillata. Revista Brasileira de Biologia 47: 465-470. [ Links ]

Bandelt HJ, Forster P, Röhl A (1999) Median-joining networks for inferring intraspecific phylogenies. Molecular Biology Evolution 16: 37-48. https://doi.org/10.1093/oxfordjournals.molbev.a026036Links ]

Bonvicino CR, Lemos B, Seuánez HN (2001) Molecular phylogenetics of howler monkeys (Alouatta, Platyrrhini): A comparison with karyotypic data. Chromosoma 110: 241-246. https://doi.org/10.1007/s004120000128Links ]

Buckner JC, Alfaro JL, Rylands AB, Alfaro ME (2014) Biogeography of the marmosets and tamarins (Callitrichidae). Molecular Phylogenetic Evolution 82: 413-425. https://doi.org/10.1016/j.ympev.2014.04.031Links ]

Casado F, Bonvicino CR, Nagle C, Comas B, Manzur TD, Lahoz MM, Seuánez HN (2010) Mitochondrial Divergence Between 2 Populations of the Hooded Capuchin, Cebus (Sapajus) cay (Platyrrhini, Primates). Journal of Heredity 101: 261-269. https://doi.org/10.1093/jhered/esp119Links ]

Cassens I, Vicario S, Waddell VG, Balchowsky H, van Belle D, Ding W, Fan C, Lal Mohan RS, Simões-Lopes PC, Bastida R, Meyer A, Stanhope MJ, Milinkovitch MC (2000) Independent adaptation to riverine habitats allowed survival of ancient cetacean lineages. Proceedings of National Academy of Sciences of the United States of America 97: 11343-11347. https://doi.org/10.1073/pnas.97.21.11343Links ]

Cavalcanti GC, Langguth A (2008) Recognition of mate and speciation in marmoset genus Callithrix (Primates, Cebidae, Callithriquinae [sic]). Revista Nordestina de Biologia 19: 59-73. [ Links ]

Coimbra-Filho AF, Pissinatti A, Rylands AB (1993) Experimental multiple hybridism among Callithrix species from eastern Brazil. In: Rylands AB (Ed.) Marmosets and tamarins: systematics, ecology and behavior. Oxford University Press, Oxford, 95-120. [ Links ]

Cortés-Ortiz L, Duda Jr TF, Canales-Espinosa D, García-Orduña F, Rodríguez-Luna E, Bermingham E (2007) Hybridization in large-bodied new world primates. Genetics 176: 2421-2425. https://doi.org/10.1534/genetics.107.074278Links ]

Cosse AA, Campbell MG, Glover TJ, Linn Jr CE, Todd JL, Baker TC, Roelofs WL (1995) Pheromone behavioral responses in unusual male European corn borer hybrid progeny not correlated to electrophysiological phenotypes of their pheromone-specific antennal neurons. Experientia 51: 809-816. https://doi.org/10.1007/BF01922435Links ]

De Morais Jr MM, Ruiz-Miranda CR, Grativol AD, Andrade CC, Lima CA (2008) Os sagüis, Callithrix jacchus e C. penicillata, como espécies invasoras na região de ocorrência do mico-leão dourado. In: Oliveira PP, Grativol AD, Ruiz-Miranda CR (Eds) Conservação do mico-leão-dourado: enfrentando os desafios de uma paisagem fragmentada. Campos dos Goytacazes, Universidade Estadual do Norte Fluminense, vol. 1, 86-117. [ Links ]

De Vivo M (1991) Taxonomia de Callithrix Erxleben, 1777 (Callitrichidae, Primates). Belo Horizonte, Fundação Biodiversitas. [ Links ]

Edward DA, Chapman T (2011) The evolution and significance of male mate choice. Trends in Ecology and Evolution 26: 647-654. https://doi.org/10.1016/j.tree.2011.07.012Links ]

Fuzessy LF, de Oliveira SI, Malukiewicz J, Silva FFR, do Carmo PM, Boere V, Ackermann RR (2014) Morphological Variation in Wild Marmosets (Callithrix penicillata and C. geoffroyi) and Their Hybrids. Evolutionary Biology 41: 480-493. https://doi.org/10.1007/s11692-014-9284-5Links ]

Garbino GST (2015) How many marmoset (Primates: Cebidae: Callitrichinae) genera are there? A phylogenetic analysis based on multiple morphological systems. Cladistics 31: 652-678. https://doi.org/10.1111/cla.12106Links ]

Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O (2010) New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Systematic Biology 59: 307-21. https://doi.org/10.1093/sysbio/syq010Links ]

Hershkovitz P (1977) Living New World Monkeys (Platyrrhini). With an Introduction to Primates. The University of Chicago Press, Chicago, vol. 1. [ Links ]

Huelsenbeck JP, Ronquist F (2001) MrBayes: Bayesian inference of phylogenetics trees. Bioinformatics 17: 754-755. https://doi.org/10.1093/bioinformatics/17.8.754Links ]

Irwin DM, Kocher TD, Wilson AC (1991) Evolution of the cytochrome b gene of mammals. Journal of Molecular Evolution 32: 128-144. https://doi.org/10.1007/BF02515385Links ]

Keane TM, Creevey CJ, Pentony MM, Naughton TJ, Mcinerney JO (2006) Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evolution Biology 6: 29. https://doi.org/10.1186/1471-2148-6-29Links ]

Librado P, Rozas J (2009) DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451-1452. https://doi.org/10.1093/bioinformatics/btp187Links ]

McCarthy C (1998) Chromas. Queensland, School of Health Science, Griffith University, version 1.45. [ Links ]

Mallet J (2005) Hybridization as an invasion of the genome. Trends in Ecology and Evolution 20: 229-237. https://doi.org/10.1016/j.tree.2005.02.010Links ]

Malukiewicz J, Boere V, Fuzessy LF, Grativol AD, French JA, Silva I, Pereira LCM, Ruiz-Miranda CR, Valença Y, Stone AC (2014) Hybridization effects and genetic diversity of the common and black-tufted marmoset (Callithrix jacchus and Callithrix penicillata) mitochondrial control region. American Journal of Physical and Anthropology 155: 522-536. https://doi.org/10.1002/ajpa.22605Links ]

Malukiewicz J, Boere V, Fuzessy LF, Grativol AD, de Oliveira e Silva I, Pereira LCM, Ruiz-Miranda CR, Valença Y, Stone AC (2015) Natural and Anthropogenic Hybridization in Two Species of Eastern Brazilian Marmosets (Callithrix jacchus and C. penicillata). PLoS ONE 10(6): e0127268. https://doi.org/10.1371/journal.pone.0127268Links ]

Menezes AN, Bonvicino CR, Seuánez HN (2010) Identification, classification and evolution of owl monkeys (Aotus, Illiger 1811). BMC Evolution Biology 10: 248. https://doi.org/10.1186/1471-2148-10-248Links ]

Natori M (1994) Craniometrical variations among eastern Brazilian marmosets and their systematic relationships. Primates 35: 167-176. https://doi.org/10.1007/BF02382052Links ]

Perelman P, Johnson WE, Roos C, Seuánez HN, Horvath JE, Moreira MA, Kessing B, Pontius J, Roelke M, Rumpler Y, Schneider MP, Silva A, O’Brien SJ, Pecon-Slattery J (2011) A Molecular Phylogeny of Living Primates. PLoS Genetic 7: e1001342. https://doi.org/10.1371/journal.pgen.1001342Links ]

Rieseberg LH, Archer MA, Wayne RK (1999) Transgressive segregation, adaptation and speciation. Heredity 83: 363-372. https://doi.org/10.1038/sj.hdy.6886170Links ]

Rieseberg LH, Ellstrand NC (1993) What can morphological and molecular markers tell us about plant hybridization? Critical Review Plant Science 12: 213-241. https://doi.org/10.1080/07352689309701902Links ]

Ruiz-Miranda CR, Affonso AG, Martins A, Beck BB (2000) Distribuição do sagui (Callithrix jacchus) nas áreas de ocorrência do mico leão dourado no Estado de Rio de Janeiro. Neotropical Primates 8: 98-101. [ Links ]

Rylands AB, Coimbra-Filho AF, Mittermeier RA (1993) Systematics, geographic distribution and some notes on the conservation status of the Callitrichidae. In: Rylands AB (Ed.) Marmosets and tamarins: systematics, ecology and behaviour. Oxford Science Publications, Oxford, 11-77. [ Links ]

Rylands AB, da Fonseca GAB, Leite YLR, Mittermeier RA (1996) Primates of the Atlantic forest: origin, distributions, endemism, and communities. In: Norconk MA, Rosenberger AL, Garber PA (Eds) Adaptive radiations of neotropical primates. Plenum, New York, 21-51. https://doi.org/10.1007/978-1-4419-8770-9_2Links ]

Rylands AB, Coimbra-Filho AF, Mittermeier RA (2009) The systematics and distributions of the marmosets (Callithrix, Callibella, Cebuella, and Mico) and Callimico (Callimico) (Callitrichidae, Primates. In: Ford SM, Porter LM, Davis LC (Eds) The smallest anthropoids: The marmoset/Callimico radiation. Springer, New York, 25-61. https://doi.org/10.1007/978-1-4419-0293-1_2Links ]

Sambrook J, Russel DW (2001) Molecular cloning. CSHL Press, Cold Spring Harbor, 3rd ed. [ Links ]

Shurtliff QR (2013) Mammalian hybrid zones: a review. Mammal Review 43: 1-21. https://doi.org/10.1111/j.1365-2907.2011.00205.xLinks ]

Tamura K, Nei M (1993) Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Molecular Biology and Evolution 10: 512-526. [ Links ]

Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution 28: 2731-2739. https://doi.org/10.1093/molbev/msr121Links ]

Tagliaro CH, Schneider MPC, Schneider H, Sampaio I, Stanhope M (1997) Marmoset phylogenetics, conservation perspectives, and evolution of the mtDNA control region. Molecular Biology and Evolution 14: 674-684. https://doi.org/10.1093/oxfordjournals.molbev.a025807Links ]

Editorial responsibility: Maria Wilhelmina Pil

ZooBank: http://zoobank.org/D0A6F9E3-E613-40CC-B180-170074DFCF39

APPENDIX 1

Specimens provenance are summarized below.

The C. penicillata and C. jacchus specimens are deposited in the National Museum, Federal University of Rio de Janeiro. Their voucher numbers are: Callithrix penicillata (MN4260-4262, MN4264-4266, MN4268-4270, MN10681, MN11334, MN23798, MN23800 and MN30549) and C. jacchus (MN3953, MN5521, MN5528, MN5535, MN5546, MN5551, MN5573, MN17274-17276, MN17291, MN23772, MN23774, MN30544, MN30548).

The hybrid (Callithrix sp.) specimens are deposited at the Center for Ecology and Socio-Environmental Development of Macaé, Federal University of Rio de Janeiro (PRG1415-1417, PRG1454, PRG1456, PRG1702, PRG1703, PRG1706, PRG1708, TXD005).

Received: July 18, 2016; Revised: December 26, 2016; Accepted: January 13, 2017

* Corresponding author: Leila M. Pessôa (pessoa@acd.ufrj.brr)

Author Contributions:

AMC, LMP and CRB participated equally in the preparation of this article.

Competing Interests:

The authors have declared that no competing interests exist.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License