Abstract
The relationship between sublethal pesticide exposure and oxidative stress in an ecologically relevant field setting is relatively unknown for reptiles. Oxidative stress is a multi-faceted concept that dictates key survival and fitness parameters in any organism. Fipronil and fenitrothion are two pesticides widely used globally for agricultural pest management. Using a field-based, BACI designed experiment we investigated the impact of sublethal pesticide exposure on oxidative stress biomarkers protein carbonyl and DNA damage (8-OHdG), in an arid-zone lizard species, Pogona vitticeps. A single ecologically relevant dose of pesticide was applied via oral gavage to treatment animals. Lizard condition, activity measures, and blood biomarkers were measured at relevant sampling intervals. Cholinesterase (ChE) and acetylcholinesterase (AChE) enzymatic biomarkers were measured in response to fenitrothion, and fipronil blood residues were measured for fipronil-treated lizards. Results suggested no significant treatment effect of either pesticide on parameters measured, however, 8-OHdG levels decreased by ≥ 45% for both pesticide treatment groups and not controls. Protein carbonyl levels showed a high degree of individual variation that proved more influential than pesticide exposure. Building our understanding of the macromolecular impacts of sublethal pesticide exposure on wild lizard populations is an integral step in addressing the current gap in literature and management practices. Our study has also highlighted the complex nature of studying oxidative stress in the field and the sheer necessity of future study.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
Oxidative stress is integral in key biological processes across all taxa (Costantini 2019). Investigating and understanding oxidative stress in living organisms is complex and multi-faceted. It is well understood that the oxidative stress system is based on a careful balance between eustress and damage, the status of this balance determining the level of oxidative stress present at any given time (Finkel and Holbrook 2000). Reactive oxygen species (ROS) are in a continuous cycle over the life of an organism and are kept under threshold through endogenous and exogenous antioxidant defence systems (Sies 1997, Sies et al. 2017). Once this balance is broken there is a potential for detrimental oxidative stress which targets various areas of the body. Variation in individual personality and physiology, disease status, reproductive status, age, sex, external environment and climate, xenobiotics and anthropogenic changes all influence oxidative stress levels and must be accounted for (Isaksson et al. 2011, Beaulieu and Costantini 2014). Furthermore, as oxidative status is reflective of life history traits and the environment an organism inhabits, we can use tailored biomarkers to inform us about the status of organisms in that environment. Oxidative stress biomarkers have the potential to be extremely useful in identifying the impact that anthropogenic threats, such as xenobiotic compounds, have on a population. This information can be used to formulate relevant guidelines for the conservation of threatened species (Beaulieu and Costantini 2014).
Pesticides are a key anthropogenic stressor that affect populations of wild animals worldwide (Gibbons et al. 2000, Todd et al. 2010, Fasola et al. 2021). They are used in many aspects of human life, from pest control in the home to large scale spray operations to protect agricultural assets (Simon-Delso et al. 2015). Both fipronil (a phenyl-pyrazole insecticide) and fenitrothion (an organophosphate pesticide) are commonly and widely used pesticides globally: fipronil in particular contributes significantly to the world insecticide market (Simon-Delso et al. 2015). Fipronil (5-amino-3-cyano-1-(2,6-dichloro-4-trifluoromethylphenyl)-4-trifluoromethylsulfinyl pyrazole, CAS 120068-37-3) targets critical processes in the central nervous system (CNS) and disrupts the action of gamma-amino butyric acid (GABA). This is done via binding to post-synaptic receptors that block essential chloride channels necessary for GABA innervation (Hainzl et al. 1998, Narahashi 2010, Story et al. 2021). Fenitrothion (O,O-dimethyl-O-(3-methyl-4-nitrophenol)-phosphorothioate, CAS 122-14-5) similarly targets the CNS, but also the peripheral nervous system (PNS). Specifically, fenitrothion targets cholinergic innervation through the inhibition of cholinesterase (ChE) activity most importantly impacting acetylcholinesterase (AChE) activity (Story and Cox 2001, Weir et al. 2010). Both pesticides are expected to cause oxidative stress (Banerjee et al. 2001, Gardner and Oberdorster 2016, Shah and Parveen 2022), by inactivation or down-regulation of antioxidant defence systems decreasing the antioxidant potential of cells or increasing ROS levels through biotransformation to electrophilic or free-radical intermediates (Banerjee et al. 2001, Lushchak 2011). It is possible to detect these macromolecular changes through the use of oxidative stress biomarkers, which have the ability to detect early-stage impacts to environmental contaminants, such as pesticides (Mitchelmore et al. 2005).
Proteins, lipids and DNA are the key targets for oxidative stress. These macromolecule systems are important markers of oxidative stress and are used widely in oxidative stress research (Kaur et al. 2014). Although lipid peroxidation is widespread, it is highly susceptible to free radical attack of any kind (Butler et al. 2016). Eating increases circulating triglyceride levels in the blood which can be reflected in higher levels of lipid peroxidation (Pérez-Rodríguez et al. 2015) which is termed postprandial oxidative stress (Sies et al. 2005, Wallace et al. 2010). Given unknown eating habits of wild organisms, lipid peroxidation is not a particularly useful measure for measuring oxidative stress when sampling wild populations. In contrast, protein carbonyl has been suggested as a good proxy for generalised oxidative stress (Reznick et al. 1992, Shah and Parveen 2022). As proteins are highly abundant and present at high concentration in cells, they are a major targets and likely sites for oxidative damage (Davies 2005, Davies 2016, Hawkins and Davies 2019). Reactive oxygen species target proteins rapidly resulting in early formation of carbonyls throughout the oxidative stress timeline. Biochemical modifications cause post-translational changes that alter amino acids causing changes in protein structure and function (Dalle-Donne et al. 2003, Gianazza et al. 2007). Potential oxidative damage is skewed toward protein modification in many situations, and thus is a more specific marker for measuring oxidative stress and damage (Davies 2005, Davies 2016). DNA damage caused by oxidative stress is also a widely used and reliable biomarker of interest due to its links to genetic mechanisms related to aging, as well as many biological disorders (Nikitaki et al. 2015). Damage to DNA involves the formation of DNA crosslinks via oxidation, methylation, depurination and deamination and basic structural changes as a result of single and double-stranded DNA breaks (Sies 1997, Kaur et al. 2014). It is widely accepted that the most important DNA lesions and thus biomarkers are those associated with hydroxylation and modification of purine and pyrimidine bases (Slaninova et al. 2009). This is generated as ROS attacks DNA at the C-8 position of 2-deoxyguanosine in turn forming the DNA lesion 8-hydroxy-2-deoxyguanosine or 8-OHdG (Kaur et al. 2014). Furthermore, 8-OHdG impacts mechanisms used for the regulation of gene expression by altering the methylation of cytosines (via enzyme-catalysation) (Thompson 2004, Valko et al. 2007). As it is the most critical DNA lesion across taxa with high prevalence, 8-OHdG is an ideal biomarker of oxidative stress to assess individual and population changes in response to sublethal pesticide exposure (Thompson 2004, Slaninova et al. 2009). Protein carbonyl and 8-OHdG are ideal markers to select and build on previous studies on reptiles (Amaral et al. 2012, Mingo et al. 2017, Isaksson 2020, Fasola et al. 2021).
Oxidative stress research is generally lacking across reptilian taxa. With global declines in reptile species, it is important to explore whether xenobiotics or environmental contamination, particularly pesticides, are a key factor (Gibbons et al. 2000, Todd et al. 2010). Reptiles, specifically squamates are susceptible to both indirect and direct impacts of pesticide exposure given their life history traits and ecosystem role as secondary consumers (Sparling et al. 2010, Gardner and Oberdorster 2016, Contador-Kelsall et al. 2022). Lizards are excellent bioindicators of environmental pollution and reveal valuable information about the ecosystems they inhabit (Campbell and Campbell 2002, Silva et al. 2020). The handful of studies investigating the impact of pesticides on oxidative stress levels in reptiles have only begun to illuminate the complexity of this topic. The wide selection of biomarkers available in conjunction with the use of different sample types and timing creates a vast number of questions that need to be answered. To address some of these questions, we have conducted a field-based experiment to investigate the impact of sublethal pesticide exposure on a range of physiological, behavioural, and molecular biomarkers of oxidative stress in a widespread, arid-zone, native Australian lizard, the Central Bearded Dragon (Pogona vitticeps, Ahl 1926). Biomarkers of oxidative stress have been selected across all levels of biological organisation and sampled along a time series to build a comprehensive picture of such effects. We expect that sublethal pesticide exposure will alter oxidative stress status, but we will also build baseline information and validate the use of protein carbonyl and 8-OHdG as appropriate biomarkers of oxidative stress for field-based research.
Materials and methods
Study species and site
Pogona vitticeps is a widespread, medium-bodied, terrestrial arid-zone agamid species that inhabits semi-arid and arid Australia (Cogger 2018). Pogona vitticeps are semi-arboreal omnivores that inhabit various desert habitats, dry sclerophyll forests, cypress pine woodlands, Acacia scrubs, and eucalypt woodlands (Rej and Joyner 2018). This species occurs in rangelands potentially subject to pesticide applications used in locust control (Bain et al. 2004).
The study was conducted from mid-January until mid-March 2019 at Nombinnie Nature Reserve (33°04’25.34”S, 145°47’28.39”E). Nombinnie Nature Reserve is situated in central-western NSW in Australia and consists of 70,000 ha of Mallee, White Cypress Pine, Bimble and Red Box, and Belah woodlands amongst lignum and other vegetation (NSW National Parks and Wildlife Services 2022). Mean monthly temperatures from January to March in 2019 were 32.5, 25.6, and 22.7 °C respectively, with the highest temperature recorded at 46.9 °C in January 2019. Historically (1885–2009) mean spring–summer (September–March) rainfall was 224.7 mm (Australian Bureau of Meteorology 2021). In spring – summer 2017/2018 total rainfall for the area was 165.5 mm and in 2018/2019 total rainfall was 220.4 mm (Australian Bureau of Meteorology 2021).
Animal capture, handling, and tracking
Pogona vitticeps were captured, sampled, and released between 17 days – 1 day pre-exposure to pesticide, as described in Contador-Kelsall et al. (2022). Lizards were initially captured opportunistically during visual scanning surveys. Once captured, lizards were morphologically sexed by hemipenal inspection, and snout to vent length (SVL) and mass (g) were measured to determine scaled body mass indices (SBMI) (Peig and Green 2009). A blood sample of 200 µL was taken from the caudal tail vein. Whole blood was partitioned for haemoglobin (Hb) measurement, blood spot cards (PerkinElmer 226 Five Spot cards) and haematocrit tubes for later plasma collection. Further methodology detailing processing of blood samples can be found in the Appendix. Following sampling, a Sirtrack VHF tag (Sirtrack Ltd, Havelock North, New Zealand) and HOBO accelerometer (HOBO® Pendant G, Onset Computer Corporation) were fixed together and attached to each lizard’s tail dorsally using brown surgical tape. Pilot data suggested lizards were rarely active during the night (Bernich et al. 2022) so accelerometers were set to record data points every 30 s during daylight.
Following initial capture, lizards were sampled at planned sampling sessions for the duration of the experiment. Sampling sessions consisted of; post 24 h, post 7 days, post 14 days, and post 28 days after pesticide dosing. For all subsequent sampling we located each lizard using their allocated VHF tag and took body measurements (SVL, mass (g)) and a blood sample (200 μL). At the final sample session (post 28 d) VHF tags and accelerometers were removed permanently.
Pesticide dosing
To create an ecologically relevant, sublethal pesticide dose, dose was calculated based on known pesticide residue levels found on natural vegetation and insect prey (after aerial spraying) with maximum feeding rates of P. vitticeps, and the average dietary intake of an adult lizard (30% insect, 70% vegetation) from unpublished feeding data and previous pesticide deposition literature (Szabo 2005, Story et al. 2013). In field studies, mean fenitrothion residue levels were found to be 62.0 µg/g on vegetation and 39.8 µg/g on locusts (Story et al. 2013) and fipronil residue levels were between 1.4 and 5.6 µg/g on vegetation (1.4 µg/g used in this experiment) and 0.124 µg/g on locusts (Szabo 2005, Story et al. 2021). According to treatment, lizards received 16.60 mg/kg of fenitrothion or 0.51 mg/kg of fipronil in a single dose. A stock solution was formulated for each treatment based on this dosage rate, with exact individual doses calculated after the initial sampling session once body mass (g) was recorded.
Lizards were allocated to one of 3 treatments at random, with fipronil (n = 5), fenitrothion (n = 5), and technical grade corn oil (control, n = 6). Both pesticides were suspended in technical grade corn oil, to mimic the formulations used in aerial spraying (Kitulagodage et al. 2011). Fipronil was first dissolved in < 150 µL acetone given its low water solubility (acetone solubility = 54.6 g/100 mL) (Tingle et al. 2003) and then suspended in technical grade corn oil. A Covidien Kendall™ plastic feeding tube (6 x 410 mm) was used to deliver pesticides and/or corn oil via oral gavage.
Biomarker assays
Fipronil
All samples were prepared according to Raju et al. (2016) and previously described in Contador-Kelsall et al. (2022). An Agilent 6490 Triple Quad Liquid Chromatography – mass spectrometry (LCMS) was used to analyse all samples. Fipronil parent and metabolite compounds (fipronil, fipronil desulfinyl, fipronil sulfide, and fipronil sulfone) standards were purchased from Sigma-Aldrich. The R2 for each calibration curve was above 0.996. All samples were run in negative mode and clearly defined, confirmed by their most abundant product ions at optimised collision energies. To ensure there was no contamination throughout, a positive control (0.01 μg/mL fipronil) and a negative control (acetonitrile) were run every five samples.
Fenitrothion
All plasma samples were analysed using the Ellman assay (Ellman et al. 1961), as modified by Gard and Hooper (1993), and described in detail in Buttemer et al. (2008) and Contador-Kelsall et al. (2022). An optimum dilution ratio of 1:5 was applied to all samples. Samples were run in duplicates or triplicates dependent of sample availability. Blank samples without the addition of enzyme and mouse serum (Sigma-Aldrich, Australia) were used as a between-assay standard. Total plasma cholinesterase (TChE) and acetylcholinesterase (AChE) activity were directly measured from the plasma samples. Butrylcholinesterase (BChE) activity was then calculated as the difference between TChE and AChE. Exposure to fenitrothion was indicated as a suppression in TChE and AChE when compared to control samples.
Quantifying protein carbonyl
Prior to protein carbonyl analysis, we performed a bicinchoninic acid assay (BCA) to quantify protein in each frozen plasma sample. Quantification was performed using Pierce™ Rapid Gold BCA Protein Assay Kit (ThermoFisher Scientific, Sydney, Australia) following the microplate procedure using the FLUOstar Omega plate reader (BMG Labtech, Ortenberg, Germany). Each sample was run as a duplicate and absorbance was read at 450 nm, with the assay range limits between 20 to 2000 µg/mL. Each plasma sample was diluted according to BCA assay results at 10 µg/mL in 1X PBS in preparation for the protein carbonyl assay. Once diluted the samples were used within 24 h to prevent freeze-thaw cycle degradation of samples.
Oxidative stress was measured via protein carbonyl derivatives using the OxiSelect™ Protein Carbonyl ELISA kit (Cell biolabs, San Diego, CA). The ELISA was conducted according to the manufacturer’s protocol, with the final incubation time of 4–6 min, prior to adding the stop solution. Each sample was run in triplicate and an average used for final analysis. Absorbance was read at 450 nm using the FLUOstar Omega plate reader. Omega MARS data analysis software (BMG Labtech, Ortenberg, Germany) performed best-fit analysis and created a standard curve based on protocol standards, including standard 8 as the blank, where 3rd Polynomial fit based on blank corrected data was selected according to the equation: Y = offset + c1x + c2x2 + c3x3.
Quantifying DNA Damage (8-Hydroxy-2’-deoxyguanosine)
Frozen plasma samples were used for the DNA Damage Competitive ELISA kit (ThermoFisher Scientific, Sydney, Australia). Prior to running the ELISA, we created a dilution series of plasma with 1X Assay Buffer ranging from 1:4 to 1:64, based on (Olsson et al. 2012) . All plasma samples were diluted to optimal dilution (1:16) with 1X Assay Buffer prior to the assay. All samples, standards and blanks were run in triplicate for the ELISA, with controls run in duplicate. Standards 1–9 were diluted as per manufacturer’s instructions, in a serial dilution ranging from 8 000 pg/mL to 0 pg/mL. The ELISA was run according to the manufacturer’s protocol, with absorbance read at 450 nm 10 min after adding the Stop Solution using the FLUOstar Omega plate reader. Omega MARS data analysis software was used for best-fit analysis based on a 4-parameter best fit and created a standard curve based on protocol standards. Standard curve was created according to the equation:
Activity measures
Activity datapoints were measured every 30 s in g-force on an X and Y axis. To ensure no confounding impacts of human disturbance, we removed data surrounding capture and release. All data points from 10 min before capture and 20���min after release were removed from the dataset. We categorised a threshold for movement according to the smallest amount of variance that represented lizard movement, 0.0199 g2 based on Bernich et al. (2022) and assessed data accordingly. Remaining data for each sample session was grouped in 10-minute intervals and then added into daily totals of minutes moved from 6:00–19:00. Three consecutive 24 h periods were selected either before or after (depending on sample session) all set sampling times to create the dataset for activity and ensure consistency across each sample session. According to each activity dataset, a mean of total minutes moved for each session was recorded and used as the final activity measurement.
Statistical analysis
To determine which variables influence protein carbonyl and DNA damage (8-OHdG) levels in P. vitticeps, we measured 6 different explanatory variables. These variables were: treatment (control, fipronil, fenitrothion), time, mean daily temperature, scaled body mass index (SBMI), haemoglobin (Hb), and activity. A total of 13 lizards were used to analyse the protein carbonyl dataset (n = 4–5 per treatment) and 14 lizards used to analyse the DNA damage dataset (n = 4–6 per treatment). Protein carbonyl analysis was based on three sample time points: Pre, Post 24 h, Post 7 d, of which four individuals had a time point missing (total dataset = 35). DNA damage analysis was based on four sample time points: Pre, Post 7 d, Post 14 d, and Post 28 d, where four individuals had 1–2 missing time points (total dataset = 49).
We ran a Pearson correlation matrix for all continuous explanatory variables in both protein carbonyl and DNA damage datasets to assess collinearity between variables. Most explanatory variables assessed for DNA damage models had r-values < 0.5, based on the understanding that Pearson correlation coefficient > 0.5 indicated a highly correlated set of variables (Zuur et al. 2009). Time and temperature were significantly positively correlated (R = 0.50, P < 0.05) and therefore temperature was removed from the model prior to analysis. In the protein carbonyl dataset, explanatory variables time and temperature were also significantly negatively correlated (R = −0.68, P < 0.001) and thus temperature was removed as an explanatory variable prior to further analysis. Mauchly’s test of sphericity was performed for the DNA damage dataset and no assumptions were violated (P > 0.05).
Model selection via backward elimination was used to select the best-fit model to explain variation in both protein carbonyl and DNA damage levels. Linear mixed-effect models for repeated measures were used to analyse the datasets. The lme4 (Bates et al. 2015) and lmerTest (Kuznetsova et al. 2017) packages were used in R Studio (ver 1.4.1717, RStudio Team 2021). The MuMin package (Bartoń 2022) was used to obtain marginal and conditional R2 and the AICcmodavg package (Mazerolle 2020) was used to obtain AICc values. All models were run using REML and Satterthwaite’s method of approximation. Full candidate models included all explanatory variables, treatment, time, and their interaction term, activity, SBMI, and Hb as fixed effects, with lizard ID as a random effect. Candidate full models were as follows:
Protein carbonyl: lmer(Protein Carbonyl ~ Treatment*Time + Activity + SBMI + Hb + (1|Lizard ID), data = PC)
DNA damage: lmer(DNA Damage ~ Treatment*Time + Activity + SBMI + Hb + (1|Lizard ID), data = DNA)
Using guidance from AICc (Akaike’s information criterion values corrected for a small sample size) and R2 values each explanatory variable was removed until significance was found, if applicable. Multiple best-fit models were selected, including those with lowest AICc score (Tables 1, 3).
To investigate the influence cholinesterases (ChE) and acetylcholinesterase (AChE) (fenitrothion-treated group) and fipronil sulfone (fipronil-treated group) on protein carbonyl and DNA damage levels we similarly ran linear mixed-effect models. Each model was run independently with either fipronil sulfone, ChE, or AChE as explanatory variables. Models with ChE and AChE also included a treatment interaction term, as these datasets were inclusive of control lizards (Tables 2, 4).
Results
Protein Carbonyl
Individual physiological measures, SBMI and Hb differed regardless of treatment and time for lizards used in the protein carbonyl analysis (n = 11). SBMI ranged from 241.5 g to 422.6 g with a mean (±SD) of 320 g ± 42.8 g and Hb ranged from 69 g/L to 102 g/L with a mean (± SD) of 89.1 g/L ± 9.1 g/L. Change in protein carbonyl varied from −41–68% with a mean (±SD) of 7 ± 30%. Activity measures ranged from 53.3 to 246.7 mins across all individuals and all sample sessions. The mean (± SD) activity was 142.1 min ± 56.7 min.
Protein carbonyl levels (% change from pre-dose levels) did not vary significantly among treatments and across time (Table 1a). Despite backward elimination all resulting models failed to explain changes in protein carbonyl levels. Individual variation was large, as evidenced by the standard errors (Fig. 1). Furthermore, SBMI, activity and Hb. did not significantly influence protein carbonyl levels in P. vitticeps. The full model had a marginal R2 of 19.7% and a conditional R2 of 86.7% (Table 1a), suggesting individual variation across subjects is responsible for protein carbonyl level fluctuations, more so than any explanatory variable. Even in the reduced model (SBMI, Hb, and activity as explanatory variables, Table 1b) there were comparable R2 values confirming that individual variation across subjects explains a large portion of the model. Protein carbonyl levels in fenitrothion-treated lizards increased from 4% at post 24 h to 13% at post 7 d. Similarly, fipronil-treated lizards increased from 0% at post 24 h to 11% at post 7 d. Control lizards saw a decrease from 9% at Post 24 h to 4% at post 7 d.
Fipronil sulfone (metabolite of fipronil) was present in all treatment lizard blood samples and absent in all control lizards for both Post 24 h and Post 7 d. Mean values (± SD) for fipronil-treated lizards were 0.069 μg/mL ± 0.06 for Post 24 h and 0.2149 μg/mL ± 0.097 for Post 7 d. Fipronil sulfone levels increased by 212% from Post 24 h to Post 7 d. No significant relationship was detected between fipronil sulfone residue levels and protein carbonyl levels (Table 2a, Fig. 2a). Marginal R2 and conditional R2 were both very low (Table 2a).
ChE and AChE levels varied similarly in both fenitrothion treated and control groups. Fenitrothion-treated lizards saw a 12% decrease in ChE activity levels from Post 24 h to Post 7 d and a 13% decrease for AChE activity levels. Control lizard ChE activity levels decreased by 3% from Post 24 h to Post 7 d and AChE activity levels increased by 3%. As ChE and AChE was decreasing for fenitrothion-treated lizards, protein carbonyl levels were increasing by 9% over time, conversely to a decrease in protein carbonyl levels of 5% for control lizards. However, neither ChE nor AChE influenced protein carbonyl levels in P. vitticeps (Table 2b, c, Fig. 3). A marginal R2 and conditional R2 values for the ChE and protein carbonyl model were very similar to the AChE model suggesting large levels of individual variation.
DNA damage (8-OHdG)
Physiological parameters (SBMI and Hb) varied among individuals across treatments used in the DNA damage analysis (n = 14). SBMI ranged from 241.5 g to 463.3 g with a mean (± SD) of 328.9 g ± 43.5 g. Hb ranged from 69 g/L to 112 g/L ± 9.7 g/L, with a mean (±SD) of 87.4 g/L. DNA damage (8-OHdG) ranged from a −80 to 160% change from initial sample (Pre-dose), with a mean (± SD) of 54 ± 30% change. Activity measures ranged widely from 0 min (mins) – 206.6 min across all time periods with an overall mean (± SD) of 105.2 min ± 53.6 min.
While there was no significant effect of treatment or time on 8-OHdG levels in P. vitticeps (Table 3), DNA damage levels (% change from pre-dose levels) decreased across time for both fenitrothion and fipronil-treated lizards (Fig. 4). From post 7 d to post 28 d 8-OHdG decreased 60% for lizards in the fenitrothion treatment group. Fipronil-treated individuals showed a 45% reduction in 8-OHdG levels from post 7d to post 28 d. Control lizards fluctuated from an initial increase at post 14 d by 16% followed by a decrease of 8% at post 28 d. The full model prior to backward elimination resulted in a marginal R2 of 38.5% and a conditional R2 of 47.7% (Table 3a). SBMI significantly increased 8-OHdG levels in models including treatment, activity, Hb and subject as explanatory variables (See Appendix, also Table 3b). The best model to explain DNA damage levels in P.vitticeps included only SBMI and Hb (Table 3c). Both SBMI and Hb significantly influence 8-OHdG levels in P. vitticeps independent of treatment effect or time (P < 0.05, Fig. 5). There was a significant positive relationship between SBMI and 8-OHdG levels and a significant negative relationship between Hb and 8-OHdG levels (Table 3c).
Fipronil sulfone was measured in control and fipronil-treated lizards. Fipronil sulfone was present in all treatment lizards for all sample sessions, including post 28 d at the end of the experiment, and was absent from all control lizards for the duration of the experiment. Mean values (± SD) for treatment individuals include, 0.215 μg/mL ± 0.11 (Post 7 d), 0.223 μg/mL ± 0.01 (Post 14 d), and 0.239 μg/mL ± 0.06 (Post 28 d). On average sulfone levels remained high at similar levels for all time periods, with an increase of 11% at post 28 d from post 7 d. There was no significant relationship between fipronil sulfone levels and DNA damage in fipronil-treated P. vitticeps (P > 0.05, Table 4a, Fig. 2b).
Both ChE and AChE activity levels were measured in control and fenitrothion-treated lizards. ChE levels were significantly negatively correlated with 8-OHdG levels (P < 0.05, Table 4b, Fig. 6), where treatment and ChE interaction term did not (P > 0.05, Table 4b). There was no significant effect of AChE or its interaction term with treatment in explaining 8-OHdG levels (Fig. 6, Table 4c). In the fenitrothion treatment group ChE increased 7% between Post 7 d and Post 14 d and decreased slightly to a 5% change at Post 28 d. AChE activity levels increased by 16% from Post 7 d to Post 14 d, and further increased at Post 28 d at 36%. Control lizards on the other hand showed an increase of 13% in ChE between Post 7 d and Post 14 d, and a 15% increase at Post 28 d. Control lizard AChE activity levels plateaued at a 17% increase from Post 7 d to both Post 14 d and Post 28 d.
Discussion
Our results investigating the degree of oxidative stress in pesticide-treated lizards were not greater than natural background levels for both macromolecular (protein carbonyl and 8-OHdG) and enzymatic biomarkers (AChE). Individual variation accounted for a higher level of oxidative damage than sublethal exposure to pesticides, however, non-significant trends in DNA damage (8-OHdG) levels provide some evidence of a relationship with pesticide exposure. As the study was based on wild P. vitticeps individuals we were unable to account for potential drivers of these high levels of individual variation, such as biological and reproductive age, microclimate, disease, and food availability and it is integral to keep this in mind when interpreting the results (Monaghan et al. 2009, Kaur et al. 2014). There is also the potential for a compounding effect of multiple stressors resulting in the oxidative stress levels portrayed in the current study (Bertram et al. 2022). However, utilising a repeated measures time series sampling protocol in our experiment provides some accountability in investigating oxidative stress levels throughout the duration of the experiment (Fasola et al. 2021, Ritchie and Friesen 2022). Our results are therefore cautiously encouraging; ecologically relevant exposure to pesticides may not result in evidence of severe oxidative damage in P. vitticeps.
Despite no significant treatment effect on either of the macromolecular biomarkers (PC and 8-OHdG), there is evidence of a trend in relation to 8-OHdG levels. Both pesticide treatment groups show a ≥45% decrease in 8-OHdG levels from post 7 d to post 28 d, whereas control treatment group does not substantially change; evidence of pesticide exposure generating oxidative stress (Fig. 4). For example, fipronil-treated individuals show a plateau in 8-OHdG levels from post 7 d to post 14 d, followed by a further drop in 8-OhdG levels from post 14 d to post 28 d which correlates with peak fipronil sulfone levels in the blood of treated individuals. The secondary decrease in 8-OHdG levels after post 14 d in the fipronil-treated cohort could be the result increasing pesticide metabolites, as fipronil sulfone is more detrimental and toxic to body systems than fipronil itself (Hainzl et al. 1998, Peveling and Demba 2003). Any decrease seen in both treatment groups primarily suggests the mobilisation of the antioxidant defence system in response to oxidative stress (Sies 1997, Sies et al. 2017, Severo et al. 2020). As fipronil sulfone reaches its peak, it is probable there is an increase in oxidative stress, and thus we see an upregulation of the antioxidant defence system. This likely causes a decrease in 8-OHdG levels as a result of enzymatic antioxidant defences, such as superoxide dismutase (Winterbourn 2016), preventing any oxidative damage to DNA. In contrast, fenitrothion has no secondary metabolites that are more toxic (Shore and Douben 1994), so we see a faster decrease in 8-OHdG levels of 45% at Post 14 d indicating potentially a quicker mobilisation of antioxidant defence systems. The decrease in 8-OHdG levels and assumed upregulation of endogenous antioxidants can be described as eustress and may be an explanation for the results seen for fipronil and fenitrothion-exposed individuals (Sies et al. 2017, Hawkins and Davies 2019), however the potential link between decreased levels and fenitrothion exposure is less clear. The overarching enzyme group, cholinesterase’s, or B esterases indicate a negative relationship between ChE activity levels and 8-OHdG levels. There is no obvious relationship between 8-OHdG and AChE, the primary target of OPs (Thompson 1999), which suggests that this interaction is simply between pre-existing ChE activity and 8-OHdG levels, which we would expect to be present regardless of OP exposure. It remains unclear why ChE significantly correlates to 8-OHdG but we have shown there is no significant influence of sublethal fenitrothion exposure on these levels, despite a trend indicating the upregulation of antioxidant defence systems. Future investigation into the relationship between these markers, both aside from OP exposure and when subjected to higher doses of fenitrothion would yield useful information.
There are multiple pathways that pesticides, such a fipronil or fenitrothion can take to cause oxidative stress. Reactive oxygen species (ROS) created via excessive cycling of cytochrome P450 enzymes, decreased antioxidant potential or downregulation of cellular antioxidants, and alterations to the function of mitochondria (decreasing ability to metabolise and detoxify) all result in oxidative stress to lipids, proteins, and DNA (Banerjee et al. 2001, Lushchak 2011). Additionally, these xenobiotic compounds can join redox cycles and transforming cells so that they become electrophilic (Banerjee et al. 2001, Lukaszewicz-Hussain 2010, Lushchak 2011). However, it is unlikely exposure to fipronil or fenitrothion resulted in a decrease in antioxidant potential given the suggested upregulation evident in changes in 8-OHdG levels in both treatment groups (Fig. 4). Thus, as it is implied that some oxidative stress (namely 8-OHdG) has occurred in P. vitticeps due to pesticide exposure, we suggest fipronil and fenitrothion potentially disrupt cytochrome P450 enzyme systems and/or interfere in key cellular processes within mitochondria or processes associated with transcription or translation (Kitulagodage et al. 2011, Lushchak 2011, Kaur et al. 2014, Isaksson 2015). Measuring these parameters would be very interesting to further understand how ecologically relevant pesticides may impact oxidative stress levels in lizards. Laboratory methodology has been developed to assess changes in mitochondrial function and cytochrome P450/NADPH-cytochrome P450 reductase (Guengerich et al. 2009, Stier et al. 2017, Kowaltowski 2019). Measuring cytochrome P450 enzymes is relatively easy however it must be done in the laboratory with specialised equipment but is a good candidate for future research in this field (Guengerich et al. 2009). On the other hand, mitochondrial function and oxidants are more difficult to measure, and technology is still developing (Kowaltowski 2019). Stier et al. (2017) developed a method to accurately measure mitochondrial function using red blood cells, which could be applied to field studies with reptiles in the future.
In the laboratory and field, we know pesticides cause oxidative stress, and we would expect an increase as a result of exposure (Valdivia et al. 2007, Lushchak 2011, Amaral et al. 2012, Severo et al. 2020). This was evident in 8-OHdG levels in treatment lizards, yet not so for PC levels. Aside from high levels of individual variation dictating PC levels in our experiment we should also consider that the sampling protocol used in our study to detect PC levels (Pre, Post 24 h, Post 7 d) was suboptimal for detecting any potential significant changes. Significant fluctuations in PC levels can occur and return to baseline in as short as 24 h or as long as 7 days depending on the intensity of the stressor, and these temporary impacts can be easily missed (Nikolaidis et al. 2012, Costantini 2019, Isaksson 2020, Marins et al. 2021). While our sampling protocol measured the predicted critical periods, as suggested by the literature, any key changes occurring before post 24 h and between post 48 h – post 7 d may have been missed. Sample type specific to oxidative stress biomarkers is also critical for measuring correct levels of oxidative stress in response to external stressors. For example, Wu et al. (2015) showed that Zebra fish (Danio rerio) had increased levels of ROS in gills, brain and liver tissue in response to cold temperatures, however the time ROS peaked was different depending on the tissue type. Similarly, South-American Catfish (Rhamdia quelen) exposed to a mixture of pesticides at environmental concentrations saw an increase in protein carbonyl levels in liver and brain tissues but a decrease in gills and muscle (Marins et al. 2021). It is likely that in our study, alternate sample types may have shown higher levels of oxidative stress. In terms of sample type, blood is still utilised widely and accepted as a good generalised marker indicative of oxidative status throughout the body (Margaritelis et al. 2015). Margaritelis et al. (2015) concluded that majority of oxidative stress biomarkers measured in blood adequately reflected oxidative status in tissues. This was true for; reduced glutathione, SOD, CAT, Malondialdehyde, glutathione peroxidase, and vitamin C and E but not so for oxidised glutathione and reduced/oxidised glutathione ratio. Argüelles et al. (2004) found that protein carbonyl levels in blood correlated with heart and liver tissue, but not spleen and kidney, indicating once again the variable nature of oxidative stress. Additionally, as blood is highly integrated throughout the body, it has the potential to have the most variable response to oxidative stress (Isaksson 2015, Speakman et al. 2015, Ouyang et al. 2018, Costantini 2019). At present, using blood to measure oxidative stress changes exhibits methodology that is relatively non-invasive and most importantly non-lethal. Any alternate analysis of various tissues (for example, liver and spleen) would require the subject to be euthanised, which was not the intention in this study. Using blood samples also allow the analysis of a time series, which is incredibly informative when investigating oxidative stress.
This study is one of only a handful to investigate a collection of biomarkers of oxidative stress in a remote field setting (Amaral et al. 2012, Mingo et al. 2017, Fasola et al. 2021, Simbula et al. 2021) and the first to investigate ecologically relevant impacts of pesticides on an important Australian reptile. There is a necessity to continue growing this data because measuring oxidative stress status is incredibly complex. While we found indications of an upregulation of antioxidant defences, evidenced namely in 8-OHdG results, we were unable to use antioxidant (non-enzymatic or enzymatic) biomarkers in our study as these require quick access to laboratory conditions, a situation not available in the field (Selman et al. 2012). The stability of these biomarkers decreases over time depending on storage protocols in conjunction with intrinsic reactivity (Jansen et al. 2013, Pérez-Rodríguez et al. 2015). Additionally, the levels portrayed in these assays may not in fact reflect accurate antioxidant responses based on the complex interactions existing among antioxidants (Pérez-Rodríguez et al. 2015), a further reason for not including such assays in our study. Instead, we selected feasible biomarkers that would not degrade in the field and chose times relevant to oxidative stress time points.
Individual physiological traits of scaled body mass index (SBMI) and haemoglobin (Hb) significantly influenced DNA damage (8-OHdG) levels. Our results suggest that as body condition increases so does DNA damage (Fig. 5a). Page and Stuart (2012) demonstrated that body mass was the key variable that explained variation in DNA repair mechanisms in endotherms. Oligonucleotide ligation activities and nucleotide incorporation was found to be negatively correlated with body mass. Our results show a positive relationship between the two variables, the converse of Page and Stuart (2012) however we see similarity in the significant correlations between body condition and oxidative stress in DNA. A candidate hypothesis as to why SBMI and DNA damage may be related is age. Age is a key factor that must be investigated when examining body condition and oxidative stress levels (Hoekstra et al. 2020). There is conflicting evidence for ectotherms as to whether age is correlated with oxidative stress levels, as they do not consistently fit into the oxidative stress theory of ageing (Buttemer et al. 2010). Hoekstra et al. (2020) summarises all available studies that present correlations or absence of between macromolecular damage and age in reptiles. Specifically for DNA damage, a handful of studies have found a positive correlation with age (Bronikowski 2008, Robert and Bronikowski 2010, Schwartz and Bronikowski 2013), while some have shown no association with age (Woodhead et al. 1980, Regan et al. 1982, Schwartz and Bronikowski 2011). For example, Robert and Bronikowski (2010) found that long-lived ecotypes of Western Terrestrial Garter Snakes (Thamnophis elegans) were in fact smaller than the short-lived ecotypes, and experienced higher DNA damage, but the damage was repaired more efficiently. Additionally, the smaller, long-lived ecotypes had more efficient antioxidant defence systems and mitochondria. Evident here, and in our study, is a complex relationship between body size, condition, age, and oxidative stress, damage, and repair, where a multitude of future studies are needed to help clarify such relationships. It has been demonstrated that reptiles in particular pose more complexity in understanding the correlations between physiology and oxidative stress given the nature of their life histories. External environment, in particular temperature influences the physiology and behaviour of reptiles and in turn metabolic rate and growth, adding a level of intricacy to understanding SBMI and oxidative stress levels (Buttemer et al. 2010).
As haemoglobin levels (Hb) increased 8-OHdG levels decreased (Fig. 5b). Hb can be used in part to describe the fitness and health of a reptile, tending to increase with better fitness (Zayas et al. 2011, Minias 2015, Johnstone et al. 2017). Oxidative stress causes damage to red blood cells (RBC) and the denaturation of haemoglobin molecules (Mohanty et al. 2014). We hypothesise that the healthier the lizard the lower the levels of oxidative stress. Bauerová et al. (2017) and Herrera-Dueñas et al. (2014) both found that oxidative stress was detected in RBCs, in conjunction with lower Hb and packed cell volume (PCV) concentrations in birds that inhabited urban areas opposed to rural areas. This is consistent with Goodchild et al. (2022) where Song Sparrows (Melospiza melodia) inhabiting suburban areas displayed lower Hb concentrations compared to individuals inhabiting rural areas. However, the study did not find any significant differences in any of the oxidative stress markers tested (GSH, d-ROMs, HOCI neutralisation) between sites, which could suggest Hb concentrations are not dictated directly by urbanisation. Cid et al. (2018) also detected significant decreases in Hb concentrations in House Sparrows (Passer domesticus) that were exposed to sublethal doses of lead (Pb). It is probable that our results are also demonstrating a degree of impact on Hb concentrations because of a sublethal pesticide dose, or potentially alternate stressors, given there is no obvious treatment effect (Figs. 7, 8).
We can equate the level of oxidative stress witnessed in our study to that described by Lushchak (2014), as basal oxidative stress (BOS) witnessed in our protein carbonyl results and low intensity oxidative stress (LOS) witnessed in our DNA damage results. Based on this proposed classification system for oxidative stress (based on intensity), BOS represents levels that have no observable impacts caused by the applied stressor and cannot be detected by current methods, a description suited to our protein carbonyl findings. On the other hand, DNA damage clearly showed molecular changes associated with probable increases in antioxidant defence systems (Keap1/nrf2 system upregulation) but with no endpoint parameter changes (increases in 8-OHdG), placing the level of oxidative stress experienced between BOS and LOS (Lushchak 2014). It is evident that a sublethal pesticide dose has likely caused a small increase in oxidative stress, resulting in the upregulation of endogenous antioxidant defence systems. Additionally, there is strong evidence that individual physiological factors and life history traits influence oxidative stress levels, potentially to a larger extent, in wild P. vitticeps. Wild individuals and populations are inherently difficult to study due to the many external factors that cannot be accounted for (Beaulieu and Costantini 2014, Kaur et al. 2014). Despite this, there is urgent need for field-based, ecologically relevant studies so that a relevant understanding of basic oxidative stress responses to not only anthropogenic stressors, such as sublethal pesticide exposure, but also the extent external environment and individual biology influences these levels. Oxidative stress influences critical biological functioning related to all aspects of life and can predict the survivability and reproductive capacity of a species (Isaksson et al. 2011, Beaulieu and Costantini 2014). It is imperative to understand this relationship when assessing impacts on wild reptile populations. The lack of ecologically relevant field studies is impeding formulation of guidelines that can be used for reptile conservation (Ritchie and Friesen 2022). Our study has successfully investigated various biomarkers of oxidative stress ranging from cellular level to whole body physiological changes in response to sublethal exposure to pesticides. Importantly, the study has validated protein carbonyl and 8-OHdG as biomarkers of oxidative stress in Central Bearded Dragons (P. vitticeps) in a wild setting. We have built upon similar previous studies, further investigating potential informative markers, and most significantly demonstrating that the time at which a sample is taken is a crucial element of biomarker selection that must be included in all future studies.
References
Amaral MJ, Bicho RC, Carretero MA, Sanchez-Hernandez JC, Faustino AMR, Soares AMVM, Mann RM (2012) The use of a lacertid lizard as a model for reptile ecotoxicology studies: Part 2 – Biomarkers of exposure and toxicity among pesticide exposed lizards. Chemosphere 87(7):765–774. https://doi.org/10.1016/j.chemosphere.2012.01.048l
Argüelles S, García S, Maldonado M, Machado A, Ayala A (2004) Do the serum oxidative stress biomarkers provide a reasonable index of the general oxidative stress status? Biochimica et Biophysica Acta (BBA) - General Subjects 1674(3):251–259. https://doi.org/10.1016/j.bbagen.2004.06.023l
Australian Bureau of Meteorology. (2021). Monthly climate statistics. Retrieved 26 March, 2021, from http://www.bom.gov.au/climate/averages/tables/cw_049032.shtml.
Australian Bureau of Meteorology. (2021). Monthly rainfall Euabalong (Mount Hope AWS). Retrieved 26 March, 2021, from http://www.bom.gov.au/jsp/ncc/cdio/weatherData/av?p_nccObsCode=139&p_display_type=dataFile&p_startYear=&p_c=&p_stn_num=049136.
Bain D, Buttemer WA, Astheimer L, Fildes K, Hooper MJ (2004) Effects of sublethal fenitrothion ingestion on cholinesterase inhibition, standard metabolism, thermal preference, and prey-capture ability in the Australian central bearded dragon (Pogona vitticeps, agamidae). Environ Toxicol Chem 23(1):109–116. https://doi.org/10.1897/02-555l
Banerjee BD, Seth V, Ahmed RS (2001) Pesticide-Induced Oxidative Stress : Perspective and Trends. Revs Environ Health 16(1):1–40. https://doi.org/10.1515/REVEH.2001.16.1.1l
Bartoń K (2022). MuMIn: Multi-Model Inference.
Bates D, Mächler M, Bolker B, Walker S (2015) Fitting Linear Mixed-Effects Models Using lme4. J Statistical Softw 67(1):1–48. https://doi.org/10.18637/jss.v067.i01l
Bauerová P, Vinklerová J, Hraníček J, Čorba V, Vojtek L, Svobodová J, Vinkler M (2017) Associations of urban environmental pollution with health-related physiological traits in a free-living bird species. Sci Total Environ 601-602:1556–1565. https://doi.org/10.1016/j.scitotenv.2017.05.276l
Beaulieu M and Costantini D (2014) Biomarkers of oxidative status: missing tools in conservation physiology. Conservation physiol 2(1). https://doi.org/10.1093/conphys/cou014l
Bernich A, Maute K, Contador-Kelsall IC, Story PG, Hose GC, French K (2022) Space use and daily movement patterns in an arid zone agamid lizard. Wildl Res 49(6):557–570. https://doi.org/10.1071/WR20152l
Bertram MG, Martin JM, McCallum ES, Alton LA, Brand JA, Brooks BW, Cerveny D, Fick J, Ford AT, Hellström G, Michelangeli M, Nakagawa S, Polverino G, Saaristo M, Sih A, Tan H, Tyler CR, Wong BBM, Brodin T (2022) Frontiers in quantifying wildlife behavioural responses to chemical pollution. Biol Rev 97(4):1346–1364. https://doi.org/10.1111/brv.12844l
Bronikowski AM (2008) The evolution of aging phenotypes in snakes: a review and synthesis with new data. AGE 30(2):169. https://doi.org/10.1007/s11357-008-9060-5l
Butler MW, Lutz TJ, Fokidis HB, Stahlschmidt ZR (2016) Eating increases oxidative damage in a reptile. J Exp Biol 219(13):1969–1973. https://doi.org/10.1242/jeb.138875l
Buttemer WA, Abele D, Costantini D (2010) From bivalves to birds: oxidative stress and longevity. Funct Ecol 24(5):971–983. https://doi.org/10.1111/j.1365-2435.2010.01740.xl
Buttemer WA, Story PG, Fildes KJ, Baudinette RV, Astheimer LB (2008) Fenitrothion, an organophosphate, affects running endurance but not aerobic capacity in fat-tailed dunnarts (Sminthopsis crassicaudata). Chemosphere 72(9):1315–1320. https://doi.org/10.1016/j.chemosphere.2008.04.054l
Campbell KR, Campbell TS (2002) A logical starting point for developing priorities for lizard and snake ecotoxicology: A review of available data. Environ Toxicol Chem 21(5):894–898. https://doi.org/10.1002/etc.5620210502l
Cid FD, Fernández NC, Pérez-Chaca MV, Pardo R, Caviedes-Vidal E, Chediack JG (2018) House sparrow biomarkers as lead pollution bioindicators. Evaluation of dose and exposition length on hematological and oxidative stress parameters. Ecotoxicol Environ Saf 154:154–161. https://doi.org/10.1016/j.ecoenv.2018.02.040l
Cogger HG (2018) Reptiles & amphibians of Australia/Harold G Cogger, CSIRO Publishing, Australia
Contador-Kelsall I, Maute K, Story P, Hose GC, French K (2022) Sublethal pesticide exposure influences behaviour, but not condition in a widespread Australian lizard. Conserv Physiol 10(1):coac024. https://doi.org/10.1093/conphys/coac024l
Costantini D (2019) Understanding diversity in oxidative status and oxidative stress: the opportunities and challenges ahead. J of Exp Biol 222(13). https://doi.org/10.1242/jeb.194688l
Dalle-Donne I, Rossi R, Giustarini D, Milzani A, Colombo R (2003) Protein carbonyl groups as biomarkers of oxidative stress. Clinica chimica acta 329(1-2):23–38
Davies MJ (2005) The oxidative environment and protein damage. Biochimica et Biophysica Acta (BBA) - Proteins and Proteomics 1703(2):93–109. https://doi.org/10.1016/j.bbapap.2004.08.007l
Davies MJ (2016) Protein oxidation and peroxidation. Biochem J 473(7):805–825. https://doi.org/10.1042/bj20151227l
Ellman GL, Courtney KD, Andres V, Featherstone RM (1961) A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem Pharmacol 7(2):88–95. https://doi.org/10.1016/0006-2952(61)90145-9l
Fasola E, Biaggini M, Ortiz-Santaliestra ME, Costa S, Santos B, Lopes I, Corti C (2021) Assessing Stress Response in Lizards from Agroecosystems with Different Management Practices. Bull Environ Contam Toxicol https://doi.org/10.1007/s00128-021-03404-3l
Finkel T, Holbrook NJ (2000) Oxidants, oxidative stress and the biology of ageing. Nat J 408(6809):239. https://doi.org/10.1038/35041687l
Gard NW, Hooper MJ (1993) Age-dependent changes in plasma and brain cholinesterase activities of Eastern Bluebirds and European Starlings. J Wildl Dis 29(1):1–7. https://doi.org/10.7589/0090-3558-29.1.1l
Gardner SC, Oberdorster E (2016) Toxicology of reptiles, CRC Press, United States of America
Gianazza E, Crawford J, Miller I (2007) Detecting oxidative post-translational modifications in proteins. Amino Acids 33(1):51–56. https://doi.org/10.1007/s00726-006-0410-2l
Gibbons JW, Scott DE, Ryan TJ, Buhlmann KA, Tuberville TD, Metts BS, Greene JL, Mills T, Leiden Y, Poppy S, Winne CT (2000) The Global Decline of Reptiles, Déjà Vu Amphibians: Reptile species are declining on a global scale. Six significant threats to reptile populations are habitat loss and degradation, introduced invasive species, environmental pollution, disease, unsustainable use, and global climate change. BioSci 50(8):653–666. https://doi.org/10.1641/0006-3568(2000)050[0653:TGDORD]2.0.CO;2l
Goodchild CG, VanDiest I, Lane SJ, Beck M, Ewbank H, Sewall KB (2022) Variation in hematological indices, oxidative stress, and immune function among male song sparrows from rural and low-density urban habitats. Front Ecol Evol 97:1–12
Guengerich FP, Martin MV, Sohl CD, Cheng Q (2009) Measurement of cytochrome P450 and NADPH–cytochrome P450 reductase. Nat Protoc 4(9):1245–1251. https://doi.org/10.1038/nprot.2009.121l
Hainzl D, Cole LM, Casida JE (1998) Mechanisms for selective toxicity of fipronil insecticide and its sulfone metabolite and Desulfinyl Photoproduct. Chem Res Toxicol 11(12):1529–1535. https://doi.org/10.1021/tx980157tl
Hawkins CL, Davies MJ (2019) Detection, identification, and quantification of oxidative protein modifications. J Biol Chem 294(51):19683–19708. https://doi.org/10.1074/jbc.REV119.006217l
Herrera-Dueñas A, Pineda J, Antonio MT, Aguirre JI (2014) Oxidative stress of House Sparrow as bioindicator of urban pollution. Ecol Indicators 42:6–9. https://doi.org/10.1016/j.ecolind.2013.08.014l
Hoekstra LA, Schwartz TS, Sparkman AM, Miller DAW, Bronikowski AM (2020) The untapped potential of reptile biodiversity for understanding how and why animals age. Funct Ecol 34(1):38–54. https://doi.org/10.1111/1365-2435.13450l
Isaksson C (2015) Urbanization, oxidative stress and inflammation: A question of evolving, acclimatizing or coping with urban environmental stress. Funct Ecol 29(7):913–923
Isaksson C (2020) Urban ecophysiology: beyond costs, stress and biomarkers. J Exp Biol 223(22):jeb203794
Isaksson C, While GM, Olsson M, Komdeur J, Wapstra E (2011) Oxidative stress physiology in relation to life history traits of a free-living vertebrate: the spotted snow skink, Niveoscincus ocellatus. Integr. Zoology 6(2):140–149. https://doi.org/10.1111/j.1749-4877.2011.00237.xl
Jansen EH, Beekhof PK, Cremers JW, Viezeliene D, Muzakova V, Skalicky J (2013) Long-term stability of parameters of antioxidant status in human serum. Free Radic Res 47(6-7):535–540
Johnstone CP, Lill A, Reina RD (2017) Use of erythrocyte indicators of health and condition in vertebrate ecophysiology: A review and appraisal. Biol Rev 92(1):150–168. https://doi.org/10.1111/brv.12219l
Kaur R, Kaur J, Mahajan J, Kumar R, Arora S (2014) Oxidative stress—implications, source and its prevention. Environ Sci Pollut Res 21(3):1599–1613. https://doi.org/10.1007/s11356-013-2251-3l
Kitulagodage M, Buttemer WA, Astheimer LB (2011) Adverse effects of fipronil on avian reproduction and development: maternal transfer of fipronil to eggs in zebra finch Taeniopygia guttata and in ovo exposure in chickens Gallus domesticus. Ecotoxicology 20(4):653–660. https://doi.org/10.1007/s10646-011-0605-5l
Kitulagodage M, Isanhart J, Buttemer WA, Hooper MJ, Astheimer LB (2011) Fipronil toxicity in northern bobwhite quail Colinus virginianus: Reduced feeding behaviour and sulfone metabolite formation. Chemosphere 83(4):524–530. https://doi.org/10.1016/j.chemosphere.2010.12.057l
Kowaltowski AJ (2019) Strategies to detect mitochondrial oxidants. Redox Biol 21:101065. https://doi.org/10.1016/j.redox.2018.101065l
Kuznetsova A, Brockhoff PB, Christensen RHB (2017) lmerTest Package: Tests in Linear Mixed Effects Models. J Statistical Softw 82(13):1–26. https://doi.org/10.18637/jss.v082.i13l
Lukaszewicz-Hussain A (2010) Role of oxidative stress in organophosphate insecticide toxicity – Short review. Pesticide Biochem Physiol 98(2):145–150. https://doi.org/10.1016/j.pestbp.2010.07.006l
Lushchak VI (2011) Environmentally induced oxidative stress in aquatic animals. Aquat Toxicol 101(1):13–30. https://doi.org/10.1016/j.aquatox.2010.10.006l
Lushchak VI (2014) Free radicals, reactive oxygen species, oxidative stress and its classification. Chemico-Biological Interacts 224:164–175. https://doi.org/10.1016/j.cbi.2014.10.016l
Margaritelis NV, Veskoukis AS, Paschalis V, Vrabas IS, Dipla K, Zafeiridis A, Kyparos A, Nikolaidis MG (2015) Blood reflects tissue oxidative stress: a systematic review. Biomarks 20(2):97–108. https://doi.org/10.3109/1354750X.2014.1002807l
Marins AT, Cerezer C, Leitemperger JW, Severo ES, Costa MD, Fontoura DO, Nunes MEM, Ribeiro LC, Zanella R, Loro VL (2021) A mixture of pesticides at environmental concentrations induces oxidative stress and cholinergic effects in the neotropical fish Rhamdia quelen. Ecotoxicology 30(1):164–174. https://doi.org/10.1007/s10646-020-02300-6l
Mazerolle MJ (2020). AICcmodavg: Model selection and multimodel inference based on (Q)AIC(c).
Mingo V, Lötters S, Wagner N (2017) The impact of land use intensity and associated pesticide applications on fitness and enzymatic activity in reptiles—A field study. Sci Total Environ 590-591:114–124. https://doi.org/10.1016/j.scitotenv.2017.02.178l
Minias P (2015) The use of haemoglobin concentrations to assess physiological condition in birds: a review. Conserv Physiol 3(1):cov007–cov007. https://doi.org/10.1093/conphys/cov007l
Mitchelmore CL, Rose CL, Place AR (2005) Tools for Assessing Comtaminant Exposure and Effects in Reptiles. Toxicology of Reptiles. E. O. Susan C. Gardner. United States of America Taylor & Francis Group, United States of America
Mohanty J, Nagababu E, Rifkind J (2014) Red blood cell oxidative stress impairs oxygen delivery and induces red blood cell aging. Fronts in Physiol 5. https://doi.org/10.3389/fphys.2014.00084l
Monaghan P, Metcalfe NB, Torres R (2009) Oxidative stress as a mediator of life history trade-offs: mechanisms, measurements and interpretation. Ecol Lett 12(1):75–92. https://doi.org/10.1111/j.1461-0248.2008.01258.xl
Narahashi T (2010) Chapter 31 - Neurophysiological Effects of Insecticides. Hayes’ Handbook of Pesticide Toxicology (Third Edition). R. Krieger. Academic Press, New York, p 799–817
Nikitaki Z, Hellweg CE, Georgakilas AG, Ravanat J-L (2015) Stress-induced DNA damage biomarkers: applications and limitations. Fronts in Chem 3. https://doi.org/10.3389/fchem.2015.00035l
Nikolaidis MG, Kyparos A, Spanou C, Paschalis V, Theodorou AA, Vrabas IS (2012) Redox biology of exercise: an integrative and comparative consideration of some overlooked issues. J Exp Biol 215(10):1615–1625. https://doi.org/10.1242/jeb.067470l
NSW National Parks and Wildlife Services. (2022). Nombinnie Nature Reserve. Retrieved 2/2/2022, 2022, from https://www.nationalparks.nsw.gov.au/visit-a-park/parks/nombinnie-nature-reserve.
Olsson M, Tobler M, Healey M, Perrin C, Wilson M, Memeg, Lund U, Lunds U (2012) A significant component of ageing (DNA damage) is reflected in fading breeding colors: an experimental test using innate antioxidant mimetics in Painted Dragon lizards. Evolution 66(8):2475–2483. https://doi.org/10.1111/j.1558-5646.2012.01617.xl
Ouyang JQ, Isaksson C, Schmidt C, Hutton P, Bonier F, Dominoni D (2018) A New Framework for Urban Ecology: An Integration of Proximate and Ultimate Responses to Anthropogenic Change. Integr Comparative Biol 58(5):915–928. https://doi.org/10.1093/icb/icy110l
Page MM, Stuart JA (2012) Activities of DNA base excision repair enzymes in liver and brain correlate with body mass, but not lifespan. Age 34(5):1195–1209
Peig J, Green AJ (2009) New perspectives for estimating body condition from mass/length data: the scaled mass index as an alternative method. Oikos 118(12):1883–1891. https://doi.org/10.1111/j.1600-0706.2009.17643.xl
Pérez-Rodríguez L, Romero-Haro AA, Sternalski A, Muriel J, Mougeot F, Gil D, Alonso-Alvarez C (2015) Measuring Oxidative Stress: The Confounding Effect of Lipid Concentration in Measures of Lipid Peroxidation. Physiol Biochem Zool 88(3):345–351. https://doi.org/10.1086/680688l
Peveling R, Demba SA (2003) Toxicity and pathogenicity of Metarhizium anisopliae var. acridum (Deuteromycotina, Hyphomycetes) and fipronil to the fringe-toed lizard Acanthodactylus dumerili (Squamata: Lacertidae). Environmental Toxicol Chem 22(7):1437–1447. https://doi.org/10.1002/etc.5620220704l
Raju KSR, Taneja I, Rashid M, Sonkar AK, Wahajuddin M, Singh SP (2016) DBS-platform for biomonitoring and toxicokinetics of toxicants: proof of concept using LC-MS/MS analysis of fipronil and its metabolites in blood. Scientific Rep 6(1):22447. https://doi.org/10.1038/srep22447l
Regan JD, Carrier WL, Samet C, Olla BL (1982) Photoreactivation in two closely related marine fishes having different longevities. Mech Ageing Development 18(1):59–66. https://doi.org/10.1016/0047-6374(82)90030-6l
Rej JE, Joyner TA (2018) Niche modeling for the genus Pogona (Squamata: Agamidae) in Australia: predicting past (late Quaternary) and future (2070) areas of suitable habitat. PeerJ 6:e6128. https://doi.org/10.7717/peerj.6128l
Reznick AZ, Witt E, Matsumoto M, Packer L (1992) Vitamin E inhibits protein oxidation in skeletal muscle of resting and exercised rats. Biochemical and Biophysical Res. Biochem Biophys Res Commun 189(2):801–806. https://doi.org/10.1016/0006-291X(92)92273-Zl
Ritchie DJ, Friesen CR (2022) Invited review: Thermal effects on oxidative stress in vertebrate ectotherms. Comparative Biochem Physiol Part A: Mol Integrative Physiol 263:111082
Robert Kylie A, Bronikowski AM (2010) Evolution of Senescence in Nature: Physiological Evolution in Populations of Garter Snake with Divergent Life Histories. Am Naturalist 175(2):147–159. https://doi.org/10.1086/649595l
Schwartz TS, Bronikowski AM (2011) Molecular stress pathways and the evolution of life histories in reptiles. Mechanisms of Life History Evolution: The Genetics and Physiology of Life History Traits and Trade-Offs, Oxford University Press: 0, United States of America
Schwartz TS, Bronikowski AM (2013) Dissecting molecular stress networks: identifying nodes of divergence between life-history phenotypes. Mol Ecol 22(3):739–756. https://doi.org/10.1111/j.1365-294X.2012.05750.xl
Selman C, Blount JD, Nussey DH, Speakman JR (2012) Oxidative damage, ageing, and life-history evolution: where now? Trends Ecol Evol 27(10):570–577. https://doi.org/10.1016/j.tree.2012.06.006l
Severo ES, Marins AT, Cerezer C, Costa D, Nunes M, Prestes OD, Zanella R, Loro VL (2020) Ecological risk of pesticide contamination in a Brazilian river located near a rural area: A study of biomarkers using zebrafish embryos. Ecotoxicol Environmental Saf 190:110071. https://doi.org/10.1016/j.ecoenv.2019.110071l
Shah ZU, Parveen S (2022) Oxidative, biochemical and histopathological alterations in fishes from pesticide contaminated river Ganga, India. Scientific Rep 12(1):3628. https://doi.org/10.1038/s41598-022-07506-8l
Shore RF, Douben PE (1994) Predicting ecotoxicological impacts of environmental contaminants on terrestrial small mammals. Revs of environmental contam and toxicol 49–89
Sies H (1997) Oxidative stress: Oxidants and antioxidants. Exp Physiol 82(2):291–295
Sies H, Berndt C, Jones DP (2017) Oxidative Stress. Annu Rev Biochem 86(1):715–748. https://doi.org/10.1146/annurev-biochem-061516-045037l
Sies H, Stahl W, Sevanian A (2005) Nutritional, diet and postprandial oxidative stress. J Nutr 135(5):969–972
Silva JM, Navoni JA, Freire EMX (2020) Lizards as model organisms to evaluate environmental contamination and biomonitoring. Environmental Monit Assess 192(7):454. https://doi.org/10.1007/s10661-020-08435-7l
Simbula G, Moltedo G, Catalano B, Martuccio G, Sebbio C, Onorati F, Stellati L, Bissattini AM, Vignoli L (2021) Biological responses in pesticide exposed lizards (Podarcis siculus). Ecotoxicology (London, England) 30(6):1017–1028. https://doi.org/10.1007/s10646-021-02440-3l
Simon-Delso N, Amaral-Rogers V, Belzunces LP, Bonmatin JM, Chagnon M, Downs C, Furlan L, Gibbons DW, Giorio C, Girolami V, Goulson D, Kreutzweiser DP, Krupke CH, Liess M, Long E, McField M, Mineau P, Mitchell EAD, Morrissey CA, Noome DA, Pisa L, Settele J, Stark JD, Tapparo A, Van Dyck H, Van Praagh J, Van der Sluijs JP, Whitehorn PR, Wiemers M (2015) Systemic insecticides (neonicotinoids and fipronil): trends, uses, mode of action and metabolites. Environmental Sci Pollut Res 22(1):5–34. https://doi.org/10.1007/s11356-014-3470-yl
Slaninova A, Smutna M, Modra H, Svobodova Z (2009) A review: Oxidative stress in fish induced by pesticides. Neuro Endocrinol Lett 30(Suppl 1):2–12
Sparling DW, Linder G, Bishop CA, Krest S (2010) Ecotoxicology of Amphibians Reptiles. CRC Press, United States of America
Speakman JR, Blount JD, Bronikowski AM, Buffenstein R, Isaksson C, Kirkwood TBL, Monaghan P, Ozanne SE, Beaulieu M, Briga M, Carr SK, Christensen LL, Cochemé HM, Cram DL, Dantzer B, Harper JM, Jurk D, King A, Noguera JC, Salin K, Sild E, Simons MJP, Smith S, Stier A, Tobler M, Vitikainen E, Peaker M, Selman C (2015) Oxidative stress and life histories: unresolved issues and current needs. Ecol Evol 5(24):5745–5757. https://doi.org/10.1002/ece3.1790l
Stier A, Romestaing C, Schull Q, Lefol E, Robin J-P, Roussel D, Bize P (2017) How to measure mitochondrial function in birds using red blood cells: a case study in the king penguin and perspectives in ecology and evolution. Methods Ecol Evol 8(10):1172–1182. https://doi.org/10.1111/2041-210X.12724l
Story P, Cox M (2001) Review of the effects of organophosphorus and carbamate insecticides on vertebrates. Are there implications for locust management in Australia? Wildl Res (East Melbourne) 28(2):179–193. https://doi.org/10.1071/WR99060l
Story P, French K, Maute K, Warden AC, Dojchinov G, Hose GC (2021) A method for topical dosing of invertebrates with pesticide for use in feeding experiments. Ecotoxicology 30:381–386. https://doi.org/10.1007/s10646-020-02324-yl
Story PG, Mineau P, Mullié WC (2013) Insecticide residues in Australian plague locusts (Chortoicetes terminifera Walker) after ultra-low volume aerial application of the organophosphorus insecticide fenitrothion. Environ Toxicol Chem 32(12):2792–2799. https://doi.org/10.1002/etc.2366l
Szabo J (2005) Avian-locust interactions in eastern Australia and the exposure of birds to locust control pesticides. Doctor of Philosophy, Texas Tech University, United States of America
Thompson HJ (2004) DNA oxidation products, antioxidant status, and cancer prevention. J Nutr 134(11):3186S–3187S
Thompson HM (1999) Esterases as Markers of Exposure to Organophosphates and Carbamates. Ecotoxicology 8(5):369–384. https://doi.org/10.1023/A:1008934505370l
Tingle CCD, Rother JA, Dewhurst CF, Lauer S, King WJ (2003) Fipronil: Environmental Fate, Ecotoxicology, and Human Health Concerns. Reviews of Environmental Contamination and Toxicology: Continuation of Residue Reviews. G. W. Ware. Springer New York, New York, NY, p 1–66
Todd BD, Willson JD, Gibbons JW (2010) The global status of reptiles and causes of their decline. Ecotoxicology of Amphibians and Reptiles, second ed. CRC Press, Boca Raton, USA
Valdivia PA, Zenteno-Savín T, Gardner SC, Alonso AA (2007) Basic oxidative stress metabolites in eastern Pacific green turtles (Chelonia mydas agassizii). Comparative Biochem Physiol Part C: Toxicol Pharmacol 146(1):111–117. https://doi.org/10.1016/j.cbpc.2006.06.008l
Valko M, Leibfritz D, Moncol J, Cronin MTD, Mazur M, Telser J (2007) Free radicals and antioxidants in normal physiological functions and human disease. Int J Biochem Cell Biol 39(1):44–84. https://doi.org/10.1016/j.biocel.2006.07.001l
Wallace J, Johnson B, Padilla J, Mather K (2010) Postprandial lipaemia, oxidative stress and endothelial function: a review. Int J Clin Pract 64(3):389–403
Weir SM, Suski JG, Salice CJ (2010) Ecological risk of anthropogenic pollutants to reptiles: Evaluating assumptions of sensitivity and exposure. Environmental Pollut 158(12):3596–3606. https://doi.org/10.1016/j.envpol.2010.08.011l
Winterbourn CC (2016) Revisiting the reactions of superoxide with glutathione and other thiols. Arch Biochem Biophys 595:68–71. https://doi.org/10.1016/j.abb.2015.11.028l
Woodhead AD, Setlow RB, Grist E (1980) DNA repair and longevity in three species of cold-blooded vertebrates. Exp Gerontol 15(4):301–304. https://doi.org/10.1016/0531-5565(80)90034-0l
Wu SM, Liu J-H, Shu L-H, Chen CH (2015) Anti-oxidative responses of zebrafish (Danio rerio) gill, liver and brain tissues upon acute cold shock. Comparative Biochem Physiol Part A: Mol Integrative Physiol 187:202–213. https://doi.org/10.1016/j.cbpa.2015.05.016l
Zayas MA, Rodríguez HA, Galoppo GH, Stoker C, Durando M, Luque EH, Muñoz-de-Toro M (2011) Hematology and blood biochemistry of young healthy broad-snouted Caimans (Caiman latirostris). J Herpetol 45(4):516–524. https://doi.org/10.1670/10-158.1l
Zuur A, Ieno EN, Walker N, Saveliev AA, Smith GM (2009) Mixed effects models and extensions in ecology with R. Springer Science & Business Media, United States of America
Acknowledgements
We would like to thank our field research assistants Bre Potter and Adam Bernich, who were integral in all aspects of fieldwork. Also integral in fieldwork were all the volunteers who helped throughout the field season. We would like to thank NPWS and Daniel Turacto for providing access to sites and accommodation at Hunthawang homestead and Nombinnie Nature Reserve. Additionally, we would like to thank Paul Story for integral ChE/AChE analysis on all samples as well as Greg Dojchinov and Andrew Warden (Black Mountain Laboratories, CSIRO) for performing fipronil blood residue assays. Finally, thank you to Chris Friesen for initial advice on oxidative stress assays.
Author contributions
ICK, KM, and KF conceived the ideas and designed the methodology. ICK and KM led field work and ICK collected all samples. ICK and MDB led oxidative stress assay lab work. ICK led writing of manuscript and performed all data analysis. All authors contributed to editing the manuscript, KF, and KM reviewed and supervised.
Funding
This work was supported by an ARC Linkage grant (LP160100686), funding was provided by Australian Research Council (ARC) and the Australian Plague Locust Commission (APLC). The grant was awarded to the University of Wollongong in collaboration with Macquarie University. Additional funding was provided through the Holsworth Wildlife Research Endowment Grant awarded to Isabella Contador-Kelsall. Open Access funding enabled and organized by CAUL and its Member Institutions.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Ethics approval
All samples and protocols were performed and collected under Animal Ethics approval through the University of Wollongong under AE17/19 and permitted through NSW National Parks and Wildlife Service Scientific License (SL100109).
Supplementary information
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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Contador-Kelsall, I., Maute, K., de Beer, M. et al. Individual variation within wild populations of an arid-zone lizard dictates oxidative stress levels despite exposure to sublethal pesticides. Ecotoxicology 32, 470–486 (2023). https://doi.org/10.1007/s10646-023-02653-8
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10646-023-02653-8