Effect of different anaesthetic techniques on gene expression profiles in patients who underwent hip arthroplasty

Abstract

Objectives

To investigate the modulation of genes whose expression level is indicative of stress and toxicity following exposure to three anaesthesia techniques, general anaesthesia (GA), regional anaesthesia (RA), or integrated anaesthesia (IA).

Methods

Patients scheduled for hip arthroplasty receiving GA, RA and IA were enrolled at Rizzoli Orthopaedic Institute of Bologna, Italy and the expression of genes involved in toxicology were evaluated in peripheral blood mononuclear cells (PBMCs) collected before (T0), immediately after surgery (T1), and on the third day (T2) after surgery in association with biochemical parameters.

Results

All three anaesthesia methods proved safe and reliable in terms of pain relief and patient recovery. Gene ontology analysis revealed that GA and mainly IA were associated with deregulation of DNA repair system and stress-responsive genes, which was observed even after 3-days from anaesthesia. Conversely, RA was not associated with substantial changes in gene expression.

Conclusions

Based on the gene expression analysis, RA technique showed the smallest toxicological effect in hip arthroplasty.

Trial registration

ClinicalTrials.gov number NCT03585647.

About the Authors:

Renata Alleva

Roles Conceptualization, Data curation, Writing – original draft

* E-mail: [email protected]

Affiliation: Department of Biomedical and Neuromotor Sciences, University of Bologna, Bologna, Italy

ORCID logo http://orcid.org/0000-0001-7327-3833

Andrea Tognù

Roles Methodology

Affiliation: Department of Anaesthesia and Postoperative Intensive Care, Rizzoli Orthopaedic Institute, Bologna, Italy

Marco Tomasetti

Roles Formal analysis, Investigation, Methodology, Visualization

Affiliation: Department of Clinical and Molecular Sciences, Polytechnic University of Marche, Ancona, Italy

Maria Serena Benassi

Roles Formal analysis, Investigation

Affiliation: Laboratory of Experimental Oncology, Rizzoli Orthopaedic Institute, Bologna, Italy

Laura Pazzaglia

Roles Formal analysis, Methodology

Affiliation: Laboratory of Experimental Oncology, Rizzoli Orthopaedic Institute, Bologna, Italy

Hanna van Oven

Roles Investigation, Validation

Affiliation: Department of Anaesthesia and Postoperative Intensive Care, Rizzoli Orthopaedic Institute, Bologna, Italy

Ettore Viganò

Roles Investigation, Validation

Affiliation: Department of Anaesthesia and Postoperative Intensive Care, Rizzoli Orthopaedic Institute, Bologna, Italy

Nicola De Simone

Roles Formal analysis, Investigation

Affiliation: Department of Anaesthesia and Postoperative Intensive Care, Rizzoli Orthopaedic Institute, Bologna, Italy

Ilaria Pacini

Roles Formal analysis, Methodology

Affiliation: Department of Anaesthesia and Postoperative Intensive Care, Rizzoli Orthopaedic Institute, Bologna, Italy

Sandra Giannone

Roles Investigation, Validation

Affiliation: Department of Anaesthesia and Postoperative Intensive Care, Rizzoli Orthopaedic Institute, Bologna, Italy

Sanjin Gagic

Roles Data curation, Formal analysis

Affiliations Department of Biomedical and Neuromotor Sciences, University of Bologna, Bologna, Italy, Research Unit of Anaesthesia and Pain Therapy, Rizzoli Orthopaedic Institute, Bologna, Italy

Raffaele Borghi

Roles Investigation, Methodology

Affiliation: Department of Anaesthesia and Postoperative Intensive Care, Rizzoli Orthopaedic Institute, Bologna, Italy

Sara Picone

Roles Data curation, Validation

Affiliations Department of Biomedical and Neuromotor Sciences, University of Bologna, Bologna, Italy, Research Unit of Anaesthesia and Pain Therapy, Rizzoli Orthopaedic Institute, Bologna, Italy

Battista Borghi

Roles Project administration, Resources, Writing – review & editing

Affiliations Department of Biomedical and Neuromotor Sciences, University of Bologna, Bologna, Italy, Research Unit of Anaesthesia and Pain Therapy, Rizzoli Orthopaedic Institute, Bologna, Italy

Introduction

Arthroplasty surgery is nowadays the technique of choice to treat different articular conditions, such as knee or hip arthritis. In hip arthroplasty, regional anaesthesia (RA) represents one of the preferred options for anaesthesia, due to its safety and convenience for both patients and surgeons [1]. Peripheral nerve block is frequently used, reducing the dose and thus limiting the adverse effects of opioid and non-opioid analgesics [2]. Use of RA is known to prevent or attenuate excessive stress response during and after surgery. RA markedly decreased recovery time and postoperative pain, with similar levels of patient satisfaction when compared to GA [3,4].

