Breeding biology of a relictual Maghreb Magpie (Pica mauritanica) population in Tunisia

The Maghreb Magpie (Pica mauritanica) is an endemic North African species. Available knowledge on this species is limited to historic descriptive data with no ecological information provided. Populations continue to dramatically decline in Tunisia, where only one relic population survives. Investigating the breeding biology of this species is essential for conservation purposes. The purpose of this study was to increase our understanding of the Tunisian relic population and provide detailed data on breeding biology over two breeding seasons (2017 and 2018). This study occurred on a private farm of 650 ha, located 10 km from Dhorbania village at Kairouan Governorate, in central Tunisia. Active nests were monitored weekly during egg laying period and twice a week during hatching period. The Ivlev’s electivity index was used to assess whether the frequency of use of nesting trees and bushes matched their availability in the study area. We recorded nest measurements and positions, and compared them using Wilcoxon signed-rank test. Variations of breeding parameters as number of eggs laid, hatchlings, and fledglings over years were performed using Mann–Whitney U-test and χ2 tests. We used a Generalized Linear Mixed Model (GLMM) to investigate how egg volume varied with clutch size and laying date. We investigated clutch size, egg size, hatching and fledging success, and evaluated how these parameters varied according to laying date and nest characteristics. Clutch size averaged 5.00 ± 0.19 but was significantly greater in 2017. Hatching success was 2.78 ± 0.34 eggs hatched per nest and fledging success reached 1.69 ± 0.30 young/nest. Causes of nest failure included the depredation of nestlings by shrikes, cobras and rats (e.g. Lanius meridionalis, Naja haje and Rattus rattus), death of parents by the Black-shouldered Kite (Elanus caeruleus) and nest parasitism by the Great Spotted Cuckoo (Clamator glandarius). Clutch size, brood size and fledgling success were unaffected by laying date, nest volume and nest elevation. Egg volume decreased with laying date but was unaffected by clutch. Our study provides the first and only detailed data on reproductive parameters of the Maghreb Magpie in its entire geographic range (North Africa). Information gleaned from this study provides valuable information for monitoring and long-term conservation plans of the endangered Tunisian Magpie population. Additionally, our data provide an avenue of large-scale comparative studies of the reproductive ecology of the magpie complex.


Background
Investigating the breeding biology of bird species is an essential step for understanding their population dynamics and for focusing on conservation efforts. This is particularly important for poorly known species, especially those showing signs of declining populations.
The Eurasian Magpie (Pica pica) is a common bird species in the Palearctic area. Until recently, the North African population of this species was considered as a subspecies (Pica pica mauritanica). However, Ebels (2003) highlighted morphological and acoustic differences between the Maghreb population and those of Europe. Later, taxonomic and phylogenetic studies (Kryukov et al. 2017;Song et al. 2018) illustrated that the North African clade is the most divergent compared to other lines. Consequently the mauritanica subspecies is now considered as a separate species, the Maghreb Magpie (Pica mauritanica), endemic to North Africa (del Hoyo et al. 2018).
At the beginning of the twentieth century, the Maghreb Magpie was locally widespread in North Africa. In the past century, 12 isolated local populations were recorded in Tunisia (Lavauden 1924;Heim de Balzac and Mayaud 1962;Etchécopar and Hue 1964;Isenmann et al. 2005), 16 in Algeria (Etchécopar and Hue 1964;Isenmann and Moali 2000) and approximately 30 in Morocco (Etchécopar and Hue 1964;Thévenot et al. 2003;Bergier et al. 2017). Nowadays, the distribution of this species in North Africa is more fragmented and limited. For instance, only one relic population still survives in Tunisia, at the Sbikha region (Kairouan Governorate) Nefla et al. 2020). However, despite this critical situation and the urgency of adequate conservation measures, the ecology and dynamics of this relic population remain poorly understood due to the lack of detailed studies. Indeed, available knowledge on this species is limited to old descriptive data on its morphology and distribution (Lavauden 1924;Heim de Balzac and Mayaud 1962;Etchécopar and Hue 1964;Isenmann et al. 2005).
The breeding biology of the Eurasian Magpie has been well studied throughout its whole geographic range, particularly in the Palearctic region (Sachteleben et al. 1992;Eguchi 1995;Jerzak 1995;Dolenec 2000;Soler et al. 2001;Antonov and Atanasova 2003;Ponz and Gil-Delgado 2004;Tucakov and Kucsera 2008;Wang et al. 2008). Overall, these studies showed that egg laying generally occurs from mid-March to mid-April. They also showed that clutch size varied among populations, with an average value ranging from five to six eggs per nest, but fledging success was always less than two fledglings per nest. However, detailed data on the breeding parameters of the endemic Maghreb Magpie are still lacking.
The objective of this study was to increase our understanding of the ecology of this relic population and provide detailed data on its breeding biology over two consecutive breeding seasons. Specifically, we aimed to: (1) describe nesting sites and nest characteristics; (2) determine key breeding parameters, notably clutch size, egg volume, hatching success and fledging success, and; (3) investigate how these parameters varied according to laying date. We expected clutch size, egg volume and chick survival to decrease with laying date, in line with the general trend known in birds (Bengtson 1972;Parsons 1972;Newton and Campbell 1975;Hill 1984;Sedinger 1992;Christians 2002;Verhulst and Nilsson 2008). This trend is more likely to occur in our arid study area because food availability deteriorates rapidly as the breeding season progresses. Moreover, we expected that breeding success would be strongly related to nest position, especially nest elevation, as this parameter is supposed to influence nest conspicuousness and accessibility to predators.

