Sex-Biased Habitat Use by Phyllostomid Bats on Riparian Corridors in a Human Dominated Tropical Landscape

Some animal species exhibit sex-specific patterns as an adaptation to their habitats, however, adaptability to a human-dominated landscape is commonly explored without considering intraspecific sexual differences. Differences between males and females lead to a sexual segregation in habitat use. In southern Mexico, we explored sex-specific responses to landscape modification of six common species of phyllostomid bats: Artibeus jamaicensis, A. lituratus, Sturnira lilium, Carollia perspicillata, Glossophaga soricina, and Platyrrhinus helleri using riparian corridors within continuous forest and cattle pastures. Furthermore, we explored sex related responses to vegetation attributes (i.e., tree height and basal area) and seasonality (i.e., wet and dry seasons). Overall, capture rates were significantly skewed toward females and riparian corridors in pastures. Females of G. soricina exhibited a strong positive relationship with greater tree height and basal area. Seasonality was important for A. lituratus and S. lilium females, only. The results indicate a sexual driven response of bats to habitat modification. The high energetic demands of females associated to reproduction could lead to foraging into riparian corridors in pastures. The presence of large trees along riparian corridors in pastures may help maintaining a diverse and dynamic bat community in modified tropical landscapes.


INTRODUCTION
The configuration of tropical landscapes is highly dynamic as a consequence of changes in landuse and cover (Mayaux et al., 2005;Fagan et al., 2013). A common feature in tropical agricultural landscapes is the presence of relict natural vegetation along streams which persists even when being exposed to long-term agricultural practices and small-scale land use changes (Lundy and Montgomery, 2010). Habitat disturbance can alter the spatial arrangement of critical resources for animals within a given landscape, potentially resulting in disrupted demographic patterns (e.g., sexratio, abundance, age classes, etc.) among habitat patches. The degree to which habitat disturbance alters the demography of animal populations, however, should vary as a function of the behavioral plasticity of individual species (Bender et al., 1998;Nupp and Swihart, 2000).
Males and females of the same species may differ in several aspects of their biology (e.g., sexual dimorphism, different thermoregulatory strategies), which may result in sexual segregations of habitat use (Lintott et al., 2014) and/or between sex competition (Lemaître et al., 2014;Benitez-Malvido et al., 2016). Sexual segregation can be broadly categorized into the following types: habitat segregation and social segregation. Habitat segregation occurs where the sexes differ in the use of the physical environment, whilst social segregation when a species tends to form single-sex groups (Wearmouth and Sims, 2008). Habitat segregation hypothesis suggests that inherent sexual differences in reproductive strategies i.e., reproductive energy demands, breeding period and predation risk (Dietz et al., 2006;Nardone et al., 2015;Beerman et al., 2016;Benitez-Malvido et al., 2016) result in females trading off habitat quality in favor of offspring safety (Wearmouth and Sims, 2008).
Most studies on the effects of human-modified landscapes on bat communities have focused at the species level responses, while the potential importance of intra-specific differences are often ignored (Lintott et al., 2014). Bats are an ideal taxon for studying sexual segregation, since sexual dimorphism in bats is rare but sexual segregation is widespread (Senior et al., 2005). Seasonal and maternal sexual segregation have been documented for many bat species (Sgroi and Wilkins, 2010;Encarnação, 2012;Diamond and Diamond, 2014). In many tropical bat species, females in resource-rich habitats roost in groups with few or no males present (harem groups) (Ortega and Arita, 1999;McCracken and Wilkinson, 2000;Altringham, 2011). Particularly due to reproductive and parental costs, females have higher energy requirements, so they are less abundant in habitats with limited food resources Ramos Pereira et al., 2010). Sex should be considered separately whenever possible in the study of bats because males and females of the same bat species may have different seasonal distribution and roosts with different characteristics (Broders et al., 2006;Safi et al., 2007;Weller et al., 2009). For instance, the response of two Neotropical frugivorous bats to local and landscape scale attributes were sex and seasonally specific; females were more abundant than males in edge and matrix habitats, and females seem to increase their foraging movements during pregnancy and low fruit availability (Rocha et al., 2017). In another study, bats showed sexual differences in the habitat use within urban landscapes, with males being more widely distributed and females more abundant in highly connected areas. Moreover, access to water was a limiting factor in determining female distribution (Lintott et al., 2014;Patriquin et al., 2019). The importance of fine-scale spatiotemporal and demographically precise data is essential for effective conservation strategies (Hutson et al., 2001;Russo et al., 2010;van Toor et al., 2011).
Habitat loss and fragmentation are important threats to bat populations as they eliminate or reduce suitable foraging habitats and forest structures for roost (Kingston, 2010). Information on the abundance and sex ratios of bat populations throughout the year is important for understanding their ecology in periods of resource scarcity (Perry et al., 2010). Therefore, obtaining sex-specific information on the behavior and habitat requirements of bats should be one of the primary goals in conservation efforts (Weller et al., 2009;Perry et al., 2010). In order to understand the mechanisms by which some bat species are affected by habitat loss, it is necessary to determine not only if habitat disturbance affects life-history parameters, but also if habitat loss generates changes in their social structure. In this study we assessed six common species of phyllostomid bats to identify sex-related patterns within a human-dominated landscape. Our study provides insights into the importance of habitat type on sex ratio and on sex distribution throughout the year. The objective of this study was to determine if habitat affects bat sex ratio in a human dominated landscape in Southern Mexico. For this, we sampled individuals from six abundant bat species in conserved continuous forests and cattle pastures along and away from riparian corridors. We expected that because riparian corridors provide food, water and roosts, capture rates of females would be greater along them (Naiman et al., 2000). We hypothesized that sex specific differences in habitat use will be caused by the reproductive energetic demands in females (i.e., pregnancy and lactation). At the local scale, vegetation attributes such as tree height may affect female abundance, because females are frequently restricted to high-quality habitats for foraging (Lintott et al., 2014).

