Habitat quality affects the social organization in mouse lemurs (Microcebus griseorufus)

Species can either maintain a certain social organization in different habitats or show different social organizations in similar habitats. The reasons underlying this variability are not always clear but might have consequences for population dynamics, especially under changing environmental conditions. Among mammals, the primate genus Microcebus lives in small groups of closely related females, derived from female philopatry and dispersed males, as illustrated by the well-studied Microcebus murinus. Here, we studied the genetic structure of a population of the congeneric Microcebus griseorufus, inhabiting three adjacent habitats with different resource availabilities. In order to learn more about the plasticity of the species’ social organization under these different conditions, we analyzed the spatial arrangement of mitochondrial haplotypes of 122 individuals. The study revealed high haplotype diversity and a pronounced difference in spatial distribution between the sexes. Females exhibited spatial aggregation of haplotypes, suggesting a system of female philopatry and matrilines, similar to M. murinus. Male haplotypes were dispersed, and males were more likely to carry rare haplotypes, indicating higher dispersal activity. These findings hint towards the unity of the social organization across the genus Microcebus, suggesting a phylogenetic origin of the social organization. Yet, with decreasing resources, the clustering of female haplotypes declined and approached a random distribution in the marginal habitat, with cluster sizes correlating with resource availability as predicted by the socioecological model. Our study supports the notion that social organization is shaped by both phylogenetic origin and ecological conditions, at least in these small primates. Impacts of habitat degradation are mostly described in terms of changes in population densities in relation to the reduction of resources. This neglects the possible effects of altered social organizations due to declining resources or population densities. Using a genetic sampling of three subpopulations of mouse lemurs in Madagascar along a gradient of food availability up to the limit of the species’ ecological tolerance, we show that their social organization consisting of spatial clusters of closely related females and overdispersed males converges towards random spatial distributions of both sexes with declining food availability.


