Margay (Leopardus wiedii) in the southernmost Atlantic Forest: Density and activity patterns under different levels of anthropogenic disturbance

The margay (Leopardus wiedii) is a small Neotropical arboreal wild cat. This species is thought to be forest-dependent, although few studies so far have directly evaluated the relationships between spatiotemporal aspects of its ecology and landscape characteristics. The aim of this study was to estimate margay population density and activity patterns in six areas with different habitat types and levels of anthropogenic disturbance in the southernmost Atlantic Forest of Brazil. Our working hypothesis was that density and activity patterns differed between areas in response to differences in forest cover and anthropogenic disturbance. Margay records were obtained using camera trapping, during spring and summer from 2017 to 2019. In all areas, the sampling scheme consisted of 20 un-baited stations, set 1km apart, each containing two paired cameras. We assessed the potential effects of environmental variables, including anthropogenic factors, on margay density, rate of detection and space use by comparing nine spatial capture-recapture (SCR) models. Activity patterns of the margay, its potential prey, and competitors were described and compared using the date and time of the records. We obtained 66 records of margay. Two of the six sampled areas were excluded from subsequent analyses due to the small number of records. The density estimated by the top-ranked model varied from 9.6±6.4 individuals/100km2 in the area with the highest human disturbance to 37.4±15.1 individuals/100km2 in a less disturbed area. Margay densities responded positively to vegetation cover, supporting the hypothesis of forest dependence by the species. Both the margay and their potential prey (small rodents and marsupials) were found to be mostly nocturnal. Margay activity also overlapped with that of the ocelot, Leopardus pardalis, and with mammals associated with human presence (wild boar, cattle, domestic dogs and cats). This is the first multi-area study on patterns of density and activity of the margay in the Brazilian Atlantic Forest. We concluded that the margay is mostly nocturnal, and while its densities are positively influenced by forest cover and negatively influenced by human disturbance, the activity pattern of the species does not seem to change across landscapes with distinct levels of human modification. Margay populations seem to be able to persist under moderate levels of habitat modification, highlighting the importance of preserving even small native forest remnants in the highly fragmented Atlantic Forest.

Introduction spatial structure of the ecological processes [1]. Spatial capture-recapture (SCR) models should be able to overcome this limitation by incorporating spatial information from the individual detections [1,34,35]. In contrast to density estimates, a much greater number of studies has described activity patterns of margay across its range, suggesting that, overall, the species is nocturnal [6,8,17,30,31,[36][37][38][39], though in some areas it may show a cathemeral pattern [40].
The first cycle of the Brazilian National Action Plan for the Conservation of Small Cats (CENAP/ICMBio) defined seven specific objectives, one of which is to assess how different natural and anthropogenic processes influence the populations of small felids [41]. To address this objective, here we aimed to estimate population density of margays using spatial capturerecapture models, and to describe its activity patterns across a range of areas with different levels of human disturbance in the southernmost limit of the Atlantic Forest. We hypothesize: i) density differs across study areas; ii) activity pattern is mostly crepuscular/nocturnal but may change across areas. We expect margay densities to respond positively to forest cover and negatively to human disturbance and ocelot (Leopardus pardalis) presence, as margays tend to respond to the presence of other felids. In fact, Oliveira et al. (2010) reported an 'ocelot effect', suggesting that the presence of the ocelot, negatively impacts densities of margay and other small cats [32]. Additionally, we expect L. wiedii to present more nocturnal activity in areas with higher levels of anthropogenic modification (with more occurrence of domestic and exotic species) [15,42]. Ultimately, we aim to generate baseline information for the definition of management actions towards margay conservation at the southern extreme of its distribution, where it may be particularly sensitive to population fluctuations [28,29] and is categorized as "Vulnerable" [43,44].

