Skip to main content

Vaginal temperature and collar activity variations in relation to fecal progestagens of captive moose

Abstract

Background

Timing of reproductive events can be crucial for a species’ population growth and stability. Accurate detection of reproductive phenology presents a challenge to scientists studying wild species, including moose (Alces alces). Currently, there are several established methods for monitoring reproductive activity and events in domestic ruminants, including the use of biologging devices. The main objective of this study was to determine whether female moose display a distinct thermal and activity pattern associated with luteal activity during the estrous cycle, which could be used to determine the onset of their breeding season. We deployed biologging devices and collected fecal samples from 12 captive female moose on the Kenai Peninsula, Alaska, USA to explore variation in vaginal temperature and collar activity and the relationship between these variables and fecal progestagen concentrations. Fecal samples were collected from mid-August to mid-October and analyzed using radioimmunoassay to determine the concentration of fecal progestagens to classify luteal activity.

Results

Captive female moose displayed an identifiable thermal pattern during the onset of luteal activity from mid-September to mid-October, associated with the initial estrous cycle of their breeding season. In contrast, we did not observe a distinct pattern in activity during this period. Recurring patterns in both vaginal temperature and activity were identified between mid-October and mid-November, however, which were likely associated with subsequent estrous cycles but not included in our fecal sampling period.

Conclusions

This study supports that female moose display an identifiable pattern in vaginal temperature which is associated with luteal activity of the initial estrous cycle of the breeding season. An identifiable pattern was observed for both vaginal temperature and activity registrations at the presumed timing of subsequent estrous cycle of the breeding season.

Background

Understanding how environmental factors affect populations of wild species and predicting the expected ecological impacts is crucial to consider in management and conservation of wildlife [1]. Environmental fluctuations caused by climate change may shift the timing of phenological events and consequently impact the general reproductive success and rate of survival of wild animals and plants [2]. Different moose (Alces alces) populations are experiencing fluctuations, and recent studies show how populations of moose living in the southern ranges of their habitat are declining, with a reduced calf survival rate [3, 4], while other populations appear to be relatively stable [5], expanding to new areas [6], or increasing in size [7, 8]. Studies in Scandinavian [9,10,11] and North American [12,13,14] moose have shown that seasonal climatic changes, physical condition, fertility problems, predation, forage nutritional quality and availability impact reproductive success and calf survival. It has been hypothesized that calves which are conceived later in the breeding season (i.e., not during the first estrous cycle) may be born when forage quality and availability are suboptimal, resulting in reduced calf growth and probability of survival [15]. Furthermore, calves born later are likely more vulnerable to predation [16]. Thus, a deeper understanding of the reproductive phenology and physiology of moose would benefit research aimed at studying population fluctuations.

The estrous cycle and breeding season in Scandinavian and North American moose is reported to occur from early September and into November [17,18,19]. Moose are seasonal breeders [3], meaning that they only breed during a specific time of the year. Changes in day length (photoperiod) affect endocrine signaling pathways which are known to regulate the transition from reproductive inactivity (anestrus) to reproductive activity (estrus) in seasonal breeders [20]. Sexual maturation in female moose is speculated to occur after 1.5 years of age but varies between different cohorts and their respective environmental conditions [3, 17, 21,22,23]. Peak ovulation dates for Scandinavian moose [21] are estimated to be synchronized and occur in late September and early October with regional variations, which align well with when the majority of mating events occur during the rut in North American moose [24]. Schwartz et al. [17, 25] estimated each estrous cycle to last about 24 (range = 22–28) days and a possible recurrence of 4–7 cycles during the breeding season. Duration of sexual receptivity, or estrus, is limited to 1–36 h at the beginning of each cycle. Malmsten et al. [18] found that for a limited number of female moose mating had likely occurred 2 weeks into pregnancy, meaning that copulation may not be a reliable indicator of estrus. In addition, the timing and relation between observed estrus and ovulation are not yet known for moose, so it is not possible to estimate a day of ovulation solely based on external signs of estrus. It is, however, likely that ovulation occurs during or shortly after estrus [3].

Today, we have well-established methods for detecting and monitoring reproduction in domestic ruminants, including behavior or physiological observations, measuring progesterone hormone, or monitoring variables, such as body temperature and activity [26, 27]. In contrast, we have less experience and few established methods for monitoring the estrous cycle in wild animals, including moose. As behavioral cues in wild ungulates may be less prominent, the reliability of behavioral observations to determine or validate reproductive events in wild species are considered to be less accurate [28]. This further supports the need for more accurate reproduction monitoring techniques in wildlife. Captive wildlife is more habituated to human presence, which facilitates the feasibility and accuracy of reproduction monitoring and decreases capture-related stress. Therefore, reproduction studies in wild species have initially been conducted in captive animals, before being extended to free-ranging ones [29].

To overcome practical and ethical challenges associated with studying non-domesticated animals, less-invasive methods are essential. Measuring fecal steroid hormones associated with reproductive physiology offers an alternative to the more invasive method of acquiring a blood sample to monitor reproduction in wildlife species. Progesterone is a steroid hormone associated with the estrous cycle and is produced by the corpus luteum in the ovary before it is released into the blood, where it circulates either freely or bound to plasma proteins. Further metabolization occurs in the liver, before excretion as conjugates in urine, bile (and subsequently feces), and saliva [30]. Measuring progesterone metabolites in feces (i.e., fecal progestagens) is a widely used method to monitor reproductive physiology in several wild ruminant species [31,32,33,34,35], including moose [25]. As an illustration, estrus in captive moose has been suggested to occur on average 0.6 days prior to the nadir (95% confidence interval (95% CI) = 4 days before to 3 days after the nadir) in fecal progestagen levels, by visual confirmation of estrus such as copulation with a bull or ruffled rump hairs as an indicator of having been mounted [17, 25].

Previous monitoring studies have shown that moose display seasonal patterns of body temperature and activity [36, 37]. In general, body temperature decreases from July to September in Scandinavian moose [36], and from July to October in Alaskan moose [37], and increases from April to June [36, 37]. In Scandinavian moose, activity is highest in June and lowest in February [36]. These patterns are likely associated with seasonal hypometabolism. Moreover, biologging devices have recently been used to investigate patterns in body temperature and activity during pregnancy and parturition in wild Swedish moose [38], revealing that body temperature is higher in pregnant moose compared to non-pregnant moose [37]. Furthermore, monitoring the combination of body temperature and activity changes can help detect the time of calving, an event associated with a drop in both variables [39]. These findings illustrate that moose show identifiable patterns in body temperature and activity related to reproductive events, similar to what has been reported in domestic cattle [40,41,42]. However, a precise method to accurately determine the onset of estrous cyclicity in free-ranging moose is yet to be established and would increase researchers’ abilities to evaluate consequences of a shift in reproductive phenology, such as estrous cyclicity, breeding, and conception.

This study describes the application of novel reproductive monitoring techniques in female moose, with a general aim to broaden our knowledge about their reproductive phenology and physiology.

The main aim of this study was to explore whether vaginal temperature and collar activity data from captive female moose could be used to determine the onset of their breeding season. In dairy cattle it has been documented that body temperature decreases just prior to estrus, and then sharply increases on the day of estrus [40, 43]. Following this peak there is a decrease around the expected time of ovulation before a steady increase post-ovulation, during the luteal phase. The rise in body temperature around estrus has been attributed to the thermogenic effect of progesterone [42], increased vaginal blood flow [44], and the accompanying increase in activity [40]. An increase in activity is associated with the timing of estrus in dairy cattle [45,46,47,48], and in wild ungulate species, such as white-tailed deer (Odocoileus virginianus), where peak daytime movement usually occurs on the day of mating or the day before [49]. Based on these findings and use of similar methodology, we made the following predictions related to the estrous cycle in moose. Since it has been found that fecal progestagen values can be utilized to classify luteal activity during the estrous cycle in moose [25], and the study area and animals enabled daily fecal collection, we predicted that this variable could be used to detect and confirm the initial luteal activity of the breeding season (P1). In addition, we predicted that there would be an identifiable pattern in vaginal temperature (P2) and collar activity (P3) associated with luteal activity in moose. Our expectations for P1 through P3 were that the three variables would share similarities with physiological changes which are associated with specific parts of the estrous cycle and ovarian activity in dairy cattle; a low basal level of plasma progesterone preceding the onset of luteal activity, combined with an increasing trend in vaginal temperature and a peak in activity [50]. Dairy cattle were selected as a reference species for the physiological changes that occur during the estrous cycle due to the amount of rigorous research on validation of activity and body temperature to detect ovarian activity which is available compared to other ungulates that may share more similarities with moose.

