Thermal Response of Laboratory Rats (Rattus norvegicus) during the Application of Six Methods of Euthanasia Assessed by Infrared Thermography

Simple Summary The present study aimed to assess the infrared thermal response of laboratory rats (Rattus norvegicus) during the application of six euthanasia methods (injectable, inhalational, and physical) to determine the method that prevents or diminishes the stress response. The surface temperature was assessed in four thermal windows: ocular (T°ocu), auricular (T°ear), interscapular (T°dor), and caudal (T°tai). The results showed that inhalant methods (CO2 and isoflurane) had temperature alterations that could be suggestive of a marked stress response, in contrast to other methods such as pentobarbital, decapitation, and xylazine + ketamine. In conclusion, according to the thermal response of the rats, it is suggested that CO2 and isoflurane might cause distress and this needs to be considered when selecting these techniques as the method of euthanasia for laboratory rats. Abstract Refinement is one of the principles aiming to promote welfare in research animals. The techniques used during an experimental protocol, including euthanasia selection, must prevent and minimize suffering. Although the current euthanasia methods applied to laboratory rodents are accepted, the controversial findings regarding the potential stress/distress they can cause is a field of research. The objective was to assess the thermal response of Wistar rats during various euthanasia methods using infrared thermography (IRT) to determine the method that prevents or diminishes the stress response and prolonged suffering. Pentobarbital (G1), CO2 (G2), decapitation (G3), isoflurane (G4), ketamine + xylazine (G5), and ketamine + CO2 (G6) were evaluated at five evaluation times with IRT to identify changes in the surface temperature of four anatomical regions: ocular (T°ocu), auricular (T°ear), interscapular (T°dor), and caudal (T°tai). Significant differences (p < 0.05) were found in G2 and G4, registering temperature increases from the administration of the drug to the cessation of respiratory rate and heart rate. Particularly, isoflurane showed a marked thermal response in T°ocu, T°ear, T°dor, and T°tai, suggesting that, in general, inhalant euthanasia methods induce stress in rats and that isoflurane might potentially cause distress, an effect that must be considered when deciding humane euthanasia methods in laboratory rodents.


Introduction
The use of animals is a key element for improvements in biomedical science [1,2], where rats and mice represent 87-98% of the total of species used in the scientific community [3,4].The potential pain and stress that laboratory animals might experience is highly controversial [5].Ethical animal research necessitates the selection of suitable euthanasia methods to minimize pain and distress, as proposed by the National Centre for the Replacement, Refinement, and Reduction of Animals in Research [6], not only to provide welfare but also to ensure the quality of results.These initiatives need be applied not only during the life of research animals, but also during the application of euthanasia methods with the aim of providing humane endpoints [7].
Currently, there is a debate around the euthanasia methods that are approved by the American Veterinary Medicine Association (AVMA).Injectable drugs (e.g., barbiturates and general anesthetics), inhalant agents (e.g., CO 2 and isoflurane), and physical methods (e.g., decapitation and cervical dislocation) are recognized as acceptable techniques to induce a humane death without suffering [8].However, some methods are under discussion since studies have shown potential adverse effects during their application.For example, inhalation of CO 2 is aversive for rats [9] and induces bradycardia and potential anxiety due to hypoxia before loss of consciousness [10,11].Moreover, CO 2 forms carbonic acid and induces the activation of pain receptors [12].Nonetheless, systematic reviews have contrasting results regarding the suitability of CO 2 and its potential distress [13].Isoflurane is considered an alternative to CO 2 euthanasia.However, it is known as an aversive agent to rats for its mild pungency [14,15].
On the other hand, the administration of injectable pentobarbital has been associated with pain-related behaviors (e.g., writhing and back arching) due to intraperitoneal (IP) irritation [16][17][18], while the combination of xylazine and ketamine, although a common anesthetic protocol, demonstrates limited action as a euthanasic agent.However, high Mouse Grimace Scale scores and anxiety-related behaviors were found after repeated doses of the combination [19].Moreover, Wellington et al. [20] found that IP administration of ketamine + xylazine to rats caused acute muscle and tissue necrosis, poor tolerance, and pain/discomfort behavioral reactions.In the case of decapitation, this procedure leads to the question of whether brain activity is present immediately after the procedure or not, as well as whether changes in the electroencephalogram (EEG) are associated with nociception during the first 15 s (s) following decapitation [21], as determined in rats by Derr [22] who reported that EEGs during the first 2.7 s after decapitation might indicate conscious awareness of pain and distress.
The refinement of procedures performed in research animals includes the implementation of non-invasive tools to assess their welfare without causing additional stress.Infrared thermography (IRT) is a technique that detects surface temperature changes as a neuroendocrine response of the Sympathetic Nervous System (SNS) after stressful/distressful and painful events [23,24].Stress-known as the reaction of the organism when its homeostasis or psychological well-being is perturbated-and distress-a negative and aversive state when the organism cannot adapt or return to homeostasis [25]-activates two main systems: the hypothalamic-pituitary-adrenal (HPA) and the locus coeruleus sympathetic adrenomedullary (SAM) axes [26].Both axes lead to the release of glucocorticoids and catecholamines, as well as the physiological changes required to adjust homeostasis [27], including alterations in body temperature and microcirculation.Therefore, temperature variations have been used as a stress-related marker in animals, as stress may cause central hyperthermia and peripheral reduction of the temperature due to vasoconstriction [28].
IRT detects these vasomotor changes as a difference in the amount of dissipated heat in different anatomical regions, where heat exchange is facilitated through the arteriovenous anastomosis and peripheral blood vessels, also called thermal windows [29].In laboratory rodents, thermal windows such as the ocular, auricular, dorsal or interscapular, and tail region have been used to assess stress [28,30] or pain [31].For example, Lecorps et al. [32] found that eye temperature increased in mice undergoing an elevated plus maze test, while tail temperature diminished as a physiological response to stress (a result that was associated with anxiety-related behaviors).The ocular surface has great vascular sensitivity because the two main arteries responsible for its irrigation (the arteria supraorbitalis and angularis occuli), as well as the innervation through the facial nerve, rapidly respond to autonomous tone changes and endogenous catecholamines [29].Likewise, Zevgolis et al. [33] reported that ocular IRT increased during the experimental handling of wild mice due to stress-induced hyperthermia (SIH).A similar response is observed in the auricular window, as shown by Wokke [34] in mice.In this study, restraining methods to administer IP drugs increased ear temperature and corticosterone levels.These temperature variations are mediated by sympathetic activity and vasodilation in the main blood vessels supplying irrigation (external jugular vein, external carotid artery, and its branches into marginal ear arteries) [29,35].Hutu et al. [36] determined that IRT measured in the ear is positively correlated with rectal temperature in rabbits.Therefore, considering that SIH also causes changes in the amount of dissipated heat in thermal windows, ear temperature could be a way to assess acute stress.
For laboratory rodents, a thermal window that is closely related to sympathetic activation and norepinephrine (NE) release after the activation of the SAM axis is the interscapular region.In this zone, small mammals have large deposits of brown adipose tissue (BAT), a thermogenic structure whose activity depends on NE binding to β 3 -adrenoreceptors located in BAT [37].The increased thermogenic activity of this tissue has been associated with corticosterone secretion and with the administration of β 3 -adrenoreceptor agonists and NE [38,39].Furthermore, SIH is also related to BAT thermogenesis in rats and humans after excessive stress [40].
Lastly, the tail of rats is considered an important thermal window because it contributes to up to 25% of heat dissipation (by vasoconstriction) due to arteriovenous anastomosis (from the coccygeal artery) [29,41].Vasoconstriction of peripheral regions such as the tail and paws is mediated by the sympathetic redistribution of blood flow to key organs (e.g., the heart and brain).In mice exposed to acute stressors, the superficial tail temperature decreased during different handling procedures, while the surface temperature of the body (assessed in the dorsal region of the mice) increased as a response to SIH [42].Gjendal et al. [30] also reported a decrease in tail temperature (up to 3.5 • C) in mice exposed to three stressors (a maze test, IP injection, and isoflurane anesthesia) as a result of vasoconstriction of tail blood vessels.
The current literature suggests IRT as a non-invasive complementary tool to assess well-being in animals, including stress-related responses [24,43].Although there are some studies regarding IRT and the pre-slaughter or antemortem period in domestic species such as pigs [44,45], there is no study up to now where IRT has been used to evaluate the effect of euthanasia methods on laboratory species.Therefore, the present study aimed to assess the infrared thermal response of laboratory rats (Rattus norvegicus) during the application of six euthanasia methods to determine the method that prevents or diminishes the stress response and prolonged suffering.Injectable drugs (pentobarbital, ketamine + xylazine), inhalant agents (CO 2 , isoflurane), physical methods (decapitation), and the combination of inhaled and injectable anesthetics (CO 2 + ketamine) were evaluated with IRT to identify changes in the surface temperature of four anatomical regions (ocular, auricular, interscapular, and caudal).Also, differences by sex according to the thermal window and the euthanasia method will be studied.The hypotheses of the study were as follows: (i) the use of IRT during different euthanasia methods will help to recognize changes in the surface temperature of laboratory rats (R. norvegicus) in response to stress perception related to the method and (ii) the combination of an injectable anesthetic overdose (ketamine) with CO 2 exposure as the euthanasia method will reduce thermal alterations associated with stress.