Study area
We sampled six areas in Rio Grande do Sul state, southernmost Brazil (Fig 1). This region comprises the southern portion of the margay's distribution [22] and the southernmost limit of the Atlantic Forest biome. The Atlantic Forest extends beyond tropical climates, with semideciduous and ombrophylus mixed forests (Araucaria forests) at higher elevations gradually replacing the dense ombrophylus forests that are typical of the lower altitudes [45][46][47]. This region also includes the ecotone between the Atlantic Forest and the Brazilian Pampas, which increases the structural and compositional landscape complexity, with changes in climate, elevation, vegetation and, consequently, beta diversity [48][49][50]. Only 12% of the original area of the Atlantic Forest persists, with the remaining area having been almost completely replaced by croplands and other human-modified landscapes [27]. Indeed, due to habitat loss and fragmentation [22,43] the margay is categorized as "Vulnerable" in Rio Grande do Sul state [43,44].
We collected data from 120 sites in six areas spanning a range of habitat types and human land-use intensity; the six areas were: 1. Serra Geral National Park (SGNP) (29˚08'2"S, 49˚59'40"W), a federal conservation unit of strict protection, constituted by a mosaic of natural landscapes of different vegetation formations, such as Araucaria forests [51]. The park portion sampled in this study was in the process of expropriation, with cattle still roaming in some patches; 2. Banhado dos Pachecos Wildlife Refuge (BPWR) (30˚06'22"S, 50˚52'11"W), a protected area located close to large urban centers (ca. 28km from the state capital, Porto Alegre); this region is composed of large areas of plains and grasslands [22,52] and forest fragments among roads and residences; two camera-trap stations ended up set outside the protected area and close to human settlements, due the standard distance between the stations in the sampling scheme. This study area is located in the ecotone between the Atlantic Forest and the Pampa biomes [53].
3. Pró-Mata Centre for Research and Conservation of Nature (PROMATA) (29˚28'54"S, 501 0'35"W) is a private natural heritage reserve, the largest private protected area in the state https://doi.org/10.1371/journal.pone.0232013.g001 [54]; previously occupied by small farms and agropastoral activities, this area has been under natural regeneration for over 20 years [55]; 4. Turvo State Park (TUSP) (27˚08'44"S, 53˚53'10"W), a strict protected area created in 1954, this park is one of the largest conservation units in the state [56] and the most pristine area in this study; 5. A rural area close to the city of Teutônia (TEUT) in the central region of the state (292 6'36"S, 51˚47'57"W); this is a non-protected area comprising a matrix of small private properties with different agricultural activities (corn and soy crops) and livestock (chicken, pigs, dairy cattle); in high hills of the region, there are fragments of native and exotic forests that are functionally connected to each other; 6. Passo Fundo National Forest (PFNF), in the northern region of the state (28˚18'47"S, 521 0'55"W), is a protected area with sustainable land use, resulting from the restoration of an agricultural area, predominantly planted with native trees [57,58].
These study areas were categorized with respect to different anthropogenic disturbance levels in a parallel study performed simultaneously at the same areas [59]. The authors generated the levels base on different variables (e.g. distance to forest edge, distance to nearest cell phone tower, stable light at night values, predicted abundances of domestic carnivores, etc.). Here we follow the gradient from the most disturbed area to the most preserved one: BPWR, TEUT, PFNF, SGNP, PROMATA, TUSP.

Field sampling
We collected margay records between 2017 and 2019, during summer and spring, averaging 60-sampling days (56 to 62 days) per area, to assume closed margay populations [60]. In each area we installed 20 camera-trap stations placed ca.1 km apart ( Fig 1C). This layout was defined based on the diameter of female home ranges for the species (<1km 2 ) estimated in previous studies [22,32]. Each sampling station was composed of two passive infrared digital camera-traps, one at each side of a wildlife trail or road, totaling 40 camera-traps per sampling area. Pairing the cameras allowed recording the two flanks of the detected animals for individual identification through the individual-specific pelage pattern [61]. We set up the unbaited camera traps at 30-40 cm above the ground. We programmed the cameras in video mode (10 s, with a 5 s motion triggered delay), to remain active 24 h per day, recording date and time. We used several camera models and brands: Bushnell Trophy Cam™ and Nature View 1 (Bushnell Outdoor Products, Overland Park, Kansas), Digital Game Camera Moultrie (Moultrie Products, LLC, Birmingham, Alabama), Browning 1 Trail Cameras (Prometheus Group, LLC, Birmingham, Alabama) and Scout Guard Infrared Digital Scouting Camera (Boly Media Communications Co. Ltd., Shenzhen, China).

Density and detection covariates
Based on the literature and previous knowledge on margay's biology, we defined a set of environmental covariates, including anthropogenic-related ones, as predictor variables for the species' spatial scale, density and rate of detections. We used the individuals' sex as a covariate to investigate the spatial scale, eight covariates for fitting models for rate of detection, and five covariates for density modeling; we evaluated if the three parameters differed between areas (Table 1). We excluded highly collinear predictors using the variance inflation factor (VIF), excluding variables with VIF > 6 [62](S1 Table).

Density modeling
We identified individual margays from the videos based on their unique spot patterns (Fig 2) and determined the sex of each individual, whenever possible, through the visualization of presence/absence of male gonads. We used spatial capture-recapture (SCR) models [1,67] following the workflow of package oSCR 0.42.0 [66,68] in R 3.6.0 [69] with multi-session sexstructure models to investigate the population density (D), the rate of detection (p) and the space use (σ) for margay in the study areas.
Each sampling area was considered as harboring a distinct margay population, thus each one represents an independent "session" for the analyses. We removed from the analyses the areas of SGNP and TUSP due the low number of records. To create the oSCR data object, a single EDF file (encounter data file) with information on the individual encounter history data diswater Euclidean distance raster created in ArcGis based on shapefile of water bodies of the Regional Executive Organization for Environmental Protection [64] Values range from 0 to 10km Density will be higher closer to natural water bodies Human population density (1km 2 ) popdens Estimate of human population density (ranging from 0 to 10000 people per square kilometer) Socioeconomic Data and Applications Center; Gridded Population of the World, Version 4 (GPWv4) for2015 [65] Density will decrease with increasing human densities Distance to roads (1km 2 ) disroads Euclidean distance raster created in ArcGis based on shapefile of roads of the Regional Executive Organization for Environmental Protection [64]; missing roads included manually through own observations Values range from 0 to 10km Density will be higher as one moves further from roads (paved or unpaved roads) Higher PIR detection range values, i.e., larger detection zone will increase the rate of detection of margay

Spatial-scale
Sex sex Sex of the individual (female, male or undetermined) [66] Spatial use will differ between the sexes Sessions b session Data from groups which can be sampled in different spatial or temporal independent studies, these data groups are called "sessions" [66]. We used the sampling areas as our sessions The density will differ between the areas a we used the setting best camera-trap model to represent the camera-trap station (site). b we used this covariate for all parameters tested. https://doi.org/10.1371/journal.pone.0232013.t001

PLOS ONE
for all areas (session, individual ID, trap, occasion, sex) was linked with the TDF files (trap deployment file) of each session. The TDF files contained name and coordinates of each site (camera trap station), the trap-by-occasion binary operation data (1 = operational, 0 = not operational) and trap-specific covariates that either varied with occasion (time varying) or were the fixed site covariates ( Table 1). The number of occasions was the same among sessions (60 occasions), except for PFNF (56 occasions). An important component of the SCR analysis is the state space (S), that in oSCR is required for each session [66]. The state space data object (ssDF-state space data frame) was created by defining a buffer distance around the camera traps and a specific resolution defining the state-space centroids, based on the session-specific trap coordinates. We used a buffer distance of ca. four times the space use parameter (σ) estimated (2000m) and a resolution value of half σ [66], (250 x 250 m) (Fig 3). We clipped out non-habitat points (e.g. water bodies), categorized using the MapBiomas [70] raster, from the buffers to avoid bias in the density estimates [1,71]. We used a three-step approach for modelling margay densities: 1) first we analysed the space used by the individuals, fitting the σ parameter models with covariates, such as the sex (female or male) and the session. In these models we set the D and p parameters as constants: D~1; p~1, σ [covariate]; 2) in the second step we investigated potential effects of covariates in the rate of detection (p) of margay, including sex, session, and constant or time-varying traplevel covariates [66]. In this step, the σ parameter was set according to the best model resulting from step 1, while the D parameter continued to be set as a constant (D~1) (D~1, p [covariate], σ[first step]); 3) in this final step we allowed the D parameter to vary as a function of a single covariate or of an additive combination of two parameters (D~[covariate], p [second step], σ[first step]). The models built in all steps represent our biological hypotheses regarding the effects of covariates on margay density (D), rate of detection (p) and spatial scale (σ). We ranked the models using the Akaike Information Criterion (AIC) [72], considering equally fitted models those with ΔAIC � 2 [72]. The covariates presented in the top model or models were considered as possible determinants of species density, rate of detection and spatial use.