Methods

Study area and animals

We collected data registrations from 12 captive female moose, ages 2–19 years, at the Kenai Moose Research Center, on the Kenai Peninsula, Alaska, USA. The Kenai Moose Research Center is operated by the Alaska Department of Fish and Game and is located on the Kenai National Wildlife Refuge (60°43’ N, 150° 26’ W). The study area is characterized as boreal forest, including dense forests, wetlands, and open meadows. Present tree species are dominated by white spruce (Picea glauca), black spruce (Picea mariana), Alaska birch (Betula neoalaskana), Scouler’s willow (Salix scouleriana), and quaking aspen (Populus tremuloides). The study animals were kept outdoors in two 2.6 km2 enclosures with free access to water and natural forage, which enabled the desired frequency of individual-specific fecal sample collection. Each individual had a specific color combination of duct tape on their GPS collar for identification and fecal pellets were only collected after confidently identifying the animal and defecation had been observed. The fence kept the study animals physically separated from wildlife in the surrounding area. Occasionally wild bulls were observed along the fence lines of both enclosures and one of the enclosures also shared a small section of fence with a pen that housed three adult bulls, both of which allowed a certain degree of indirect contact between females and males.

In May and July 2021, the captive female moose were immobilized following an earlier described protocol [37] with a combination of Thiafentanil oxalate (0.001–0.004 mg/kg estimated body mass; 10 mg/mL; Wildlife Pharmaceuticals Inc., Windsor, CO, USA) and Xylazine (0.03–0.05 mg/kg estimated body mass; 100 mg/mL; Lloyd Laboratories, Shenandoah, IA, USA) by intramuscular hand-injection and equipped with Vertex Plus Global Positioning System (GPS) collars (Vectronic Aerospace GmbH, Berlin, Germany) and Vaginal Implant Transmitters (VITs; size = big; Vectronic Aerospace GmbH; Berlin, Germany). Immobilization was reversed by a combination of Atipamezole HCl (0.005 mg/kg estimated body mass; ¼ dose intravenous, ¾ dose intramuscular; 5 g/mL; Zoetis, Parsippany, NJ, USA) and intramuscularly administered Naltrexone HCl (100 mg/mg Thiafentanil oxalate intramuscular; 50 mg/mL; ZooPharm LLC, Laramie, WY, USA). GPS collars and VITs were manually removed from all animals without sedation on November 23rd, 2021.

Vaginal temperature and collar activity

The VITs transmitted vaginal temperature (Tv; °C) every 5 min to the GPS collar and were additionally programmed to record vaginal temperature and store this data on-board every 17:10 min as a backup should transmission fail. During the aforementioned immobilization procedure, a VIT was inserted into the vaginal canal of each individual using a lubricated (OB Lube; Jorgensen Laboratories Inc., Loveland; CO, USA), sterilized speculum (Sterile Disposable Vaginal Speculum; Jorgensen Laboratories Inc.) following previously established procedures [51]. Tv was successfully downloaded from 11 individuals (from GPS-collar n = 10 and from the VIT n = 1). One VIT stopped recording on the 12th of August 2021 (i.e., before the start of fecal sample collection), and so this individual was excluded from further Tv analysis. Over 150,000 registrations of Tv were made during the dates with corresponding fecal progestagen values (i.e., August 23rd to October 15th, 2021). Daily mean Tv was calculated, resulting in 594 registrations during this period.

The GPS collars included a triaxial accelerometer which registered forward–backward, left–right, and up-down movements, stored as X, Y, and Z, respectively. Movement was recorded in 5-min intervals as average values of each axis ranging between 0 and 255 at 6–8 Hz and stored in the collar. Overall activity was calculated by summing the X, Y and Z axes, giving values ranging from 0 to 765 every 5 min, with 0 representing no or low activity and 765 the highest activity level. Activity registrations were successfully stored for all 12 individuals resulting in over 185,000 registrations for each axis (X, Y, and Z) during the period of fecal sampling (i.e., August 23rd to October 15th, 2021). Daily summed activity during the fecal sampling period was calculated, resulting in 234 to 288 registrations a day.

VIT validation

We validated the accuracy and precision of the VIT temperature logger within the range of moose body temperature (i.e., within the range of vaginal temperature in this study and prior research [37]). We placed the VITs in a warm water bath heated to 37.00 °C, and we then increased the temperature of the water bath by 0.50 °C until 40.50 °C. The water bath temperature was measured with a National Institute of Standards and Technology-certified thermometer (Traceable® Digital Thermometer Model 90080–09, accuracy ± 0.05 °C, 0 to 100 °C, resolution 0.001 °C (Webster, TX, USA)). We also tested the accuracy and precision of the VIT in an ice bath, which has been traditionally used to calibrate temperature logger VITs [52]. We used linear mixed model regression, with individual VIT as a random intercept, to determine accuracy of the loggers.

Fecal samples