Study Area
The study was carried out in the tropical region of Lacandona, south of the state of Chiapas, Mexico. The original vegetation consists mainly on lowland tropical rain forests. Deforestation of the region began in the 1970's, resulting in the reduction of oldgrowth continuous forest from 95% in 1976to 56% in 1996(De Jong et al., 2000; only 36% of the original old-growth continuous forest remains today (Carabias et al., 2012). Currently, the main land-use practices in the region consist of cattle pastures, the cultivation of maize and other crops (De Jong et al., 2000;Zermeño-Hernández et al., 2015). The resulting landscape comprise a mosaic of human-modified habitats that include semiurban settlements, agricultural land, open pastures for cattle, riparian zones, patches of secondary and old-growth forests of various sizes. The region has a mean annual temperature of 24 • C; average annual rainfall is 3,000 mm with June to October as the wettest months (551 mm month −1 ) and February to April as the driest months (<100 mm month −1 ) (Comisión Federal de Electricidad, 2006;van Breugel et al., 2006).
Four different habitat types were selected for this study including the following: (i) riparian habitat within old-growth continuous forest (RM); (ii) riparian habitat in active cattle pastures (RP); (iii) old-growth continuous forest 1,000 m away from riparian vegetation (MF); and (iv) active cattle pastures (P) 1,000 m away from riparian vegetation. For a total of 12 sampling sites. Each habitat type was replicated three times and study sites were at least 1.5 km away from each other (Figure 1). Streams were all permanent (although with variable amounts of running water throughout the year) while stream width varied from 2 to 8 m. Study sites in pastures were active cattle pastures, although these sites were predominantly devoid of a tree cover, all sites had isolated trees and a few trees that serve as live fences. It is common that isolated trees are left standing to provide shade for cattle, firewood for cooking, or for aesthetic reasons (Galindo-González et al., 2000). The cattle pastures were located in the Marqués de Comillas municipality, on the south side of the Lacantún River. Old-growth continuous forest sites were located in the 330,000 ha Montes Azules Biosphere Reserve (MABR) on the north side of the river (16 • 04 N to 90 • 45 W; INE, 2000, Figure 1).