Activity patterns
Activity patterns of the margay were evaluated using date and time of the camera-trap records. To maintain temporal independence between the records and avoid autocorrelation, we only considered those with at least a 1-hour gap [10,73]. We applied the same approach to all other species expected to affect the activity pattern of the margay, including possible prey items (small mammals and small birds), a possible wild competitor (the ocelot), and exotic/domestic species such as domestic dog, cat, cattle and boar.
We tested the uniformity in the activity pattern over the 24 hours of the day for margay, small mammals, and small birds using Rao's Spacing Test and we applied Watson's Two-Sample Test of Homogeneity to compare the pairwise distribution of activity of using the package circular 0.4-93 [74] in R. If the sample group had homogeneity, we used records of the group from all areas together and applied Rao's Spacing test to measure the uniformity of activity over 24 hours for the entire region of southernmost Atlantic Forest.
We compared the activity of margay with that of the other species by estimating the coefficient of overlap (Δ) per area. We used the package overlap 0.3.2 [75] I n R [76]. The coefficient of overlap (Δ) ranges from 0 (no overlap) to 1 (complete overlap), and is described graphically through a kernel density curve [75,77]. We used the threshold proposed by [7] to classify the degree of overlap between species. We adopted the estimator Δ 1 , for small samples (n>75) [75,77] and performed smoothed bootstraps, with 1000 resamples, to obtain confidence intervals for the Δ 1 estimator [77].