A total of 468 fecal samples were collected from mid-August to mid-October 2021. The study period was selected based on previous literature on the timing of the species’ reproductive season, personal observations of the study population’s estrous behavior in previous years, and a desire to collect samples during the transition between anestrus to estrus [17, 21, 24]. Individual samples were collected every other day from the 23rd of August to the 21st of September 2021. After this, samples were collected daily until the 16th of October 2021. All sampling occurred between 07:00 and 19:00 Alaska Daylight Time. Each sample weighed between 100 and 200 g and was collected in pre-labeled Whirl–Pak bags (Nasco Whirl–Pak®, Fort Atkinson, WI, USA) during or shortly after defecation and stored on ice in a cooler bag, before being transferred to a −18 °C propane freezer. The frozen samples were freeze dried (Labconco model 7752020, Kansas City, MO, USA) and homogenized to a powder to further distribute metabolites evenly before randomly subsampling 5 g [53]. The subsamples were shipped to an endocrine laboratory (Applied Biosciences, TX, USA) for further extraction and analysis, following a previously described protocol [54]. Progestagen concentration was determined by radioimmunoassay (RIA; Catalog #07-270102; ImmuChem Double Antibody, 125I RIA Kit, MP Biomedicals, Costa Mesa, CA, USA), following the manufacturer’s protocol. Intra-assay variation was 10.6% on a low pool, medium pool at 7.93%, and 5.21% on a high pool. The extraction efficiency average on spiked samples of fecal material was 83.23% (n = 10). Fecal progestagen concentration is expressed as nanogram (ng) of immunoreactive fecal progestagen hormone metabolites per gram (g) dry fecal weight.

Environmental variables

Within the Kenai Moose Research Center there is a National Oceanic and Atmospheric Administration (NOAA), U.S. Climate Reference Network weather station which recorded environmental variables, including ambient air temperature (Ta, °C), relative humidity (%), solar radiation (W/m2), windspeed (m/s), and precipitation (mm) every 5 min [55]. Ambient air temperature and relative humidity recorded by the NOAA weather station were used to calculate dew point temperature (°C), which was further used to calculate actual vapor pressure (hPa) [56]. The NOAA weather station had 1 day of missing data for ambient temperature, relative humidity, and solar radiation (25th of August, 2021). HOBO dataloggers (HOBO U23 Pro v2 Temperature/Relative Humidity datalogger; HOBO Pendant Temperature/Light 64 K datalogger; Onset Computer Corp., Pocasset, MA, USA) were deployed adjacent to the NOAA weather station for the duration of the study. We used simple linear regression to estimate the missing NOAA values from the HOBO logger data for daily mean ambient air temperature (y = 0.96*Ta + 0.50, r2 = 0.99), daily range in ambient air temperature (0.88005*TaR + 0.49969, r2 = 0.98) relative humidity (y = 1.05859*RH −9.79112, r2 = 0.98), and solar radiation (y = 0.0004*lux—0.1266, r2 = 0.98).

Data preparation and analysis

All data handling and statistical analysis was performed in R version 4.1.0 (R Core Team, 2021), using RStudio version 1.4.1717. Means are reported with ± SD in parenthesis. Tv and collar activity data were filtered for the 3 days following capture to exclude values which would likely be elevated due to lingering effects of the immobilizing drugs or capture-related stress [37]. Since fecal collections occurred on alternate days for the first half of the sampling period, fecal progestagen values were linearly interpolated using the “approx” function from the “zoo” package [57] in R to generate daily values to classify luteal activity. Similarly, daily values for Tv and collar activity were generated by calculating the daily mean Tv and the total daily sum of activity values (referred to as the daily sum of activity).

To determine a baseline level of fecal progestagen considered as indicative non-luteal activity, we calculated the mean of the 2 lowest values for each individual. Furthermore, a threshold was created by multiplying the baseline level by 2, to differentiate between non-luteal activity and luteal activity. If fecal progestagen levels were below this threshold for at least 1 day, and then increased, and stayed above it for at least 14 days, we defined the day on which the threshold was exceeded as an onset of luteal activity (OLA), and the period during which fecal progestagen levels stayed above the threshold as luteal activity. Because the duration of the luteal phase is currently unknown for moose we based our OLA classification method on existing approaches for classifying normal luteal activity in cattle [58], along with the average luteal phase length in other seasonal breeding ruminants, such as marsh deer (Blastocerus dichotomus), domestic sheep (Ovis aries), and domestic goats (Capra hircus; 15, 14 and 16 days, respectively [34, 59, 60]). The date of the first OLA was assigned as day 0, which was used as a scale to explore variations in Tv and collar activity for that individual in relation to luteal activity. The above calculations and classifications were done by applying a rolling function using the “zoo” R package [57].

The four moose which we were not able to determine an OLA for either had erratic fecal progestagen profiles (age = 2 years), fecal progestagen levels which did not stay above the set threshold for the minimum required days, or an erratic fecal progestagen profile which stayed above the set threshold but was not preceded by a day in which the level was below the set threshold (age = 19 years). Tv and activity registrations from these individuals were, therefore, excluded from further statistical analyses. These individuals were either below or above prime age (2.5–10.5 years [21]), indicating that our method is more successful when applied on prime aged moose.

Relationship between luteal activity and vaginal temperature

To determine when changes in daily mean Tv occur in relation to the onset of luteal activity, and to account for non-linear relationships between response and explanatory variables, we used Generalized Additive Mixed Models (GAMMs; using the “mcgv” R package from Wood [61]). All explanatory variables included in the model were adjusted to daily values to match up with the daily frequency of fecal progestagen values. Daily mean Tv was selected as the response variable and OLA was included as the main explanatory variable to explore the relationship to luteal activity. In addition, explanatory variables that are known to influence core body temperature in moose were included: daily sum of activity, ordinal day (to account for seasonal variation in Tv [36]), daily mean Ta (°C), daily range of Ta (°C), daily mean vapor pressure (hPa), daily mean wind speed (m/s), total daily precipitation (mm), and total daily solar radiation (W/m2) [37]. Correlation and structure of the explanatory variables were checked using the function “ggpairs” from the “GGally” [62] extension to the “ggplot2” R package [63]. If the Pearson’s correlation coefficient between OLA and another explanatory variable was above 0.50 or below -0.50, it was considered a high correlation [64] and this variable was not included together with OLA in the same model (this was the case for ordinal day, vapor pressure and daily ambient temperature, see Additional file 1). In addition, if two explanatory variables other than OLA had a correlation with each other above or below this threshold, the combination of these variables was not included in the same model (this was the case for daily ambient temperature and solar radiation, and wind speed and relative humidity, see Additional file 2).

These steps left us with the following explanatory variables: OLA, daily range of Ta, daily mean relative humidity, total daily precipitation, daily mean wind speed, total daily solar radiation, and daily sum of activity (Additional file 9: Table S1). In addition, to test if the trend in Tv was attributed to seasonal variation rather than OLA, we included ordinal day in one of the models. We included the individual moose (“CollarID”) as a random intercept in all models to control for inter-individual variability and repeated measurements. A gaussian distribution with an identity link function and the maximum likelihood was used. We then ran a model selection based on Akaike’s Information Criterion, adjusted for small sample sizes (AICc, function “AICctab” from the “bbmle” R package [65]), by selecting the most parsimonious model within ΔAICc ≤ 2 [66]. Residuals were assessed and basis dimensions for the parameter k were checked before we inspected diagnostic plots to validate model assumptions. Finally, the selected model was refitted with restricted maximum likelihood.

Relationship between luteal activity and collar activity

GAMMs were also used to determine when changes in daily sum of activity occur in relation to the onset of luteal activity. Daily sum of activity was selected as the response variable and OLA was included as the main explanatory variable to explore the relationship between luteal activity and activity values. The same environmental variables mentioned previously were included as they have been shown to affect behavior and habitat selection in moose [67]. In addition, we included daily mean Tv as an explanatory variable since thermoregulation influences activity levels in moose [67], and the same random intercept (individual moose, “CollarID”). Due to positive-skewed data we used a gamma distribution with a log link function and applied the same model selection process as described above.

Results

Classification of luteal activity

Fecal progestagen metabolite concentrations observed in the 468 fecal samples ranged from 2.20 to 68.30 ng/g fecal matter (15.30 ng/g ± 9.30 ng/g; Additional file 9: Table S2). The classification method used to differ between luteal activity and non-luteal activity was able to identify a date for the onset of luteal activity (OLA) for 8 out of 12 individuals (Fig. 1). The OLA had a mean date of occurrence on the 23rd of September (± 4 days) and showed no significant interindividual variability (ANOVA, P = 0.59).

Fig. 1
figure 1

Fecal progestagen profiles. Fecal progestagen (ng/g) profiles from captive female moose (n = 12) from August 23rd to October 15th, 2021, at the Kenai Moose Research Center, AK, USA. The bold outlined graphs indicate the seven individuals with fecal progestagen profiles which fulfill the criteria to determine a date for the onset of luteal activity (OLA). The remaining five graphs are for moose with fecal progestagen profiles which did not fulfill the criteria to determine a date for OLA or did not have vaginal temperature data for the whole duration of the study period (moose 30680) and was, therefore, excluded from further data analysis. Dashed red horizontal line = individual progestagen threshold. Dashed blue vertical line = date for OLA

Vaginal temperature and collar activity

Vaginal implant transmitters (VIT) placed in a warm water bath (37.00–40.50 °C) had an accuracy of 0.25 °C (95% CI 0.23–0.27 °C), and the variance associated with the random intercept of VIT was insignificant (σ2 < 0.01). When the VITs were placed in an ice bath, we found that the accuracy was much lower at 1.23 °C (95%CI 0.33–2.15 °C), and the variance of the random intercept of VIT was significant (σ2 = 1.26). In both water baths, the VITs registered temperatures that were colder than the actual temperature of the water bath. Across all individuals, the registered daily mean vaginal temperature (Tv) declined from late August to late September. From late September to early October periodic variations in daily mean Tv were observed (i.e., a discernible increase over a two to 3-week period followed by a 1-day decrease). Each moose had at least one such characteristic Tv pattern during the fecal sampling period. Additional patterns in Tv were observed for eight individuals on dates exceeding the fecal sampling period. The interval between the first and second observed pattern in Tv was significantly shorter than the interval between the second and third observed pattern (paired T test, P = 0.02, 24 ± 2 days vs. 27 ± 2 days, respectively).

In general, collar activity (referred to as the daily sum of activity, see Methods) declined from August towards November. Two moose had discernible but relatively small peaks of activity during the fecal sampling period. However, a total of nine moose had one or two more obvious activity peaks from mid-October until the end of activity registrations in November, which were found 27 ± 1 day apart from each other.

For the seven moose which had Tv, collar activity, and an OLA available during the fecal sampling period, the general observation was that Tv started to increase a few days before or after day 0 of OLA (range = 7 days before to 3 days after OLA), whereas activity values did not appear to show as much of a distinct pattern around this time (Fig. 2, for the remaining individuals see Additional file 3, Additional file 4, Additional file 5, Additional file 6, Additional file 7, and Additional file 8). Activity peaks were 1.40 to 2.70 times greater (1.80 ± 0.40) than the individual mean of activity values during the fecal sampling period.

Fig. 2
figure 2

Vaginal temperature, collar activity, and fecal progestagen. Graphs from two captive female moose that were 9 years (IDs = 30681 and 30684) containing daily values for mean vaginal temperature (blue line) and daily sum of activity (red line) recorded between July 23rd and November 22nd, 2021, and fecal progestagen concentrations recorded between August 23rd to October 15th, 2021. Dashed vertical blue line = onset of luteal activity, which is on the 23rd of September for both graphs. Dashed horizontal red line = individual fecal progestagen threshold

On dates exceeding the period of fecal sample collection (i.e., after October 15th), we observed that the drop in Tv would periodically align with a peak in activity (a total of 15 times, two times for six individuals and one time for three individuals). The peak in activity was generally seen on the same day as the drop in Tv, or 1–2 days after, with an interval of 27 ± 2 days. Activity peaks were 2.30–6.80 times greater (4.00 ± 1.40) than the individual mean of activity values for this period, which was significantly higher than for the previous period that included corresponding fecal progestagen values (paired T test, P < 0.01).

Relationship between luteal activity and physiological changes

Registrations from a total of seven moose which had Tv, activity and an OLA available were included in the statistical models. The highest ranked model for predicting Tv included only the explanatory variable OLA which explained 72.30% of the deviance (Additional file 9: Tables S3 and S4). This model predicted that Tv dropped by 0.18 °C over a 17-day period, before stabilizing at 37.77 °C (95% CI 37.74–37.80 °C) 5 days before OLA (Fig. 3). Subsequently, Tv increased by 0.17 °C over a second 17-day period before entering a new decreasing phase after reaching 37.94 °C (95% CI 37.91–37.97 °C) 12 days after OLA.

Fig. 3
figure 3

Model predictions. Predicted values of a generalized additive mixed model predicting daily mean vaginal temperature from captive female moose plotted against days in relation to the onset of luteal activity (OLA). Black line = the predicted daily mean vaginal temperature over time. Shaded gray ribbon = 95% confidence interval. The colored lines represent the raw data that the model based its predictions on, and each color indicates each moose that contributed to the data (n = 7). Dashed vertical blue line = day 0 for OLA

The highest ranked model for predicting activity included the explanatory variables OLA and total daily precipitation (Additional file 9: Tables S5 and S6) which explained 39.30% of the deviance. This model predicted that activity decreased by 19.87% from 15 days before OLA (95% CI 6666.09–7910.89) to 12 days after OLA (95% CI 5340.82–6339.85), and then increased by 13.44% from 12 days after OLA to 20 days after OLA (95% CI = 5959.89–7582.65). In addition, when daily mean precipitation increased from 0.00 mm to 17.80 mm, activity increased by 15.04% (0.00 mm 95% CI 5902.4–6979.4, 17.80 mm 95% CI 6579.27–8673.83).

Discussion

The goal of this study was to gain better insights about the moose estrous cycle and investigate the timing of luteal activity using biologging technology. We observed periodic patterns in Tv from mid-September and into November in prime-aged moose (2.5–10.5 years [21]). Patterns in collar activity were visible in October and November (supporting P3) but were less obvious or absent in September (not supporting P3), during the time when fecal progestagen data were collected. There continue to be several aspects regarding the estrous cycle of moose which are yet to be described but our results align well with existing literature on the reproductive characteristics of moose [3, 17, 18, 21] and provide new perspectives on physiological variables during the breeding season and their relation to the estrous cycle.

As anticipated, we were able to classify the initial luteal activity of the estrous cycle in captive moose based on fecal progestagen values (P1). Since the study period lasted for 54 days, combined with cycle length in moose (22–28 days [17]), we predicted that we would observe at least one estrous cycle per moose, as well as the transition from seasonal anestrus to estrus. Based on the early timing of fecal sample collection (i.e., starting date on the 23rd of August), we likely captured the first estrous cycles of the season, as cycling in moose has been reported to occur between September and November [17,18,19, 21, 24]. This assumption is strengthened by the concentrations of fecal progestagen during the first half of the sampling period (i.e., from the 23rd of August until mid-September). During this time, progestagen levels stayed below or near the individually calculated threshold and thus deviated from the expected periodic pattern that is associated with multiple estrous cycles in captive moose [25].

We found that the timing of OLA was not significantly different between individuals, which supports the theory that female moose are synchronized in respect to estrous cyclicity [21, 25]. In Scandinavia, 95% of ovulations in wild moose have been reported to occur within an interval of less than 10 days from late September to early October [21]. On a global scale, most mating events in moose have been reported to occur over a 15-day period, from the 23rd of September to the 8th of October [24]. The mean date for OLA in the present study (the 23rd of September) also corresponds well with previously reported dates for the estrous cycle and breeding season in moose [17,18,19, 21, 24].

Monitoring progesterone concentrations alone is not considered sufficient to predict the specific timing of ovulation in cattle [68]. Therefore, without a method to specifically determine if the study animals were in estrus, such as observing behavioral cues like copulation with a bull [25], measuring luteinizing hormone (LH), or examining the ovaries through transrectal palpation and/or ultrasonography [69, 70], we have taken care not to draw any assumptions regarding estrus or ovulation. Our focus was to use the data to identify when progestagen concentrations were indicative of luteal activity. Future studies should include additional methods to estimate the timing of estrus and ovulation in captive moose, such as registration of behavioral cues, directly through observation or indirectly using devices that detect when a female is mounted [71]. Another strategy could be the application of remote blood collection technology [72,73,74,75] in captive moose to directly measure levels of hormones associated with the estrous cycle (i.e., LH, progesterone, estradiol, and follicle-stimulating hormone [76]).

The VITs used in this study had a lower accuracy than what was reported by the manufacturer (0.25 °C versus 0.10 °C, respectively). Furthermore, although calibration in ice water has previously shown no need for adjustment for this VIT model [52], we found this was not a good method for evaluating the accuracy of the VITs and so we recommend that temperature loggers are validated within the range of body temperatures being measured in the species of study. The observed inaccuracy should be taken into account when interpreting our results. When exploring patterns in Tv and daily sum of activity in relation to fecal progestagen and luteal activity, our expectations were met to a certain degree (P2 and P3). An increase in daily mean Tv was observed around the same time as fecal progestagen concentration increased (in line with P2). There were a few identifiable peaks in the daily sum of activity during the fecal sampling period, though larger peaks were observed afterwards (P3). However, when exploring both Tv and activity registrations which occurred after the fecal sampling period, a periodic peak in activity was found in close relation to when Tv increased (P2 and P3). The interval of this pattern (about 27 days) aligns well with the expected duration of the estrous cycle in moose [17]. Moreover, we found that the interval length between the first two Tv patterns (between September and October) was significantly shorter than the following two (between October and November). If the interval between observed Tv patterns reflects the length of the estrous cycle in moose, as it is reported for cattle [77], it is possible that the initial estrous cycle of a breeding season in moose may be shorter than the following ones, as has been seen in other seasonal breeders [78,79,80].

A possible explanation behind fewer observed patterns in activity is that accelerometry values as a measure of estrous behavior is presumably less appropriate for a free-ranging species like moose (captive or wild) compared to domestic cattle which are often kept in a more confined setting during estrous monitoring. Furthermore, the activity values included in our model were summed to daily values, meaning that higher resolution variations within a day were lost, and since estrus is limited to a short interval (1–36 h [17]), using a mean daily value could reduce the representability of estrus-associated variations. These uncertainties made it challenging to include collar activity data in a statistical model and interpret the results. To increase the interpretability and application of collar activity, future studies should investigate the relationship between accelerometry registrations and specific behavioral observations in moose, similar to what has been done for other species of both domestic [81, 82] and wild ruminants [83, 84].

The high explained deviance (72.30%) of the model with daily mean Tv as a response variable supports that there is a relationship between luteal activity and the identifiable pattern which was observed for Tv (P2). According to our model predictions, the rise in collar activity occurs 18 days after the predicted rise in Tv, and on an individual level each observed activity peak tended to occur following an observed increase in Tv. Although the model with activity as a response variable had a lower explained deviance, this indicates that the elevation in Tv is likely not associated with an increase in activity but rather an increase in progesterone hormone, which is known for its thermogenic effects [85].

Despite a high explained deviance, it appeared difficult to generalize the predictions of our model as some moose exhibited a rise in Tv either before or after OLA (range = 7 days before to 3 days after OLA), whereas our final model predicted that Tv would rise 5 days before OLA. Delay between production and excretion of fecal progestagen [86], sample collection, subsampling, and assay type are all factors that contributed to this inaccuracy. In summary, even though our study demonstrates patterns in Tv which are associated with luteal activity, this variable should not be utilized as a general predictor for luteal activity in moose without the support of other methods to monitor reproductive physiology and behavior.

In contrast, the explained deviance (39.30%) of the model with collar activity as a response variable along with OLA and total daily precipitation as explanatory variables indicates that activity is a weaker predictor for the initial onset of luteal activity in captive moose during the first estrous cycle of the season (P3). Redden et al. [47] found that Tv and activity monitoring, using pedometers, yielded similar estrus detection rates in cattle (81% and 80%, respectively) and some authors have even found that monitoring mounting activity can identify low-intensity estrous behavior which would otherwise be labeled as “silent ovulation” [87]. However, others have found that monitoring temperature can yield better detection rates than monitoring activity [88, 89], which also seems to be the case in the present study. Perhaps activity would have been a good predictor for the subsequent estrous cycles if we had been able to include data which exceeded the fecal collection period, when significantly higher activity peaks were observed.

Conclusion

This study supports the presence of an identifiable pattern in Tv associated with luteal activity in moose, and that succeeding Tv and activity patterns may be indicative of luteal activity in succeeding estrous cycles. These findings offer a new perspective on physiological variables that may be further explored to answer questions about reproductive phenology in moose and indicate that methods which are applied in domestic ruminants may be applicable to wild ruminants as well. Without further research involving more specific and/or additional combinations of methods to monitor reproductive physiology and behavior, these variables should not yet be utilized as general predictors for luteal activity in moose. Additional studies which include a method for confirmation of estrus, ovulation and/or copulation should be conducted to investigate how these events align with our findings before developing methods to monitor reproduction in wild moose. Including this information would enhance our current knowledge about reproductive characteristics of moose and serve a pivotal role in making well-informed management and conservation decisions for the species.

Availability of data and materials

The data sets used and analyzed for this study are available on reasonable request from the Alaska Department of Fish and Game and subject to a data sharing agreement.

Abbreviations

GAMM:

Generalized additive mixed model

GPS:

Global positioning system

LH:

Luteinizing hormone

OLA:

Onset of luteal activity

Ta :

Ambient temperature

Tv :

Vaginal temperature

VIT:

Vaginal implant transmitter

References

  1. Stenseth NC, Mysterud A, Ottersen G, Hurrell JW, Chan K-S, Lima M. Ecological effects of climate fluctuations. Science. 2002;297(5585):1292.

    CAS  PubMed  Google Scholar 

  2. Walther G-R, Post E, Convey P, Menzel A, Parmesan C, Beebee TJC, et al. Ecological responses to recent climate change. Nature. 2002;416(6879):389–95.

    CAS  PubMed  Google Scholar 

  3. Malmsten J. Reproduction and health of moose in southern Sweden [Doctoral thesis]. Uppsala: Swedish University of Agricultural Sciences; 2014.

  4. Monteith KL, Klaver RW, Hersey KR, Holland AA, Thomas TP, Kauffman MJ. Effects of climate and plant phenology on recruitment of moose at the southern extent of their range. Oecologia. 2015;178(4):1137–48.

    PubMed  Google Scholar 

  5. Wattles D, De Stefano S. Status and management of moose in the northeastern United States. Alces J Devoted Biol Manag Moose. 2011;47:53–68.

    Google Scholar 

  6. Darimont CT, Paquet PC, Reimchen TE, Crichton V. Range expansion by moose into coastal temperate rainforests of British Columbia, Canada. Divers Distrib. 2005;11(3):235–9.

    Google Scholar 

  7. Murray D, Hussey K, Finnegan L, Lowe S, Price G, Benson J, et al. Assessment of the status and viability of a population of moose (Alces alces) at its southern range limit in Ontario. Can J Zool. 2012;90:422–34.

    Google Scholar 

  8. Harris RB, Atamian M, Ferguson H, Keren I. Estimating moose abundance and trends in northeastern Washington State: index counts, sightability models, and reducing uncertainty. Alces J Devoted Biol Manag Moose. 2015;51:57–69.

    Google Scholar 

  9. Cederlund GN, Håkan KGS, Pehrson Å. Body mass dynamics of moose calves in relation to winter severity. J Wildl Manag. 1991;55(4):675–81.

    Google Scholar 

  10. Malmsten J, Dalin A-M. Reproductive failure in moose (Alces alces) due to embryonic mortality and unfertilized oocytes. Acta Theriol. 2014;59(3):449–55.

    Google Scholar 

  11. Sand H, Cederlund G. Individual and geographical variation in age at maturity in female moose (Alces alces). Can J Zool. 1996;74(5):954–64.

    Google Scholar 

  12. Bowyer R, Ballenberghe V, Kie J. Timing and synchrony of parturition in Alaskan moose: long-term versus proximal effects of climate. J Mammal. 1998;79:1332–44.

    Google Scholar 

  13. Cook RC, Oyster J, Mansfield K, Harris RB. Evidence of summer nutritional limitations in a Northeastern Washington moose population. Alces J Devoted Biol Manag Moose. 2021;57:23–46.

    Google Scholar 

  14. Keech MA, Lindberg MS, Boertje RD, Valkenburg P, Taras BD, Boudreau TA, et al. Effects of predator treatments, individual traits, and environment on moose survival in Alaska. J Wildl Manag. 2011;75(6):1361–80.

    Google Scholar 

  15. Schwartz CC, Hundertmark KJ, Becker EF. Growth of moose calves conceived during the first versus second estrus. Alces J Devoted Biol Manag Moose. 1994;30:91–100.

    Google Scholar 

  16. Schwartz MD. Phenology: an integrative environmental science. 2nd ed. Dordrecht: Springer; 2013.

    Google Scholar 

  17. Schwartz CC, Hundertmark KJ. Reproductive characteristics of Alaskan moose. J Wildl Manag. 1993;57(3):454–68.

    Google Scholar 

  18. Malmsten J, Söderquist L, Thulin C-G, Gavier Widén D, Yon L, Hutchings MR, et al. Reproductive characteristics in female Swedish moose (Alces alces), with emphasis on puberty, timing of oestrus, and mating. Acta Vet Scand. 2014;56(1):23.

    PubMed  PubMed Central  Google Scholar 

  19. Ballenberghe VV, Miquelle DG. Mating in moose: timing, behavior, and male access patterns. Can J Zool. 1993;71(8):1687–90.

    Google Scholar 

  20. Shinomiya A, Shimmura T, Nishiwaki-Ohkawa T, Yoshimura T. Regulation of seasonal reproduction by hypothalamic activation of thyroid hormone. Front Endocrinol. 2014. https://doi.org/10.3389/fendo.2014.00012.

    Article  Google Scholar 

  21. Garel M, Solberg EJ, Sæther B-E, Grøtan V, Tufto J, Heim M. Age, size, and spatiotemporal variation in ovulation patterns of a seasonal breeder, the Norwegian moose (Alces alces). Am Nat. 2009;173(1):89–104.

    PubMed  Google Scholar 

  22. Sand H, Cederlund G, Danell K. Geographical and latitudinal variation in growth patterns and adult body size of Swedish moose (Alces alces). Oecologia. 1995;102(4):433–42.

    PubMed  Google Scholar 

  23. Sand H. Life history patterns in female moose (Alces alces): the relationship between age, body size, fecundity and environmental conditions. Oecologia. 1996;106(2):212–20.

    PubMed  Google Scholar 

  24. Sigouin D, Ouellet JP, Courtois R. Geographical variation in the mating and calving periods of moose. Alces J Devoted Biol Manag Moose. 1997;33:85–95.

    Google Scholar 

  25. Schwartz CCC, Monfort SL, Dennis PH, Hundertmark KJ. Fecal progestagen concentration as an indicator of the estrous cycle and pregnancy in moose. J Wildl Manag. 1995;59(3):580–3.

    Google Scholar 

  26. Roelofs J, López-Gatius F, Hunter RHF, van Eerdenburg FJCM, Hanzen C. When is a cow in estrus? Clinical and practical aspects. Theriogenology. 2010;74(3):327–44.

    CAS  PubMed  Google Scholar 

  27. Sveberg G, Refsdal AO, Erhard HW, Kommisrud E, Aldrin M, Tvete IF, et al. Behavior of lactating Holstein-Friesian cows during spontaneous cycles of estrus. J Dairy Sci. 2011;94(3):1289–301.

    CAS  PubMed  Google Scholar 

  28. Sontakke SD. Monitoring and controlling ovarian activities in wild ungulates. Theriogenology. 2017;109:31–401.

    PubMed  Google Scholar 

  29. Lasley BL, Kirkpatrick JF. Monitoring ovarian function in captive and free-ranging wildlife by means of urinary and fecal steroids. J Zoo Wildl Med. 1991;22(1):23–31.

    Google Scholar 

  30. Schiffer L, Barnard L, Baranowski ES, Gilligan LC, Taylor AE, Arlt W, et al. Human steroid biosynthesis, metabolism and excretion are differentially reflected by serum and urine steroid metabolomes: a comprehensive review. J Steroid Biochem Mol Biol. 2019;194:105439.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Kirkpatrick JF, Bancroft K, Kincy V. Pregnancy and ovulation detection in bison (Bison bison) assessed by means of urinary and fecal steroids. J Wildl Dis. 1992;28(4):590–7.

    CAS  PubMed  Google Scholar 

  32. Pereira RJG, Polegato BF, de Souza S, Negrão JA, Duarte JMB. Monitoring ovarian cycles and pregnancy in brown brocket deer (Mazama gouazoubira) by measurement of fecal progesterone metabolites. Theriogenology. 2006;65(2):387–99.

    CAS  PubMed  Google Scholar 

  33. Matsuura Y, Sasamoto Y, Sato K, Takahashi Y, Suzuki M, Ohtaishi N. Monitoring ovarian cycle and conception by fecal progesterone analysis in sika deer. Ecol Res. 2004;19(4):397–404.

    CAS  Google Scholar 

  34. Polegato BF, Zanetti EDS, Duarte JMB. Monitoring ovarian cycles, pregnancy and post-partum in captive marsh deer (Blastocerus dichotomus) by measuring fecal steroids. Conserv Physiol. 2018;6(1):cox073.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Dumonceaux GA, Bauman JE, Camilo GR. Evaluation of progesterone levels in feces of captive reticulated giraffe (Giraffa camelopardalis reticulata). J Zoo Wildl Med. 2006;37(3):255–61.

    PubMed  Google Scholar 

  36. Græsli AR, Thiel A, Fuchs B, Singh NJ, Stenbacka F, Ericsson G, et al. Seasonal hypometabolism in female moose. Front Ecol Evol. 2020. https://doi.org/10.3389/fevo.2020.00107.

    Article  Google Scholar 

  37. Thompson DP, Barboza PS, Crouse JA, McDonough TJ, Badajos OH, Herberg AM. Body temperature patterns vary with day, season, and body condition of moose (Alces alces). J Mammal. 2019;100(5):1466–78.

    Google Scholar 

  38. Græsli AR, Thiel A, Fuchs B, Stenbacka F, Neumann W, Malmsten J, et al. Body temperature patterns during pregnancy and parturition in moose. J Therm Biol. 2022;109:103334.

    PubMed  Google Scholar 

  39. Nicholson KL, Warren MJ, Rostan C, Månsson J, Paragi TF, Sand H. Using fine-scale movement patterns to infer ungulate parturition. Ecol Ind. 2019;101:22–30.

    Google Scholar 

  40. Wrenn TR, Bitman J, Sykes JF. Body temperature variations in dairy cattle during the estrous cycle and pregnancy. J Dairy Sci. 1958;41(8):1071–6.

    CAS  Google Scholar 

  41. Huzzey JM, von Keyserlingk MAG, Weary DM. Changes in feeding, drinking, and standing behavior of dairy cows during the transition period. J Dairy Sci. 2005;88(7):2454–61.

    CAS  PubMed  Google Scholar 

  42. Suthar VS, Burfeind O, Bonk S, Dhami AJ, Heuwieser W. Endogenous and exogenous progesterone influence body temperature in dairy cows. J Dairy Sci. 2012;95(5):2381–9.

    CAS  PubMed  Google Scholar 

  43. Suthar VS, Burfeind O, Patel JS, Dhami AJ, Heuwieser W. Body temperature around induced estrus in dairy cows. J Dairy Sci. 2011;94(5):2368–73.

    CAS  PubMed  Google Scholar 

  44. Abrams RM, Thatcher WW, Bazer FW, Wilcox CJ. Effect of estradiol-17β on vaginal thermal conductance in cattle. J Dairy Sci. 1973;56(8):1058–62.

    CAS  PubMed  Google Scholar 

  45. Brehme U, Stollberg U, Holz R, Schleusener T. ALT pedometer – a new sensor-aided measurement system for improvement in oestrus detection. Res Agric Eng. 2018;52:1–10.

    Google Scholar 

  46. Løvendahl P, Chagunda MGG. On the use of physical activity monitoring for estrus detection in dairy cows. J Dairy Sci. 2010;93(1):249–59.

    PubMed  Google Scholar 

  47. Redden KD, Kennedy AD, Ingalls JR, Gilson TL. Detection of estrus by radiotelemetric monitoring of vaginal and ear skin temperature and pedometer measurements of activity. J Dairy Sci. 1993;76(3):713–21.

    CAS  PubMed  Google Scholar 

  48. Shahriar MS, Smith D, Rahman A, Freeman M, Hills J, Rawnsley R, et al. Detecting heat events in dairy cows using accelerometers and unsupervised learning. Comput Electron Agric. 2016;128:20–6.

    Google Scholar 

  49. Ozoga JJ, Verme LJ. Activity patterns of white-tailed deer during estrus. J Wildl Manag. 1975;39(4):679–83.

    Google Scholar 

  50. Lewis GS, Newman SK. Changes throughout estrous cycles of variables that might indicate estrus in dairy cows. J Dairy Sci. 1984;67(1):146–52.

    CAS  PubMed  Google Scholar 

  51. Patterson BR, Benson JF, Middel KR, Mills KJ, Silver A, Obbard ME. Moose calf mortality in central Ontario, Canada. J Wildl Manag. 2013;77(4):832–41.

    Google Scholar 

  52. Schmidt NM, Grøndahl C, Evans AL, Desforges J-P, Blake J, Hansen LH, et al. On the interplay between hypothermia and reproduction in a high arctic ungulate. Sci Rep. 2020;10(1):1514.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Tanaka Y, Sandoval EDP, Duarte JMB. Non-homogeneous distribution of steroids in fecal pellets: An example in brown brocket deer (Mazama gouazoubira) with progesterone metabolites. Gen Comp Endocrinol. 2019;282:113206.

    CAS  PubMed  Google Scholar 

  54. Daly-Crews K, Edell RH, Metrione LC. Fecal progesterone and estrogen metabolite monitoring for cyclicity and pregnancy in southern tamandua (Tamandua tetradactyla). Zoo Biol. 2020;39(4):239–45.

    CAS  PubMed  Google Scholar 

  55. Diamond HJ, Karl TR, Palecki MA, Baker CB, Bell JE, Leeper RD, et al. U.S. Climate reference network after one decade of operations: status and assessment. Bull Am Meteorol Soc. 2013;94(4):485–98.

    Google Scholar 

  56. Alduchov OA, Eskridge RE. Improved magnus form approximation of saturation vapor pressure. J Appl Meteorol Climatol. 1996;35(4):601–9.

    Google Scholar 

  57. Zeileis A, Grothendieck G. zoo: S3 Infrastructure for regular and irregular time series. J Stat Softw. 2005;14(6):1–27.

    Google Scholar 

  58. Martin AD, Lystad ML, Reksen O, Ropstad E, Waldmann A, Nafstad O, et al. Assessment of progesterone profiles and postpartum onset of luteal activity in spring calving Hereford beef suckler cattle. Acta Vet Scand. 2010;52(1):42.

    PubMed  PubMed Central  Google Scholar 

  59. Fatet A, Pellicer-Rubio M-T, Leboeuf B. Reproductive cycle of goats. Anim Reprod Sci. 2011;124(3):211–9.

    CAS  PubMed  Google Scholar 

  60. Brown HM, Fabre Nys C, Cognié J, Scaramuzzi RJ. Short oestrous cycles in sheep during anoestrus involve defects in progesterone biosynthesis and luteal neovascularisation. Reproduction. 2014;147(3):357–67.

    CAS  PubMed  Google Scholar 

  61. Wood SN. Generalized additive models: an introduction with R. 2nd ed. Boca Raton: Chapman and Hall/CRC; 2017. p. 496.

    Google Scholar 

  62. Schloerke B, Crowley J, Cook D. Package "GGally". Extension to "gpplot2". 2018;713.

  63. Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer-Verlag New York; 2016.

    Google Scholar 

  64. Mukaka MM. Statistics corner: a guide to appropriate use of correlation coefficient in medical research. Malawi Med J. 2012;24(3):69–71.

    CAS  PubMed  PubMed Central  Google Scholar 

  65. Bolker B. Package "bbmle". Tools for General Maximum Likelihood Estimation 2017; 641.

  66. Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. 2nd ed. New York, USA: Springer-Verlag; 2002.

    Google Scholar 

  67. Thompson DP, Crouse JA, Barboza PS, Spathelf MO, Herberg AM, Parker SD, et al. Behaviour influences thermoregulation of boreal moose during the warmseason. Conserv Physiol. 2021;9(1):coaa130.

    PubMed  PubMed Central  Google Scholar 

  68. Roelofs JB, Van Eerdenburg FJCM, Hazeleger W, Soede NM, Kemp B. Relationship between progesterone concentrations in milk and blood and time of ovulation in dairy cattle. Anim Reprod Sci. 2006;91(3):337–43.

    CAS  PubMed  Google Scholar 

  69. Bartlewski PM, Beard AP, Rawlings NC. Ovarian function in ewes at the onset of the breeding season. Anim Reprod Sci. 1999;57(1):67–88.

    CAS  PubMed  Google Scholar 

  70. Hanzen CH, Pieterse M, Scenczi O, Drost M. Relative accuracy of the identification of ovarian structures in the cow by ultrasonography and palpation per rectum. Vet J. 2000;159(2):161–70.

    CAS  PubMed  Google Scholar 

  71. Alhamada M, Debus N, Lurette A, Bocquier F. Validation of automated electronic oestrus detection in sheep as an alternative to visual observation. Small Rumin Res. 2016;134:97–104.

    Google Scholar 

  72. Säkkinen H, Tornbeg J, Goddard P, Eloranta E, Ropstad E, Seppo S. The effect of blood sampling method on indicators of physiological stress in reindeer (Rangifer tarandus tarandus). Domest Anim Endocrinol. 2004;26:87–98.

    PubMed  Google Scholar 

  73. Voigt CC, Faßbender M, Dehnhard M, Wibbelt G, Jewgenow K, Hofer H, et al. Validation of a minimally invasive blood-sampling technique for the analysis of hormones in domestic rabbits, Oryctolagus cuniculus (Lagomorpha). Gen Comp Endocrinol. 2004;135(1):100–7.

    CAS  PubMed  Google Scholar 

  74. Ingram JR, Matthews LR, McDonald RM. Remote blood sampling device a stress free blood sampling technique for free ranging animals. Proc N Z Soc Anim Prod. 1994;54:39–42.

    Google Scholar 

  75. Fønss A, Munksgaard L. Automatic blood sampling in dairy cows. Comput Electron Agric. 2008;64(1):27–33.

    Google Scholar 

  76. Ravindra JP, Rawlings NC. Ovarian follicular dynamics in ewes during the transition from anoestrus to the breeding season. Reproduction. 1997;110(2):279–89.

    CAS  Google Scholar 

  77. Piccione G, Caola G, Refinetti R. Daily and estrous rhythmicity of body temperature in domestic cattle. BMC Physiol. 2003;3:7.

    PubMed  PubMed Central  Google Scholar 

  78. Gonzalez A, Murphy BD, de Alba MJ, Manns JG. Endocrinology of the postpartum period in the pelibuey ewe. J Anim Sci. 1987;64(6):1717–24.

    CAS  PubMed  Google Scholar 

  79. Camp JC, Wildt DE, Howard PK, Stuart LD, Chakraborty PK. Ovarian activity during normal and abnormal length estrous cycles in the goat. Biol Reprod. 1983;28(3):673–81.

    CAS  PubMed  Google Scholar 

  80. I’Anson H, Legan SJ. Does the first LH surge of the breeding season initiate the first full-length cycle in the ewe? Reproduction. 1988;82(2):761–7.

    CAS  Google Scholar 

  81. Alvarenga FAP, Borges I, Palkovič L, Rodina J, Oddy VH, Dobos RC. Using a three-axis accelerometer to identify and classify sheep behaviour at pasture. Appl Anim Behav Sci. 2016;181:91–9.

    Google Scholar 

  82. Vázquez Diosdado JA, Barker ZE, Hodges HR, Amory JR, Croft DP, Bell NJ, et al. Classification of behaviour in housed dairy cows using an accelerometer-based activity monitoring system. Animal Biotelemetry. 2015;3(1):15.

    Google Scholar 

  83. Heurich M, Traube M, Stache A, Löttker P. Calibration of remotely collected acceleration data with behavioral observations of roe deer (Capreolus capreolus L.). Acta Theriol. 2012;57(3):251–5.

    Google Scholar 

  84. Löttker P, Rummel A, Traube M, Stache A, Šustr P, Müller J, et al. New possibilities of observing animal behaviour from a distance using activity sensors in GPS-collars: an attempt to calibrate remotely collected activity data with direct behavioural observations in red deer Cervus elaphus. Wildl Biol. 2009;15(4):425–34.

    Google Scholar 

  85. de Mouzon J, Testart J, Lefevre B, Pouly JL, Frydman R. Time relationships between basal body temperature and ovulation or plasma progestins. Fertil Steril. 1984;41(2):254–9.

    PubMed  Google Scholar 

  86. Schwarzenberger F, Tomášová K, Holečková D, Matern B, Möstl E. Measurement of fecal steroids in the black rhinoceros (Diceros bicornis) using group-specific enzyme immunoassays for 20-oxo-pregnanes. Zoo Biol. 1996;15(2):159–71.

    CAS  Google Scholar 

  87. Shipka MP. A note on silent ovulation identified by using radiotelemetry for estrous detection. Appl Anim Behav Sci. 2000;66(1):153–9.

    Google Scholar 

  88. Walton JS, King GJ. Indicators of estrus in Holstein cows housed in tie stalls. J Dairy Sci. 1986;69(11):2966–73.

    CAS  PubMed  Google Scholar 

  89. Sakatani M, Takahashi M, Takenouchi N. The efficiency of vaginal temperature measurement for detection of estrus in Japanese Black cows. J Reprod Dev. 2016;62(2):201–7.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Additional thanks to Jason Rupp, Sarah Newberry, Theresa M. Kirchner, Jeffrey B. Stetz, Kristin Denryter and Nick Fowler for assisting with sample collection and preparation.

Funding

This work was supported by the Alaska Department of Fish and Game Federal Wildlife Restoration Grant (Grant Number AKW-4 Project No. 1.63), the Ecophysiology and One Health Research Group at Inland Norway University of Applied Sciences, and the Faculty of Veterinary Medicine at the Norwegian University of Life Sciences.

Author information

Authors and Affiliations

Authors

Contributions

AE, DT, JC, ER and JH-P contributed to the study design. DT and JC immobilized the study animals, in addition to deploying and retrieving biologging devices. JH-P, DT and JC performed data collection, with data management and statistical analysis done by JH-P, LL, and AT. JH-P prepared the figures and drafted the manuscript. All authors participated in reviewing and approving the submitted version.

Corresponding author

Correspondence to Jennifer Høy-Petersen.

Ethics declarations

Ethics approval and consent to participate

All animal care, handling and experimentation procedures at the Kenai Moose Research Center were approved by the Animal Care and Use Committee, Alaska Department of Fish and Game, Division of Wildlife Conservation (protocol number 0086-2020-40).

Consent for publication

Not applicable.

Competing interests

The authors declare that there is no competing interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Correlations between onset of luteal activity (OLA) and other explanatory variables. Correlations between the explanatory variable onset of luteal activity (OLA) and ordinal day, daily ambient temperature (°C) and daily vapor pressure (hPa) 0.96, -0.77 and -0.71, respectively). Blue dots and line = locally estimated scatterplot smoothing line. Gray region around the blue line = 95% confidence interval.