Besides, multi-modal analgesia is increasingly used to reduce postoperative pain and accelerate rehabilitation.3 Integrated anaesthesia (IA), which consists of lumbar plexus block plus spinal anaesthesia, integrated with GA, represents another option, and has been shown to reduce intra-operative anaesthetic drugs requirement and to provide better post-operative pain relief in shoulder arthroscopic surgeries [5]. As of today, the cytotoxic effects of different anaesthesia techniques have not been clearly elucidated. A few studies have evaluated the effect of anaesthesia on gene expression (transcriptome) in cells. Lowes et al [6], reported that isoflurane anaesthesia affects differential gene expression in the brain of anesthetized rats compared to control. Isoflurane alters several genes involved with neurotransmitter transport, signalling and cellular structure [7], as well as genes involved in drug metabolism and clock gene expression [8].

In this study, we analysed differentially regulated genes involved in oxidative stress and toxicology in peripheral blood mononuclear cells (PBMCs) of patients who underwent arthroplasty under three different anaesthetic methods. We hypothesized that anaesthesia procedures trigger toxicity, thus inducing changes in the mRNA profile. The results may provide a more profound understanding of the molecular mechanism of anaesthesia and in overcoming the adverse effects arising from their use.

Methods

Ethics statement

The study was carried out according to the Helsinki Declaration and written informed consent was obtained from all participants. Ethical approval for this study (Ethical Committee N° 12630) was provided by the Regional Ethical Committee of the Rizzoli Orthopaedic Institute, via pupilli 1–40136, Bologna, Italy, (Chairperson Prof. Gian Paolo Salvioli) on 12 June 2014 (S1 File).

Study population

Patients undergoing elective hip arthroplasty were recruited and followed-up from September 2014 to November 2016 at the ‘Rizzoli’ Orthopaedic Institute, Bologna, Italy. The study design is shown in Fig 1, according to CONSORT Checklist (S2 File). The sample size was estimated according to Lee-Whitmore and G-Power F-test for ANOVA fixed effects, omnibus, one-way [9]. Taking in account a power of 0.8 and a Bonferroni correction with a significance of 0.05, the minimum sample size was established for 30 patients of each group. By considering a 10% of drop-out, the chosen sample size was 33 patients per group which leads to a total sample size of 99 patients.

[Figure omitted. See PDF.]

Fig 1. Overview of study design according to CONSORT format.

Accordantly to eligibility, 130 patients were enrolled. Of them, 31 patients were excluded from the study, then 99 patients were randomised in the three anaesthesia groups.

https://doi.org/10.1371/journal.pone.0219113.g001

By using a computer-generated randomization table, patients were randomly consecutively allocated to receive general (GA group, n = 30), regional (RA group, n = 30), or integrated (IA group, n = 30) anaesthesia. After arrival in the operating theatre, an 18-gauge intravenous cannula was placed at the forearm, and then all patients were pre-medicated with 0.05 mg/kg intravenous (IV) midazolam. A 5-mL/kg/hour intravenous infusion with ringer lactate was then started without prophylactic volume expansion. The duration of the procedure was 110 ± 47 min. in group GA, 108 ± 31 min. in group RA, and 107 ± 41 min. in group IA (p = 0.68). Patients with contraindication to spinal anaesthesia or lumbar catheter placement, as well as obese patients, with arterial hypertension not controlled by oral medication, severe pulmonary, cardiovascular, renal, hepatic, cerebrovascular, or psychiatric diseases were excluded from the study. Accounting for a 30% drop-out rate, a total of 130 patients were enrolled. The study protocol (S3 File) was registered at www.clinicaltrials.gov (NCT03585647) after enrolment of participants as the study later was defined as a Clinical Trial and all related clinical trials have been registered. All participants provided written informed consent based on documents approved by Institute Institutional Review Board.

Anaesthesia procedures

GA was induced by intravenous fentanyl (1 mcg/kg) and propofol (2 mg/kg), followed by vecuronium bromide (0.1 mg/kg) to facilitate tracheal intubation, then GA was maintained using a 50% air/oxygen mixture and sevoflurane. The end-tidal concentration of sevoflurane was adjusted to maintain heart rate and blood pressure values within 20% of baseline. Mechanical ventilation was regulated to maintain the end-tidal CO2 partial pressure ranging between 4.3 and 5.1 Kilopascal (kPa).

RA: included continuous lumbar plexus block, performed by or under supervision of an experienced operator using a nerve stimulator (Stimuplex-HNS 11, B. Braun Melsungen, Germany) and Braun Contiplex-Tuohy Continued Peripheral Nerve Block Set.

The landmark used to locate the needle insertion point was the soft tissue depression at the iliac crest prominence, the most anterior and cranial aspect of the crest. A total dose of 20 ml of 0.5% Levobupivacaine was administered at the time of catheter placement. Then, dural puncture was performed at the L3-L4 interspace using a 25-Gauge Whitacre spinal needle (Becton-Dickinson, New Jersey, USA) with the midline approach using 3 ml of 0.5% Levobupivacaine.