Ethics statement
Research permits for all protected areas were obtained from the Brazilian Ministry of the Environment (permits SISBIO-49050 and SISBIO-64647-1) and from the Environment Secretariat of Rio Grande do Sul (permit SEMARS-588), as well as directly from the landowners in the case of private lands. During this study, no animals were captured or handled, and thus no additional permits or protocols were required by the Brazilian law.

Field sampling
Our sampling effort resulted in 7220 camera-trap-nights and 66 independent margay records across the six sampling areas; six of these were in BPWR, 27 in PROMATA, 17 in TEUT, 12 in PFFN, two in SGNP and two in TUSP. The two latter sites were thus removed from subsequent statistical analyses due their small sample size.
Also, we obtained 3089 independent records of small birds (Columbidae and Passeriformes), 289 of small mammals (small rodents and small marsupials) and seven of ocelot (Leopardus pardalis). In addition we obtained the following exotic and domestic species independent records: 51 of domestic cat (Felis catus), 75 of domestic dog (Canis familiaris), 29 of wild boar (Sus scrofa), 51 of cattle (Bos taurus, Equus caballus); some of these were only detected in some of the sampled areas (S2 Table).

Population density estimates
We were able to individually identify 23 margays, in the four areas considered for analysis: two in BPWR, 10 in PROMATA, six in TEUT, and five in PFNF. We identified the sex of 19 individuals and in four individuals we were unable to determine the sex. We discarded from the analyses seven records for which individual identifications were not possible, due the lowquality videos, totalling 55 available records for the density analysis. The individual encounter frequency ranged from six captures of one individual in the same trap, to only a single capture event (for eight of the individuals). The four sampling areas considered for the analyses and the respective state space are shown in Fig 3. For the first step, regarding models created to test covariates on the spatial scale, only one model returned ΔAIC �2, representing 79% of the weight of the models ( Table 2). In the spatial scale of detection, sex influenced significantly the σ movement parameter (β = 0.66, 0.27-1.06 CI), which was larger for males (1.19 km) than for females (0.59 km) (Fig 4). In the second step, regarding the 13 alternative models testing covariate effect on rate detection, the two-top ranked (ΔAIC <2) represented 54% of the weight of all models ( Table 3). The first model estimated a small but positive influence of small bird detection on the rate of detection of the margay (β = 0.23 ± 0.09 SE [0.05-0.41 CI], P = 0.02, weight = 39%), while the second best estimated similar influence of this covariate, and an additional effect of small mammal detection, though not significant (β = 0.12± 0.42 SE [-0.67-0.93 CI], P = 0.794) (Fig 4). Because biologically it makes sense that detections of margay respond to the presence of potential prey, we choose the second ranked model as best fitted, and kept it for the spatial variation density models.
From the nine candidate models in the third modelling step, the top model included only vegetation cover as covariate for the density of the margay, with 33% of the weight of all models (Table 4), with a marginally positive influence on the density (β = 5.81±3.01, SE [0.14-11.4 CI], P = 0.05), (Figs 4 and 5). We disregarded the second and third models, because none of the covariates (distance of water and population density) were statistically significant. Density estimates varied between the studied areas ( Fig 6); moreover, the model including the session covariate (comparison between areas) indicated significant differences between the densities in BPWR and PROMATA (P = 0.04); however, this model was not top ranked ( Table 4). The density estimates from the top-ranked oSCR model varied between the areas: 9.6±6.4 individuals/100km 2 in BPWR; 37.4±15.1 individuals/100km 2 in PROMATA; 29.6±11.4 individuals/ 100km 2 in TEUT; 28.4±12.5 individuals/100km 2 in PFNF (Figs 6 and 7).