Additional file 2: Figure S2.

Correlations between explanatory variables. Correlations between the explanatory variables daily range of ambient temperature (Ta) and solar radiation, and between wind speed and relative humidity (0.74 and -0.63, respectively). Blue line = locally estimated scatterplot smoothing line. Gray region around the blue line = 95% confidence interval.

Additional file 3: Figure S3.

Vaginal temperature, collar activity, and fecal progestagen (2). Graphs from two captive female moose that were 9 years (IDs = 44494 and 44495) containing daily values for mean vaginal temperature (blue line) and daily sum of activity (red line) recorded between July 23rd and November 22nd, 2021, and fecal progestagen concentrations recorded between August 23rd to October 15th, 2021. Dashed vertical blue line = onset of luteal activity, which is on the 25th of September for the top graph and on the 24th of September on the bottom graph. Dashed horizontal red line = individual fecal progestagen threshold.

Additional file 4: Figure S4.

Vaginal temperature, collar activity, and fecal progestagen (3). Graphs from two captive female moose that were 12 years (ID = 30683) and 13 years (ID = 30679) containing daily values for mean vaginal temperature (blue line) and daily sum of activity (red line) recorded between July 23rd and November 22nd, 2021, and fecal progestagen concentrations recorded between August 23rd to October 15th, 2021. Dashed vertical blue line = onset of luteal activity, which is on the 25th of September for the top graph and on the 23rd of September on the bottom graph. Dashed horizontal red line = individual fecal progestagen threshold.

