- Open Access
Modeling the origins of mammalian sociality: moderate evidence for matrilineal signatures in mouse lemur vocalizations
Frontiers in Zoology volume 11, Article number: 14 (2014)
Maternal kin selection is a driving force in the evolution of mammalian social complexity and it requires that kin are distinctive from nonkin. The transition from the ancestral state of asociality to the derived state of complex social groups is thought to have occurred via solitary foraging, in which individuals forage alone, but, unlike the asocial ancestors, maintain dispersed social networks via scent-marks and vocalizations. We hypothesize that matrilineal signatures in vocalizations were an important part of these networks. We used the solitary foraging gray mouse lemur (Microcebus murinus) as a model for ancestral solitary foragers and tested for matrilineal signatures in their calls, thus investigating whether such signatures are already present in solitary foragers and could have facilitated the kin selection thought to have driven the evolution of increased social complexity in mammals. Because agonism can be very costly, selection for matrilineal signatures in agonistic calls should help reduce agonism between unfamiliar matrilineal kin. We conducted this study on a well-studied population of wild mouse lemurs at Ankarafantsika National Park, Madagascar. We determined pairwise relatedness using seven microsatellite loci, matrilineal relatedness by sequencing the mitrochondrial D-loop, and sleeping group associations using radio-telemetry. We recorded agonistic calls during controlled social encounters and conducted a multi-parametric acoustic analysis to determine the spectral and temporal structure of the agonistic calls. We measured 10 calls for each of 16 females from six different matrilineal kin groups.
Calls were assigned to their matriline at a rate significantly higher than chance (pDFA: correct = 47.1%, chance = 26.7%, p = 0.03). There was a statistical trend for a negative correlation between acoustic distance and relatedness (Mantel Test: g = -1.61, Z = 4.61, r = -0.13, p = 0.058).
Mouse lemur agonistic calls are moderately distinctive by matriline. Because sleeping groups consisted of close maternal kin, both genetics and social learning may have generated these acoustic signatures. As mouse lemurs are models for solitary foragers, we recommend further studies testing whether the lemurs use these calls to recognize kin. This would enable further modeling of how kin recognition in ancestral species could have shaped the evolution of complex sociality.
Maternal kin selection (the preferential treatment of matrilineal relatives [1, 2]) has been argued to be one of the driving forces in the evolution of mammalian sociality, underpinning some of the most complex and intriguing social behaviors including communal infant rearing and socialization, the evolution of group-living, alliance formation and cooperation [1, 3–5]. While such manifestations of kin selection are well documented in gregarious species that live in complex social groups [3, 4], its evolutionary foundations are likely to have emerged in less complex, ancestral species (, but see ). Given that ancestral mammals are believed to have been asocial with no social relationships maintained outside of mating and rearing infants , tracing how maternal kin selection may have formed the backbone for this transition is likely to be crucial to understanding how social complexity evolves.
A prerequisite of maternal kin selection in any mammalian social system is that maternal kin must be sufficiently distinctive from nonkin that they can be recognized and thus receive preferential treatment [1, 2]. For the asocial and nocturnal ancestral mammals , this would have also meant being distinctive over distances, through darkness, and dense foliage where visual and olfactory cues would have been inefficient. Mammals under these conditions would be expected to benefit from having matrilineal signatures in their vocalizations.
To date, much of the attention that has been given to investigating matrilineal signatures in mammalian vocalizations has focused on social species (ie. goats [8, 9], meerkats , marmots , sperm whales and killer whales [12–14], bats [15–17] and the socially variable house mouse [18, 19]). Much less has been done on solitary species (i.e., pandas ). In the solitary pandas, individual signatures were found, but there was no correlation between overall acoustic distance between individuals and their relatedness, and only a few individual parameters correlated with relatedness . Though the authors did not clarify whether relatedness was matrilineal, patrilineal, or both, the lack of stronger results may still indicate that pressure to encode kinship within vocalizations may not be as strong as in the more social species . Each of these studies that investigated kin signatures, either exclusively focused on matrilineal relatedness or had a high likelihood of relatedness from both patrilineal and matrilineal relationships, thus suggesting that matrilineal relatedness had a strong role in the signatures found. We differentiate between individual signatures that may be recognized by kin (i.e., primates: [21–24], pinnipeds , elephants , dolphins ) and matrilineal signatures. Matrilineal signatures have the important distinction that they may enable the recognition of unfamiliar maternal kin via the similarity to known maternal kin, thus facilitating the preferential treatment of unfamiliar maternal kin.
In order to better understand the evolutionary transition from asociality to social complexity, we focus on primates, an order in which some lineages have evolved highly complex, cohesive social groups while other lineages are believed to have retained the social system that is believed to be ancestral to primates: solitary foraging . In the lineages that evolved social systems with cohesive social groups, the ancestral solitary foragers are believed to have been a transition phase between asociality and group-living . Solitary foragers forage alone, but maintain a dispersed social network of relationships with conspecifics communicating through vocalizations and scent-marks, and often have consistent co-sleeping associations . It is these dispersed social networks in ancestral primates that are thought to have been the foundation for the evolution of more complex primate social systems , thus they are likely to have been crucial for kin networks and a likely pathway for kin selection [28, 29].
In order to determine whether matrilineal signatures in vocalizations may have facilitated matrilineal kin selection in solitarily foraging ancestral primates, we use the gray mouse lemur (Microcebus murinus) to model ancestral primates. Mouse lemurs are frequently used as ancestral primate models because their socioecology is thought to be similar to that of the last common ancestor of the primate order [6, 30–46]. Like mouse lemurs today, ancestral primates are thought to have been small-bodied, small-brained nocturnal solitary foragers that forage for fruits and insects in the thin, terminal ends of branches [6, 30–46]. Therefore, we use the dispersed social networks of living mouse lemurs to model ancestral primate social organization and to reconstruct the social behavior patterns from which present-day primate diversity evolved .
Our current knowledge of the gray mouse lemur’s dispersed social networks makes it an excellent model species in which to test for matrilineal signatures. Male dispersal and female philopatry are common [47, 48]. Both sexes forage solitarily in home ranges that overlap with those of other individuals of both sexes [43, 49]. During the day, adult males sleep alone [44, 50]. Females form sleeping groups with female kin and cooperatively raise their young in tree holes [42, 50]. Immature males and females are socialized within these groups [42, 50] and thus have ample opportunity to hear and learn the calls of their matrilineal kin. However, given that larger nest groups may split, it is also possible for subsequent generations to encounter matrilineal kin with whom they personally did not share a nest . It is also possible that inherited vocal tract morphology (see source-filter theory: i.e., [51–53]) could cause related individuals to produce similar calls. Thus, both genetic factors and social learning could contribute to the development of matrilineal signatures in this species.
Mouse lemurs have an elaborate vocal repertoire and use vocalizations in a diversity of social interactions (e.g., mating contexts [54, 55], mother-infant communications , emotional state [5, 57], paternal kin recognition ). We chose to investigate the individually distinctive agonistic call . It is a short, frequency modulated vocalization with an upward and downward sweep (Figure 1) containing harmonics in both the audible and ultrasonic range . Because aggressive/defensive encounters have the potential to be very costly due to injuries sustained, we predicted that it would be advantageous for agonistic calls to contain matrilineal signatures so that aggression amongst matrilineal relatives could be minimized. We hypothesized that these agonistic calls will be distinctive by matrilineal kin group and that the genetic relatedness of female dyads will negatively correlate with their acoustic distance. We found moderate evidence for matrilineal signatures and a trend suggesting that increasing relatedness is associated with decreasing acoustic distance. Further studies are needed to determine whether mouse lemurs use these signatures to recognize kin.
Acoustic differences in agonistic calls between matrilines
Qualitative differences are visible between matrilines (Additional file 1) in both frequency and temporal parameters. For example, several individuals in matriline 6 gave calls with an unusually high peak frequency, matriline 5 gave calls of longer duration, and matrilines 1 and 3 typically gave lower frequency calls, with the calls of matriline 1 being generally qualitatively shorter than those of 3.
The principal component analysis produced two components which together explained 66.2% of the variation in the original dataset. The first component was highly correlated (>0.4 or < -0.4) with all original acoustic parameters, but correlated most strongly (>0.7) with the frequency parameters, and thus, it is referred to as the frequency component (48.6% of the total variation). The second component correlated highly (>0.4 or < -0.4) with call duration and inter-call interval and is thus called the time component (17.5% of the total variation). Table 1 shows the 25% quartile, median, and 75% quartile for each of the original acoustic parameters and their loadings on the frequency component and the time component. Table 2 shows the matrix of pairwise acoustic distances calculated for each dyad of females.
Median pairwise relatedness for all dyads in the population is r = -0.02 (n = 107 individuals, min = -0.38, max = 0.91). Median pairwise relatedness for the females within the kin groups was r = 0.41 (n = 16 females, min = 0.30, max = 0.52, Table 3), whereas the between kin group median relatedness was r = -0.02 (n = 16 females, min = -0.12, max = 0.06). Table 2 shows the matrix of pairwise relatedness values of all the females in the kin groups. Within the females in the population, we found seven mitochondrial haplotypes (Figure 2). The kin groups in this study belonged to the three most frequent haplotypes (H3, H4, H6).
Kin group signatures and correlation between acoustic distance and genetic relatedness
The pDFA correctly classified 47.1% of the 160 calls by kin group (pDFA, chance level = 26.7%, p = 0.03). Figure 3 shows the separation of the kin groups produced by the frequency and time components (classification table produced by a non-permutated DFA is presented in Additional file 2). In addition, we found a statistical trend for a weak, negative correlation between genetic relatedness and acoustic distance among the 16 females (Mantel Test, g = -1.61, Z = 4.61, r = -0.13, p = 0.058, Figure 4). Thus, an increase in relatedness was associated with a tendency towards a decrease in acoustic distance.
We found moderate evidence for matrilineal signatures in mouse lemur agonistic calls. While the calls were classified to the correct matriline at a rate significantly higher than chance (47% correct vs. 26.7% chance), the false classification was still made more than half the time. In addition, while we found a statistical trend for a negative relationship between genetic relatedness and acoustic distance, the correlation coefficient was relatively low.
Given that the dispersed matrilineal social system of mouse lemurs provides the opportunity for matrilineal signatures to occur through both inherited traits in vocal morphology (see source filter theory, i.e., [51–53]) and through social learning [42, 50], both may have been important proximate mechanisms for the moderate signatures found in this study. Offspring may inherit vocal tract morphology affecting vocal cord length and thickness which would in turn affect the fundamental frequency of the calls [51–53]. And indeed, fundamental frequency parameters were important in this analysis, loading highly on principle component 1. Offspring could also inherit traits having to do with lung capacity which could influence breathing rate and thus affect acoustic parameters such as call duration and inter-call interval (see source filter theory, i.e. [51–53]), both of which loaded highly on component 2. Unfortunately, it was not possible to collect measurements of vocal tract morphology (ie, length and thickness of vocal folds, length of vocal tract) as this would be highly invasive and, in the field, very complicated. Thus, it is not possible for us to test for a relationship between vocal tract morphology and acoustics. Furthermore, we do not expect less targeted morphological measures (ie. body mass, body length, head size, etc) to be useful proxies for heritability of vocal tract morphology, because they will often vary with pregnancy status, age, season, and the availability of sufficient nutrition for the developing females, none of which we could control for in this population of wild mouse lemurs. Additionally, previous reviews have shown that body size tends to correlate with acoustic differences across age and sex classes in monkeys and humans, but within those classes the relationship is less clear [52, 61].
In addition to genetic mechanisms, it is possible that offspring may also learn to produce calls similar to the calls of the matrilineal relatives from the same nest, which they hear during socialization [42, 50]. Prior research has shown that infant mouse lemurs produce highly variable infant calls that then stabilize into the adult form around the time of weaning . This could mean that social learning during development may be crucial for the development of kin signatures. Similar findings have been found in birds  and other mammals (i.e., ). Unfortunately our data do not allow us to separate the effects of social learning and heritability. While co-sleeping promotes familiarity and thus generates opportunities for social learning, we cannot reliably compare co-sleepers with non-co-sleepers. We cannot exclude the possibility that the non-co-sleeping dyads may have co-slept when they were younger, but no longer did during our study. This is particularly likely for 28-09 and 36-11 and for 28-09 and 46-11 who were not observed to co-sleep. However, since 28-09 was not recaptured in the second year of the study, we do not know if she was still alive. If she was still alive when 36-11 and 46-11 were born, she may have co-slept with them until her death. Similarly, 06-09 was not captured the second year of the study when 41-11 was first caught, thus we do not know if both members of this dyad were alive at the same time. In addition, 19-10 and 06-09 were both at least one year old at the start of the study. Thus, they could be a sibling dyad or mother-daughter dyad which was part of a sleeping group which split as the lemurs aged. However, while it is not possible to distinguish between the two mechanisms here (genetics and social learning), we suspect that the two proximate mechanisms are not mutually exclusive and may even have additive effects (though additional interaction effects could also be possible). Thus, we expect that both mechanisms are likely to have contributed to the evolution of the moderate matrilineal signatures present in the calls.
While these matrilineal signatures are statistically present, their weakness brings up several intriguing questions. The first is whether the lemurs would be expected to use them to recognize kin. Prior work on mouse lemurs showed that females responded differently to calls from paternal kin and nonkin when the classification rate was 79% (mating calls), but not when it was 45% (alarm calls) . Though our current study focuses on only the agonistic calls, future work testing other call types for matrilineal signatures would be very valuable. It would contribute to an increasing body of literature which suggests that the strength of acoustic signatures and the corresponding responses from conspecifics varies greatly by call type (i.e., [21, 23, 58, 64, 65]). However, while playback studies like those above focus exclusively on acoustic cues, in the wild kin recognition is a multi-modal process and the relative strength of each of the cue types may vary depending upon the context in which selection is expected to occur. It is possible that weaker signatures might be present in calls typically used at short distances when visual and olfactory cues would also be available [58, 65]. Given that agonistic calls are frequently given during close-range conflicts, and mouse lemurs have not been documented to recruit kin for alliances, it is possible that kin signatures in agonistic calls may not be under strong selection (though see  for a case of nonkin recruitment).
As solitary foragers are thought to be the intermediary link between the solitary ancestral mammalian condition and the more complex, derived forms of gregarious primate sociality , our results suggest that ancestral solitary foragers may also have had moderate matrilineal signatures in their vocalizations. Such signatures, if used for kin recognition, may have been a crucial element of the dispersed social networks from which more complex, gregarious sociality is thought to have evolved in primates. However, if these signatures are not recognized, then it would lend support to an alternative theory of cryptic kin selection  in which kin-based sociality is thought to have evolved from the spatial proximity of kin alone. As mouse lemur females (as well as females of many other species ) are philopatric, they could interact preferentially with kin simply because kin are there, rather than because they discriminate kin and nonkin, and this could then be the foundation from which more complex forms of kin-based sociality evolved . Future work is underway to test for the vocal recognition of matrilineal kin in this population and is expected to facilitate testing hypotheses about the possible influences of kin-biased behavior on the evolution of complex sociality (ie., ).
The presence of both matrilineal and individual signatures in several species (i.e., the gray mouse lemur , bats ) brings up the question of whether different pressures select for kin vs. individual signatures [2, 69, 70]. It could be that there is an optimal amount of divergence between individuals within a kin group which optimizes both types of signatures, thus enabling an individual to be categorized by kin group and be simultaneously individually distinctive within its kin group [12, 71]. Or, it is possible that one of the two levels of distinctiveness is more strongly targeted by selection and that the other is merely a side effect of inherited vocal tract morphology and social learning of call production . For example, if individual distinctiveness is highly selected for, how distinctive an individual could be might be constrained by inherited vocal tract morphology and socially learned call production . Alternatively, if kin group distinctiveness is highly selected for, within kin group similarity might be constrained by their individual genetics and individual learning experiences. To tease the two apart, future work should compare the acoustic distances between individuals within kin groups across species with different social systems. Future work should also examine whether females use these moderate signatures to discriminate familiar kin, unfamiliar kin, and familiar nonkin. If only familiar kin are recognized, then it would suggest that the recognition of kin occurs primarily through familiarity with individuals who happen to be kin (see  for work on cooperatively breeding birds). Such future work, testing whether the lemurs actually recognize matrilineal kin will be highly important to determining the biological and evolutionary significance of these signatures.
We found moderate evidence for matrilineal signatures in mouse lemur agonistic calls. In addition, there was a tendency for acoustic distance between individuals to decrease as relatedness increased. We expect that both inherited morphological traits and social learning are proximate mechanisms for these signatures. Given that mouse lemurs are solitary foragers, they serve as models for the ancestral solitary foragers that are believed to have been the link between ancestral solitary mammals and derived, more complex forms of sociality in primates . Thus, our results suggest that the ancestral solitary foragers might have had similar, moderate, matrilineal signatures in their calls and we recommend further studies testing whether the lemurs use these calls recognize kin. Such studies would enable further modeling of how kin recognition in ancestral primates might have impacted the evolution of more complex forms of sociality in primates.
Field site and animal housing
This study was conducted at the Ankarafantsika National Park near the Ampijoroa forestry station (16○19’S, 46○48’E) in northwestern Madagascar during the dry seasons (May through November) of 2010 and 2011 in the designated research area of the park called Jardin Botanique A (JBA). Mouse lemurs were trapped in Sherman Live Traps baited with banana, marked with subcutaneously injected, individually distinctive transponders (ID-100, Trovan Small Animal Marking System, Telinject®, Römberg, Germany), and had small (1-2 mm2) ear biopsies taken as genetic samples. Previous generations were already marked (for methodological details regarding trapping and sampling techniques see: [42, 72]). Tissue samples were stored in approximately 1 ml of Queen’s lysis buffer  for up to 7 months at ambient temperature in the field and then at 4°C until extraction (up to 6 years for archived samples collected in prior years) .
A subset of the trapped lemurs were temporarily kept in cages in the forest near the campsite to facilitate recording of vocalizations (total number trapped lemurs =107, total caged lemurs = 45). No lactating females were kept in the cages. Animals were kept either singly in cages of 0.5 m by 0.5 m by 1 m (width x depth x height) or in small groups (two to four animals) in sets of two adjoining cages, each approximately 1 m wide by 1.2 m high by 0.5 m deep. Each set of cages had two passages (0.3 m × 0.2 m × 0.2 m) connecting them. Cardboard cans were provided as nest boxes (one for each lemur) and the cages were furnished with branches for climbing. The lemurs were fed fresh fruit and could be observed catching insects that flew into the cages. They were provided with additional insects as often as possible. Water was available ad libitum. These housing conditions are comparable to those in captive colonies  and no lemurs were injured by the cages or by a cage-mate. Lemurs were released at their capture site after recording was completed (1 night – approximately 2 weeks, mean = 5 nights). Methods were approved by Madagascar National Parks (2010 permits: N102/ 10/MEF/SG/DGF/DCB.SAP/SCBSE, N103/10/MEF/SG/DGF/DCB.SAP/SCBSE, 2011 permits: N101/11/MEF/SG/DGF/DCB.SAP/SCB, N102/11/MEF/SG/DGF/DCB.SAP/SCB) and the Arizona State University Institutional Animal Care and Use Committee (Protocol: 10-1077R).
Before release, 25 adult female gray mouse lemurs (2010: n = 13, 2011: n = 15, three collared in both years, adult = 50 g) were fitted with a radio-collar (either a PicoPip or a Pip3 collar from BioTrack Ltd., United Kingdom, weight 2.3-3.1 g). We then used radio-telemetry to locate the females’ daytime sleeping sites using a TR-4 receiver (Telonics, Mesa, AZ, USA). We read the transponders of lemurs sleeping inside the nests with a handheld microchip reader (Trovan Small Animal Marking System, Telinject®, Römberg, Germany). We checked the sleeping sites on a total of 118 days (65 days in 2010, 53 days in 2011), which resulted in a range of 11-74 days of data per collared female (mean = 29 days), depending on the lifespan of the radio-collar and survival of the female.
Recording methods and acoustic measurements
We recorded all calls given during controlled social encounters when two lemurs were introduced within the cages (or, occasionally during coincidental encounters when a free-ranging lemur outside the cage approached). The introductions inside the cages were observed and the elicited calls were considered to be agonistic when they were associated with aggressive/defensive behaviors such as fighting, chasing, fleeing, etc. When the lemurs were first introduced, the experimenter remained present during the entire night so that she could separate the lemurs if necessary. However, this was rarely necessary, and no lemurs were injured during the introductions.
We recorded the calls with a D1000X Bat Detector (flat frequency response: 5-235 kHz, sampling frequency 200 kHz, 16-bit resolution, Pettersson Elektronik, Upsala, Sweden) from a distance of approximately 2-4 meters from the inside of an observation tent. Under these conditions, agonistic calls were recorded from 15 female gray mouse lemurs. None of the lemurs were caged with female kin group members when the calls were recorded. For one additional female, calls were recorded at a distance of approximately 3 meters while she ate at a feeding platform in the forest after a conflict with another lemur.
Calls were measured in Signal 4.0 (Engineering Design) using the macro written by M. Scheumann for agonistic gray mouse lemur calls and previously used in Leliveld et al. . Ten high quality calls were selected from each female. High quality calls were those that had a clearly visible fundamental frequency, low background noise, and no overlaps with other sound-producing organisms. As the calls are typically given in series, we selected 2-3 series per lemur. Each series consisted of 2-7 calls for a total of 10 calls for each of the 16 lemurs. Figure 1 and Table 4 provide a description of the acoustic parameters that were measured or calculated.
Genetic analyses were conducted at the University of Veterinary Medicine Hannover in the Institute of Zoology. Extractions were performed with a proteinase K digestion and a phenol / chloroform extraction. Eight microsatellite loci (Table 5) were successfully amplified using one of three methods: 1) We used a Qiagen Multiplex PCR Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions, but reduced to the final reaction volume to 10 μl. Ratios followed the instructions with the exception that only 1 μl of Q Solution was used. Cycling conditions followed the provided protocol with annealing temperatures of 48-58°C and up to 48 cycles. 2) We used a MyTaq DNA Polymerase kit (Bioline GmbH, Luckenwalde, Germany) following the manufacturer’s instructions and concentrations, but reducing the reaction volume to 10 μl and using 0.15 μM of each primer and 0.05 μl MyTaq. 3) We performed PCR reactions with final concentrations of 1.5-2.0 mM MgCl2 Solution (Invitek , Berlin, German), 1 x NH4-reaction buffer (50 mM Tris-HCL (pH8.8), 16 mM (NH4)2SO4, 0.1% Tween ©20, Invitek, Berlin, Germany) or 1 x PARR buffer (Cambio, Cambridge, UK), 225 μM of each dNTP (Fermentas Life Sciences), 0.15-0.5 μM of each primer, and 0.025 U of Taq DNA Polymerase. Cycling conditions for this procedure and the MyTaq kit consisted of an initial denaturation phase of 2-4 min at 92-94°C, denaturing for 20-60s at 92-94°C, annealing for 20-60s at 48-58°C, extension for 30-90s at 72°C, and a final extension phase of 5-7 min at 72°C. We used up to 48 cycles. For one marker, M3, the cycling conditions were as follows: initial denaturing at 94°C for 4 min, denaturing at 94°C for 30s, annealing at 55°C for 20s, extension at 72°C for 30s (6-7 cycles), denaturing at 94°C for 30s, annealing at 53°C for 20s, extension at 72°C for 30s (6-7 cycles), denaturing at 94°C for 30s, annealing at 50°C for 20s, extension at 72°C for 30s (25-30 cycles), and a final extension phase at 72°C for 7 min.
The length of the resulting PCR products were determined on an Applied Biosystems 3500 capillary sequencing machine (Applied Biosystems, Life Technologies, GmbH, Darmstadt, Germany). Alleles were scored in Genemapper 4.1 (Applied Biosystems, Life Technologies, GmbH, Darmstadt, Germany) and checked by eye. All homozygous samples were amplified at least twice, following the procedures in prior studies (i.e., [42, 78]).
The mitochondrial D-loop was sequenced using the universal mammalian control region primers H16498 and L15997  for all captured females. PCR was conducted in a 25 μl reaction volume with the following concentrations: 3 mM MgCl2, 1 × NH4-reaction buffer (50 mM Tris-HCL (pH8.8), 16 mM (NH4)2SO4, 0.1% Tween©20, Invitek, Berlin, Germany), 400 μM of each dNTP (Fermentas Life Sciences), 0.8 μM of each primer, 0.125 U of Taq DNA Polymerase. We used an initial denaturation phase of 3 min at 94°C, a denaturing phase of 1 min at 94°C, an annealing phase of 1 min at 50°C, an extension phase of 1 min at 72°C (35-50 cycles), and a final extension phase of 5 min at 72°C. For samples that did not amplify well and were weak when visualized on a 1.5% agarose gel (containing 1.3 × 10-4 mg/ml ethidium bromide), we used the MyTaq DNA Polymerase kit (Bioline GmbH, Luckenwalde, Germany). We followed the manufacturer’s instructions and concentrations, but reduced the reaction volume to 25 μl and used 1 μl of each primer (10 pM/μl) and 0.1 μl MyTaq. Cycling conditions were the same as above. PCR products were then cleaned using the MSB Spin PCRapace kit (Stratec Molecular GmbH, Berlin, Germany). Sequencing followed one of two procedures. We either mailed the samples to Macrogen Ltd. (http://dna.macrogen.com) where they were sequenced using an ABI 3730XL automatic DNA sequencer or we performed the sequencing reactions ourselves using the ABI Prism BigDye Terminator v. 3.1 Cycle Sequencing Kit (Applied Biosystems, Life Technologies, GmbH, Darmstadt, Germany). We used 10 μl reactions consisting of 6.5 μl cleaned PCR product, 1 μl ABI Prism BigDye Terminator Ready Reaction Mix, 2 μl 5× Sequencing Buffer, and 0.5 μl primer (10 mM) and performed 25 cycles of 96°C for 10s, 57°C for 5 s, and 60°C for 3 min. After a final cleaning step with an ethanol precipitation, subsequent sequencing was performed on an Applied Biosystems 3500 capillary sequencer.
Sequences of 446-563 bp (mean = 531.6) were edited, analyzed and aligned in SeqMan 7.0 (DNASTAR Inc., Madison, WI, USA). The final alignment and a matrix of the number of pairwise differences was calculated in Mega 5 , and a haplotype network was produced in Network 184.108.40.206 (Fluxus Technology Ltd., Suffolk, UK).
Sequences have been deposited in GenBank (Accession numbers: KJ183142-KJ183177).
Of the eight microsatellite markers, one (M21) was not in Hardy-Weinberg equilibrium and displayed a significant deficit in heterozygotes (Fis = 0.155, P = 0.0003, calculated in Fstat 220.127.116.11 ). Because this could influence the relatedness calculations, this marker was dropped from the analysis. The remaining markers and the calculations over all loci were in Hardy-Weinberg equilibrium (Table 5) and were therefore included in the relatedness calculations. While we acknowledge that increasing the number of markers improves the resolution of the kinship relationships , using 7 microsatellites is within the range used in similar studies on mouse lemurs (7 microsatellites in Radespiel et al. , 6 in Wimmer et al. , Radespiel et al. ). In addition, we maximized the genetic information obtained from these microsatellites we used by selecting markers that are highly polymorphic (9-25 alleles, see Table 2) as advocated by Harrison et al. .
We calculated pairwise relatedness in Kinship 1.3.1  according to Queller and Goodnight  based on the genotypes of 107 individuals (72 males, 35 females) that were captured during the study period. We chose this relatedness estimator for two reasons. First, it has been shown to perform well on samples with a high percentage of highly related pairs , which we expected to have, given that we were focusing on co-sleeping females. (Prior research has shown that co-sleeping females are typically closely related , and indeed, in this study, all co-sleeping dyads were closely related). Second, it will allow for comparisons with previous studies on mouse lemur relatedness using this estimator (i.e., [42, 84]). Using Kinship we used a simulation procedure which uses the allele frequencies within the population to test the likelihood that the r-value between each dyad was produced by a relationship of rmaternal = 0.5 and rpaternal = 0 against a null hypothesis of rmaternal = 0 and rpaternal = 0. This was performed for all possible dyads among the 107 individuals. By doing so, we distinguished between dyads with a close matrilineal relatedness and dyads that were matrilineally unrelated. This procedure is based upon Van Horn et al. ’s findings showing that though pairwise relatedness may not be precise enough to distinguish small differences in relatedness (e.g.., full- and half-siblings), unrelated dyads can be accurately distinguished from closely related dyads and vice versa. Van Horn et al.  showed that that closely related dyads are rarely misclassified as unrelated and unrelated dyads are rarely misclassified as closely related (Van Horn et al. , page 1177, Table 1). Kinship’s pairwise relatedness values have been shown to correlate with known pedigree relationships , and negatively with allelic exclusions in this population .
Distinguishing rmaternal = 0.5 and rpaternal = 0 from rmaternal = 0.0 and rpaternal = 0.5 was possible because we integrated the pairwise relatedness data with the mitochondrial haplotype data and co-sleeping data (discussed in greater detail below). The mitochondrial data enabled us to exclude closely related pairs with no matrilineal relationship. While we acknowledge that it could be possible for closely related paternal relatives to have the same mitochondrial haplotype, it is very unlikely that such dyads would also co-sleep as prior research on this population  and others  showed that sleeping groups consist of close matrilineal, not patrilineal, relatives.
The probability of identity between two individuals in the population was <1-6 , calculated according to Botstein et al.  in PopAssign 3.9e (written by S.M. Funk). The probabilities of exclusion, according to Jamieson and Taylor , calculated in PopAssign 3.9e, were 0.999941 for one parent, 1.000000 for the second parent, and 0.998505 in the case of a missing parent. Alpha was set at 0.05 for all statistical tests in this study unless otherwise specified.
Marker characteristics are shown in Table 5. Expected and observed heterozygosity (He and Ho) for each locus and over all loci were calculated in PopAssign 3.9e. The observed Fis for each locus and over all loci and the associated P values testing for a deficit in heterozygotes were calculated in Fstat 18.104.22.168 . The statistical error p was calculated by randomizing alleles among individuals over 7000 randomizations. P values are the proportion of randomizations that gave a larger Fis than the observed. The Bonferroni corrected alpha was set at <0.00714.
Kin group selection
In order to minimize the confounding effects of paternal relatedness when testing for matrilineal signatures, we selected dyads of females within kin groups that had high pairwise relatedness and strong genetic and behavioral evidence of matrilineal relationships. While we realize that the inclusion of full sister dyads (and therefore some cases of paternal relatedness) cannot be excluded with certainty, we assume that due to the promiscuous mating system, possible multiple paternities within litters, and the high turn-over rate of mouse lemurs across field seasons [72, 88, 89], most of our dyads are likely to consist of mother-daughter pairs or half sisters. This would mean that on average, barring severe inbreeding, matrilineal relatedness should be much higher than patrilineal relatedness within the dyads. We grouped the dyads into matrilineal kin groups based upon three criteria. Within a kin group: 1) females had the same mitochondrial haplotype, 2) behavioral evidence showed that they co-sleep, and 3) females had a Queller and Goodnight relatedness value  that is significantly likely to result from a maternal relatedness of 0.5 with all other individuals in the group. In three out of six groups all dyads met all three criteria for kin groups. Within the remaining three groups (groups 1, 5, and 6) not all of the criteria were fulfilled for all dyads (Table 3). Within group 1, female 06-09 was not observed to share a sleeping site with the other females in her group. However, because she shared her mitochondrial haplotype with one of the other females in the group (the third could not be determined), and was closely related to both of the other two females (r = 0.61, P < 0.001 and r = 0.29, P < 0.01), 06-09 is included in the kin group. The mitochondrial haplotype of a second female (19-10) from kin group 1 was unknown, but she shared a nest and had an r-value likely to result from a maternal relatedness of 0.5 (r = 0.44, P < 0.001) with one of the other females in her group. As sleeping groups in this population have been shown to typically consist of close matrilineal relatives , 19-10 is also included in this matrilineal kin group. Within group 5, the two females were not observed to co-sleep, but they fulfilled the other two criteria, including having a significant r value (r = 0.36, P < 0.001) and thus are still considered a kin group. Within sleeping group 6, co-sleeping data is unavailable for two dyads. However, 17-10 co-slept with 28-09 in 2010 and with 36-11 and 46-11 in 2011. It is unknown whether 28-09 lived long enough to have the opportunity to share a sleeping site with 36-11 and 46-11 because she was not recaptured in 2011. In total, we divided the 16 females into 6 kin groups: one group of four females, two groups of three females, and three groups of two females (Table 3).
Test of kin group signatures
In order to test whether agonistic calls are distinctive by kin group, we conducted a discriminant function analysis. We performed a principal component analysis with no rotation on the correlation matrix conducted in SPSS 21 to reduce the dimensionality of the dataset. Then, because we have a nested design (individuals are nested within kin groups), we conducted a permutated linear discriminant function analysis (pDFA) in R 2.14.0 (The R Foundation for Statistical Computing, 2011) with kin group as the test factor and individual as the control factor  and 10,000 permutations. As the maximum number of input parameters is one less than the number of objects in the smallest class (two individuals in some of the kin groups), we could only include one principal component in the analysis . We included the first principal component because it accounted for the greatest amount of variation in the original dataset relative to the other components. Cross-validation was performed using the leave-one-out method (Mundry, R., personal communication). Because the pDFA does not produce a classification table, we present the table produced by a nonpermutated discriminant function analysis conducted in SPSS 21.
Correlation between acoustic distance and genetic distance
We used the first principal component to calculate an acoustic distance for all dyads. First we calculated a mean value for each individual for PCA1. We then calculated the Euclidean distances between each pair of individuals producing a matrix of acoustic distances between the individuals. We conducted a Mantel test in Mantel 2.0  using 1000 permutations to test for a correlation between acoustic distance and genetic relatedness.
Hamilton WD: The genetical evolution of social behavior. I and II. J Theor Biol. 1964, 7: 1-52. 10.1016/0022-5193(64)90038-4.
Rendall D: Recognizing kin: mechanisms, media, minds, modules, and muddles. Kinship and Behavior in Primates. Edited by: Chapais B, Berman CM. 2004, New York: Oxford University Press, 295-316.
Kinship and Behavior in Primates. Edited by: Chapais B, Berman CM. 2004, New York: Oxford University Press
Animal Social Complexity: Intelligence, Culture, and Individualized Societies. Edited by: de Waal F, Tyack PL. 2003, Cambridge: Harvard University Press
The Evolution of Emotional Communication: From Sounds in Nonhuman Mammals to Speech and Music in Man. 2013. Edited by: Altenmüller E, Schmidt S, Zimmermann E. 2013, Oxford: Oxford University Press
Müller AE, Thalmann U: Origin and evolution of primate social organisation: a reconstruction. Biol Rev. 2000, 75: 405-435. 10.1017/S0006323100005533.
Shultz S, Opie C, Atkinson QD: Stepwise evolution of stable sociality in primates. Nature. 2011, 479: 219-224. 10.1038/nature10601.
Briefer EF, McElligott AG: Social effects on vocal ontogeny in an ungulate, the goat, Capra hircus. Anim Behav. 2012, 83: 991-1000. 10.1016/j.anbehav.2012.01.020.
Briefer EF, de la Torre MP, McElligott AG: Mother goats do not forget their kids' calls. P Roy Soc B-Biol Sci. 2012, 279: 3749-3755. 10.1098/rspb.2012.0986.
Townsend SW, Hollen LI, Manser MB: Meerkat close calls encode group-specific signatures, but receivers fail to discriminate. Anim Behav. 2010, 80: 133-138. 10.1016/j.anbehav.2010.04.010.
Blumstein DT, Nguyen KT, Martin JGA: Ontogenetic variation of heritability and maternal effects in yellow-bellied marmot alarm calls. P Roy Soc B-Biol Sci. 2013, 280: 7-
Tyack PL: Convergence of calls as animals form social bonds, active compensation for noisy communication channels, and the evolution of vocal learning in mammals. J Comp Psychol. 2008, 122: 319-331.
Yurk H, Barrett-Lennard L, Ford JKB, Matkin CO: Cultural transmission within maternal lineages: vocal clans in resident killer whales in southern Alaska. Anim Behav. 2002, 63: 1103-1119. 10.1006/anbe.2002.3012.
Whitehead H, Dillon M, Dufault S, Weilgart L, Wright J: Non-geographically based population structure of south Pacific sperm whales: dialects, fluke-markings and genetics. J Anim Ecol. 1998, 67: 253-262. 10.1046/j.1365-2656.1998.00187.x.
Gillam EH, Chaverri G: Strong individual signatures and weaker group signatures in contact calls of Spix's disc-winged bat, Thyroptera tricolor. Anim Behav. 2012, 83: 269-276. 10.1016/j.anbehav.2011.11.002.
Chaverri G, Kunz TH: All-offspring natal philopatry in a Neotropical bat. Anim Behav. 2011, 82: 1127-1133. 10.1016/j.anbehav.2011.08.007.
Scherrer JA, Wilkinson GS: Evening bat isolation calls provide evidence for heritable signatures. Anim Behav. 1993, 46: 847-860. 10.1006/anbe.1993.1270.
Hoffmann F, Musolf K, Penn DJ: Spectrographic analyses reveal signals of individuality and kinship in the ultrasonic courtship vocalizations of wild house mice. Physiol Behav. 2012, 105: 766-771. 10.1016/j.physbeh.2011.10.011.
Latham N, Mason G: From house mouse to mouse house: the behavioural biology of free-living Mus musculus and its implications in the laboratory. Appl Anim Behav Sci. 2004, 86: 261-289. 10.1016/j.applanim.2004.02.006.
Charlton BD, Zhang ZH, Snyder RJ: Vocal cues to identity and relatedness in giant pandas (Ailuropoda melanoleuca). J Acoust Soc Am. 2009, 126: 2721-2732. 10.1121/1.3224720.
Rendall D, Notman H, Owren MJ: Asymmetries in the individual distinctiveness and maternal recognition of infant contact calls and distress screams in baboons. J Acoust Soc Am. 2009, 125: 1792-1805. 10.1121/1.3068453.
Rendall D, Cheney DL, Seyfarth RM: Proximate factors mediating "contact" calls in adult female baboons (Papio cynocephalus ursinus) and their infants. J Comp Psychol. 2000, 114: 36-46.
Rendall D, Rodman PS, Emond RE: Vocal recognition of individuals and kin in free-ranging rhesus monkeys. Anim Behav. 1996, 51: 1007-1015. 10.1006/anbe.1996.0103.
Nunn CL: Maternal recognition of infant calls in ring-tailed lemurs. Folia Primatol (Basel). 2000, 71: 142-146. 10.1159/000021742.
Insley SJ, Phillips AV, Charrier I: A review of social recognition in pinnipeds. Aquat Mamm. 2003, 29: 181-201. 10.1578/016754203101024149.
McComb K, Moss C, Sayialel S, Baker L: Unusually extensive networks of vocal recognition in African elephants. Anim Behav. 2000, 59: 1103-1109. 10.1006/anbe.2000.1406.
Sayigh LS, Tyack PL, Wells RS, Solow AR, Scott MD, Irvine AB: Individual recognition in wild bottlenose dolphins: a field test using playback experiments. Anim Behav. 1998, 57: 41-50.
Kessler SE, Scheumann M, Nash LT, Zimmermann E: Paternal kin recognition in the high frequency / ultrasonic range in a solitary foraging mammal. BMC Ecol. 2012, 12: 26-10.1186/1472-6785-12-26.
Nash LT: Kinship and behavior among nongregarious nocturnal prosimians: What do we really know?. Kinship and Behavior in Primates. Edited by: Chapais B, Berman CM. 2004, New Yor: Oxford University Press, 200-222.
Cartmill M: New views on primate origins. Evol Anthropol. 1992, 1: 105-111.
Cartmill M: Rethinking primate origins. Science. 1974, 184: 436-443. 10.1126/science.184.4135.436.
Cartmill M: Arboreal adaptations and the origin of the order Primates. The functional and evolutionary biology of primates. Edited by: Tuttle RH. 1972, Chicago: Aldine-Atherton, 97-122.
Gebo DL: A shrew-sized origin for primates. Yearb Phys Anthropol. 2004, 47: 40-62. Yearbook of Physical Anthropology]
Sussman RW: Primate origins and the evolution of angiosperms. Am J Primatol. 1991, 23: 209-223. 10.1002/ajp.1350230402.
Rasmussen DT: Primate origins - Lessons from a neotropical marsupial. Am J Primatol. 1990, 22: 263-277. 10.1002/ajp.1350220406.
Rasmussen DT: The origin of primates. The Primate Fossil Record. Edited by: Hartwig W. 2002, Cambridge: Cambridge University Press, 5-9.
Martin RD, Soligo C, Tavare S: Primate origins: Implications of a cretaceous ancestry. Folia Primatol (Basel). 2007, 78: 277-296. 10.1159/000105145.
Charles-Dominique P, Martin RD: Evolution of lorises and lemurs. Nature. 1970, 227: 257-10.1038/227257a0.
Piep M, Radespiel U, Zimmermann E, Schmidt S, Siemers BM: The sensory basis of prey detection in captive-born grey mouse lemurs, Microcebus murinus. Anim Behav. 2008, 75: 871-878. 10.1016/j.anbehav.2007.07.008.
Radespiel U: Ecological diversity and seasonal adaptations of mouse lemurs (Microcebus spp. Lemurs: Ecology and Adaptation. Edited by: Gould L, Sauther ML. 2006, New York: Springer, 211-234.
Primate Origins: Adaptations and Evolution. Edited by: Ravosa MJ, Dagosto M. 2007, Chicago: Springer
Radespiel U, Sarikaya Z, Zimmermann E, Bruford MW: Sociogenetic structure in a free-living nocturnal primate population: sex-specific differences in the grey mouse lemur (Microcebus murinus). Behav Ecol Sociobiol. 2001, 50: 493-502. 10.1007/s002650100402.
Radespiel U: Sociality in the gray mouse lemur (Microcebus murinus) in northwestern Madagascar. Am J Primatol. 2000, 51: 21-40.
Radespiel U, Cepok S, Zietemann V, Zimmermann E: Sex-specific usage patterns of sleeping sites in grey mouse lemurs (Microcebus murinus) in northwestern Madagascar. Am J Primatol. 1998, 46: 77-84. 10.1002/(SICI)1098-2345(1998)46:1<77::AID-AJP6>3.0.CO;2-S.
Kappeler PM: Nests, tree holes, and the evolution of primate life histories. Am J Primatol. 1998, 46: 7-33. 10.1002/(SICI)1098-2345(1998)46:1<7::AID-AJP3>3.0.CO;2-#.
Kappeler PM: Determinants of primate social organization: Comparative evidence and new insights from Malagasy lemurs. Biol Rev Camb Philos Soc. 1997, 72: 111-151. 10.1017/S0006323196004999.
Radespiel L, Lutermann H, Schmelting B, Bruford MW, Zimmermann E: Patterns and dynamics of sex-biased dispersal in a nocturnal primate, the grey mouse lemur, Microcebus murinus. Anim Behav. 2003, 65: 709-719. 10.1006/anbe.2003.2121.
Schliehe-Diecks S, Eberle M, Kappeler PM: Walk the line-dispersal movements of gray mouse lemurs (Microcebus murinus). Behav Ecol Sociobiol. 2012, 66: 1175-1185. 10.1007/s00265-012-1371-y.
Eberle M, Kappeler PM: Mouse lemurs in space and time: a test of the socioecological model. Behav Ecol Sociobiol. 2002, 51: 131-139. 10.1007/s002650100409.
Eberle M, Kappeler PM: Family insurance: kin selection and cooperative breeding in a solitary primate (Microcebus murinus). Behav Ecol Sociobiol. 2006, 60: 582-588. 10.1007/s00265-006-0203-3.
Fitch WT, Hauser MD: Vocal production in nonhuman primates: acoustics, physiology, and function constraints on "honest" advertisement. Am J Primatol. 1995, 37: 191-219. 10.1002/ajp.1350370303.
Ey E, Pfefferle D, Fischer J: Do age- and sex-related variations reliably reflect body size in non-human primate vocalizations? A review. Primates. 2007, 48: 253-267. 10.1007/s10329-006-0033-y.
Owren MJ, Rendall D: An affect-conditioning model of nonhuman primate vocal signaling. Perspect Ethol. Edited by: Owings DH, Beecher MD, Thompson NS. 1997, New York: Plenum Press, 12: 299-346. 10.1007/978-1-4899-1745-4_10.
Braune P, Schmidt S, Zimmermann E: Acoustic divergence in the communication of cryptic species of nocturnal primates (Microcebus ssp.). BMC Biol. 2008, 6: 19-10.1186/1741-7007-6-19.
Buesching CD, Heistermann M, Hodges JK, Zimmermann E: Multimodal oestrus advertisement in a small nocturnal prosimian. Microcebus murinus. Folia Primatol (Basel). 1998, 69: 295-308. 10.1159/000052718.
Scheumann M, Zimmermann E, Deichsel G: Context-specific calls signal infants' needs in a strepsirrhine primate, the gray mouse lemur (Microcebus murinus). Dev Psychobiol. 2007, 49: 708-718. 10.1002/dev.20234.
Zimmermann E: Vocal expression of emotion in a nocturnal prosimian primate group, mouse lemurs. Handbook of Mammalian Vocalizations: An Integrative Neuroscience Approach. Edited by: Brudzynski SM. 2009, Oxford: Academic Press, 215-225.
Leliveld LMC, Scheumann M, Zimmermann E: Acoustic correlates of individuality in the vocal repertoire of a nocturnal primate (Microcebus murinus). J Acoust Soc Am. 2011, 129: 2278-2288. 10.1121/1.3559680.
Goodnight KF, Queller DC: Computer software for performing likelihood tests of pedigree relationship using genetic markers. Mol Ecol. 1999, 8: 1231-1234. 10.1046/j.1365-294x.1999.00664.x.
Queller DC, Goodnight KF: Estimating relatedness using genetic markers. Evolution. 1989, 43: 258-275. 10.2307/2409206.
Rendall D, Kollias S, Ney C, Lloyd P: Pitch (F-0) and formant profiles of human vowels and vowel-like baboon grunts: The role of vocalizer body size and voice-acoustic allometry. J Acoust Soc Am. 2005, 117: 944-955. 10.1121/1.1848011.
Zimmermann E: Ontogeny of acoustic communication in prosimian primates. Primatology Today. Edited by: Ehara A, Kimura T, Takenaka O, Iwamoto M. 1991, Amsterdam: Elsevier Science, 337-340.
Sharp SP, McGowan A, Wood MJ, Hatchwell BJ: Learned kin recognition cues in a social bird. Nature. 2005, 434: 1127-1130. 10.1038/nature03522.
Rendall D, Owren MJ, Rodman PS: The role of vocal tract filtering in identity cueing in rhesus monkey (Macaca mulatta) vocalizations. J Acoust Soc Am. 1998, 103: 602-614. 10.1121/1.421104.
Mitani JC, GrosLouis J, Macedonia JM: Selection for acoustic individuality within the vocal repertoire of wild chimpanzees. Int J Primatol. 1996, 17: 569-583. 10.1007/BF02735192.
Eberle M, Kappeler PM: Mutualism, reciprocity, or kin selection? Cooperative rescue of a conspecific from a boa in a nocturnal solitary forager the gray mouse lemur. Am J Primatol. 2008, 70: 410-414. 10.1002/ajp.20496.
Hatchwell BJ: Cryptic kin selection: Kin structure in vertebrate populations and opportunities for kin-directed cooperation. Ethology. 2010, 116: 203-216. 10.1111/j.1439-0310.2009.01732.x.
Widdig A: The impact of male reproductive skew on kin structure and sociality in multi-male groups. Evol Anthropol. 2013, 22: 239-250. 10.1002/evan.21366.
Grafen A: Do animals really recognize kin?. Anim Behav. 1990, 39: 42-54. 10.1016/S0003-3472(05)80724-9.
Tang-Martinez Z: The mechanisms of kin discrimination and the evolution of kin recognition in vertebrates: a critical re-evaluation. Behav Process. 2001, 53: 21-40. 10.1016/S0376-6357(00)00148-0.
Fripp D, Owen C, Quintana-Rizzo E, Shapiro A, Buckstaff K, Jankowski K, Wells R, Tyack P: Bottlenose dolphin (Tursiops truncatus) calves appear to model their signature whistles on the signature whistles of community members. Anim Cogn. 2005, 8: 17-26. 10.1007/s10071-004-0225-z.
Radespiel U, Dal Secco V, Drogemuller C, Braune P, Labes E, Zimmermann E: Sexual selection, multiple mating and paternity in grey mouse lemurs, Microcebus murinus. Anim Behav. 2002, 63: 259-268. 10.1006/anbe.2001.1924.
Seutin G, White BN, Boag PT: Preservation of avian blood and tissue samples for DNA analyses. Can J Zool. 1991, 69: 82-90. 10.1139/z91-013.
Wrogemann D, Radespiel U, Zimmermann E: Comparison of reproductive characteristics and changes in body weight between captive populations of rufous and gray mouse lemurs. Int J Primatol. 2001, 22: 91-108. 10.1023/A:1026418132281.
Radespiel U, Funk SM, Zimmermann E, Bruford MW: Isolation and characterization of microsatellite loci in the grey mouse lemur (Microcebus murinus) and their amplification in the family Cheirogaleidae. Mol Ecol Notes. 2001, 1: 16-18.
Hapke A, Eberle M, Zischler H: Isolation of new microsatellite markers and application in four species of mouse lemurs (Microcebus sp.). Mol Ecol Notes. 2003, 3: 205-208. 10.1046/j.1471-8286.2003.00398.x.
Wimmer B, Tautz D, Kappeler PM: The genetic population structure of the gray mouse lemur (Microcebus murinus), a basal primate from Madagascar. Behav Ecol Sociobiol. 2002, 52: 166-175. 10.1007/s00265-002-0497-8.
Radespiel U, Juric M, Zimmermann E: Sociogenetic structures, dispersal and the risk of inbreeding in a small nocturnal lemur, the golden-brown mouse lemur (Microcebus ravelobensis). Behaviour. 2009, 146: 607-628. 10.1163/156853909X426372.
Guschanski K, Olivieri G, Funk SM, Radespiel U: MtDNA reveals strong genetic differentiation among geographically isolated populations of the golden brown mouse lemur, Microcebus ravelobensis. Conserv Genet. 2007, 8: 809-821. 10.1007/s10592-006-9228-4.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Goudet J: Fstat version 1.2: A computer program to calculate Fstatistics. J Hered. 1995, 86: 485-486.
Harrison HB, Saenz-Agudelo P, Planes S, Jones GP, Berumen ML: Relative accuracy of three common methods of parentage analysis in natural populations. Mol Ecol. 2013, 22: 1158-1170. 10.1111/mec.12138.
Csilléry K, Johnson T, Beraldi D, Clutton-Brock T, Coltman D, Hansson B, Spong G, Pemberton JM: Performance of marker-based relatedness estimators in natural populations of outbred vertebrates. Genetics. 2006, 173: 2091-2101. 10.1534/genetics.106.057331.
Rode JE, Nekaris A, Markolf M, Schliehe-Diecks S, Seiler M, Radespiel U, Schwitzer C: Social organisation of the northern giant mouse lemur Mirza zaza in Sahamalaza, north western Madagascar, inferred from nest group composition and genetic relatedness. Contrib Zool. 2013, 82: 71-83.
Van Horn RC, Altmann J, Alberts SC: Can't get there from here: inferring kinship from pairwise genetic relatedness. Anim Behav. 2008, 75: 1173-1180. 10.1016/j.anbehav.2007.08.027.
Botstein D, White RL, Skolnick M, Davis RW: Construction of a genetic-linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet. 1980, 32: 314-331.
Jamieson A, Taylor SS: Comparisons of three probability formulae for parentage exclusion. Anim Genet. 1997, 28: 397-400. 10.1111/j.1365-2052.1997.00186.x.
Eberle M, Kappeler PM: Sex in the dark: determinants and consequences of mixed male mating tactics in Microcebus murinus, a small solitary nocturnal primate. Behav Ecol Sociobiol. 2004, 57: 77-90. 10.1007/s00265-004-0826-1.
Eberle M, Kappeler PM: Selected polyandry: female choice and inter-sexual conflict in a small nocturnal solitary primate (Microcebus murinus). Behav Ecol Sociobiol. 2004, 57: 91-100. 10.1007/s00265-004-0823-4.
Mundry R, Sommer C: Discriminant function analysis with nonindependent data: consequences and an alternative. Anim Behav. 2007, 74: 965-976. 10.1016/j.anbehav.2006.12.028.
Liedloff AC: Mantel Nonparametric Test Calculator Version 2.0. School of Natural Resource Sciences. 1999, Australia: Queensland University of Technology
We would like to thank the Malagasy government, Ministère De L’Environnement et des Forêts, and Madagascar National Parks (MNP) for their permission to work in the Ankarafantsika National Park and for their support during our fieldwork. We are grateful to Solofonirina Rasoloharijaona, Blanchard Randrianambinina and the Faculté des Sciences, Université de Mahajanga for valuable training and logistical support in the field, Jhonny Kennedy and Jean de la Croix for guide services, Hella Breitrueck for assistance in the genetics lab, Dr. Marina Scheumann for sharing the macro for the acoustic analysis, Soenke von den Berg and Philipp Hohenbrink for technical support. This research was funded by an NSF Dissertation Improvement Grant (#0961779), PEO Scholar Award, the Animal Behavior Society, Lewis and Clark Fund of the American Philosophical Society, German Academic Exchange Service (#A/09/81743), American Society of Primatologists, Sigma Xi (National Chapter, grant #G2009101504), Sigma Xi (Arizona State University chapter), Arizona State University Graduate and Professional Student Association, Arizona State University School of Human Evolution and Social Change, an Arizona State University Graduate College Dissertation Writing Fellowship, and the Institute of Zoology at The University of Veterinary Medicine Hannover.
The authors declare that they have no competing interests.
SK designed the study, collected and analyzed the data and drafted the manuscript. UR participated in study design, supervised the genetic analyses, and assisted with data analysis and manuscript preparation. AH and LL assisted with data collection in the field. LN participated in study design, data analysis and manuscript preparation. EZ participated in study design, supervised the acoustic analyses, and assisted with data analysis and manuscript preparation. UR and EZ also provided logistical support for the fieldwork, and equipment and materials for field and lab work. All authors read and approved the final manuscript.