Location and Ethical Statement
The present study was conducted at the Animal Facility and Experimental Surgery Facility from the Biotechnological Research Sub-Department of the Instituto Nacional de Rehabilitación Luis Guillermo Ibarra Ibarra, Mexico City, Mexico.All procedures were approved by the Committee for the Care and Use of Laboratory Animals (INRLGII/CICUAL/014/2021) at the National Institute of Rehabilitation Luis Guillermo Ibarra-Ibarra.
The handling and care of the laboratory animals was in accordance with the Mexican norm for laboratory animals NOM-062-ZOO-1999, published by the Department of Agriculture, Rural Development, Fisheries and Alimentation [46].All dead animals were disposed of by incineration following NOM-062-ZOO-1999.

Animals and Housing Conditions
A total of 60 adult Wistar rats (R. norvegicus), 30 male and 30 female, were purchased from the Center for Research and Advanced Studies at the National Polytechnic Institute (CINVESTAV-IPN).The animals had an average weight of 311 ± 62 g at 8-10 weeks old (in puberty) and were obtained with an animal health certificate to ensure they were free of infectious pathogens (bacteria, viruses, and parasites).The sample size was calculated using G*Power 3.1.9.7 (Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany).The total sample size was 48 animals, considering an α error probability of 0.05, confidence level of 95%, power (1-β error probability) of 0.90, and correction among repeated measures of 0.5 for six experimental groups with five measurements.
According to the principles of the 3Rs [47], reduction was applied to the total number of animals used by reusing animals from finished protocols related to behavioral tests (e.g., balance beam or maze tests).Rats that were part of the control groups were selected to avoid the inclusion of animals undergoing invasive procedures or those with residual drug levels.Through a general physical examination, the animals were classified as healthy without signs of disease, stress, or pain-related behaviors.The physical exam considered body weight, posture, level of consciousness, secretions, the color of the mucosa, sneezing, and a species-specific behavioral repertoire.Rats showing signs of disease, injury, or pain were eliminated.Pregnant females were excluded.
Rats were housed in separate rooms according to sex.They were placed in groups of five animals per cage in standard polycarbonate cages for rats (47 × 36 × 21 cm) with wood shavings as bedding (Aspen, Nepco, Riverside, RI, USA) and without enrichment.The rats were maintained under a 12 h day-night cycle with lights on between 0500 h and 1700 h.The controlled temperature inside the housing room and the testing room was set at an average of 23.2 ± 0.5 • C and 22.9 ± 0.5 • C, respectively, with respective humidity levels of 48% and 52%.The rats had ad libitum access to food (LabDiet 5010, LabDiet, St. Louis, MI, USA) and purified water (in 500 mL drinking water bottles), and the cages were cleaned once a week.Visual health inspection was performed twice daily.