IA: First, the patients received regional anaesthesia (lumbar plexus block + spinal anaesthesia) as described above. Then GA was induced by propofol 1% and a laryngeal mask airway of appropriate size was inserted. GA and mechanical ventilation were maintained as described above.

Postoperative analgesia: patients in the RA and IA groups received a continuous lumbar plexus infusion of levobupivacaine 0.25% at 7 mL/h injected through an elastomeric pump, whereas the GA group was treated with a continuous infusion of 1 mg/h morphine after an intravenous bolus of 0.1 mg/kg. If such analgesic treatment was not sufficient to relieve the postoperative pain, a rescue dose of 30 mg ketorolac (Recordati, Italy) was available for the patients. The pain intensity level was evaluated by Numerical Rating Scale (NRS) values. The prediction of the operative risk has been evaluated by American Society of Anaesthesiologists (ASA) physical status classification [10].

Blood samples

Whole blood samples (10 mL) were obtained from all enrolled patients at three time points: early morning on the operation day (T0), immediately after surgery (T1) and third day (T2) after surgery. The samples were collected in heparin tubes and PBMCs were isolated as previously described for gene expression evaluation [11]. Briefly, blood samples (6 ml) were layered onto 4 ml of Lympholyte-H (Cederlane, Hornby, Ontario, Canada) and centrifuged at 1000 g (20°C, 30 min). After centrifugation, the cloudy layer was collected and placed in a 15 ml Falcon tube, filled with PBS, pH 7.4, and centrifuged at 1000 g (20°C, 15 min). After removing the supernatant, the pellet of PBMCs containing 80–90% of lymphocytes was collected and stored at -80°C for RNA extraction. Whole blood and serum obtained after centrifugation were used for the analysis haematological and biochemical, such as Glutamate Oxaloacetate Transaminase (GOT), Glutamate-Pyruvate Transaminase (GPT), Bilirubin (BIL), Creatinine (CREA), Creatine Phosphokinase (CPK), and Haemoglobin (HB).

Microarray analysis

Total RNA was extracted from PBMCs using PerfectPure RNA Kit (5Prime, Hamburg, Germany) according to the manufacturer’s instructions. The cDNA was synthesized using RT2-First Strand Kit (SA Biosciences, Frederick, MD, USA) according to the manufacturer’s instructions. Human Stress and Toxicity PathwayFinder RT2 Profiler PCR Expression Array (PAHS-003 SABiosciences) was used for gene expression profiling. The expression of 84 stress- and toxicity-related genes was assessed by qRT- PCR (Mastercycler EP Realplex, Eppendorf, Milano, Italy) using RT2 SYBR Green qPCR Master Mix (SABiosciences). A set of controls (6 housekeeping genes, 3 control genes for quality of the retro-transcription and 3 control genes for qPCR) were included in the array. Among the tested housekeeping genes, GAPDH showed smaller changes in its expression across different sample groups (differences in Ct values < 1). Gene-specific products were then normalized to GAPDH and expressed as fold change (2-ΔΔCt). The PCR array raw dataset are shown in theS4 File.

Gene ontology analysis

Gene ontology analysis with Benjamini and Adchberg false discovery rate (FDR) correction, (cut-off<0.05), was performed to identify the sets of genes grouped to biological process and molecular function (nodes) that are significantly different as expressed by the anaesthetic procedures. Of all genes being significantly expressed, the fold change versus preoperative control was determined. The data were analysed by Cytoscape software.

Quantitative RT-PCR analysis

Total RNA from PBMCs was obtained using the RNeasy Mini Kit (Qiagen), according to the manufacturer's instructions. The concentration of RNA was determined spectrophotometrically at 260 nm (NanoDrop ND-1000, Thermo Scientific, Wilmington, USA). The cDNA was synthesized using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems). The expression of selected genes was quantified using the TaqMan system (Applied Biosystems, Foster City, CA, USA), using the Mastercycler EP Realplex instruments (Eppendorf). The GAPDH was used for normalization and the results were expressed as ΔCT, and fold changes in relative mRNA expression were calculated using the equation 2-ΔΔCT.

Statistical analysis

Results were expressed as mean ± S.D. Normalization and microarray analysis were performed by RT2 Profiler PCR Array Analysis software version 3.5 (SABiosciences), and two-tailed Student’s t-test was used to compare gene expression at T0 versus T1 and T2 time points. Comparisons among groups of data (age and BMI) were made using one-way analysis of variance (ANOVA) with Tukey post-hoc analysis. The Chi-square test was used to compare categorical variables (gender and ASA). ANOVA repeated measure with the Sidak post-hoc test was used to evaluate the changes of biochemical parameters (GOT, GPT, BIL, CREA, CPK, HB) in each group over time (T0, T1, T2). Generalized linear model (GLM) multivariate regression analysis was used to model the changes (ΔT = T0-T2) of biochemical parameters (GOT, GPT, BIL, CREA, CPK, HB) using continuous variables (age and body mass index, BMI) as covariates, and categorical variables (anaesthetic methods, gender and ASA) as fixed factors. The data were analysed by the Statistical Package Social Sciences (version 19) software (SPSS, Chicago, IL, USA) and p-values less than 0.05 were considered significant.