Activity patterns
To evaluate the activity pattern of the margay we considered 62 independent records, excluding the records from SGNP and TUSP (n = 2 for each). Because there were no significant differences in margay and small mammal activity patterns between the sampled areas (Table 5), we merged the records for each taxon to tested for uniformity in the activity pattern over the 24 hours of the day. Margay showed statistically significantly non-uniform activity pattern (Rao's Spacing Test, r = 261.29; P< 0.001), with average between 01:00 and 02:00am (Fig 8A).
Regarding the activity patterns of the potential prey, small mammals also presented non-uniform activity pattern throughout the daily cycle (Rao's Spacing Test r = 341.31, P< 0.001) with average at 00:00 ( Fig 8B). Temporal coefficient overlap between margay and small mammals was 0.73 (0.62-0.87 CI) (Fig 8C). The activity pattern of small birds was dissimilar between the sampled areas (Table 5), so activity overlap between margay and small birds was tested separately for each area; temporal coefficient overlap between margay and small birds ranged from 0.15 to 0.39 (Fig 9). Overlap coefficients between margay and dogs ranged from 0.24 to 0.56 (Fig 9). The estimation of the temporal overlap between margay and domestic cats was only possible for BPWR and TEUT, with 0.48 and 0.73 overlap coefficients, respectively (Fig 9). The lowest temporal overlap estimated was with humans, cattle (horse and cow) and small birds (Fig 9). Temporal overlap between margay and ocelot, and margay and wild boar were estimated only for PRO-MATA, resulting in 0.45 and 0.44 overlap coefficients, respectively (Fig 9).

Discussion
Few studies estimated the density of margay [30,36,78] and compared temporal activities of any Neotropical small cats with human related species or human related impacts [17,79]. Our study represents the first large-scale evaluation of density and activity patterns of the margay for the southernmost limit of the Atlantic Forest and may be used as model for other studies seeking the management and conservation of small felids.

Margay: Undoubtedly a forest cat
In the southern limit of the Brazilian Atlantic Forest, the margay seems to occur at higher densities in areas with higher vegetation cover, as hypothesized, probably resulting from the known arboreal habit of the species [20,23]. This arboreal behavior also seems to contribute to the potential relevance of small birds in the diet of margays [40,80]; indeed, the rate of

PLOS ONE
detection of small birds positively impacted rate of detection of the margay. Though we did not evaluate diet composition, this pattern was expected as small birds and small mammals constitute important food resources for the species [40,80]. The effect of the detection of small mammals on margay's detection was not statistically significant as occurred with the detection of small birds. Indeed, we obtained a low number of small mammal records, which may have resulted from the camera-trap position (angle), as distinct camera-trap position is necessary for the maximum detection of small mammals [81,82]. Margay males seem to walk longer distances (1.19 km) than females (0.59 km). A similar pattern was already reported in Mexican margay populations, who estimated average home ranges of 4.1 km 2 for males and of 0.72 km 2 for females [52]. Similar differences have been reported for other felids in the Atlantic Forest, including the jaguar and the ocelot [83,84]. Indeed, larger home ranges of solitary males of the Carnivora seem to respond to resource availability and as response to reproductive opportunities [85]. We found a lower use of space by females than reported by IUCN for the Mexican margay populations [22,52]. Reported influence of roads and human occupation on jaguar female use of space [86] may perhaps explain our results with the margay.
The highest densities of margay occurred in PROMATA, the largest privately protected area of Rio Grande do Sul, and where ocelots seem to be rare (unpublished data; (S2 Table). This area presents low level of human disturbance [59] and is in large part occupied by primary forests, and forests and natural fields under natural regeneration from agricultural use for the last two decades [87]. The existence of high margay densities indicates that, unless a landscape is profoundly modified, active management practices towards forest regeneration should allow the persistence of a significant number of individuals in smaller fragments, most likely because arboreal shelters and prey exist in sufficient numbers to reduce intraspecific competition [88,89]. Consistently, the lowest density estimate was obtained for BPWR. In this area, forest coverage is small, as the area is dominated by open physiognomies such as permanent wetlands, marshlands, restingas and dry grasslands [90] and, indeed, suitability models suggest the species to be negatively correlated to flooded grasslands and savannas [24]. From our sampled areas, the BPWR is that in the southernmost latitude, and situated in the ecotone between the Atlantic Forest and the Brazilian Pampa. In Pampean landscapes, the marginal occurrence of the margay should be tied with the natural forested fragments of the biome, occurring mostly bordering water courses [24]. Besides that, the BPWP, although a strictly protected conservation unit, is located only 28 km away from the largest urban center of Rio Grande do Sul (Porto Alegre, with ca.1,409 million habitants) [59], inevitably suffering from the expected

PLOS ONE
associated impacts: illegal hunting, water drainage, and intensive agricultural schemes, particularly rice production [91].
TEUT and PFNF presented intermediate densities suggesting similar responses of margay populations in these two areas. PFNF is a relevant representative of the Araucaria Forest physiognomy in the southern Atlantic Forest, showing a gradient of regeneration after somewhat intensive human use in the last decades [57,58]. Although a conservation area of sustainable use, it is located in a matrix of monocultures and livestock production, representing an isolated fragment with some, while low, connectivity with other forest fragments; agriculture and illegal hunting represent the major negative impacts on the native fauna. These characteristics may explain why PFNF, being a protected area with somewhat similar features to PROMATA, presents densities of margay more similar to TEUT an altogether non-protected area. TEUT is basically an area of small properties with different agricultural uses, harboring forested fragments in those areas considered unsuitable for crop production. The estimated margay density here is higher than that estimated for other neotropical cats, including the southern tiger cat (Leopardus guttulus) and the jaguarundi (Herpailurus yagouaroundi) (8 individuals/100km 2    Table). https://doi.org/10.1371/journal.pone.0232013.g009

PLOS ONE
and 4 individuals/100km 2 , respectively [92], though these studies used different methods for estimating population densities). We were unable to estimate densities for TUSP and SGNP, due to the small number of records obtained for these areas. In TUSP, low margay detection contrasts with high ocelot densities (14 (14-24 CI)-66  individuals/100km 2 [93] and 15.5 (6.4-36 CI) individuals/ 100km 2 [94]). The suspected 'ocelot effect' may be affecting margay; still, because this is quite a pristine area, with high native tree cover, we expected to collect enough data to at least estimate the density in the area, which was not the case. In SGNP, we suspect that the somewhat intensive presence of exotic/domestic fauna within the park (cattle, boars and buffalos), as several areas are still ongoing private expropriation, could explain the low number of margay records [24]. Additionally, the area is a mosaic of forest and open fields and, within our sampled sites, is that with the largest proportion of grassland habitats, which may explain the scarce number of records of the margay, positively influenced by vegetation cover and seemingly forestdependent.
To our knowledge there is no published information on density patterns of the margay in the Atlantic Forest. Our estimates are, however, lower than those reported for southwest and southeastern Mexico, which varied from 12 individuals/100 km 2 [30] to 81 individuals/100 km 2 [31], respectively. Overall, margay population densities seem to respond positively to vegetation cover. Higher estimated values of density were found in more preserved areas and the lower values in human-altered landscapes [59]. The intermediate values of density were found in areas with moderate human land use and reduced natural vegetation cover [59], suggesting that conservation strategies focusing on less pristine or small forest fragments may have positive effects on the density of margays across its distribution range.

Margays prefer the night: Avoiding antagonists or simply following prey?
Margay was nocturnal across our studied areas, consistent with what has been found in other studies [17,30,31,39,95]. In fact, our results did not support our prediction of changes in the activity pattern of the margay in response to human disturbance: margay was strictly nocturnal across a range of human-altered landscapes. This pattern probably co-evolved as a response to the activity pattern of their preferential prey-small mammals [80]-but may have also been intensified in the last centuries, by human disturbance [15], including hunting pressure and the presence of antagonists-potential predators such as domestic dogs or simply antagonists in the sense of profound habitat modifiers such as wild boars [79,96,97].
Small mammals are recognized as the most important dietary item of margays [40,98,99]. Similar nocturnal activity patterns probably co-evolved over millions of years increasing on the one hand the chances of prey capture by the margay, but on the other, increasing possibilities of the potential prey to escape. However, small birds, also an important item in diet of margay [40,80], are mostly diurnal, suggesting that if eaten in our sampled areas, they are probably preyed upon while resting [100].
Ocelot activity did not overlap significantly with that of the margay, and that is a pattern already described for southeast remnants of the Brazilian Atlantic Forest [40]. Temporal segregation facilitates species coexistence, by reducing competition for space and prey [6,40,101]. However, it is important to notice that the number of ocelot records in our study was small, and more data will be necessary to support this conclusion. Apart from temporal segregation, segregation by vertical stratification between the ocelot and the margay still needs to be tested, as this has been suggested as the factor behind the higher activity overlap between the two species in the Argentine Atlantic Forest [17]. Moreover, the margay nocturnal activity reported in this study can be associated with a temporal segregation between this species and another felid, as diurnal jaguarundi (Herpailurus yagouaroundi) [6,17]. These two species seem to have the highest similarity in morphological traits associated with trophic ecology among Neotropical cats, and the contrasting activity pattern may be relevant to allow their coexistence [6,102].
Temporal overlap between the margay and domestic dogs was low in general, with the highest overlap occurring in TEUT. Dogs presented diurnal activity, consistently with the pattern found in other Atlantic Forest areas [103,104]. Theoretically, such low overlap suggests that domestic dogs probably do not prey upon margays. Still, they may have negative indirect effects, representing an important cause of decreased prey populations [105] or as reservoirs for pathogenic agents [106]. Domestic cats, on the other hand, revealed to be mostly nocturnal, and this was particularly evident in TEUT where a considerable higher number of domestic cat records was obtained, probably related to the higher number of human habitations [107]. Domestic cats, therefore, may represent a potentially important competitor for food regarding the margay and other felids (e.g. southern tiger cat, jaguarondi), as they show similar opportunistic predatory behavior [97,108]. Disease transmission is a very likely indirect negative impact resulting from the occurrence of both species in the wild, representing an additional threat to native felids [18,[107][108][109][110][111].
We also found cattle and wild boar in some of the studied areas, although both were mainly diurnal. Potential impacts of both species on margay populations are thus probably related to changes in the landscape [112][113][114][115][116], especially by the wild boar, known for intense nest predation and destruction of native wildlife habitats [96,[117][118][119].

Conclusion
To our knowledge, this is the first study specifically attempting to estimate multi-area density and activity patterns of the margay, and also to compare the activity pattern of the margay with those of human related species. Our results supported our hypothesis that densities differ across the sampled areas, reflecting differences in the composition of the landscape and in the levels of human disturbance. On the other hand, we found no significant changes in the activity pattern of the species between the sampled areas. Indeed, the margay seems to be mostly a forest nocturnal cat, whose densities are positively influenced by forest cover and negatively influenced by intense human-related disturbances, while not changing its activity pattern across landscapes with distinct intensities of human use. Undoubtedly, large pristine forest areas, with high prey-small birds and mammals-availability, are critical for the persistence of dense populations of the species. However, under moderate levels of habitat modification and human disturbance, the margay is still able to persist-at least at intermediate densitiessuggesting that the conservation of even small native forest remnants, especially those showing some degree of connectivity among them, is key for margay population management and conservation at the southern limit of the Brazilian Atlantic Forest.
Supporting information S1