Experimental Design
This was an experimental prospective-comparative study.All measurements were performed by a single trained and unblinded evaluator.Once the rats were selected for the study, they underwent habituation for 15 days to the customized euthanasia chamber, handling, and the evaluator's presence.The animals were randomly divided into six groups by number generation (Microsoft Excel; Microsoft 365).A total of 10 rats were assigned in each group (5 males and 5 females) as follows: G 1 : Pentobarbital (Pentobarbital, Aranda ® , Mexico City, Mexico) overdose at 400 mg/kg performed via IP injection with a 3 mL sterile syringe (Ambiderm ® , Baja California, Mexico) following Lofgren et al. [48]'s procedure.The dose was calculated through a pilot study (no published data) using the minimal and maximal doses that appear in Reimer et al. [17]'s study.The selected dose resulted in rapid loss of consciousness, cardiorespiratory depression without excitation, loss of reflexes, and clinical death.G 2 : CO 2 overdose administered inside a customized acrylic euthanasia chamber (Acrifactory, Mexico City, Mexico) (32.5 × 42 × 21 cm).The chamber had five gates with hermetically sealed doors that used neodymium magnets to avoid gas leaks.Each door was fitted to the size of the thermal camera lens to allow for thermal imaging (Figure 1).According to the AVMA [8], the flow rate was set at 30% of the chamber volume/min.G 3 : Decapitation using a rodent guillotine (51330, Senna, Mexico) [48].G 4 : Inhalation of isoflurane (Fluriso, VET ONE ® , Delhi, India) using the open-drop exposure method (two cotton swabs soaked with 2 mL of isoflurane each).The dose was calculated using the studies of Risling et al. [49] and de Brito [50] as a basis.The cotton swabs were placed where animals could not have direct contact with the inhalant anesthetic drug.G 5 : Ketamine (Ketamin-Pet, Aranda ® , Mexico) + xylazine (Procin, Pisa Agropecuaria ® , Nuevo México, Mexico) overdose administered at doses of 450 mg/kg IP and 45 mg/kg IP, respectively [8].G 6 : Combination of ketamine (100 mg/kg IP) + CO 2 (after 5-10 min of ketamine administration) [51].
Animals 2023, 13, x FOR PEER REVIEW 5 of 21 al. [17]'s study.The selected dose resulted in rapid loss of consciousness, cardiorespiratory depression without excitation, loss of reflexes, and clinical death.G2: CO2 overdose administered inside a customized acrylic euthanasia chamber (Acrifactory, Mexico City, Mexico) (32.5 × 42 × 21 cm).The chamber had five gates with hermetically sealed doors that used neodymium magnets to avoid gas leaks.Each door was fitted to the size of the thermal camera lens to allow for thermal imaging (Figure 1).According to the AVMA [8], the flow rate was set at 30% of the chamber volume/min.G3: Decapitation using a rodent guillotine (51330, Senna, Mexico) [48].G4: Inhalation of isoflurane (Fluriso, VET ONE ® , Delhi, India) using the open-drop exposure method (two cotton swabs soaked with 2 mL of isoflurane each).The dose was calculated using the studies of Risling et al.
[49] and de Brito [50] as a basis.The cotton swabs were placed where animals could not have direct contact with the inhalant anesthetic drug.G5: Ketamine (Ketamin-Pet, Aranda ® , Mexico) + xylazine (Procin, Pisa Agropecuaria ® , Nuevo México, Mexico) overdose administered at doses of 450 mg/kg IP and 45 mg/kg IP, respectively [8].G6: Combination of ketamine (100 mg/kg IP) + CO2 (after 5-10 min of ketamine administration) [51].Rats from all groups were assessed at five evaluation times.The basal time point represents assessment that was performed 24 h before the euthanasia method inside the housing room, and Ti1 represents three minutes before the application of the euthanasia in the test room.On the trial day, the rats were moved from the housing room to the test room, allowing for 30 min of rest and room acclimatization before starting the trial at Ti2, the time during the application of the method (e.g., while the animal received the IP dose of pentobarbital, ketamine + xylazine, or while it was inside the induction chamber or placed in the guillotine).Ti3 represents the time immediately after the application of the euthanasia method until loss of the righting reflex (LORR) as a sign of unconsciousness, and Ti4 represents the time until the cessation of breathing and heartbeat by visual assessment and assessment using a stethoscope (3M™ Littmann ® Classic III™, 3M, Saint Paul, MN, USA).The absence of palpebral, interdigital, and righting reflex was also used to confirm the euthanasia method.It is noteworthy to mention that all groups, except G3 Rats from all groups were assessed at five evaluation times.The basal time point represents assessment that was performed 24 h before the euthanasia method inside the housing room, and Ti 1 represents three minutes before the application of the euthanasia in the test room.On the trial day, the rats were moved from the housing room to the test room, allowing for 30 min of rest and room acclimatization before starting the trial at Ti 2 , the time during the application of the method (e.g., while the animal received the IP dose of pentobarbital, ketamine + xylazine, or while it was inside the induction chamber or placed in the guillotine).Ti 3 represents the time immediately after the application of the euthanasia method until loss of the righting reflex (LORR) as a sign of unconsciousness, and Ti 4 represents the time until the cessation of breathing and heartbeat by visual assessment and assessment using a stethoscope (3M™ Littmann ® Classic III™, 3M, Saint Paul, MN, USA).The absence of palpebral, interdigital, and righting reflex was also used to confirm the euthanasia method.It is noteworthy to mention that all groups, except G 3 (decapitation), had the same five evaluation times.In G 3 , the separation of the head from the body and unconsciousness was considered the same event; therefore, only four evaluation times were considered for this group (basal, Ti 1 , Ti 2 , and Ti 3 ).

Assessed Parameters 2.4.1. Infrared Thermography (IRT)
Thermal imaging was performed using an FLIR™ E60 camera (FLIR Systems, Orlando, FL, USA) positioned 1 m from the rats while maintaining a perpendicular angle to the subject (90 • ).Radiometric images were taken with an emissivity of 0.95, an IR resolution of 320 × 240 pixels, thermal sensitivity of <0.05 • C, and accuracy of ±2%.To prevent reflective heat affecting the acrylic cages placed in the housing room and the test room, the walls were covered with kraft paper.The handler used latex gloves when restraining the rats for IP injection and decapitation.Moreover, thermal imaging was performed at the same time of the day in all experimental groups (between 0800 h and 1500 h).
During basal, Ti 1 , Ti 2 , Ti 3 , and Ti 4 , the four evaluated body regions (or thermal windows) were the ocular (T • ocu ), auricular (T • ear ), interscapular (T • dor ), and tail (T • tai ) regions.Thermal imaging for T • ocu and T • ear was taken from the right side of the animals.The delimitation of these regions of interest (ROIs) is shown in Figure 2. The thermal images were processed with FLIR Tools software (FLIR Systems, USA) to obtain the maximum, minimum, and average temperatures for T • ocu , T • ear , and T • dor .For T • tai , only the average temperature in the proximal, medial, and distal parts of the tail was recorded.This is due to the delimitation of the ROI with a spot, which only provides the average value.
Animals 2023, 13, x FOR PEER REVIEW 6 of 21 (decapitation), had the same five evaluation times.In G3, the separation of the head from the body and unconsciousness was considered the same event; therefore, only four evaluation times were considered for this group (basal, Ti1, Ti2, and Ti3).