Results

There were no significant differences in age, BMI, gender, clinic-pathologic parameters and ASA scores among GA, RA and IA groups. However, an increase of transaminases (GOT and GPT) was observed on the third postoperative day in patients underwent to GA as result of a hepatic damage (Table 1).

[Figure omitted. See PDF.]

Table 1. Demographic characteristics and biochemical parameters of the enrolled subjects.

https://doi.org/10.1371/journal.pone.0219113.t001

To uncover differentially expressed genes in relation to anaesthesia procedures, we performed a screening phase where 84 genes involved in cell stress and toxicity were evaluated in PBMCs from a sub-population of 9 patients. Patients undergoing to arthroplasty surgery under GA (n = 3, age = 61±6; M/F = 1/2; BMI = 27±3, ASAI/II = 0/3), receiving RA (n = 3, age = 63±5; M/F = 2/1; BMI = 29±1, ASAI/II = 1/2), or receiving IA (n = 3, age = 64±1; M/F = 1/2; BMI = 29±4, ASAI/II = 0/3), before (T0) and after operation (T1) and on postoperative day three (T2) were randomly selected for gene expression profiling. The complete list of genes used to identify deregulated pathways is provided on S1 Table. Differentially expressed genes within different anaesthesia techniques at T1 and T2 time points respect to T0 were identified as having a ± 2-fold expression change and a p-value less than 0.05. As shown in Fig 2A, GA modulates stress-responsive genes such as RAD50, a protein involved in DNA double-strand break repair, and the small heat shock protein (HSP) family. Most genes were down regulated; the fold changes of the deregulated genes at time points T1 and T2 are reported in Table 2. Gene ontology analysis (biological and molecular process) revealed significantly altered genes involved in age-dependent response to reactive oxygen species (ROS) and in the negative regulation of nuclease and endonuclease activity by affecting DNA binding (Fig 2B). Few genes were induced in the RA group, mostly involved in the age-dependent response to metabolic and oxidative stress. This is represented graphically in Fig 3A and 3B and fold changes showed in Table 3, where most of the genes were under-expressed or not affected by the anaesthesia. Conversely, induction of gene expression was observed in the IA group. The gene expression profile drastically changes at T2 time point with the 42% of up-regulated genes and 19% of down-regulated genes (Table 4). Gene ontology analysis identified nodes involved in the age-dependent response to metabolic decline and oxidative stress (Fig 4A and 4B). Fig 5A and 5B shows the Venn diagram of the distribution of the co-deregulated genes across the anaesthesia groups and clustering analysis at T2 time point.

[Figure omitted. See PDF.]

Fig 2. Gene expression analysis in patients receiving general anaesthesia (GA).

(A) Heatmap of significant gene expression in PBMCs of patients (n = 3) undergoing elective hip arthroplasty receiving GA, immediately after operation (T1) and on the third postoperative day (T2), adjusted at p-value less than 0.05. Genes with greater and lower abundance after anaesthesia (FC, fold change) are shown in red and green, respectively, with significance highlighted in yellow. Normalization and microarray analysis were performed by RT2 Profiler PCR Array Analysis software version 3.5 (SABiosciences). (B) Gene-gene interaction analysis of the significantly differentially expressed genes showing the networks of deregulated pathways at T2 time point.

https://doi.org/10.1371/journal.pone.0219113.g002

[Figure omitted. See PDF.]

Fig 3. Gene expression analysis in patients receiving regional anaesthesia (RA).

(A) Heatmap of significant gene expression in PBMCs of patients (n = 3) undergoing elective hip arthroplasty receiving RA, immediately after operation (T1) and on the third postoperative day (T2), adjusted at p-value less than 0.05. Genes with greater and lower abundance after anaesthesia (FC, fold change) are shown in red and green, respectively, with significance highlighted in yellow. Normalization and microarray analysis were performed by RT2 Profiler PCR Array Analysis software version 3.5 (SABiosciences). (B) Gene-gene interaction analysis of the significantly differentially expressed genes showing the networks of deregulated pathways at T2 time point.

https://doi.org/10.1371/journal.pone.0219113.g003

[Figure omitted. See PDF.]

Fig 4. Gene expression analysis in patients receiving integrated anaesthesia (IA).

(A) Heatmap of significant gene expression in PBMCs of patients (n = 3) undergoing elective hip arthroplasty receiving IA, immediately after operation (T1) and on the third postoperative day (T2), adjusted at p-value less than 0.05. Genes with greater and lower abundance after anaesthesia (FC, fold change) are shown in red and green, respectively, with significance highlighted in yellow. Normalization and microarray analysis were performed by RT2 Profiler PCR Array Analysis software version 3.5 (SABiosciences). (B) Gene-gene interaction analysis of the significantly differentially expressed genes showing the networks of deregulated pathways at T2 time point.