Study area
The relic population inhabits a private, 650 ha, farm located 10 km from Dhorbania village (35° 58′ 16″ N-10° 01′ 15″ E) at Kairouan Governorate, in central Tunisia (Fig. 1). The climate is semi-arid with annual rainfall of 290 mm and an average summer (June, July and August) temperature of 38 °C (Mougou et al. 2011

Nest monitoring
From March to late June in 2017 and 2018, we visited the study site to locate nests of breeding magpies. Located nests were visited weekly during the egg laying period (in order to minimize disturbance) and twice a week during hatching until 25 days post-hatch. All active and unoccupied nests were geo-referenced using GPS (Garmin Etrex 75049).
In 2017, we recorded the height (m), diameter (m) and tree species the nest was in. For each nest, we also recorded nest elevation (ground to nest distance) (m) and nest to tree top distance (m). Nest volume (in liters) was also calculated as: V = (4/3) × (π × a × b 2 )/1000, where a is the largest radius of the ellipsoid nest (cm) and b half nest width (cm) (Soler et al. 1995).

Reproductive data collection
We used an endoscopic camera to check nest contents. For each nest we recorded the date when the first egg was laid, the number of eggs in each nest (i.e. clutch size) and the number of nestlings that survived to 25 days posthatch (age at which nestlings scrambling made further visits impractical); we considered 25 days post-hatching as our operational criterion for fledging. The clutch was considered complete when there were no additional eggs after two consecutive visits. We considered a nest unsuccessful if it was found empty before the estimated hatching/fledging date or with damaged eggs/nestlings, since the incubation period was estimated to 15 days (Ouni 2018). We identified predators based on the nest condition and evidence in the nest. For example, rats typically leave eggshell fragments, while snakes feed on nests without leaving a trace (Pietz and Granfors 2000;Klug et al. 2010). In 2018, we used vernier calipers to measure (± 0.01 mm) maximum length (L) and breadth (B; mm) of 35 eggs from seven clutches. We calculated egg volume (mm 3 ) following Hoyt (1979)

Data analyses
We assessed whether the frequency of use of nesting trees and bushes matched their availability in the study area, using the Ivlev's electivity index: where r i is the relative abundance of the tree/shrub species i in the subset of trees/bushes used as nesting supports, while p i denotes the relative abundance of this tree/shrub species in the study area (Ivlev 1961). Values of E range from − 1 to + 1, with negative values indicating avoidance, positive values suggesting active selection, and values close to zero indicating random use. Moreover, as an investigation of nest position within the nest tree, we compared nest elevation to nest distance from tree top using the Wilcoxon signed-rank test.
We used data collected at active nests (during the two breeding seasons of 2017 and 2018), to determine clutch size (number of eggs laid), hatching success (number of eggs hatched), and fledging success (number of fledglings produced). We assessed whether these parameters varied between the two years of study using Mann-Whitney U-test. Moreover, a series of colony-level parameters were also determined: egg hatching rate (total number of eggs hatched/total number of eggs laid in the colony), nestling survival rate (total number of fledglings/total number of eggs hatched) and breeding success (number of fledglings/number of eggs laid). We also calculated nest survival (number of nests with at least one hatched egg/total number of nests), fledging survival (number of nests raising at least one fledgling/number of nests with at least one hatched egg), and the rate of successful nests (number of nests raising at least one fledgling/ total number of nests). χ 2 tests were then conducted to check whether these rates varied between the two years of study. Linear regressions were carried out to assess whether clutch size, brood size and the number of fledglings varied among nests according to laying date (Julian days, day 1 = January 1st), nest elevation and nest volume.
Finally, we used a Generalized Linear Mixed Model (GLMM) to investigate how egg volume varied with clutch size (three categories: four-egg clutches, five-egg clutches and six-egg clutches) and laying date (days after March 1st), while accounting for clutch identity as a random factor. All statistical tests and analyses were carried out using Statistica 12 software (StatSoft Inc 2013).

Nesting sites
In 2017 we recorded 119 nests: 13 (11%) active nest and 106 old unoccupied nests. One breeding pair has laid two clutches (initial and replacement) in two different nests. In 2018, we recorded 20 (15%) active nests with no replacement clutch. Among these 20 nests we distinguished 4 (20%) restored old nests, 3 (15%) reused nests (from active nests in 2017) and 13 (65%) new built nests. Thus the total number of nests recorded in the study area in 2018 was 132, showing an increase of 11% from 2017 to 2018.
All nests were located in the uncultivated part of the study site and none were found in fruit trees. Nests were mainly built in Searsia tripartita shrubs (121 nests; 91.6%), which is the most available nesting support in the area (76.5% of nesting supports), corresponding to an Ivlev's index value close to zero and suggesting a random use of this shrub species (Fig. 2). Remaining nests were located in Ziziphus lotus (4 nests), Olea oleaster (3 nests), Lycium europaeum (2 nests), Acacia cyclops (1 nest) and Acacia horrida (1 nest). Together, these plant species accounted for 9.6% of potential nesting sites and supported 8.3% of magpie nests. By contrast, no nest was built in Periploca angustifolia, which is the second most abundant shrub species in the study area (10.5% of potential nesting vegetation). This gives a negative Ivlev's index value and indicates avoidance of this shrub species (Fig. 2).
Mean height and diameter of nesting shrubs (± SE) were 2.87 ± 0.43 m (range 2.29-3.90 m) and 4.19 ± 1.17 m (range 2.50-6.80 m) respectively. Nest volume varied between 0.08 and 0.52 m 3 and averaged 0.25 ± 0.04 m 3 . Nest elevation ranged from 1.10 to 2.09 m, and averaged 1.54 ± 0.26 m. The height of the nest in the shrub was significantly higher than the distance from the nest to the shrub top (Z = 3.17, P < 0.01), showing that nests were built above the middle of the shrubs.

Reproductive parameters
Egg laying occurred between mid-March and mid-April, with 69% of nests initiated during the second half of March and 31% in the first half of April.
The breeding parameters remained consistent between years, except for clutch size which was greater in 2017 than 2018 (Table 1).
Over the two years of study, 44% of eggs laid in the colony were lost during incubation and 39% of hatchlings were lost before fledging. The main causes of nest failure were depredation of nestlings by Southern Grey Shrike (Lanius meridionalis) (34% of losses), Egyptian Cobra (Naja haje) (24%) and rats (Rattus rattus) (18%). The depredation of parents by raptors, particularly Blackshouldered Kite (Elanus caeruleus) and nest parasitism by the Great Spotted Cuckoo (Clamator glandarius) also accounted for 6% and 9% of losses respectively. However, we were unable to identify the cause of 9% of failed nests.

Discussion
To our knowledge, this is the first investigation of the breeding ecology of the Maghreb Magpie. Our main objective was to provide detailed data on the nesting habits and breeding success of this declining relictual population. We obtained new and valuable data describing the nesting habitat use and the breeding performance of the Maghreb Magpie. These results are the first and only provided data related to Maghreb Magpie in Tunisia and entire North African region, which is of great interest for the implementation of a possible monitoring and longterm conservation plan of the endangered Tunisian magpie population.
Our results showed that magpie nests were built on thorny bushes, in particular Searsia tripartita which is the most available nest habitat in the area. Non-thorny shrubs, such as Periploca angustifolia, were not used despite their wide availability. In addition, nests were often built just above the middle of the bushes, which is similar to what has been reported in the Eurasian Magpie in temperate rural environments (e.g. Dolenec 2000;Tucakov and Kucsera 2008). This tendency to place nests in the middle of thorny shrubs might be a defense mechanism to limit the accessibility of mammalian and avian nest predators.
We found that clutch size and nest survival in 2018 were lower than those in 2017. Difference seem to be related to weather condition mainly rainfall. Firstly, winter rainfall (September to February) in 2017 (109 mm) was higher than in 2018 (46 mm). According Rodríguez and Bustamante (2003) nest success rate is influenced positively by rainfall in winter. In fact, winter rainfall may increase food availability during incubation since dry winters will force birds to quit incubation to forage when food abundance is low. Secondly, we found that spring rainfall (March to May) in 2017 (16 mm) was lower than recorded in 2018 (34 mm). This is in accordance with Rodríguez and Bustamante (2003) demonstrating that rainfall during the nestling period has a negative effect on nest success rate. Thus, spring rainfall at the time of hatching cause temporal declines in foraging activity (Dawson and Bortolotti 2000;Radford et al. 2001) or in food availability (Avery and Krebs 1984), causing mortality of complete clutches of recently hatched chicks. According Senapathi et al. (2011), spring rainfall affects birds to breed later in the season. Delays in breeding cause reduction of reproductive success as birds get exposed to risks associated with adverse climatic conditions later in the breeding season, which reduce nesting success.
We showed that the depredation of Maghreb Magpie nestlings by Southern Grey Shrike was the main cause of nest failure. This species shared the same breeding habitat with Maghreb Magpie. Contrary to previous reports that state Magpies prey on shrikes nests (Cramp and Perrins 1994;Strnad et al. 2012), we found that shrikes were the most common magpie nest predators. Southern Grey Shrike nests inside thorny bushes showing a great ability to squeeze between thorns. Given the two bird species share the same breeding habitat, we believe shrikes kill the nestlings of neighboring magpie pairs to cause nest abandonment. This appears to be an active defense strategy different from the known behavior of shrikes to prevent predation of their nests by magpies: avoid attacking predators that approach the nest so as not to show its location (Syrová et al. 2016). A more detailed study of the predatory interactions between these two species is warranted.
Our preliminaries results showed that the breeding biology of Maghreb Magpie slightly differs from the Eurasian Magpie by nesting earlier and laying smaller clutches and eggs (Högstedt 1980;Tatner 1982;Arias de Reyna et al. 1984;Redondo and Carranza 1989;Ponz and Gil-Delgado 2004;Birkhead 1991). Differences in laying date may obey to a latitudinal gradient effect on birds breeding phenology (Lack 1947(Lack , 1968Cardillo 2002;Jetz et al. 2008). More data and further analysis are needed to verify this hypothesis. The earlier spring in the latitudes of North Africa in comparison with the temperate zones of Tunisia would explain the early entry into reproduction of the Maghreb Magpie compared to the populations of Eurasian Magpie. According de Neve et al. (2004) food availability directly influences laying date in many species, including the Eurasian Magpie, where food supplementation significantly advanced the onset of laying (Meijer and Drent 1999).
It is known that egg size of magpies significantly increases in response to food supplementation (Högstedt 1981;Hochachka and Boag 1987;Ramsay and Houston 1997;Meijer and Drent 1999;de Neve and Soler 2002;de Neve et al. 2004) despite the relatively inflexibility within individual female birds showed by Christians (2002). The smaller clutches and eggs recorded in our population in comparison with those known in the Eurasian Magpie in temperate zones may be due to possible food scarcity in our arid study area and/or the shorter days at our latitudes. Unfortunately, we did not assess food availability in our study area to test this hypothesis and to investigate more deeply the possible role of food availability in shaping the reproductive investment of the studied birds.
Egg volume varied significantly with laying date, being larger earlier in the season. This pattern matches the general trend known in birds (Bengtson 1972;Parsons 1972;Newton and Campbell 1975;Hill 1984;Sedinger 1992;Christians 2002). This might be due to poorer quality of late breeders compared to early ones. Indeed, some young and less experienced birds often start breeding later and produce less-quality eggs compared to more experienced birds (Christians 2002;Williams 2012;Verhoeven et al. 2019). Alternatively, deteriorating breeding conditions, notably food availability, as the season progress may also result in decreasing egg volume. These hypotheses were not investigated since no data on monitored breeding pairs ages and food availability seasonal changes were available.
Because the number of eggs laid by a female could be traded-off their sizes, a negative relationship between egg size and clutch size could be expected (Blackburn 1991;de Neve and Soler 2002;Ręk 2010). The opposite trend could also be expected as high-quality females may lay larger clutches of larger eggs compared to lowquality females (de Neve and Soler 2002;de Neve et al. 2004;Lifjeld et al. 2005). Given that no significant relationship between egg volume and clutch size was found, our results did not give support to any of the previous hypotheses.

Conclusion
In conclusion, we believe that without claiming to have completely filled a lack of information, our study provides the first and only detailed data on the reproductive parameters of the Maghreb Magpie in Tunisia and whole North African region. In addition to their potential importance for a possible large-scale comparative study of the reproductive ecology of the magpie complex, these data could be of great interest for the implementation of a possible monitoring and long-term conservation plan of the endangered Tunisian magpie population. In this context, the results of our study would suggest that the preservation of thornier bushes in protected areas, such as the private farm where the study population currently live, could ensure the increase of population size by providing it with more suitable and safe nesting habitat. This would require a short and medium-term conservation program, involving the responsible Tunisian authorities (i.e. general directorate of forests) but also local environmental NGOs, and aiming at sensitizing local farmers and helping them to dedicate parts of their lands to thorny bush plantations.