Infrared Thermography (IRT)
Thermal imaging was performed using an FLIR™ E60 camera (FLIR Systems, Orlando, FL, USA) positioned 1 m from the rats while maintaining a perpendicular angle to the subject (90°).Radiometric images were taken with an emissivity of 0.95, an IR resolution of 320 × 240 pixels, thermal sensitivity of <0.05 °C, and accuracy of ±2%.To prevent reflective heat affecting the acrylic cages placed in the housing room and the test room, the walls were covered with kraft paper.The handler used latex gloves when restraining the rats for IP injection and decapitation.Moreover, thermal imaging was performed at the same time of the day in all experimental groups (between 0800 h and 1500 h).
During basal, Ti1, Ti2, Ti3, and Ti4, the four evaluated body regions (or thermal windows) were the ocular (T°ocu), auricular (T°ear), interscapular (T°dor), and tail (T°tai) regions.Thermal imaging for T°ocu and T°ear was taken from the right side of the animals.The delimitation of these regions of interest (ROIs) is shown in Figure 2. The thermal images were processed with FLIR Tools software (FLIR Systems, USA) to obtain the maximum, minimum, and average temperatures for T°ocu, T°ear, and T°dor.For T°tai, only the average temperature in the proximal, medial, and distal parts of the tail was recorded.This is due to the delimitation of the ROI with a spot, which only provides the average value.T°ear was evaluated using a circle (El2) in the external ear canal to assess the irradiated temperature of the tympanic membrane and inner ear.(B).For T°dor, a rectangle (Bx1) was placed in the dorsal area over the interscapular space.(C) For T°tai, three spots (Sp1, Sp2, y Sp3) were placed at the proximal (T°prox), medial (T°medial), and distal (T°distal) segments of the tail.

Time to Death
To record the duration of each euthanasia method, after Ti2, the evaluator started a timer to register the time of death, the time to LORR, and the time to the cessation of breathing (visual assessment) and heartbeat (thoracic auscultation).The results were expressed as seconds, and the average value ± standard deviation (SD) of the 10 animals per group was recorded on an evaluation sheet.

Procedures
After the 15-day habituation period and the random assignment of the rats into the six experimental groups, the procedures were performed as shown in Figure 3.

Time to Death
To record the duration of each euthanasia method, after Ti 2 , the evaluator started a timer to register the time of death, the time to LORR, and the time to the cessation of breathing (visual assessment) and heartbeat (thoracic auscultation).The results were expressed as seconds, and the average value ± standard deviation (SD) of the 10 animals per group was recorded on an evaluation sheet.

Procedures
After the 15-day habituation period and the random assignment of the rats into the six experimental groups, the procedures were performed as shown in Figure 3.During the basal time point, IRT, as well as room temperature and relative humidity (% RH), was recorded inside the housing room of the selected experimental group with a wireless indoor and outdoor weather station with a hygrometer (Taylor ® , Oak Brook, IL, USA).The equipment and surfaces were conditioned to IRT readings 24 h after by covering the walls of the polycarbonate cages with kraft paper and the use of wood shavings as bedding.The surfaces to place the guillotine, cages, and induction chamber were also covered by either sterile drapes or cork pads to avoid reflective heat.The rats from the corresponding group were moved from the housing room to the testing room so that euthanasia was not performed where the rest of the animals were housed.A period of 30 min was given to the rats to acclimatize them to the controlled temperature in the testing room and avoid stress related to transportation.Following this 30 min, the euthanasia method started.Room temperature, % RH, IRT, and time of death were recorded for each individual in all experimental groups.

Statistical Analyses
All analyses were performed using the GraphPad Prism 10.0.0 (San Jose, Ca, USA) statistical package.The Shapiro-Wilk test was performed to establish data normality in the data set collected from T°ocu, T°ear, T°dor, and T°tai.Descriptive statistics were obtained and results were expressed as mean ± standard deviation (SD).A linear mixed model for repeated measures was used to evaluate the effect of the six euthanasia methods (treatments G1, G2, G3, G4, G5, and G6) at the five time points (basal, Ti1, Ti2, Ti3, and Ti4) for each of the four thermal windows.Multiple comparison of means was performed with the post-hoc Tukey test.In every case, the significance level was set at p < 0.05.The following statistical model was used: Yijk= µ + τi + τj+ τiτj + βk +еij where the symbols indicate the following: Y = variable response (IRT); τi = fixed effect (G1, G2, G3, G4, G5, G6); τj = evaluation times (basal, Ti1, Ti2, Ti3, and Ti4); β = aleatory effect (rat); µ = general mean; е = error.
To determine if there were differences between males and females from each experimental group, repeated measure ANOVA was performed with a Greenhouse-Geisser correction and a post-hoc Tukey test for multiple comparisons.Time of death, time to LORR, and time to the cessation of breathing (visual assessment) and heartbeat (thoracic auscultation) were expressed as mean ± SD.To establish the correlation between During the basal time point, IRT, as well as room temperature and relative humidity (% RH), was recorded inside the housing room of the selected experimental group with a wireless indoor and outdoor weather station with a hygrometer (Taylor ® , Oak Brook, IL, USA).The equipment and surfaces were conditioned to IRT readings 24 h after by covering the walls of the polycarbonate cages with kraft paper and the use of wood shavings as bedding.The surfaces to place the guillotine, cages, and induction chamber were also covered by either sterile drapes or cork pads to avoid reflective heat.The rats from the corresponding group were moved from the housing room to the testing room so that euthanasia was not performed where the rest of the animals were housed.A period of 30 min was given to the rats to acclimatize them to the controlled temperature in the testing room and avoid stress related to transportation.Following this 30 min, the euthanasia method started.Room temperature, % RH, IRT, and time of death were recorded for each individual in all experimental groups.

Statistical Analyses
All analyses were performed using the GraphPad Prism 10.0.0 (San Jose, Ca, USA) statistical package.The Shapiro-Wilk test was performed to establish data normality in the data set collected from T • ocu , T • ear , T • dor , and T • tai .Descriptive statistics were obtained and results were expressed as mean ± standard deviation (SD).A linear mixed model for repeated measures was used to evaluate the effect of the six euthanasia methods (treatments G 1 , G 2 , G 3 , G 4 , G 5 , and G 6 ) at the five time points (basal, Ti 1 , Ti 2 , Ti 3 , and Ti 4 ) for each of the four thermal windows.Multiple comparison of means was performed with the post-hoc Tukey test.In every case, the significance level was set at p < 0. To determine if there were differences between males and females from each experimental group, repeated measure ANOVA was performed with a Greenhouse-Geisser correction and a post-hoc Tukey test for multiple comparisons.Time of death, time to LORR, and time to the cessation of breathing (visual assessment) and heartbeat (thoracic auscultation) were expressed as mean ± SD.To establish the correlation between the ther-mal windows, Pearson correlation coefficients were calculated.All values with p < 0.05 were considered significant.