https://doi.org/10.1371/journal.pone.0219113.g004

[Figure omitted. See PDF.]

Fig 5. Venn diagram and gene clustering.

(A) Venn diagram that shows differentially expressed genes that were shared among three anaesthesia methods. Cluster analysis of differentially expressed genes in general anaesthesia (Group-1), regional anaesthesia (Group-2), and integrated anaesthesia (Group-3) at third postoperative day (T2), respect to control Group (before anaesthesia, T0). (B) Gene clustering by their expression levels.

https://doi.org/10.1371/journal.pone.0219113.g005

[Figure omitted. See PDF.]

Table 2. Deregulated gene expression in PBMCs of patients (n = 5) undergoing elective hip arthroplasty receiving general anaesthesia (GA), immediately after operation (T1) and on the third postoperative day (T2).

https://doi.org/10.1371/journal.pone.0219113.t002

[Figure omitted. See PDF.]

Table 3. Deregulated gene expression in PBMCs of patients (n = 5) undergoing elective hip arthroplasty receiving regional anaesthesia (RA), immediately after operation (T1) and on the third postoperative day (T2).

https://doi.org/10.1371/journal.pone.0219113.t003

[Figure omitted. See PDF.]

Table 4. Deregulated gene expression in PBMCs of patients (n = 5) undergoing elective hip arthroplasty receiving integrated anaesthesia (IA), immediately after operation (T1) and on the third postoperative day (T2).

https://doi.org/10.1371/journal.pone.0219113.t004

To evaluate the effects of anaesthesia methods on interactions between genes, gene ontology analysis of microarray data was performed. According to the gene ontology analysis, genes involved in two principal nodes, such as DNA repair (ERCC1, ERCC3, RAD50, XRCC1, XRCC2), and stress response (EGR1, HSF1, SOD, GPX) were selected and evaluated in all enrolled population by qRT-PCR.

As showed in Fig 6, there was agreement between the microarray and qRT-PCR differential expression levels, confirming that IA resulted in a postoperative increase of gene expression. GLM multivariate regression analysis performed to evaluate the association between hepatic and renal marker changes and the anaesthetic method adjusted by age, BMI, gender, and ASA showed that the anaesthesia methods, ASA and sex affected bilirubin changes (Table 5).

[Figure omitted. See PDF.]

Fig 6. Expression of significantly deregulated stress-responsive genes.

Expression of DNA repair proteins (ERCC1, ERCC3, RAD50, XRCC1, XRCC2), early growth response protein 1 (EGR1), heat shock factor 1 (HSF1), glutathione peroxidase (GPX), and superoxide dismutase (SOD2) in patients undergoing arthroplasty before anaesthesia (T0), immediately after operation (T1), and third postoperative day (T2). General anaesthesia (GA, n = 25), regional anaesthesia (RA, n = 25), integrated anaesthesia (IA, n = 25). The data shown are expressed as fold change at T1 and T2 with respect to T0 time points. The symbol ‘*’ denotes significant differences in data at T1 and T2 versus T0 time points, with symbol ‘°’ significance between data at T1 and T2 time points with p<0.05.

https://doi.org/10.1371/journal.pone.0219113.g006

[Figure omitted. See PDF.]

Table 5. GLM-multivariate regression analyses showing associations between biochemical parameters, anaesthesia and demographic variables.

https://doi.org/10.1371/journal.pone.0219113.t005

Discussion

In the present study, we evaluated the influence of three anaesthetic techniques on gene expression in circulating cells at early time point and 3 days after anaesthesia. From the array analysis, we found that among the three anaesthesia methods, RA did not significantly affect gene expression. Deregulated genes were observed in patients under GA, and to a greater extent in subjects undergoing IA, even 3-days after anaesthesia induction, thus implying that multiple molecular processes remain altered over time. The affected genes are known to participate in a wide array of processes including cell cycle and DNA repair, signal transduction, transcriptional regulation, stress response proteins.