Introduction
The vast majority of nonhuman primate populations are declining due to the degradation and destruction of their habitats, hunting, and climate change (Estrada et al. 2017). Studies on habitat degradation focus on identifying resources that are important for the persistence of species in degraded habitats. In most cases, "importance" is defined through correlations of specific resource densities with population densities of primates as well as of other taxa (e.g., Johns and Skorupa 1987;Chapman et al. 2000;Irwin et al. 2010;Peres et al. 2010;Sodhi et al. 2010;Steffens et al. 2023). Conclusions about the viability of populations based on the focus on resources neglect behavioral responses and consequences of habitat change on social systems and life histories that could either provide solutions or reflect constraints, leading to intra-specific variation in social organization due to social flexibility (Lott 1984;Schradin 2013;Schradin et al. 2018;Ozgul et al. 2023). The socioecological model for the evolution of different social systems, made popular by Krebs and Davies (1978) and refined continuously thereafter, has provided a useful framework for the study of the social and mating systems of animals. The basic assumption is that females and males evolved social systems that maximize their reproductive success, sometimes with diverging sex-specific interests (Trivers 1972;Cassini 2021). Due to their impressive diversity of social systems and a large number of case studies, primates have often served as examples to test predictions of the model (Crook 1964;Emlen and Oring 1977;Wrangham 1980;Terborgh and Janson 1986;Dunbar 1988;van Schaik 1989;Goldizen 1990;Sterck et al. 1997;Kappeler and van Schaik 2002;Koenig and Borries 2009;Clutton-Brock and Janson 2012;Koenig et al. 2013). The socioecological model uses ecological parameters such as resource distribution in space and time to explain the social organization such as spatial distributions and group size. Simplified, the spatial distribution and relationships of females are determined by food availability with associated within-group and between-group feeding competition and risk avoidance, while male distribution and relationships are mainly based on the distribution of mating opportunities (van Schaik 1983;Dunbar 1988;Janson 1992). Where resources are not worth defending, females can be expected to disperse. Males are expected to compete among themselves for receptive females and should adjust their mating strategy to optimize access to receptive females. If females are spatially clumped and exhibit non-synchronized estrus cycles, the chances of monopolization of females are high. When females are dispersed, and their estrus cycles are synchronized, the potential to monopolize females is reduced, leading to dispersed males (Ims 1988;Schwagmeyer 1988). While the socioecological model has provided a useful framework for studies on the driving forces of sexual selection, many findings are inconsistent with the model's prediction, arguing that other factors or phylogenetic constraints might be more important than ecological ones and that optimization has to be considered a concept rather than a state that can be achieved (Pulliam and Caraco 1984;Brockman 2001;Koenig and Borries 2009;Clutton-Brock and Janson 2012;Koenig et al. 2013;Strier 2017;Fischer et al. 2019).
Evolutionary processes and aspects of social organization, such as dispersal and breeding strategies, affect the genetic structure of a population (Hamilton 1964;Wright 1965;Wilson 1975;Ross 2001). Because of this association between social organization and genetic structure, information on the genetic structure is essential in helping us understand social behavior and analyze social systems (Kappeler and van Schaik 2002). Genetic data such as measurements of genetic similarity can provide insights into social systems that are not attainable through observational methods alone , especially in solitary and nocturnal species, whose social structures are difficult to observe.
The different species of mouse lemurs (genus Microcebus spp.) are nocturnal solitary foragers (Kappeler et al. 2022), whose females are hypothesized to exhibit spatial philopatry with the males dispersing. Female philopatry results in spatial clusters of closely related females with highly overlapping home ranges (Radespiel et al. 1998(Radespiel et al. , 2001Wimmer et al. 2002;Génin 2008Génin , 2010Schliehe-Diecks et al. 2012;Génin and Rambeloarivony 2018). Yet, the genus consists of 20+ species, some with very small distributional ranges. These small ranges suggest that the different species differ in specific traits, and question whether or not all species can be assigned to the same ecological and behavioral categories. Studies of different aspects of their social organization are restricted to a handful of species and the genetic structure of their social organization has been described in detail only for Microcebus murinus (Radespiel et al. 2001;Wimmer et al. 2002;reviewed by Kappeler et al. 2022). Different species of Microcebus living in sympatry with M. murinus, share some of the characteristics of the social organization of M. murinus but deviate in others (Microcebus berthae: Dammhahn and Kappeler 2005; Microcebus ravelobensis: Weidt et al. 2004). The differences between sympatric species and the small geographic ranges of others raise the question of whether the genetic correlates described for the social organization of M. murinus can be extrapolated to other species of Microcebus and whether or not the organization is altered intra-specifically by different environmental conditions. The Grey-brown Mouse Lemur, Microcebus griseorufus (Kollman, 1910), occurs in the spiny forest in the south and southwest of Madagascar (Rasoazanabary 2004;Mittermeier et al. 2010). The vegetation types inhabited by the species range from dry forest resembling a real forest with a closed canopy during the wet season, to the spiny forest and xerophytic bush with few trees and discontinuous vegetation cover (Moat and Smith 2007;Ratovonamana et al., 2011), the latter seemingly representing the dry limit for the species (Bohr et al. 2011;Steffens et al. 2017). In the dry and spiny forest, the species forms population nuclei where related females live in close proximity and defend resources that provide food reliably against other animals (Génin 2008(Génin , 2010Génin and Rambeloarivony 2018). These population nuclei might be equivalent to the spatial clusters of closely related females in M. murinus. According to the socioecological model, the number of females defending these food resources should decline with the declining size of food patches due to increased intra-group competition. Eventually, the patches should become too small to warrant defense. The question is then, whether or not the social organization of clustered females collapses or is maintained, and might have come about in the first place, for other reasons, such as phylogenetic inertia. For the present study, we chose a population of M. griseorufus inhabiting three distinct adjacent vegetation types with different food and tree availability, ranging from "lush" dry forest to xerophytic bush that represents the environmental limit to the distribution of M. griseorufus (Fig. 1). The region is characterized by very high unpredictability of rainfall with severe effects on resource availability (e.g., flowers and fruits) (Dewar and Richard 2007;Ratovonamana et al. 2011;Kasola et al. 2020). The unreliability of flowering and fruiting might make it unprofitable to defend specific parts of an area year-round because it is uncertain whether or not the plants in question will actually provide food. Apart from adding to the understanding of the socioecological model, the study also offers the possibility to understand effects of future habitat changes on the social organization with possible consequences on population dynamics (desiccation and habitat degradation due to climatic or anthropogenic impacts; (Hannah et al. 2008;Tadross et al. 2008;Ozgul et al. 2023).
In this study, we investigate the genetic structure of M. griseorufus in three distinct, adjacent habitats with different food availability. Towards this end, we supplemented the existing database by Scheel et al. (2015) by sequencing the mitochondrial D-Loop of another 74 individuals. Following the published approach to infer the social organization of Microcebus murinus from genetic information (Radespiel et al. 2001;Wimmer et al. 2002), we use the spatial distribution of mitochondrial haplotypes as a proxy for the social organization. First, we wanted to know whether the social organization described for M. murinus also applies to M. griseorufus. To test this, we analyzed and compared the published population genetic data on M. murinus at the center of its geographic range (Radespiel et al. 2001;Wimmer et al. 2002) with our new genetic data from M. griseorufus, living at the dry environmental limits of the genus (Bohr et al. 2011;Steffens et al. 2017).
Second, we tested whether the social organization remains invariant across the environmental gradient from the dry forest to the even dryer xerophytic bush of the species' occurrence or whether it changes in relation to resource abundance.

Study site
This study was conducted in the northwestern part of Tsimanampetsotsa National Park in southwestern Madagascar, ca. 85 km south of Toliara (S 24°01′; E 43°44′). Annual rainfall averages around 400 mm with recurrent droughts without rain for several years. The study area is characterized by two different seasons: eight dry months (April-November) and four wet months (December-March) .
The study site can be divided into three distinct habitat types ( Fig. 1; Table 1), as a result of the topography and edaphic differences (Mamokatra 1999;Andriatsimietry et al. 2009;Rakotondranary et al. 2010;Ratovonamana et al. 2011). These are with declining suitability for M. griseorufus: (1) dry forest on sandy soil (Habitat 1), a 500 m wide area situated between a soda lake to the west and the 40-100 m cliff of the Mahafaly plateau to the east; (2) dry forest on ferruginous soil (Habitat 2), in depressions of the plateau, filled with red sand, and (3) xerophytic bush (spiny bush) on calcareous soil (Habitat 3), covering the slope from the soda lake to the limestone plateau and extending on the plateau towards the east.
Microcebus griseorufus is omnivorous, feeding on gum, fruits, and insects (Génin 2008). Gum is available year-round, provided mostly by large trees, and consumed mainly during the dry season. Fruits are available and consumed especially during the wet season. Abundance data on insects are not available. Food plants have been identified during tracking studies (Bohr et al. 2011;PG unpubl. data). The general vegetation structure of the three vegetation types was described in the study area by a total of 17 30 × 30 m 2 plots, recording all plants (herbs to trees; Ratovonamana et al. 2011). Food plants (all plants > 1 m in height) were counted in each habitat in two 5 × 200 m plots, also used to monitor plant phenology (Bohr et al. 2011;Ratovonamana et al. 2011). The xerophytic bush on calcareous soil differs floristically from the other two habitats on the sand. The three habitats differ in the densities of total food plants [ind./ha] and the density of trees with a diameter at breast height (DBH) ≥ 10 cm [ind./ha]. In the dry forest on sandy soil, the total food plant density is 6032 plants/ha, almost twice as high as that of dry forest on ferruginous soil (3121 plants/ha). 1889 food plants/ha are found in the xerophytic bush on calcareous soil. The largest numbers of trees with DBH ≥ 10 cm were found in Habitat 1 (mean ± standard error: 1029 ± 692) and Habitat 2 (1176 ± 823), while in Habitat 3 only roughly one-third of that was recorded (346 ± 280 trees with DBH ≥ 10 cm). Based on the vegetation structure and the number of food plants, we consider the dry forest on sandy soil (Habitat 1) the most suitable habitat for M. griseorufus, followed by the dry forest on ferruginous soil (Habitat 2) and the xerophytic bush on calcareous soil (Habitat 3).

Microcebus captures
We used capture sites as proxies for the spatial association of individuals. Home ranges measure less than 1 ha with males having larger home ranges than females. Males seem to increase their home ranges during the breeding season, though this observation is anecdotal (Bohr et al. 2011;Génin 2008;PG unpubl. data). Female and male home ranges overlap. Since there is no study on M. griseorufus yet that radio-tracked all individuals within a site simultaneously, we do not know the spatial arrangements of the population members. We use the Euclidian distance between capture sites of individuals to interpret possible grouping patterns as detailed below.
Microcebus griseorufus were captured between April 2007 and March 2009 in 6 ha trapping grids (150 m × 400 m) in each of the three habitats. Traps were spaced regularly at 25-m intervals. The trapping grids were ca. 500 m apart. M. griseorufus lemurs were trapped using 119 Sherman live traps (H. B. Sherman Traps, Tallahassee, FL: 7.5 × 7.5 × 30.5 cm) per grid, equipped with ripe bananas as bait and set at a height of 1-2 m shortly before sunset. Traps were checked before sunrise or in the lactation and weaning season at midnight. Traps were set for four nights per trapping session. Captured mouse lemurs were anesthetized with Ketaminhydrochlorid (Ketamin® 100 mg/ml, Parke-Davis, Berlin) and marked individually either by coded ear clipping or a subcutaneous transponder (Trovan® Passive Transponder System, EURO ID, Identifikationssysteme GmbH and CoKG, Weilerswist, Germany). Small tissue samples were taken from the ear and preserved in 70% ethanol for genetic analyses. Betadine was used to disinfect the mouse lemur skin. None of the recaptured individuals showed signs of infection as a result of the treatment. After examination, animals were kept in their Sherman traps to allow recovery from anesthesia. During that time, they were provided with bananas and water. Animals were released at dusk or at the dawn of the trapping day at their capture sites (Scheel et al. 2015). It was not possible to record data blind because our study involved focal animals in the field. Yet, samples for genetic analyses were only numbered and analyzed by people not involved in the fieldwork (fieldwork by PG, YRR, and SJR; genetic analyses by CA, BMS, and TLL). The two datasets were combined after the completion of the genetic analyses.
Age is difficult to estimate accurately in mouse lemurs since young and older adults cannot be reliably distinguished by their outer appearance outside the mating season (Radespiel et al. 2019). Offspring are born between January and February and are sexually mature after the dry season (ca. July-September). Males in the congeneric M. murinus disperse when they are a few months old (Schliehe-Diecks et al. 2012). Under the assumption that life histories are similar in M. griseorufus, all animals caught after August should have been sexually mature. The body mass of the lightest individual caught after August was 27.5 g; thus, we assume that individuals with a body mass > 27 g were adults. Under this assumption, 13 of the 122 individuals had a body mass between 22 and 27 g and thus were potentially juveniles. Yet, the majority of them were females caught in February and April which are more likely to represent adult females after having given birth or being lactating than being juveniles as they could not have achieved a body mass of > 22 g between birth and the capture date. If they would have been juveniles, their inclusion would make the pattern of diverging haplotype distribution between males and females weaker and thus would represent a conservative error.

Molecular methods and data processing
Genomic DNA was extracted prior to this study from the collected tissue samples with a DNeasy Blood and Tissue Kit (Qiagen) following the standard protocol. In addition to the samples already extracted by Scheel et al. (2015), DNA extraction was performed on 74 samples (NCBI accession numbers: OQ605019 -OQ605092). Polymerase chain reaction (PCR) was carried out according to Scheel et al. (2015) using a standard PCR protocol in a 25 μl reaction mixture containing 17.4 μl of nuclease-free water, 2.5 μl of 10× DreamTaq Green Buffer (Thermo Fisher Scientific, Waltham, MA), 1 μl of dNTP mix (5 mM each; Thermo Fisher Scientific, Waltham, MA), 1 μl of each primer (10 μM), 0.1 μl of DreamTaq Green DNA Polymerase (5U/μl; Thermo Fisher Scientific, Waltham, MA) and 2 μl of template DNA.
A fragment of the mitochondrial D-loop containing the hypervariable Region 1 (HV1) was amplified using the primers TsimMgCytbfw2 (forward; 5′-TCG GAC AAG TGG C-CTC TAT -3′) and mih1coau (reverse; 5′-GTT ATA GTT TCA GGT TAG TCA-3′) (Hapke et al. 2011). DNA amplification of the D-Loop was conducted in a thermocycler (Analytik Jena, Jena, Germany) with the following settings: Initialization at 92.0 °C for 2 min was followed by 35 cycles of denaturation at 92 °C for 40 s, annealing at 55 °C for 60 s and extension at 72 °C for 60 s, with one final elongation step at 72 °C for 5 min.
After PCR, 3 μl of the PCR products was verified on a 1.5% agarose gel and then purified by adding 0.5 μl Exonuclease I (Thermo Fisher Scientific, Waltham, MA) and 1 μl FastAP (Thermo Fisher Scientific, Waltham, MA) to 5 μl PCR product. The reaction mixture was put into a thermocycler with the following conditions: 37 °C for 15 min followed by 85 °C for 15 min. Sequencing of the reverse strand was done by Macrogen (Amsterdam, Netherlands).
Sequence processing and alignment were done using Bioedit (Hall 1999), the software package Geneious Version 10.2.6 (Biomatters Limited, Auckland, New Zealand) (Kearse et al. 2012), and the integrated MUSCLE algorithm (Edgar 2004). The received sequences were checked for contamination and clipped to a common length of 330 bp in order to obtain a complete alignment. Scheel et al. (2015) have published D-Loop sequences from additional samples of the same collection effort, which allowed for the inclusion of 48 additional sequences in all analyses (NCBI accession numbers: KP793607 -KP793669), creating a final alignment containing 122 samples (see Supplementary Information Table S1 for a complete list of samples).

Clustering of haplotypes
Haplotypes were inferred and a median-joining haplotype network (Bandelt et al. 1999) with ε = 0 was created using PopArt (Leigh and Bryant 2015). We then wanted to measure whether the geographic distance between capture sites among individuals with the same haplotype differed between sexes, thus indicating different degrees of clustering. For this, distances between all individuals of the same habitat and sex were calculated with Microsoft Excel using the trapping coordinates. We then calculated the mean distance between all individuals per habitat. Then, for each individual, the mean distance between the individual in question and individuals with the same haplotype occurring in the same habitat was calculated. This measure was called "Indi-vidual_Mean." If haplotypes were distributed at random, the Individual_Mean would not differ from the mean distance between all individuals per habitat. If individuals with the same haplotype would be spatially clumped, the Individual_ Mean would be smaller than the mean for all individuals; if individuals with the same haplotype were overdispersed, the Individual_Mean would be larger.
Microcebus griseorufus occurred in different densities in the three habitats, and thus, the distance measures differed between habitats due to the different densities. If there were the postulated differences in the spatial arrangement of haplotypes between habitats and we would use the "raw" means, we would not be able to distinguish between the effect of animal densities and the effect of habitat differences. Therefore, we calculated a standardized index for clustering, named Individual_Distance_by_Grand_Mean (ID). For this, we divided each individual mean distance by the mean distance per habitat. This accounted for habitat-specific differences in lemur densities. Values below 1 indicate lower distance than average (=clustering) and values above 1 indicate overdispersion. Graphics were created and statistical analyses were performed using Microsoft Excel and SPSS (Field 2013). Confidence intervals (95%) for medians were calculated using the bootstrap procedure with 10,000 runs of the Statistics Calculators by Montgomery College (https:// press books. montg omery colle ge. edu/ statc alcs/ chapt er/ boots trap-confi dence-inter vals/).

Haplotype group size
While the analysis of clustering provides a measure of spatial cohesiveness between individuals with the same haplotype, we use the number of individuals per haplotype as an indication of the number of individuals in population nuclei. Since we recorded the positions of animals during their active phase, the resulting measure of "group size" is not meant as a measure of the size of sleeping associations. Rather, it should represent the size of a group that, in the case of females, jointly defends important food resources that provide food reliably (sensu Génin 2008Génin , 2010Génin and Rambeloarivony 2018).

Haplotype diversity
Using PopArt, we inferred 24 different mtDNA haplotypes from the sequence alignment containing 122 samples/ sequences (Fig. 2, see Supplementary Information Table  S2 for haplotype sequences). Among the 64 males, all 24 mtDNA haplotypes were represented, while only 12 of the haplotypes were found among the 58 females. The number of individuals per haplotype was higher in females than in males (Mann-Whitney U test; U = 167, n female = 15, n male = 36, p = 0.034). All rare haplotypes (N=6), only represented by one individual, belonged to male mouse lemurs. The spatial arrangement of capture sites of the 122 sequenced mouse lemurs is illustrated in Fig. 3 for all three habitats.

Spatial distribution
Plotting the exact trapping coordinates for each vegetation type and sex shows clear differences in spatial distribution between the sexes (Fig. 4). While females show a clear pattern of haplotype clustering in Habitats 1 and 2 and less so in Habitat 3, males exhibit a more heterogeneous distribution of haplotypes, indicating increased dispersal. Female mouse lemurs inhabiting the dry forest on sandy soil (Habitat 1) clustered into five different haplotypes representing a total of 26 individuals, with two large central clusters containing 11 and 9 individuals, respectively. Males inhabiting Habitat 1 exhibited a total of 14 different haplotypes, represented by 26 individuals. Mouse lemurs inhabiting the dry forest on ferruginous soil (Habitat 2) exhibit a similar pattern. Female individuals in this habitat exhibited four different mtDNA haplotypes, but nearly all of them (N=15/17) belonged to only one of two haplotype clusters. The other two (haplotypes Page 7 of 16 65 5 and 17) were represented by single individuals. The 22 males observed here exhibited a total of 13 different haplotypes, which were heterogeneously dispersed across the whole transect. In the driest habitat, the xerophytic bush on calcareous soil (Habitat 3), females exhibit less clustering, constituting six different haplotypes represented by 15 individuals, with the two largest clusters representing five and four individuals, respectively. In the same habitat, nine mtDNA haplotypes were represented by 16 males.

Clustering
In order to compare the clustering of mouse lemurs between sexes and among habitats, individual aggregation was characterized by the Individual_Distance_by_Grand_Mean (ID; Fig. 5). Values below 1 indicate increased clustering while values above 1 indicate increased overdispersion. In all three habitats, females exhibited median ID values below 0.7; however, their significance differed according to habitat type. In the two forest habitats on sand (Habitat 1 and Habitat 2), the 95% confidence intervals of the median did not include 1, indicating significant clustering of haplotypes for females (C 95% of the median, lower and upper limits: females: Habitat 1: 0.54-0.69; Habitat 2: 0.43-0.60). The lower and upper limits of the 95% confidence interval of the ID for females in Habitat 3 were 0.13-1.26 and thus indicated no consistent clustering of females in Habitat 3 (xerophytic bush) (Fig. 5). The median ID of females did not differ among the three habitats (Kruskal-Wallis analysis of variance: H = 1.894, df = 2, p = 0.388).
The 95% confidence intervals for male IDs in Habitat 1 and Habitat 2 include 1 and thus indicate no difference from randomness, with a tendency to being overdispersed, i.e., ID values above 1 (C 95% of the median, lower and upper limits: males: Habitat 1: 1.00-1.13; Habitat 2: 0.97-1.15). In Habitat 3, the C 95% of the median ID of male mouse lemurs also include 1 but exhibit a tendency towards clustering with a median ID of 0.86 (C 95% 0.57-1.00). The median IDs of males do not differ between the three habitats (Kruskal-Wallis analysis of variance: H = 4.843, df = 2, p = 0.089).
The degree of clustering differs between females and males in Habitat 1 and Habitat 2 (Mann-Whitney U tests: Habitat 1: U = 624, p < 0.001; Habitat 2: U = 319, p < 0.001). The medians do not differ in Habitat 3 (Mann-Whitney U test: U = 138, p = 0.474). Thus, the distribution pattern (clustering) of females and males is distinct in the forests but converges towards randomness in the xerophytic bush (Fig. 5).

Haplotype group size
In females, the clusters of individuals belonging to the same haplotype in Habitat 1 and Habitat 2 consist of more individuals than in Habitat 3, while the number of individuals per haplotype remains constant in males (Fig. 6). Yet, since the actual number of haplotypes per habitat is small, neither the number of females nor the number of males per haplotype differs significantly among habitats (Kruskal-Wallis analysis of variance: females: H = 1.47; males: H = 0.16, df = 2, p > 0.4 in both tests).

Discussion
Overall, we recovered 24 different mtDNA haplotypes among 122 individuals, which shows a substantial level of genetic diversity in the mitochondrial DNA compared with other species. Wimmer et al. (2002) found 13 mtDNA haplotypes in 85 individuals of M. murinus,  revealed 4 haplotypes among 46 individuals (sampled in 1 year) of the sympatric Mirza coquereli, and Gerloff et al. (1999) detected 5 haplotypes in a population of 36 bonobos. The mtDNA haplotype diversity of M. griseorufus (0.921), calculated according to Saitoh (2021), is higher than that of M. murinus (0.508) and M. coquereli (0.660). While it is difficult to compare different higher taxa, social units, and sampling areas, M. grisoeurufus is very diverse and ranks among the top 15% of the 66 terrestrial mammal species compiled by Saitoh (2021).

Genetic structure and social organization of M. griseorufus
The results of this study revealed that the genetic structure of Microcebus griseorufus is characterized by spatial clustering of females with shared haplotypes and dispersed males. This study detected a pronounced sex Fig. 3 Distribution map of collected samples. All 122 sequenced mouse lemurs are displayed according to their respective capture coordinates across all three habitats. Some circles represent multiple individuals that were caught in the same trap difference in the haplotype distribution of M. griseorufus. The high general diversity in haplotypes is mostly driven by males, which represent the majority of unique and rare haplotypes. Females (within a population) express fewer haplotypes that are generally shared with other females in proximity. There is a significant difference in the number of individuals per haplotype between males and females. This genetic structure is consistent with the results on M. griseorufus in Berenty where the population consists of population nuclei where related females live in close proximity and defend resources that provide food reliably (Génin 2008(Génin , 2010Génin and Rambeloarivony 2018) and the findings on M. murinus with small family units of closely related females (Radespiel et al. 2001;Wimmer et al. 2002) and female philopatry, while high haplotype diversity and lack of spatial clustering of males imply dispersal. Males might have immigrated from other matrilineal clusters outside the study area. Since within subpopulations there are males sharing the most common haplotypes in females, it is possible that not all males disperse immediately after weaning. Alternatively, sons from distant female clusters with the same haplotype could have moved into the vicinity of females with the same haplotype. In other studies, this genetic organization has also been linked to or derived from sleeping associations (Radespiel et al. 2001;Wimmer et al. 2002;Rode (Table 1). But M. griseorufus seeks shelter in bundles of spiny branches that could compensate for the lack of hollows in large trees. While we did not analyze the genetics of sleeping associations, observations from Microcebus ganzhorni provide evidence that sleeping associations have different dynamics than the spatial associations of individuals during their nightly activity time (Lahann 2008). Since our study relied on captures of active individuals and did not reconstruct home ranges or monitored sleeping associations, we cannot resolve this issue.
Patterns of kinship are often invoked to explain the spatial aggregation of females in combination with affiliative and cooperative behavior to increase inclusive fitness (Sherman 1977;Gouzoules 1984;Packer et al. 1991;Moore 1992). The degree of kinship impacts cooperative behavior in anthropoid primates (Morin et al. 1994;Pope 2000;Chapais 2001;Silk 2002), though cooperation between unrelated individuals has also been reported for bonobos and chimpanzees (Goldberg and Wrangham 1997;Gerloff et al. 1999;Mitani et al. 2000). In M. griseorufus, the evolutionary benefit of forming groups of related females could be a better defense of rewarding food resources that result in higher reproductive success of kin (Génin and Rambeloarivony 2018), as Fig. 5 Clustering of female and male haplotypes by individ-ual_distance_by_grand_mean (ID) value, comparing all three vegetation types across both sexes. Values are medians and quartiles; whiskers indicate the 1.5*interquartile range, outliers (circles) are outside the 3rd quartile + 1.5*interquartile range, extreme values (asterisks) are outside the 3rd quartile + 3*interquartile range Fig. 6 Group size per haplotype across sexes and habitats. Shown are medians (horizontal bars), means (x) and quartiles; whiskers indicate the 1.5*interquartile range, outliers (circles) are outside the 3rd quartile + 1.5*interquartile range well as mutual nursing of offspring from related females (Eberle and Kappeler 2006). Alternatively, females avoid the risks of dispersal and increase thus their fitness even without any benefits from defending resources.
These findings on the social organization of M. griseorufus are in line with the research of the well-studied Microcebus murinus, which exhibits female sleeping associations and dispersed males (Radespiel et al. 2001(Radespiel et al. , 2003Wimmer et al. 2002;Fredsted et al. 2004Fredsted et al. , 2005. Microcebus murinus, at the center of its distribution range, shows a similar genetic structure as M. griseorufus, living at the dry environmental limits of the genus (Bohr et al. 2011;Steffens et al. 2017). This study of M. griseorufus might hint towards the unity of social organizations, across a large environmental gradient. This is consistent with behavioral observations of different species of Microcebus from all major forest types in Madagascar. In all species, males have larger home ranges than females (M. murinus : Fietz 1999;M. rufus: Atsalis 2000;M. berthae: Schwab 2000;Schwab and Ganzhorn 2004;Dammhahn and Kappeler 2005;M. ganzhorni: Lahann 2008). But this is where comparative data come to their limits. The composition of sleeping groups has been used to interpret the social organization in a similar way as the spatial clusters of related females. Data are available for M. rufus (Karanewsky and Wright 2015), M. berthae (Schwab 2000;Dammhahn and Kappeler 2005), M. griseorufus (Génin 2008), M. ganzhorni (Lahann 2008), M.ravelobensis (Radespiel et al. 2009), and M. sambiranensis (Hending et al. 2017). Most of these studies described sleeping groups as highly dynamic with changing composition, sometimes involving members of both sexes. The analogy of sleeping groups and clusters of related females may not be wrong but sleeping groups are something different than the spatial clusters documented by animals during their active phase. Though some sleeping groups allow for alloparental care (Eberle and Kappeler 2006;Génin 2008), the evolutionary relevance of sleeping groups is unclear.
It is not uncommon for members of a genus to share a social organization. In Madagascar, all nine species of wooly lemurs (Avahi) live in pairs, and all nine species of the genus Propithecus are organized in groups Lawler and Richard 2022). However, there are also studies indicating a range of different social organizations and social systems within the same species or genus. Varecia shows substantial variation in grouping patterns (Vasey et al. 2022), and most Eulemur species share a social organization across a variety of habitats, but Eulemur rubriventer and Eulemur mongoz form pair-bonded family groups while the other Eulemur species live in larger multi-male multi-female groups (Johnson et al., 2022). Baboons of the genus Papio comprise six closely related species that show ecological plasticity and different social organizations. Four species of baboons exhibit uni-level organization with big groups consisting of multiple males and females, while two other species display multi-level social organizations with one-male-units, parties, and gangs (Barton et al. 1996;Fischer et al. 2019;Zinner et al. 2021).

Plasticity of the social organization along an ecological gradient
The socioecological model predicts female spatial distribution according to ecological factors such as resource availability and risk avoidance (e.g., predation pressure, dispersal risks, disease transmission). The availability of receptive females, in contrast, affects the spatial distribution of males trying to maximize their access to said females. Analyzing M. griseorufus across their natural habitat range in three distinct vegetation types offers the advantage of a pronounced ecological gradient at a small local scale. The habitat data published by Ratovonamana et al. (2011) show a decline in the density of food plants as well as fewer large trees (with tree holes as shelters) available in the xerophytic bush compared to the forest habitats. Assuming that other ecological factors such as predation pressure and risk of dispersal should be similar because of the close proximity of the different habitats, changes in their social organization should largely reflect a response to resource availability. The capture data in these habitats show a corresponding abundance of mouse lemurs to resource availability.
The general social organization of female philopatry with dispersed males is found across all three habitats regardless of the ecological conditions. These results are in line with other mammalian studies providing evidence of fixed social systems in genera or species across different ecological conditions. Even though it was hypothesized that the social organization of hamadryas baboons is an adaptation to the harsh ecological conditions they are facing, Guinea baboons exhibit the same multi-level social organization, despite mostly living in vastly different habitats (Zinner et al. 2021). Thus, some adaptivity in social characteristics might be explained by phylogeny rather than by ecology (Thierry et al. 2000;Menard 2004;Ossi and Kamilar 2006;Thierry 2007;Shultz et al. 2011;Schradin 2013). There are studies finding contrary evidence, however, showing plasticity in social organization as a hypothesized response to ecological change. Koenig et al. (1998) found differences in female social relationships in Hanuman langurs (Semnopithecus entellus) between populations as an assumed adaptation to ecological conditions. This is paralleled in Propithecus diadema, Eulemur collaris, and other lemur species where group size and group cohesion declined within populations with increasing fragmentation and degradation of forests and reduced food patch sizes (Ganzhorn 1988;Irwin 2007;Donati et al. 2011). Eulemur exhibits a variety of behaviors in response to their environment across the genus (Kappeler and Fichtel 2016), while the social organization is more related to phylogenetic distance between species (Ossi and Kamilar 2006). Izar et al. (2012) found that social relationships between female tufted capuchin monkeys (Sapajus) changed between two species in accordance with their food sources, while the mating system seemed to be constrained by phylogenetic inertia. Similarly, macaques (Macaca), as the most geographically widespread primate taxon show the same grouping and dispersal patterns across their 22 species, while displaying unparalleled diversity in behavior and relationships (Thierry et al. 2000;Thierry 2007).
It is to be noted that while the clustering of females with shared haplotypes in M. griseorufus does not significantly differ among habitats, female aggregation is lower in the xerophytic bush and is not significantly different from a random distribution. This suggests that when resources are not defendable or not worth defending, female associations might decrease. At the same time, males tend to be overdispersed in the forest habitats, indicating competition among individuals. While the distribution of males does not differ significantly from randomness according to our measure, their organization changes drastically between forests and the xerophytic shrub (Habitat 3). It might be possible, that the genetically determined social organization of M. griseorufus starts to crumble at the dry edge of their tolerance. The reduced density of mouse lemurs inhabiting the xerophytic shrub is probably a response to the scarcity of food and nesting sites. The spatial clustering of females is more pronounced in the forest ecosystems, with higher availability of food and trees. This pattern matches the predictions of the socioecological model that anticipate a decline in group size with the declining size of available food patches.
While our findings are coherent with Thierry's and Clutton-Brock and Janson's assessment of the socioecological model (Thierry 2008; Clutton-Brock and Janson 2012), we cannot rule out the possibility that M. griseorufus has different population densities in the three habitats that are unrelated to resource availabilities. If so, the observed differences in social organization could be a consequence of differing population densities. Changes of this kind have been reviewed by Schradin et al. (2020). By using several populations of African striped mice with different constellations of resource abundance and population densities, they could show that changes in the social organization were density-dependent but unrelated to food resources. We tried to account for different population densities by introducing the measure of "Individual_Distance_by Grand_Mean" that represents a standardization that accounts for the different population densities. Yet, this option cannot be ruled out.
Some primate behavior, such as feeding competition, is undoubtedly influenced by their ecological surroundings. Nevertheless, other traits might be better explained by phylogeny. This is exemplified by Guinea and hamadryas baboons that live under very different ambient conditions but show very similar social organizations (Fischer et al. 2019;Zinner et al. 2021). It is also to be considered that primate societies might not be optimal (Janson 2000), and the number of social organizations that are possible is underlain by stabilizing selection preventing adaptations needed for better survival in any given environment (Thierry 1997(Thierry , 2007(Thierry , 2008Hemelrijk 2002). This would allow for some diversity in social organization (Di Fiore and Rendall 1994) while restricting adaptations to ecological characteristics (Fleagle and Reed 1996;Thierry 2008).
The present study offers a perspective on how the social organization of a mouse lemur can withstand the ecological pressures of degrading habitats to a certain extent, but eventually, the genetic structure seems to collapse. The consequences of a collapse of the social organization are hard to predict for different species of this genus. Cooperative breeding and nursing offspring of related females improves the reproductive success of M. murinus and M. griseorufus (Eberle and Kappeler 2006;Génin 2008). This option is probably reduced at the dry limits of the species ranges. But since our data are based on year-round captures and did not allow for a finer temporal resolution, we could not exclude the option of seasonal affiliations for raising offspring. At least, M. murinus seems to be able to compensate for environmental changes by accelerating their life histories (Ozgul et al. 2023). Other Microcebus species, such as the sympatric M. berthae, do not seem to have this plasticity and seem to be on their way to extinction (Kappeler et al. 2022). (In conclusion, the social organization of the gray-brown mouse lemur, Microcebus griseorufus, can be characterized by spatial clusters of females sharing the same haplotype while male haplotypes are overdispersed.) The principle of this structure persists across an environmental gradient towards the environmental limit for the species, but the number of individuals per cluster with the same haplotype changes in a similar way as group size is predicted to change by the socioecological model. At the dry limit of their tolerance, when less food and nesting sites are available, this structure seems to collapse, and the spatial arrangement of haplotypes does not deviate from random. For the time being, the consequences of these changes in the social organization are not known.

Declarations
Ethics approval All animal work was conducted according to Malagasy and German guidelines. Our research was conducted in collaboration with the Département Biologie Animale and the Département Biologie et Ecologie Végétale of the Université d'Antananarivo. Authorization to enter Tsimanampetsotsa National Park as well as to capture and handle small mammals were delivered by the Ministère de l'Environement, des Eaux et Forêts et du Tourisme of Madagascar in accordance of Madagascar National Parks (MNP, former ANGAP; permit n° 057/07 issued on March 12, 2007, permit n° 009/08 issued on January 15, 2008, and permit n° 261/08 issued on October 9, 2008). We hereby confirm that our study was conducted in accordance with the recommendations of the Weatherall report "The use of non-human primates in research". The research was approved by the Ethics commission of the Institute of Zoology of the Universität Hamburg.

Conflict of interest The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.