Results
Differences in the thermal response of the rats grouped in different experimental groups were obtained according to the thermal window, assessing their maximum (T • max), minimum (T • min), and mean temperature (T • mean).In general, G 2 and G 4 registered significant differences between evaluation times and between groups in three of the four thermal windows.Additionally, G 4 individuals showed a progressive increase in temperature in all thermal windows, in contrast to the other experimental groups.

Ocular Surface Temperature (T •
Ocu ) Table 1 shows the mean and standard deviation (SD) values for the temperature of T • ocu .For T • max, T • min, and T • mean, differences among evaluation times were recorded in G 2 (p = 0.0097) and G 4 (p = 0.0001).For G 2 , a decrease in T • mean of up to 2.15 • C at Ti 4 was observed when compared to basal values.Similarly, T • mean in G 4 decreased by up to 5.82 • C at Ti 2 .Regarding differences between groups, a progressive temperature decline was registered in all euthanasia methods.However, differences between the groups were observed during Ti 2 (p = 0.0007), Ti 3 (p = 0.0002), and Ti 4 (p = 0.0019), with the lowest a,b,c different literals indicate significant differences (p < 0.05) between events (basal, Ti 1 , Ti 2 , Ti 3 , Ti 4 ). 1,2,3,4different numerals indicate significant differences (p < 0.05) between treatments (G 1 , G 2 , G 3 , G 4 , G 5 , G 6 ).Bold p-values represent statistically significant differences between events and treatments.ND = not determined due to the experimental group.Treatments (G 1 : pentobarbital; G 2 : CO 2 ; G 3 : decapitation; G 4 : isoflurane; G 5 : ketamine + xylazine; G 6 : ketamine + CO 2 ).Evaluation times (basal: 24 h before the procedure; Ti 1 : three minutes before the procedure; Ti 2 : during the application of the euthanasia method; Ti 3 : immediately after the application until loss of righting reflex (LORR); Ti 4 : from LORR to cessation of heartbeat and breathing).
Between groups, significant differences were reported in Basal values of the three segments (p = 0.002, 0.001 and 0.005).Particularly, during Ti 2 of the T • prox segment, differences were observed in G 2 and G 4 (p = 0.05), with the lowest T Correlations between the four thermal windows and the experimental groups were obtained and are shown according to the experimental group in Supplementary Tables S1-S6.In all groups, significant (p < 0.001) and strong correlations (r > 0.96) were found between thermal windows.
Table 5 summarizes descriptive analysis of the recorded times (in seconds) for each group.Considering the total time of death, the longest duration was observed in G 6 , followed by G 4 (294.2± 74.3) and G 2 (390.2 ± 171.4).Similarly, cessation of RR and HR was longer in G 6 and G 4 .The groups that reached LORR faster were G 5 and G 2 (67 ± 10.3 and 78 ± 29.0 s, respectively).

Effect of Sex on the Thermal Response of the Rats
Figure 4 illustrates the impact of sex on the temperatures of each thermal window according to the euthanasia methods.In general, no marked effect by sex was found in the present study.Only four statistically significant differences were registered for T • ocu , T • ear , and T • dor , while no effect was found on T • tai .For T • ocu (Figure 4A).A significant difference (p = 0.04) was found in the Tmean of G 6 , where males had the highest temperatures (36.0 • C) in comparison to females (34.9 • C).However, a tendency to show a difference was observed in Tmax (p = 0.08) and Tmin (0.07).For T • ear (Figure 4B), temperatures between males and females significantly differed in terms of Tmin and Tmean for G 5 (p = 0.002 and p = 0.007, respectively), finding the highest temperatures in males rather than females (33.0 and 35.9 • C vs. 30.8and 34.2 • C, respectively).Figure 4C shows a statistical significance between sexes in terms of Tmax for G 1 (p < 0.0001) for T   °C) in comparison to females (34.9 °C).However, a tendency to show a difference was observed in Tmax (p = 0.08) and Tmin (0.07).For T°ear (Figure 4B), temperatures between males and females significantly differed in terms of Tmin and Tmean for G5 (p = 0.002 and p = 0.007, respectively), finding the highest temperatures in males rather than females (33.0 and 35.9 °C vs. 30.8and 34.2 °C, respectively).Figure 4C shows a statistical significance between sexes in terms of Tmax for G1 (p < 0.0001) for T°dor.G1 males registered a Tmax T°dor of 33.3 °C, while females recorded 32.2 °C.

Discussion
Based on the findings of this study, in contrast to injectable, physical, and combined euthanasia methods, inhalant agents (CO2 and isoflurane) resulted in substantial alterations in the T°ocu, T°ear, T°dor, and T°tai of Wistar rats.To date, there are no studies where IRT was used to evaluate different euthanasia methods or their effect on the thermal