Propofol and sevoflurane are widely used in GA, and both have been reported to modulate gene expression [12–14]. Culley et al. first demonstrated that altered gene expression profile in rats induced by GA was associated with persistent changes in hippocampal gene expression, suggesting that recovery of the brain from anaesthesia was considerably slower than generally recognized [15]. Anaesthesia by isoflurane had marked effects on genes in the brain with differential regulation. Gene ontology analysis showed that some genes were functionally related to the anaesthesia and were involved with neurotransmitter release, transport and secretion [16]. We found that most deregulated genes were involved in DNA repair. The gene expression profile revealed an overall down regulation of DNA damage response genes after exposure to propofol, with significant under-expression of genes associated with nucleotide excision repair (ERCC1), double strand break repair (RAD50, RAD23A) and single strand break repair (XRCC1). Since up regulation of DNA repair genes occurs in response to DNA damage, a down regulation of DNA repair genes would be expected if a protective effect is induced. The protective effect of propofol was further supported by the down regulation of stress responsive genes, such as the inducible SOD2 and inflammatory cytokines/chemokines (TNF, IL6, CCLs). The protection of propofol against inflammation and oxidative stress has been previously reported [17–20]. Both propofol and sevoflurane attenuated the extent of hepatic ischemia/reperfusion (I/R) injury by inhibiting Nuclear factor kappa B (NFkB) activation and subsequent alterations in inflammatory cytokines [21,22].

On the other hand, increased DNA damage was observed in PBMCs of subjects who underwent GA: oxidized DNA bases occurred after 15 min of isoflurane exposure associated with an enhancement of DNA repair activity. However, most DNA damage was repaired on the first postoperative day [23]. The mutagen effect of anaesthetics was also observed in occupationally exposed subjects referred as population risk [24]. The isoflurane induced DNA damage as consequence of oxidative stress and inhibition of the repair of DNA damage through the p53 signalling pathway [25]. Besides, desflurane anaesthesia induces DNA strand breaks/alkali-labile sites on the day after minimally invasive surgery in healthy patients [26]. There is controversy over the genotoxic effects of volatile anaesthetics. However, a rise in serum liver enzymes was found in patients who underwent GA, supporting its hepatic toxicity.

No significant alteration in gene expression was found in RA. A toxic effect was observed for local anaesthetics. Lidocaine, bupivacaine, and ropivacaine treatment induced a significant decrease in viability with a concomitant increase in the number of apoptotic cells [27,28]. Mitochondrial DNA damage and decreased ATP and mitochondrial protein levels were found in cells exposed to local anaesthetics [23]. Several mechanisms have been proposed to explain toxicity of local anaesthetics to cells, including the blockade of potassium channels and mitochondrial injury [28]. Rather than single anaesthesia, we reported a synergic interaction between propofol used in GA and regional anaesthetics (IA) activating a set of genes associated with the repair of impaired proteins and structures including genes related to apoptosis, and extensive stress responses probably as result of induced oxidative stress.

Among them, early growth response protein 1 (EGR1) is known as redox-sensitive factor that plays a protective role when cells suffer starvation, ultraviolet light irradiation, hypoxia, and oxidative stress. ROS induce nuclear translocation of apurinic/apyrimidinic endonuclease 1 (APE1), which in turn induces DNA binding of transcriptional regulators such as EGR1 [29]. A number of genes are regulated directly by EGR1, which are highly associated with growth, vascular cell proliferation, cell survival programs and apoptosis [30]. It has been reported that EGR1 binding activity promotes tumour progression or atherosclerosis [31].

Increased expression of the heat shock factor 1 (HSF1) was observed at day 3 of IA induction. HSF1 enable the cell to adapt to various forms of oxidative, electrophilic, thermal, and inflammatory stress by orchestrating elaborate transcriptional programs termed the heat shock response (HSR). The heat-induced genes included many genes known to encode heat shock proteins (HSPs) and other proteins involved in protein folding and degradative pathways [32].

Activation of human heat shock genes was associated by increased antioxidant defence grid that relies on endogenous enzymatic antioxidants (SOD2 and GPX). The induction of the oxidative stress response pathway would suggest a potentially deleterious effect, which affects all forms of major surgery including cardiac surgery, general surgery, trauma surgery, plastic surgery and orthopaedic surgery [33]. This is seen more frequently in arthroplasty, which involves open surgical approaches and, compared to closed surgery, lead to a higher degree of soft tissue injury. The inflammation associated with oxidative stress has been shown to contribute to postoperative complications and to slow down recovery from surgery. In this context, it was reported that patients scheduled for total hip replacement receiving RA showed lower rate of complications and faster circulating cells mass recovery compared with GA and IA [3]. Therefore, anaesthesia procedures may affect postoperative recovery by modulating stress response and RA over the other anaesthesia techniques causes less surgical stress response in terms of changes in stress gene expression. Gene enrichment analysis was performed for regulated genes from microarray; further study using transcriptome sequencing (RNA-seq) technology will quantify overall expression levels and the degree of changes.

Conclusions

Changes in gene expression in arthroplasty occur early following general anaesthesia, and to a greater extent when integrated with lumbar plexus block, spinal anaesthesia and represents mainly genes involved in stress responses. Surgery itself results in expression of genes whose modulation may prove to be beneficial for acceleration of repair process [34]. Anaesthesia procedure may alter surgery-induced gene expression, thus affecting tissue repair. RA has advantages over the other anaesthesia techniques in terms of changes in stress gene expression.

Supporting information

[Figure omitted. See PDF.]

S1 File. Raw dataset.

https://doi.org/10.1371/journal.pone.0219113.s001