Bat Sampling
Sampling of bats was conducted twice during the dry season (December to May) and twice during the wet season (June to November) for 3 consecutive years (2011, 2012, and 2013), using a standardized method of four nights per site. A previous study in the same locations showed that over 70 sampling nights, 34 during the dry season and 36 during the rainy season, a total of 1,752 individuals belonging to 28 species of Phyllostomidae were captured (de la Peña-Cuéllar et al., 2015). For this study however, we considered the six most abundant species of phyllostomid bats in the region, including the following: Artibeus jamaicensis, A. lituratus, Sturnira lilium, Carollia perspicillata, Glossophaga soricina, and Platyrrhinus helleri, the minimum number of captures needed to be considered in the study were one capture per habitat per season, species abundance cut off of n = 60.
Five nets (12 m long × 2.6 m high) were set at ground level and were opened at dusk (1800-1830) for four consecutive hours, which corresponds to the peak foraging time for most phyllostomid species (La Val, 1970). The bat sampling nets were arranged according to habitat type: (1) in the riparian habitats (including continuous forest and pastures), nets were located parallel and/or diagonally across the stream, depending on site characteristics; (2) in continuous forest, nets were positioned across natural flying corridors; (3) in active pastures, nets were located under the canopy of isolated trees. In all sites we searched for similar physical characteristics that allow the same mist net arrangement, this is one individual net and two pairs of nets in an "L" shape (two nets connected perpendicularly). Nets were located ca. 50 m apart. Nights with a full moon or heavy rain were avoided during bat sampling in order to prevent variation in capture success associated with these conditions (Morrison, 1978). Captured individuals were temporarily stored in cloth bags and identified to species following Medellín et al. (2011). For all captured bats, we determined sex by inspecting genitalia . In females, we detected pregnancy by palpation , and lactation by the occurrence of enlarged nipples surrounded by a hairless skin area and by extruding milk with a gentle finger pressure on the nipple base. Sex ratio was calculated as the ratio of males to females in each site (Russo et al., 2010).
There is a risk of exposure to some significant zoonotic agents for any person handling bats in the field, following biosecurity recommendations (Newman et al., 2011). For this, all participants that sampled bats were appropriately trained to handle bats; bites and scratches were avoided by using leather gloves, previous anti-rabies vaccination and hand washing before and after bat manipulation.

Vegetation Structure
To determine the influence of vegetation structure on the abundance of male and female bats, for each site we recorded all trees ≥10 cm diameter at breast height (dbh) within a 0.1 ha (20 × 50 m) plot (Gentry, 1982). Transects were located along streams in riparian continuous forest and riparian pasture habitats and randomly located in continuous forest and pasture habitats. We considered the following vegetation attributes: density of individuals, species richness, forest basal area, and tree height.