Discussion
Based on the findings of this study, in contrast to injectable, physical, and combined euthanasia methods, inhalant agents (CO 2 and isoflurane) resulted in substantial alterations in the T • ocu , T • ear , T • dor , and T • tai of Wistar rats.To date, there are no studies where IRT was used to evaluate different euthanasia methods or their effect on the thermal response of rats.However, studies addressing the mechanism of action of CO 2 and isoflurane, as well as its overdose to induce euthanasia, have shown that both drugs trigger physiological stress and cardiovascular alterations in laboratory rodents [26].Although both are considered safe, inexpensive, and effective methods to induce unconsciousness, the results obtained agree with previous findings suggesting that gases with anesthetic properties cause high levels of aversion and stress in mice and rats [15,52].
Before discussing the IRT results according to the thermal windows, the current literature regarding endocrine and behavioral changes after CO 2 and isoflurane euthanasia in rats and mice suggests that the thermal response due to stress/distress could be associated with both methods.CO 2 inhalation used in G 2 rats induces hypercapnia, acidosis, and suppression of the synaptic potentials [52][53][54] and activates the stress-mediated sympathetic HPA and SAM axes [55].This was corroborated by Borovsky et al. [56], who reported that exposure of rats to CO 2 increased their blood pressure (50-60 mmHg) and NE concentration up to ten times in response to hypoxia.CO 2 euthanasia at 10% is also related to distress due to a high exhibition of anxiety behavior in rats [11].Moreover, aversion to high concentrations of CO 2 (more than 40% of the induction chamber volume) is potentially associated with carbonic acid formation in the mucous membranes, causing irritation and discomfort behaviors such as spinning and pawing in rodents [57].Acute exposure of mice to CO 2 , increased NE and adrenocorticotropic hormone (ACTH) levels [55], and studies in outbred mice and rats have shown that CO 2 inhalation increases total serum protein levels, a biomarker associated with stress [58].In contrast, Hackbarth et al. [59] found no behavioral signs of distress or endocrine alterations (ACTH and corticosterone) in rats undergoing CO 2 euthanasia.
Studies comparing the effect of both CO2 and isoflurane euthanasia have also shown different results.Isoflurane administered in G 4 rats is a volatile anesthetic that causes depression in the cardiorespiratory centers, leading to hypoxemia and death [26,60].Although isoflurane has been considered as a method of CO 2 euthanasia refinement [14,15], in comparison with CO 2 and other volatile anesthetics, isoflurane has mild pungency [61], causing more aversion responses in laboratory rodents that is possibly due to airway irritation [52,62], air hunger, and dyspnea [52].
Powell et al. [57] found that the use of isoflurane during euthanasia increased anxietyrelated behaviors, agitation scores, and corticosterone concentrations in mice compared to the low CO 2 flow rate (30%), the same one used in the present research.Boivin et al. [63] compared isoflurane anesthesia followed by CO 2 , CO 2 , and barbiturates administration as euthanasia methods in mice.The authors found that, according to ACTH concentrations, barbiturates were less stressful than the other two methods.Nonetheless, cardiovascular alterations and pain/stress-related responses did not differ in the three methods, suggesting that isoflurane does not provide benefits above CO 2 euthanasia.In this sense, Valentine et al. [64] reported that a combination of isoflurane and CO 2 caused more signs of distress in mice and that CO 2 alone has less evidence of stress in the animals.
In contrast to what was mentioned, Makowska and Weary's [15] study reported that CO 2 and inhalant agents are aversive to rodents, though the aversion is lower for isoflurane.Likewise, exposure to high concentrations of CO 2 increased adrenaline and noradrenaline concentrations compared to isoflurane euthanasia.This could be indicative of a stress response, but since no stress-related behaviors (grooming, audible vocalizations) were reported in the CO 2 group, CO 2 could not be considered as more stressful than isoflurane [6].In rats, Zardooz et al. [65] found that plasma corticosterone and insulin levels increased in rats exposed to CO 2 , while isoflurane caused a contrary reaction.In Hickman et al. [55]'s research, ACTH, corticosterone, and noradrenaline levels were detected in rats anesthetized with isoflurane; however, the increase was not as significant as with CO 2 .
The present study did not assess behavioral or endocrine parameters to associate the thermal response of rats to the different euthanasia methods.However, the literature shows that both methods trigger stress-related responses that have physiological consequences for the organism, which can be associated with temperature variations according to the thermal window.