Additional file 5: Figure S5.

Vaginal temperature, collar activity, and fecal progestagen (4). Graph from an 18-year-old captive female moose (ID = 30682) containing daily values for mean vaginal temperature (blue line) and daily sum of activity (red line) recorded between July 23rd and November 22nd, 2021, and fecal progestagen concentrations recorded between August 23rd to October 15th, 2021. Dashed vertical blue line = onset of luteal activity, which is on the 15th of September for this individual. Dashed horizontal red line = individual fecal progestagen threshold.

Additional file 6: Figure S6.

Vaginal temperature, collar activity, and fecal progestagen from excluded animals (1). Graph from a 12-year-old captive female moose (ID = 30680) containing daily values for mean vaginal temperature (blue line) and daily sum of activity (red line) recorded between July 23rd and November 22nd, 2021, and fecal progestagen concentrations recorded between August 23rd to October 15th, 2021. This individual had missing vaginal temperature data from August 12th. Dashed vertical blue line = onset of luteal activity, which is on the 18th of September for this individual. Dashed horizontal red line = individual fecal progestagen threshold.

Additional file 7: Figure S7.

Vaginal temperature, collar activity, and fecal progestagen from excluded animals (2). Graphs from two 2-year-old captive female moose (IDs = 44496 and 44497) containing daily values for mean vaginal temperature (blue line) and daily sum of activity (red line) recorded between July 23rd and November 22nd, 2021, and fecal progestagen concentrations recorded between August 23rd to October 15th, 2021. No onset of luteal activity was identified for these individuals. Dashed horizontal red line = individual fecal progestagen threshold.