(PDF)

S2 File. Protocol.

https://doi.org/10.1371/journal.pone.0219113.s002

(DOC)

S3 File. CONSORT checklist.

https://doi.org/10.1371/journal.pone.0219113.s003

(DOCX)

S4 File. Study protocol.

https://doi.org/10.1371/journal.pone.0219113.s004

(XLS)

S1 Table.

https://doi.org/10.1371/journal.pone.0219113.s005

(DOCX)

Citation: Alleva R, Tognù A, Tomasetti M, Benassi MS, Pazzaglia L, van Oven H, et al. (2019) Effect of different anaesthetic techniques on gene expression profiles in patients who underwent hip arthroplasty. PLoS ONE 14(7): e0219113. https://doi.org/10.1371/journal.pone.0219113

References

1. Barroso Rosa S, James D, Matthews BD. Is knee arthroscopy under local anaesthetic a patient-friendly technique? A prospective controlled trial. Eur J Orthop Surg Traumatol. 2016; 26: 633–638. pmid:27324262

2. Kratz T, Dette F, Schmitt J, Wiesmann T, Wulf H, Zoremba M. Impact of regional femoral nerve block during general anaesthesia for hip arthroplasty on blood pressure, heart rate and pain control: A randomized controlled study. Technol Health Care. 2015; 23: 313–22. pmid:25669214

3. Borghi B, Casati A, Iuorio S, Celleno D, Michael M, Serafini PL, et al. Effect of different anaesthesia techniques on red blood cell endogenous recovery in hip arthroplasty. J Clin Anesth. 2005; 17: 96–101. pmid:15809124

4. Liszka H, Gądek A. Pre-emptive Local Anaesthesia in Ankle Arthroscopy. Foot Ankle Int. 2016; 37: 1326–1332. pmid:27621264

5. Kumara AB, Gogia AR, Bajaj JK, Agarwal N. Clinical evaluation of post-operative analgesia comparing suprascapular nerve block and interscalene brachial plexus block in patients undergoing shoulder arthroscopic surgery. J Clin Orthop Trauma. 2016;7: 34–39. pmid:26908974

6. Lowes DA, Galley HF, Moura AP, Webster NR. Brief isoflurane anaesthesia affects differential gene expression, gene ontology and gene networks in rat brain. Behav Brain Res. 2017; 317: 453–460. pmid:27659556

7. Pan JZ, Wei H, Hecker JG, Tobias JW, Eckenhoff RG, Eckenhoff MF. Rat brain DNA transcript profile of halothane and isoflurane exposure. Pharmacogenet Genomics. 2006; 16: 171–182. pmid:16495776

8. Sakamoto A, Imai J, Nishikawa A, Honma R, Ito E, Yanagisawa Y, et al. Influence of inhalation anaesthesia assessed by comprehensive gene expression profiling. Gene. 2005;356: 39–48. pmid:15967596

9. Lee ML, Whitmore GA. Power and sample size for DNA microarray studies. Stat Med. 2002;21: 3543–3570. pmid:12436455

10. Chen LH, Liang J, Chen MC, Wu CC, Cheng HS, Wang HH, et al. The relationship between preoperative American Society of Anesthesiologists Physical Status Classification scores and functional recovery following hip-fracture surgery. BMC Musculoskelet Disord. 2017;18: 410. pmid:29017476

11. Tomasetti M, Alleva R, Borghi B, Collins AR. In vivo supplementation with coenzyme Q10 enhances the recovery of human lymphocytes from oxidative DNA damage. FASEB J. 2001;15: 1425–1427. pmid:11387245

12. Yu D, Huang LJ, Chen NM. Anesthetic Propofol-Induced Gene Expression Changes in Patients Undergoing Coronary Artery Bypass Graft Surgery Based on Dynamical Differential Coexpression Network Analysis. Comput Math Methods Med. 2016;2016: 7097612.

13. Yamashita K, Matsumoto H, Saito F, Takeyoshi M. Differences in gene expression profiles in liver caused by different types of anaesthesia: cases of CO2-O2 and isoflurane. J Toxicol Sci. 2015;40: 829–836. pmid:26558464

14. Edmands SD, Ladow E, Hall AC. Microarray Analyses of Genes Regulated by Isoflurane Anaesthesia In Vivo: A Novel Approach to Identifying Potential Preconditioning Mechanisms. Anesth Analg. 2013;116: 589–595. pmid:23400992

15. Culley DJ, Yukhananov RY, Xie Z, Gali RR, Tanzi RE, Crosby G. Altered hippocampal gene expression 2 days after general anaesthesia in rats. Eur J Pharmacol. 2006;549: 71–78. pmid:16979155

16. Lowes DA, Galley HF, Moura AP, Webster NR. Brief isoflurane anaesthesia affects differential gene expression, gene ontology and gene networks in rat brain. Behav Brain Res. 2017;317: 453–460. pmid:27659556