T • ocu
A significant increase in the T • mean of T • ocu from Ti 2 to Ti 3 and from Ti 2 to Ti 4 in G 2 and G 4 , respectively, was observed in the rats.The stress-mediated thermoregulatory impairment that CO 2 and isoflurane cause on thermosensitive neurons due to the acidosis and hypoxic effect [66,67] could explain the increase in T • ocu .
Several studies have shown that epinephrine, NE, ACTH, and corticosterone levels increase after CO 2 and isoflurane exposure due to the potential stress that both drugs cause [6,[68][69][70].Although the present study did not consider these biomarkers for evaluation, their release modifies the vasomotor reaction of the microvasculature, inducing vasodilation in key organs (e.g., the eye) and an increased amount of dissipated heat, registered as higher IRT temperatures [29,71] like the ones observed in T • ocu for G 2 and G 4 rats.
Ocular surface temperature in animals has been used as a thermal window to assess acute stress and pain, indicated by a recorded increase in both cases [72,73].To the authors' knowledge, there are only two studies combining IRT and the effect of isoflurane as an anesthetic [30,74], though these studies did not compare CO 2 and isoflurane as a euthanasic.Gjendal et al. [30] determined that, from three different types of stimulus, isoflurane anesthesia in mice had a marked stress response due to the alterations in ocular temperature.Similarly, Vogel et al. [74] used isoflurane anesthesia and found that ocular temperature changed according to the isoflurane concentration and that this temperature also reflects rectal temperature in rodents.Nonetheless, no association was made with stress.
Conversely, while there is no published evidence on euthanasia and ocular IRT, an increase in ocular temperature was reported in wild rodents (Apodemus mystacinus) as a reflection of SIH during the manipulation of individuals [33] during a fear-conditioned test in rats (increasing the eye temperature by up to 1.5 • C) [75], while in mice SIH and active behaviors were positively correlated [76].Furthermore, in guinea pigs, the ocular temperature increased in relation to negative human interaction (petting) [77].Similarly, Wongsaengchan et al. [78] used eye temperature to assess acute exposure to a stressor (small cage, handling, and restraint cone).The authors found significant increases in the left ocular temperature of females during restraint, together with corticosterone increases.
The data suggest that the peripheral vasomotor changes might respond to the flightfight response when exposed to a stressor.Increases in T • ocu from Ti 2 in all experimental groups suggest that rats perceived stress regardless of the euthanasia method.Nonetheless, knowing that CO 2 and isoflurane inhalation might trigger stress-mediated pathways, this could explain the significant changes observed only in G 2 and G 4 from the application of the drug to LORR, probably due to an increased stress response.Finally, although both groups showed significant increases in T • ocu , G 2 and G 4 maintained overall lower temperatures than the rest of the groups, possibly due to heat loss facilitation due to the vasodilator properties of both drugs [53].A similar result was obtained in Gjendal et al. [30]'s study, where isoflurane anesthesia in mice led to a reduction in T • max due to the hypothermia caused by general anesthetics.

T • ear
Comparable to T • ocu , significant increases in the T • mean of T • ear from Ti 2 were observed in both inhalant groups (G 2 and G 4 ).This pattern was expected because ear temperature assessed at the external ear canal is associated with the carotid artery and hypothalamic temperature, the main structure involved in central and peripheral thermoregulatory adaptations [79], particularly when exposed to stressors.Studies conclude that CO 2 and isoflurane exert acute stress in rodents [52,55,58].
In animals, auricular temperature was associated with stress due to the administration of intraperitoneal drugs and restraining techniques in Wokke [34]'s study, as well as in rabbits during handling [80].In rats, increases ranging between 0.8 and 1.5 • C were observed during conditioned fear reactions [75].In the present study, in all experimental groups, an increase in the T • mean of T • ear was observed.However, only CO 2 and isoflurane caused significant increases.This response and its association with previously reported behavioral and endocrine responses with CO 2 and isoflurane euthanasia/anesthesia might cause SIH.Since authors such as Hutu et al. [36] have concluded that superficial ear temperature is correlated to core temperature (around 37.1 ± 0.2 • C), the increase in T • ear could be the reflection of SIH in G 2 and G 4 .
In contrast to the reported findings, some studies have not found significant changes or decreases in the ear temperature of mice and rats.This might be because of the lack of arteriovenous anastomosis present in other species, such as rabbits [81].Additionally, conflicting results can be derived from the thermal window delimitation used by other authors (e.g., external ear canal or auricular pavilion).

T • dor
An expected increase in T • mean values recorded for T • dor after the administration of the euthanasia method was found in all experimental groups.Particularly, significant differences were reported in G 2 and G 4 , maybe due to the induced acute stress that CO 2 and halogenated anesthetics induce in rats.In the anatomical region where T • dor was evaluated, large amounts of BAT can be found [37].This thermogenic tissue responds to NE release.Borovsky et al. [56] and Hickman [55] reported NE increases after exposure of rats and mice to CO 2 and isoflurane as potential stressors.
Due to these elements the interscapular region was used in the present research to determine the effect of the different euthanasia methods, finding that G 2 and G 4 had significant increases in BAT activity.Similarly, a study by Blenkuš [42] reported the highest dorsal superficial temperatures in mice exposed to stressors (daily handling) and behavioral tests (voluntary interaction and elevated plus maze).Miyazono et al. [82] found that body surface temperature (assessed in the dorsal region of mice) increased after acute stress (e.g., reaction to a predator odor), while SIH is also related to BAT thermogenesis in rats and humans after excessive stress [40].Pain perception in spinal lesion murine models has also shown increases in interscapular temperature, an effect that can be lessened with the administration of analgesic drugs [43].
Therefore, the data suggest that the significant local hyperthermia detected in T • dor of G 2 and G 4 subjects could be perceived as a negative stimulus, particularly in both experimental groups.

T • tai
A progressive reduction in T • tai was observed for T • prox, T • medial, and T • distal, regardless of the experimental group.This is due to the vasoconstrictor effect of catecholamines on the microcirculation of peripheral regions such as the tail and paws, as well as the subsequent reduction in radiated heat detected by IRT [29,41].
In different studies, the superficial temperature of the tail has been used to assess stress and the emotional responses of laboratory rodents, where the effect is a reduction from basal values after the exposure of the stressor, as found in the present research.Exposure to stressors such as handling and restraint have been shown to reduce the tail temperature of rats [83].Fear-conditioned rats registered a gradual decrease in tail surface temperature of up to 5.3 • C [75], while the temperature in the tail also decreased in mice during an elevated plus maze test as a result of stress and anxiety [32].Likewise, Blenkuš et al. [42] reported decreases in tail temperature after 30-60 s of exposure to a stressor (daily handling, voluntary interaction test, and elevated plus maze test).Furthermore, Weitkamp [83] has mentioned that T • tai can not only serve to identify acute stressors, but also provides insight into their intensity.This is relevant and consistent with the present results because, even though all euthanasia methods resulted in thermal changes associated with stress, only the inhalant agents caused significant effects in T • tai and all thermal windows.Although a reduction in T • tai was reported in all euthanasia methods (confirming that euthanasia elicits stress-related changes regardless of the method), significant decreases were observed in G 2 and G 4 , and the lowest T • tai values were recorded in both groups.This could be due to the potent vasodilation properties of CO 2 and isoflurane [53,84], triggering a circulatory shift to restrict peripheral circulation (T • tai ) and redirect the blood flow to central sites (T • ocu , T • ear , and T • dor ).

Effect of Sex on the Thermal Response of Animals during Euthanasia Methods
The results obtained when evaluating the effect that sex has on the thermal response showed, in general, no marked differences between males and females.In total, four statistical differences were found in G 6 for T • ocu , G 5 for T • ear , and G 1 for T • dor , showing an inconsistency in the results.This agrees with what was reported by Zevgolis et al. [33] regarding the eye temperature of mice, where no differences by sex were found.
In contrast, Faraji and Metz [28] reported differences between male and female mice.Females exposed to rearing deprivation as a stressor exhibited increased superficial temperature in the head and the back and a decrease in tail temperature, while males did not have differences.Another study from the same authors concluded that rats also show temperature differences when evaluated through IRT and that the stress responses of males and females differ depending on the sex of the experimenters [85].Likewise, apart from the differences reported between females and males, the temperature of a specific thermal window can also differ depending on the sex, as shown in a study where, according to IRT, females were prone to show an exacerbated stress response to restraint [78].
A possible explanation for the lack of significant differences between males and females in the present study could be due to the short period of evaluation used for each euthanasia method.Euthanasia times are (and must be) short so as to avoid high levels of stress.Although, as shown in Table 5, G 4 (294.2± 74.3 s) and G 2 (390.2 ± 171.4) were two of the three euthanasia methods with longer time of death, this time might not have permitted the finding of differences according to sex.Powell et al. [57] mention that rodents require at least two minutes of stressor exposure to increase corticosterone values in response to stress.Nonetheless, since IRT has not been previously evaluated during euthanasia methods considering both sexes, future research needs to consider these factors.

Time of Death and Additional Findings
Regarding the time of death, time of LORR, and cessation of RR and HR, the times obtained in the six experimental groups are in accordance with previous studies evaluating time of death with pentobarbital [18], CO 2 [86], decapitation [87], isoflurane [6], and ketamine + xylazine [8].
Lastly, a distinct pattern and pronounced difference between G 4 and the rest of the experimental groups should be noted.T • ocu , T • ear , T • dor , and T • tai for G 4 rats showed a progressive increase in the surface temperature from Ti 2 to the death of the animals, apart from recording the lowest temperatures from Ti 2 to Ti 4 when compared to the other five groups.In contrast, animals from the other groups, including G 2 , presented a temperature increase from Ti 2 to Ti 3 and a subsequent decrease in all thermal windows.This suggests that the anesthetic stress and physiological response triggered by isoflurane is more marked than that induced by other inhalant, injectable, and physical methods of euthanasia.The present results are in agreement with what other authors have stated regarding isoflurane as a refinement method for CO 2 [52,62] and affirm that precautions should be taken when deciding to use isoflurane as a sole method for the humane killing of research animals.

Limitations and Future Recommendations
The main limitation of the current study, and a field for complementary research using IRT aimed to evaluate euthanasia methods, is the lack of monitoring using physiological markers such as NE, ACTH, corticosterone, glucose, and other parameters (e.g., rectal temperature) that have been reported to increase their concentration during the application of different types of euthanasia [26,55].Moreover, histological analyses could also help to identify the possible tissular changes associated with an inflammatory response to different drugs, providing additional information according to the euthanasia method.Additionally, analyzing the time of death, IRT response, and other biomarkers could help to understand the influence of the application speed and the thermal response of rodents.In the present study, the novel conception of an anesthesia induction chamber designed to allow for IRT readings during inhalant euthanasia is a valuable tool that might serve to further assess how euthanasia drugs, in combination with other physiological, endocrinal, and behavioral parameters, can contribute to the refinement of animal research.
In this sense, an important finding of the present study that can be considered for future research as a refinement in euthanasia procedures in Wistar rats is the combination of injectable agents and CO 2 .As the results showed, contrary to the use of CO 2 alone, the combination administered in G 6 diminished the thermal alterations observed in G 2 .This could be due to the sedative properties of ketamine before CO 2 exposure, antagonizing NMDA receptors, modulating neuronal activity, and reducing the discomfort sensation with CO 2 [88,89].This could prevent physiological responses due to induced hypoxia, acidosis, and stress-related changes [12].
Regarding the non-significant effect of sex in the thermal response of the subjects, the contrasting information between the present results and the published literature shows the complexity of using IRT as a tool to evaluate stress.For example, this variable and the other factors mentioned by Wongsaengchan et al. [78] (e.g., period of evaluation, sex, left/right side for the taking of thermal image) are important elements that need to be considered in further studies where IRT is intended to be used as a tool to improve the welfare of laboratory rodents.Similarly, the weight of rodents should also be considered when using IRT because the thermal response of animals might differ according to their energy reserves and metabolic activity (e.g., obesity in mammals is associated with increased depots of adipose tissue) [90].Moreover, studies have shown that external traits such as coat color or type of fur can affect the amount of radiated heat [91].In the present study we only used Wistar rats (white coat); however, when using IRT in other strains or species, these traits need to be addressed to objectively interpret thermal imaging.
Considering that IRT serves as a non-invasive method to assess the thermal response that can be associated with vasomotor changes due to sympathetic activation, IRT could be implemented as a complementary tool to evaluate stress under other conditions (e.g., heat stress).Likewise, pain assessment and even disease detection can be other fields where thermal imaging could be applied together with biomarkers and other technologies with the aim of improving laboratory animal welfare [24,[29][30][31]73].

Conclusions
Based on the results obtained, it can be concluded that CO 2 and isoflurane elicit stressmediated thermal responses during rat euthanasia.In particular, isoflurane exposure might be a euthanasia method that causes potential distress, and this must be considered when deciding to use this drug as part of a euthanasic protocol.Refinement techniques such as the combination of ketamine + CO 2 were shown to minimize the alterations observed with the sole use of CO 2 , but further research is required to perform a comprehensive evaluation of this alternative.Furthermore, the present study shows the usefulness of IRT as a non-invasive tool for the evaluation of euthanasia techniques and the thermal response of laboratory rodents.In this way, thermal imaging could be recommended together with other physiological, endocrinal, and behavioral parameters to assess and improve the welfare of research animals.

Figure 1 .
Figure 1.Customized acrylic euthanasia chamber for thermal imaging.(A) shows the components of the chamber, with the five hermetic seal doors and the respective gas inlet and vent port.(B) shows the position of the thermal camera that ensured the lens fit into each gate during the evaluation of inhalant euthanasia.

Figure 1 .
Figure 1.Customized acrylic euthanasia chamber for thermal imaging.(A) shows the components of the chamber, with the five hermetic seal doors and the respective gas inlet and vent port.(B) shows the position of the thermal camera that ensured the lens fit into each gate during the evaluation of inhalant euthanasia.

Figure 2 .
Figure 2. Representation of the four evaluated thermal windows.(A).T°ocu was delimited by a circle (El1) covering the entire ocular region or ocular globe, without including the upper or lower eyelid.T°ear was evaluated using a circle (El2) in the external ear canal to assess the irradiated temperature of the tympanic membrane and inner ear.(B).For T°dor, a rectangle (Bx1) was placed in the dorsal area over the interscapular space.(C) For T°tai, three spots (Sp1, Sp2, y Sp3) were placed at the proximal (T°prox), medial (T°medial), and distal (T°distal) segments of the tail.

Figure 2 .
Figure 2. Representation of the four evaluated thermal windows.(A).T • ocu was delimited by a circle (El1) covering the entire ocular region or ocular globe, without including the upper or lower eyelid.T • ear was evaluated using a circle (El2) in the external ear canal to assess the irradiated temperature of the tympanic membrane and inner ear.(B).For T • dor , a rectangle (Bx1) was placed in the dorsal area over the interscapular space.(C) For T • tai , three spots (Sp1, Sp2, y Sp3) were placed at the proximal (T • prox), medial (T • medial), and distal (T • distal) segments of the tail.

Figure 3 .
Figure 3. Experimental timeline for the euthanasia methods applied in rats.HR: heart rate; RR: respiratory frequency.* for G3, Ti3 includes LORR and HR/RR cessation.

Figure 3 .
Figure 3. Experimental timeline for the euthanasia methods applied in rats.HR: heart rate; RR: respiratory frequency.* for G 3 , Ti 3 includes LORR and HR/RR cessation.

Table 1 .
T • mean values in G 2 (33.79 ± 0.92 • C) and G 4 (29.30± 1.23 • C) registered at Ti 2 .Effect of the six euthanasia methods, assessed at five evaluation times, on the maximum, minimum, and mean surface temperature (mean ± standard deviation, SD) of T • ocu ( • C) in Wistar rats.

Table 2 .
Mean± standard deviation (SD) of T • ear ( • C) values of the six euthanasia methods, assessed in five evaluation times, registering the maximum, minimum and mean surface temperatures.

Table 3 .
Mean ± standard deviation (SD) of T • dor ( • C) maximum, minimum and mean values of the six euthanasia methods, assessed in five evaluation times.

Table 4 .
Mean ± standard deviation of T • tai ( • C) values at the proximal (T • prox), medial (T • medial), and distal (T • distal) segments of the tail, assessed on the six euthanasia methods at five evaluation times.

Table 5 .
Comparison between recorded times (in seconds) according to the experimental group (mean ± SD).