Data Analyses
First, we compared capture rates for each sex among habitat types by using a standardized capture rate (captures/1,000 mist net hour) that compensated for differences in number of nets, size of nets and length of time nets were open (Perry et al., 2010). We compared capture rates using analysis of variance on ranks (ANOVA). Data were checked to meet assumption of homoscedasticity and normality. We performed non-metric multidimensional scaling analyses (NMDS) based on the identity and abundance of tree species occurring in the sampling sites, to obtain a continuous synthetic variable summarizing dissimilarity patterns among vegetation species composition. The matrix used in the analysis was built using the Bray-Curtis index (Magurran, 2004). This iterative method of ordination has the advantage of properly handling non-linear species response of any shape (Oksanen, 2013) and has a good performance even when beta diversity is high (McCune and Grace, 2002). It is one of the preferred ordination methods for analyzing community data (McCune and Grace, 2002). The scores of axis 1 were used as an explanatory variable for evaluating differential sex response to tree species composition. Second, we fitted a general linear mixed model (GLMM) for each species separately with binomial error distribution and logit link function to determine the influence of vegetation traits on male and female abundance. In order to assess the relative effects of the explanatory variables on males in comparison to females, the model was run with the proportion of females to males per night (n = 70) as the response variable, with "site" as a random factor (Lintott et al., 2014). We considered the following explanatory variables: habitat type; season (dry and rainy), Vba, forest basal area per site; Vab, total number of trees; Vrich, trees species richness; Vh, average tree height; Vspcomp, scores of NMDS axis 1 (see Supplementary Tables 1, 2). For each model, we calculated Akaike's information criterion (AICc) corrected for small sample size following (Burnham and Anderson, 2004). This approach allowed us to select the most plausible models from a set of models. The set of models considered for every response variable, at each scale, included the null model (without explanatory power) and other models that considered each explanatory variable independently. We compared the model using i, which is the difference of AICc between a given model and the best (lowest AICc) model. We also calculate the AIC weights (wi) for each model. The wi represents the weight of the evidence that a certain model is the best model given the data and the set of candidate models. The 95% confidence set of the best models was defined by summing the wi, from the largest to the smallest, until the sum is = 0.95. Only models with an AICc lower than the null model were considered to define the 95% confidence set of plausible models.
All analyses were performed with R v. 1.0.136 (R Core Team, 2019).

Vegetation Attributes and Seasonality Response
Riparian habitats in mature forests (RM) and pastures (RP) showed higher vegetation structural complexity than nonriparian habitats. Despite of the fact that RM contained more trees than RP, these habitats showed similar average canopy height (Table 1), whereas the lowest canopy height was recorded in pastures. On the other hand, basal area was greater in riparian habitats that in non-riparian habitats, in both old-growth forests and pastures ( Table 1). The bioplot resulting of NMDS ordination is in Supplementary Figure 1.
Seasonality and habitat type were the best explanatory variables describing the presence of A. lituratus females; while for S. lilium, tree species composition and seasonality appeared as the most important variables explaining the incidence of females. In the case of G. soricina we found that the presence of females was positively related to vegetation structure including forest basal area and average tree height ( Table 2). The results for all species are in Supplementary Table 3.

DISCUSSION
Overall, the results showed that males and females of the selected bat species cope differently with habitat disturbance. Sex ratio was skewed toward females, which is an expected pattern for harem species (A. jamaicensis and A. lituratus) (de Mello and Fernandez, 2000). Furthermore, the results provide insights into FIGURE 2 | Capture rate (bats/mist net hour) of males and females of six bat species across different habitat types at the Lacandona rain forest during the dry and rainy seasons. The values correspond to 10, 50 × 2 m transects (0.1 ha) per site per habitat type. Only trees with diameters at breast height >10 cm were considered. the relative importance of habitat type for a specific sex. Except for mature forest, overall capture rates of females were greater than those of males, for the six studied species implying that during reproductive period, females of some bat species, may increase their activity in these habitats (e.g., foraging and/or drinking). Female bats can show changes in foraging activity probably due to lactation when energy and water requirements increase (Adams and Hayes, 2008;Barclay, 2012). Unlike males, females need to return at night to the maternity roosts to nurse their young which probably limits female foraging areas and restricts foraging females to the rewarding areas located in the proximity of their roosts (van Toor et al., 2011). Moreover, males may have greater survival than females (Keen and Hitchcock, 1980;Kurta and Matson, 1980), because males are not subject to the additional energetic pressures associated with pregnancy and lactation. For instance, in arid environments the drinking passes of lactating female bats were significantly higher than those of non-reproductive adult females, thus, survival of reproductive female bats seems to be conditioned to the availability or frequent and uninterrupted access to free-standing water sources (Adams and Hayes, 2008). We found that capture rates in active pastures were predominately toward females. The prevalence of females may be related to the high energy demanding of flying in cluttered habitat than flying in more open areas (Grodzinski et al., 2009), due to the elevated energetic costs associated with higher vegetation complexity might represent a particularly high burden for females during pregnancy and while nursing, males otherwise prefer sites with greater vegetation cover possibly related to roost defense (Henry and Kalko, 2007;Rocha et al., 2017).