17. Xu Z, Lu Y, Wang J, Ding X, Chen J, Miao C. The protective effect of propofol against TNF-α-induced apoptosis was mediated via inhibiting iNOS/NO production and maintaining intracellular Ca2+ homeostasis in mouse hippocampal HT22 cells. Biomed Pharmacother. 2017;91: 664–672. pmid:28499237

18. Yu X, Sun X, Zhao M, Hou Y, Hou Y, Li J, et al. Propofol attenuates myocardial ischemia reperfusion injury partly through inhibition of resident cardiac mast cell activation. Int Immunopharmacol. 2018;54: 267–274. pmid:29172064

19. Yoon JY, Jeon HO, Kim EJ, Kim CH, Yoon JU, Park BS et al. Propofol protects human keratinocytes from oxidative stress via autophagy expression. J Dent Anesth Pain Med. 2017;17: 21–28. pmid:28879325

20. Szyfter K, Stachecki I, Kostrzewska-Poczekaj M, Szaumkessel M, Szyfter-Harris J, Sobczyński P. Exposure to volatile anaesthetics is not followed by a massive induction of single-strand DNA breaks in operation theatre personnel. J Appl Genet. 2016;57: 343–348. pmid:26685861

21. Xu Z, Yu J, Wu J, Qi F, Wang H, Wang Z, et al. The Effects of Two Anesthetics, Propofol and Sevoflurane, on Liver Ischemia/Reperfusion Injury. Cell Physiol Biochem. 2016;38: 1631–1642. pmid:27119513

22. Zhang Y, Chen Z, Feng N, Tang J, Zhao X, Liu C et al. Protective effect of propofol preconditioning on ischemia-reperfusion injury in human hepatocyte. Thorac Dis. 2017;9: 702–710.

23. Alleva R, Tomasetti M, Solenghi MD, Stagni F, Gamberini F, Bassi A, et al. Lymphocyte DNA damage precedes DNA repair or cell death after orthopaedic surgery under general anaesthesia. Mutagenesis. 2003;18: 423–428. pmid:12960410

24. Souza KM, Braz LG, Nogueira FR, Souza MB, Bincoleto LF, Aun AG, et al. Occupational exposure to anesthetics leads to genomic instability, cytotoxicity and proliferative changes. Mutat Res. 2016;791–792: 42–48. pmid:27639372

25. Ni C, Li C, Dong Y, Guo X, Zhang Y, Xie Z. Anesthetic Isoflurane Induces DNA Damage Through Oxidative Stress and p53 Pathway. Mol Neurobiol. 2017;54: 3591–3605. pmid:27194299

26. Nogueira FR, Braz LG, de Andrade LR, de Carvalho AL, Vane LA, Módolo NS, et al. Evaluation of genotoxicity of general anaesthesia maintained with desflurane in patients under minor surgery. Environ Mol Mutagen. 2016;57: 312–316. pmid:27062561

27. Grishko V, Xu M, Wilson G, Pearsall AW 4th. Apoptosis and mitochondrial dysfunction in human chondrocytes following exposure to lidocaine, bupivacaine, and ropivacaine. J Bone Joint Surg Am. 2010;92: 609–618. pmid:20194319

28. Dragoo JL, Braun HJ, Kim HJ, Phan HD, Golish SR. The in vitro chondrotoxicity of single-dose local anesthetics. Am J Sports Med. 2012;40: 794–799. pmid:22287644

29. Pines A, Bivi N, Romanello M, Damante G, Kelley MR, Adamson ED, et al. Cross-regulation between Egr-1 and ape/ref-1 during early response to oxidative stress in the human osteoblastic hobit cell line: Evidence for an autoregulatory loop. Free Radic Res. 2005;39: 269–281. pmid:15788231

30. Pagel JI, Deindl E. Disease Progression Mediated by Egr-1 Associated Signaling in Response to Oxidative Stress. Int J Mol Sci. 2012;13: 13104–13117. pmid:23202940

31. Hasan RN, Phukan S, Harada S. Differential regulation of early growth response gene-1 expression by insulin and glucose in vascular endothelial cells. Arterioscler Thromb Vasc Biol. 2003;23: 988–993. pmid:12689920

32. Dayalan Naidu S, Kostov RV, Dinkova-Kostova AT. Transcription factors Hsf1 and Nrf2 engage in crosstalk for cytoprotection. Trends in Pharmacol Sci. 2015; 36: 6–14.

33. Rosenfeldt F, Wilson M, Lee G, Kure C, Ou R, Braun L, et al. Oxidative stress in surgery in an ageing population: pathophysiology and therapy. Exp Gerontol. 2013;48: 45–54. pmid:22465624

34. Bech NH, Hulst AH, Spuijbroek JA, van Leuken LL, Haverkamp D. Perioperative pain management in hip arthroscopy; what options are there? J Hip Preserv Surg. 2016;3: 181–189. pmid:27583156