Additional file 8: Figure S8.

Vaginal temperature, collar activity, and fecal progestagen from excluded animals (3). Graphs from two captive female moose that were 2 years (ID = 44498) and 19 years (ID = 30678) containing daily values for mean vaginal temperature (blue line) and daily sum of activity (red line) recorded between July 23rd and November 22nd, 2021, and fecal progestagen concentrations recorded between August 23rd to October 15th, 2021. No onset of luteal activity was identified for these individuals. Dashed horizontal red line = individual fecal progestagen threshold.

Additional file 9: Table S1.

Summary of the response and explanatory variables included in the model selection process. Table S2. Fecal progestagen classification data from female moose (n = 12) from August to October 2021. Table S3. Estimated parameters for explanatory variables in the model with vaginal temperature as a response variable. Table S4. Model combinations and ranking with vaginal temperature as the response variable. Table S5. Estimated parameters for explanatory variables in the model with activity as a response variable. Table S6. Model combinations and ranking with activity as the response variable.

Rights and permissions

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://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Høy-Petersen, J., Lemière, L., Thompson, D.P. et al. Vaginal temperature and collar activity variations in relation to fecal progestagens of captive moose. Anim Biotelemetry 11, 36 (2023). https://doi.org/10.1186/s40317-023-00345-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40317-023-00345-y

Keywords