Sex Ratio and Vegetation Structure
Contrary to our hypothesis, we found that females of A. lituratus did not exhibit relationship with habitat quality. The largest species of Stenodermatinae is A. lituratus, and as in other mammals, larger species could be more sensitive to human habitat disturbance (Lande, 1987). Nevertheless, the high number of females in this bat species, reflects selective foraging in a resource rich environment and higher roost availability in forested habitats. Even though Artibeus species may cross inhospitable matrix areas in fragmented landscapes, covering different vegetation types and flying distances ranging from 5 to 10 km (Galindo-González, 1998), females of A. lituratus might locally depend on temporal foliar roosts, and prefer larger trees within the dense and shaded mature forest that can provide energetic and thermal requirements to leave the young while foraging (Evelyn and Stiles, 2003;Bianconi et al., 2006;Arnone et al., 2016). For the same study area, A. lituratus has shown to select roosts with high humidity located in trees with the greatest basal areas (Ortiz-Ramírez et al., 2006).
Furthermore, some studies have argued that Glossophaginae are resilient to land use change (Willig et al., 2007). Our data indicate that the presence of G. soricina females is significantly and positively correlated to vegetation attributes such as tree height and forest basal area, supporting the idea that the species is an habitat specialist (Aguiar et al., 2014). The association between females of G. soricina with large trees could be due to the fact that large trees provide greater availability of roosts and foraging opportunities (Evelyn and Stiles, 2003;Ortiz-Ramírez et al., 2006). There is evidence showing food differentiation between specimens of G. soricina in the same areas where females preferred a plant food item different from males (Alvarez and Sánchez-Casas, 1999). During pregnancy and breeding seasons females might feed on the nearest available resource, whereas males fly larger distances in search of other feeding areas (Sosa et al., 1996). This foraging behavior might reduce the activity of G. soricina to habitats with high resource availability limiting its activity to small home ranges increasing its susceptibility to local extinction (Arita and Santos-del-Prado, 1999). This increases the importance of vegetation traits associated to the reproductive energetic demand (pregnancy, lactation and roost defense) (Charles-Dominique, 1991;Klingbeil and Willig, 2010).
Tree species composition was the strongest predictor variable for S. lilium, this frugivorous bat is known for its preference for understory shrubs and pioneer tree species, females of S. lilium are able to forage among forest elements as a result of nonrandom distribution of resources across the landscape (Loayza and Loiselle, 2008). Moreover, frugivorous bats like S. lilium which can visit different vegetation types can be considered as indicator taxa of habitat change in riparian vegetation, rather than highly specialized taxa in which population decline rapidly under environmental changes (de la Peña-Cuéllar et al., 2015).

Sex Ratio and Seasonality
The importance of seasonality for A. lituratus and S. lilium females might reside on differences in abundance and diversity of food resources between the wet and dry seasons. Seasonal fluctuations in rainfall influence phenology of fruiting plants and affects productivity in tropical forests (Ramos Pereira et al., 2010). In tropical regions, usually the rainy season corresponds to greatest fruit abundance than the dry season (Smythe, 1986). Even, many tropical bats timing their reproductive phenology to match periods of peak food availability, female bats may be constrained by the energetic requirements associates with reproduction, which might force them to alter their foraging behavior (Lintott et al., 2014). Resource availability due to seasonality can result in shortened flights during the exploration for food and shelter, whereas during the dry season when food resources are often scarce, females respond to localscale vegetation structure increasing foraging movements into resource-rich pioneer fruiting plant species areas like secondary forests, while males tend to select areas close to old-growth forests (Rocha et al., 2017).

Implications for Conservation
Our results show that responses of bat to human disturbance are sex-specific. Taking into account sex ratios of bat populations may help to a better understanding of the pervasive consequences of habitat loss and fragmentation. Sex-specific studies are important for bat conservation practices in order to promote habitat conditions favorable for both, females and males (Perry et al., 2007). Vegetation attributes like three height and basal area reflect the age and vertical complexity of forest, and more varied niche opportunities for bats, enhancing greater taxonomic and phylogenetic biodiversity (Martins et al., 2017). Our results suggest that the structural complexity of the vegetation and large trees influence the presence of females of G. soricina. In this sense, management efforts should promote riparian vegetation cover with large forest basal areas, important for the conservation of the entire bat community. Even though our analysis was restricted to six common bat species, we assumed that the maintenance of habitats that favor habitat generalist species should also benefit bat species that are habitat specialists (Istvanko et al., 2016;Rocha et al., 2018). Therefore, we encourage the inclusion of species of sensitive trophic guilds (gleaning insectivores and carnivores) and particularly roosting habits emphasizing adequate protection of females in conservation plans. Conservation actions toward female protection are particularly important due to their high level of parental investment associated with rearing pups (Istvanko et al., 2016). Management decisions that do not guarantee the protection of the habitat frequently used by female bats would likely have detrimental long term consequences on their reproduction, jeopardizing the dynamics and long-term persistence (van Toor et al., 2011;Frank et al., 2016).
In human dominated landscapes, the presence of isolated trees in pasture promote bat flights across pasture and increase bat detectability (Galindo-González et al., 2000), based on foraging behavior frugivorous bats travel across pastures and visit isolated trees while foraging, and use canopies for roosts or to decrease predation risks (Galindo-González, 1998), also isolated trees in pastures act as stepping stones for traveling across fragmented landscapes to different forest remnants (Guevara et al., 1989). However, pastures cannot sustain the same species richness of bats as old-growth forest and riparian vegetation (de la Peña-Cuéllar et al., 2015). Moreover, land use change in tropical landscapes seems to have considerable effects on bat population dynamics, for instance evidence suggests that forest-adapted insectivorous species are particularly sensitive to habitat conversion (Medellín et al., 2000;Williams-Guillén and Perfecto, 2010). Furthermore, frugivorous bats respond to matrix quality in different manners, whereas studies have found that frugivorous abundance was positively associated with the proportion of high quality habitats (Pinto and Keitt, 2008;Avila-Cabadilla et al., 2012;de la Peña-Cuéllar et al., 2015), some other studies have shown that frugivorous bat richness and abundance are higher in moderately fragmented landscapes than in oldgrowth forest (Willig et al., 2007;Klingbeil and Willig, 2009). Overall agricultural intensification may cause detrimental effects on bats and thus presumably on the ecosystem services they provide (Williams-Guillén et al., 2015). Additional research is needed to directly examine the effects of pregnancy and lactation on habitat selection by bats. We encourage radio-tracking studies that can show specific habitat use of males and females (roost and foraging areas) and if there is temporal segregation between sexes; this could provide information about how different habitats have an impact in the demography of bat populations.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
Ethical review and approval was not required for the animal study because the handling was very brief and we did not perform animal sacrifices.

AUTHOR CONTRIBUTIONS
EP-C and JB-M conceived and designed the experiments, contributed reagents, materials, and analysis tools, wrote the manuscript, and revised the manuscript. EP-C performed the experiments and analyzed the data. Both authors contributed to the article and approved the submitted version.
Conservation International (BCI Student Research Scholarship Program). We are grateful to R. Lombera-Estrada, G. and I. Lombera, G. Rodríguez-Barrera, and J. L. Peña-Mondragón for their assistance in the field. We are also grateful for logistical support provided by J. M. Lobato-García.