Gene expression profiles in skeletal muscle after gene electrotransfer
© Hojman et al; licensee BioMed Central Ltd. 2007
Received: 09 May 2007
Accepted: 29 June 2007
Published: 29 June 2007
Gene transfer by electroporation (DNA electrotransfer) to muscle results in high level long term transgenic expression, showing great promise for treatment of e.g. protein deficiency syndromes. However little is known about the effects of DNA electrotransfer on muscle fibres. We have therefore investigated transcriptional changes through gene expression profile analyses, morphological changes by histological analysis, and physiological changes by force generation measurements. DNA electrotransfer was obtained using a combination of a short high voltage pulse (HV, 1000 V/cm, 100 μs) followed by a long low voltage pulse (LV, 100 V/cm, 400 ms); a pulse combination optimised for efficient and safe gene transfer. Muscles were transfected with green fluorescent protein (GFP) and excised at 4 hours, 48 hours or 3 weeks after treatment.
Differentially expressed genes were investigated by microarray analysis, and descriptive statistics were performed to evaluate the effects of 1) electroporation, 2) DNA injection, and 3) time after treatment. The biological significance of the results was assessed by gene annotation and supervised cluster analysis.
Generally, electroporation caused down-regulation of structural proteins e.g. sarcospan and catalytic enzymes. Injection of DNA induced down-regulation of intracellular transport proteins e.g. sentrin. The effects on muscle fibres were transient as the expression profiles 3 weeks after treatment were closely related with the control muscles. Most interestingly, no changes in the expression of proteins involved in inflammatory responses or muscle regeneration was detected, indicating limited muscle damage and regeneration. Histological analysis revealed structural changes with loss of cell integrity and striation pattern in some fibres after DNA+HV+LV treatment, while HV+LV pulses alone showed preservation of cell integrity. No difference in the force generation capacity was observed in the muscles 2 weeks after DNA electrotransfer.
The small and transient changes found in the gene expression profiles are of great importance, as this demonstrates that DNA electrotransfer is safe with minor effects on the muscle host cells. These findings are essential for introducing the DNA electrotransfer to muscle for clinical use. Indeed the HV+LV pulse combination used has been optimised to ensure highly efficient and safe DNA electrotransfer.
In vivo gene transfer to skeletal muscle is a promising strategy for treatment of muscular disorders and for systemic delivery of therapeutic proteins. Transgene expression has been reported following intramuscular injection of naked plasmid DNA , yet the clinical use is limited due to low efficiency and large variation. By combining intramuscular plasmid injection with local application of electric pulses a 200-fold increase in transfection efficiency with high reproducibility has been achieved in vivo [2, 3]. This proves that DNA electrotransfer can be an efficient and feasible way of introducing genes into tissues. DNA electrotransfer allows high production of systemic delivered proteins e.g. erythropoietin [4, 5] and cytokines [6, 7] with expression detected more than a year after treatment in rodents [8, 9].
The high level of gene expression in muscle after DNA electrotransfer is a consequence of plasmid distribution, membrane permeabilisation and plasmid electrophoresis [10–12]. Recently efforts have been made to optimise the electric pulses for gene transfer, resulting in a combination of a short high voltage (HV) pulse for membrane permeabilisation, followed by a long low voltage (LV) pulse for electrophoretic transport of plasmid towards the cell membrane. In skeletal muscle and skin the HV + LV pulse combination has led to increased marker gene expression .
DNA electrotransfer is moving rapidly towards clinical use. In fact, electroporation in combination with chemotherapy (electrochemotherapy) has been used in the clinic for several years now [14, 15], showing encouraging results for anti-tumour treatment with good tolerability for the patients [16–18]. However the biological response to electroporation and DNA electrotransfer in vivo remains poorly defined. Rubenstrunk et al studied the expression of stress related genes after DNA electrotransfer, and found no significant variation between treated and non-treated muscles . Furthermore, studies have shown that transient changes in force generation, muscular ion content and ATP levels were observed following DNA electrotransfer . Yet no thorough description of the transcriptional changes caused by DNA electrotransfer has been reported.
In this first comprehensive microarray analysis covering the entire murine genome, we examine the transcriptional changes underlying the muscular response to DNA electrotransfer by evaluating the gene expression profiles of mice tibialis cranialis muscles 4 hrs, 48 hrs and 3 weeks after DNA electrotransfer using oligonucleotide microarrays and quantitative PCR (Q-PCR). Furthermore we performed histology and physiological tests such as force generation measurements and reflex and motor testing to support our findings.
The effect of DNA electrotransfer on gene expression profiles
For precise description of the transcriptional changes following DNA electrotransfer, total RNA was extracted and pooled from 4 muscles excised 4 hrs, 48 hrs and 3 weeks after treatment and gene expression profile analysis was performed. Statistical analysis was employed by comparing electroporated (EP) and non-EP groups, groups either injected or not with plasmid, or groups at the different time points.
By two-way ANOVA testing with a significance level of p = 0.001, 29 genes were found to be differentially expressed between EP and non-EP groups, 38 genes were differentially expressed in groups with or without intramuscular plasmid injection, and 145 genes were differentially expressed across time. To further filter the gene lists, unpaired t-test with 50 permutations accepting a false discovery rate (FDR) less than 5% and significance levels of p = 0.05 were performed. In the EP vs. non-EP groups, the unpaired t-test revealed that 27 of the 29 genes were significantly differentially expressed (p = 0.05) of which 7 were up-regulated and 20 down-regulated. Unpaired t-test among the groups injected with plasmid or not, revealed that 36 genes were significantly differentially expressed (p = 0.05) with 9 genes up-regulated and 27 down-regulated. Unpaired t-test among genes differing across time resulted in 130 genes with a false discovery rate of 8.5%, therefore the significance level was changed to p ≤ 0.005, resulting in 45 differentially expressed genes and a FDR of 0%.
Descriptive statistics and unsupervised clustering
Electroporated versus non-electroporated muscles
Summary of transcriptional changes following gene transfer by electroporation
phosphoenolpuryvate carboxykinase 2
protein phosphatase 6
microtubule-associated protein 7
junction adhesion molecule 4
Forkhead box P1
F-box and leucine-rich repeat protein
ring finger protein 184
splicing factor (SRp75)
Transducer of ERBB2
eukaryotic translation initiation factor
host cell factor C1
CCR4-NOT transcription complex
cytoplasmic polyadenylation element
FK506 binding protein 1a
ADP-ribosylation factor-like 6 and 10
ubiquitin-like 1 (sentrin)
nuclear transport factor 2
Chloride channel Ka
Kruppel-like factor 15
K+ voltage gated channel, Shab related
Solute carrier family 6
ATP-binding cassette (ABC1)
Furthermore expression of genes involved in the metabolism e.g. phosphoenolpuryvate carboxykinase (BB024477) and dipeptidase (NM007876) were also down-regulated following electroporation, indicating that the general catabolism might be depressed. Ephrin (D38146), which is involved in neural development and cell-cell interaction among vascular endothelial cells were found to be up-regulated in all EP groups.
Plasmid vs. no plasmid injection
GO annotation showed that the expression of genes involved in intracellular trafficking e.g. sentrin (ubiquitin-like 1, NM016682), nuclear transport factor 2 (AA920031) and ADP-ribosylation factor-like 10C (NM026011) were down-regulated in the plasmid injected groups. In addition polymerase II (BC004810) was also down-regulated in plasmid injected groups, indicating the transcription also was depressed.
Furthermore the ER Ca2+ binding protein reticulocalbin (NM011992) and Ca2+ release modulator FK506 binding protein (AF483488) were both down-regulated following plasmid injection, which might be important for controlling the Ca2+ homeostasis. Indeed we have previously shown that DNA injection especially in combination with electroporation affects the intracellular level of Ca2+ (manuscript in preparation).
At 4 hrs the chloride channel Ka (C76734) is up-regulated in all groups compared to all other time point. At this initial time point after treatment and membrane permeabilisation, the ion homeostasis across the membrane might still be disturbed, stimulating the expression of the Cl-channel. At 48 hrs and in the EGT-treated group at 4 hrs, the transducer of ERBB2 (AV174616), was up-regulated. Up-regulation of this inhibitor of proliferation might indicate disturbances in the muscles, which limit the proliferation potential transiently.
Compared to the control groups, the Cl- channel activator Kruppel-like factor 15 (BC013486) and the K+ voltage-gated channel (NM011862) were down-regulated at 4 hrs, 48 hrs and 3 weeks, indicating the ion homeostasis might be disturbed longer than the initial phase. Similarly the solute carrier family 6 protein (bg069516), which transports creatine into the muscle cells and is essential for energy consumption, was also down-regulated at all time points compared to the control groups.
At 4 and 48 hrs the transcription activator host cell factor C1 (bg864583), the regulator of polymerase II CCR4-not transcription complex (BM202561) and the polyadenylation element binding protein (NM026252) were down-regulated compared to the control and 3 weeks groups; indicating that the protein synthesis apparatus might be negatively affected within the first 48 hrs after treatment.
Q-PCR validation of differentially expressed genes
Validation of differential expressed genes after DNA electrotransfer
HV + LV
DNA + HV + LV
FK506 binding pt.
Solute carrier family
Nuc. transport factor
Histological changes in tibialis cranialis muscles after DNA electrotransfer
No loss in muscle performance after DNA electrotransfer
Additional file 4: pre group 1 run. Representative film of running pattern in mice in group 1 (control) before treatment. (MOV 10 MB)
Additional file 5: pre group 1 seize. Representative film of seizing pattern in mice in group 1 (control) before treatment. (MOV 4 MB)
Additional file 6: pre group 2 run. Representative film of running pattern in mice in group 2 (DNA alone) before treatment. (MOV 10 MB)
Additional file 7: pre group 2 seize. Representative film of seizing pattern in mice in group 2 (DNA alone) before treatment. (MOV 5 MB)
Additional file 8: pre group 3 run. Representative film of running pattern in mice in group 3 (HV + LV) before treatment. (MOV 9 MB)
Additional file 9: pre group 3 seize. Representative film of seizing pattern in mice in group 3 (HV + LV) before treatment. (MOV 5 MB)
Additional file 10: pre group 4 run. Representative film of running pattern in mice in group 4 (DNA + HV + LV) before gene electrotransfer. (MOV 10 MB)
Additional file 11: pre group 4 seize. Representative film of seizing pattern in mice in group 4 (DNA + HV + LV) before gene electrotransfer. (MOV 6 MB)
Additional file 12: 48 H group 1 run. Representative film of running pattern in mice in group 1 (control) 48 hours after treatment. (MOV 10 MB)
Additional file 13: 48 H group 1 seizing. Representative film of seizing pattern in mice in group 1 (control) 48 hours after treatment. (MOV 4 MB)
Additional file 14: 48 H group 2 run. Representative film of running pattern in mice in group 2 (DNA alone) 48 hours after treatment. (MOV 10 MB)
Additional file 15: 48 H group 2 seize. Representative film of seizing pattern in mice in group 2 (DNA alone) 48 hours after treatment. (MOV 3 MB)
Additional file 16: 48 H group 3 run. Representative film of running pattern in mice in group 3 (HV + LV) 48 hours after treatment. (MOV 10 MB)
Additional file 17: 48 H group 3 seize. Representative film of seizing pattern in mice in group 3 (HV + LV) 48 hours after treatment. (MOV 3 MB)
Additional file 18: 48 H group 4 run. Representative film of running pattern in mice in group 4 (DNA + HV + LV) 48 hours after gene electrotransfer. (MOV 10 MB)
Additional file 19: 48 H group 4 seize. Representative film of seizing pattern in mice in group 4 (DNA + HV + LV) 48 hours after gene electrotransfer. (MOV 3 MB)
DNA electrotransfer to skeletal muscle is highly efficient and holds great promise for future treatment of monogenetic disorders, protein deficiency syndromes, cancers, etc. For electric pulsing of muscle in vivo trains of long low field strength pulses have typically been employed [2, 21]. However studies have shown that the electric pulse-mediated gene transfer might consist of two components; a component, which permeabilises the membrane, and an electrophoretic component, which can facilitate the transport of plasmid DNA across the membrane [10, 12]. Using combination pulses, a short high voltage permeabilising pulse followed by a long low voltage electrophoretic pulse, less intense pulsing parameters are needed for obtaining the same expression levels, which reduces the discomfort from the pulsing regiment. Indeed in this study we confirm these observations as little changes in the gene expression profiles are observed. Moreover we find less histological changes than previously described after DNA electrotransfer with trains of long, low field strength pulses [22, 23]. In addition we find no reductions in muscle performance 2 weeks after treatment, indicating that no major muscle regeneration is required after DNA electrotransfer.
Array analysis, clustering and time aspect
Due to the expenses of the study, the replica number was low in this study, however this was accommodated by designing our statistical analysis as a multiple factor analysis. During the inferential statistics we performed an ANOVA test, which gives an efficient and comprehensive evaluation despite the low replicate number in this set-up . Moreover, permutation test and unpaired T-tests were performed using stringent level of significance to correctly identify genes that were differentially expressed. Furthermore extensive descriptive statistics were performed, which supported the findings of differentially expressed genes. Unsupervised cluster analyses showed clear separations between EP and non-EP muscles and DNA injected or non-injected muscles, indicating that the differentially transcribed genes in these clusters were highly characteristic for the given treatment. Indeed some of the changes we observed were in line with previous reported changes following DNA electrotransfer to cells. Not surprisingly the gene expression profiles at 4 and 48 hrs were relatively comparable, while the gene expression profiles after 3 weeks were more similar to the control muscles. This indicates that the observed changes after DNA electrotransfer are transient, followed by a return to the naïve state of the muscle. The most striking point about the expression profiles is that no genes involved in stress, immune responses, cell death or muscle regeneration were found to be differentially expressed, demonstrating that EGT is non-toxic. Indeed guided supervised cluster analysis looking at cell homeostasis, ion homeostasis and immune response ontologies showed no significant changes in the gene expression profiles.
DNA electrotransfer-induced transcriptional changes
Based on previous in vitro studies certain effects of electroporation and plasmid injection would have been expected. Electroporation of cells is known to cause cytoskeleton rearrangements and electro-conformational changes of ion channels, while DNA injection is associated in cytotoxic effects on the cells.
Changes in cytoskeleton-associated proteins
During the application of an electric field across cells, electro-compressive forces exert pressure on the membrane leading to elongation of cells. During this elongation phase the pressure on the membrane enlarges, resulting in the development of membrane destabilisations, which allow for entry of plasmids [25, 26]. Formation of microvilli and blebs associated with deformation of the cytoskeleton has been observed immediately after pulsing [27, 28]. In fact, the same phenomenon is observed in cyclic stretch-induced gene transfer, where cytoskeleton rearrangements are involved in mediating the gene transfer . In this study we found down-regulation of cytoskeleton proteins after electro-pulsing, corresponding to the reported deformations of the cytoskeleton after pulsing. Interestingly, studies have shown that muscle contractile function returns to normal within 4 hours after EGT with the HV+LV pulse combination (manuscript in preparation), indicating that the cytoskeleton rearrangements only have a minor effect on the muscle function.
Electro-conformational changes of ion channels
During electro-pulsing, voltage-dependent ion channels are particularly sensitive to the external electric field due to their intrinsic voltage sensors. High intensity electric fields cause structural changes, also known as electro-conformational changes, which especially affect the moveable charge particles in the channel [30, 31]. Chen et al have shown that the movement of charge particles in the delayed rectifier K+ channels was reduced after applying an external transmembrane potential of 400 mV [31, 32]. This applied membrane potential is considerably larger than the 150 – 200 mV required for lipid bilayer breakdown. At lower than 400 mV field intensities no current changes demonstrating conformational changes were observed . In our study we found down-regulation of the K+ voltage-gated channel. Although our applied electric field does not reach amplitude of 400 mV at the membrane level, weaker field strengths might be sufficient to induce electro-conformational changes of ion channels in muscle fibres, as these cells are larger. This suggests that electro-conformational changes might play a role in the regulation of voltage-gated channels after DNA electrotransfer in vivo.
Effects of plasmid injection
Intramuscular injection of plasmid can be associated with toxicity in the target tissue [33, 34]. Hartikka et al (2001) observed that muscle damage after DNA electrotransfer with super-coiled endotoxin-free plasmid was DNA dependent. When introducing 50 μg of plasmid into the quadriceps muscle, they observed lesions with necrotic myofibres and infiltration of inflammatory cells. After injection of vehicle and electric pulses delivery no such lesions were observed. We have used a low concentration of plasmid, which might explain why we observed little infiltration at our histological evaluation and no effects on the transcriptional level.
Intracellular plasmid DNA transport
Unexpectedly we found that proteins involved in intracellular transport were down-regulated following DNA injection. After entry into muscle fibres plasmid DNA must be transported to the nuclei for transcription, but as the cytoplasmic latticework is quite condense, it is unlikely that this transport occur by simple diffusion . In HeLa and COS cells it has been shown that plasmid DNA has a limited half-life (1 – 2 hours) in the cytosol due to degradation by cytoplasmic nucleases [36, 37], meaning that efficient transport to the nucleus is required. In fact studies have shown that the cytoplasmic movement of DNA towards the nucleus in cells is facilitated by the microtubular network . The down-regulated intracellular transport proteins are not directly involved in these transport pathways, but they belong to a group of proteins, which might be a limiting factor in efficient gene transfer. This could therefore be of importance for the transfection efficacy.
DNA electrotransfer by HV+LV pulses demonstrates small and transient adverse effects on the muscle. This is highly important for introduction in the clinical setting, where numerous applications of this technology are envisaged. Examples of interesting pre-clinical experiments with gene transfer to muscle tissue include correction of anaemia [39–42] or β-thalassemia [43, 44] by electrotransfer of EPO. Furthermore DNA electrotransfer can also be used to target the muscle itself, either for research purposes or for correction of muscle anomalies e.g. muscular dystrophies. In this regard, it is important that the gene transfer procedure itself imposes minimal effects on the target tissue.
In conclusion DNA electrotransfer using the HV+LV pulse combination is safe with only small changes in the expression of cytoskeletal and intracellular transport proteins, while most importantly no detrimental changes in muscle performance or gene expression profiles of proteins involved in cell death, inflammation or muscle regeneration was verifiable both shortly and long after treatment. Histology did, however, revealed some cells with morphological changes after 4 and 48 hrs and few cells containing central nuclei after 3 weeks.
Six to eight weeks old female C57Black/C mice (average weight 22 g, Taconic, Denmark) or 4 – 5 weeks old male or female Wistar rats (own breed, University of Aarhus) were kept under pathogen-free conditions at 22°C in a 14/10 hrs light/dark cycle with food and water ad libitum. All animal experiments were conducted in accordance with the recommendations of the European Convention for the Protection of Vertebrate Animals used for Experimentation and after permission from the Danish Animal Welfare Committee.
Plasmid DNA and in vivo DNA electrotransfer
The plasmid phGFP-S65T, encoding the green fluorescent protein (GFP), was obtained from Clontech (Palo Alto, CA, USA). DNA preparations were performed using Qiafilter Plasmid Maxiprep kits (Qiagen, Germany), and the concentration and quality of the plasmid preparations were controlled by spectrophotometry, accepting a 260/280 ratio of 1.8–1.9 and by gel electrophoresis.
Five micrograms of plasmid was dissolved in 20 μl PBS (Phosphate-buffered saline without MgCl2, the Hospital Pharmacy at Herlev, Denmark), and the plasmid solution was injected intramuscularly along the fibres into the tibialis cranialis muscle of anaesthetised C57Bl/6 mice using an insulin syringe. For rats, 5 μg/10 μl PBS were injected into the EDL muscle. Plate electrodes with 4-mm gap were fitted around the hind legs and an electric field was applied using the Cliniporator™ (IGEA, Italy). The muscle was pulsed with a combination of a high voltage pulse (100 μsec, 1000 V/cm) followed by a long low voltage pulse (400 ms, 100 V/cm) with 1 s lag between the pulses.
RNA extraction and gene expression analysis
Four hrs, 48 hrs and 3 weeks after treatment the mice were euthanized and the muscles were removed and placed in 1 ml solution D (guandinium thiocyanat, sodium citrate, sarcosyl and mercaptoethanol) on ice. Muscle extracts were prepared by homogenising the muscles with a rotor-stator homogeniser (S8N – 5G, IKA-werke, Germany), and RNA were extracted using the Chomczynski and Sacchi method  and digested with DNase (Promega, Ramcon, Denmark). The extracted RNA amount ranged from 5 – 33 μg with a mean concentration of 14.3 +/- 7.9 μg.
Due to low yield of total RNA, mRNA were amplified yielding antisense RNA (aRNA) using the MessageAMP II aRNA kit (Ambion, Europe). Briefly, total RNA were reverse transcribed into first strand cDNA using a T7 Oligo(dT) primer and M-MLV reverse transcriptase and incubated at 42°C for 2 hrs. Second strand cDNA was synthesized using DNA polymerase and purified on cDNA filter cartridge according to the manufacturers protocol. An in vitro transcription (IVT) reaction was set up to generate multiple copies of aRNA from the double stranded cDNA through incubation at 37°C for 14 hrs.
Affymetrix microarray analyses were preformed at the Affymetrix core facility at the Microarray Center at National University Hospital. In short 5 μg aRNA was used to synthesize double stranded cDNA using Superscript® Choice System (Invitrogen, Denmark) with an oligo(dT) primer containing a T7 RNA polymerase promoter (GenSet, Evry, France). The cDNA was used as template for an IVT reaction to generate biotin-labelled antisense cRNA (BioArray™ High Yield RNA Transcript Labelling Kit; Enzo Diagnostics, Farmingdale, NY, USA). After fragmentation at 94°C for 35 min in fragmentation buffer (40 mM Tris, 30 mM MgOAc, 10 mM KOAc), the labelled cRNA were hybridised for 16 hrs to Affymetrix MOE 430 2.0 arrays (Affymetrix, Santa Clara, CA, USA). The arrays were washed and stained with phycoerythrin conjugated streptavidin using the Affymetrix Fluidics Station® 400, and the arrays were scanned in the Affymetrix GeneArray® scanner to generate fluorescent images, as described in the Affymetrix GeneChip protocol.
Microarray data analysis and statistics
The cel image files (Affymetrix) were imported, pre-processed and analysed in DNA-Chip Analyser 2006. The array files were normalised using the non-linear invariant set normalisation method, choosing array (DNA-EP-48 HRS) as baseline. Normalisation curves were visually inspected and none of the arrays were critical. In order to interpret the probe signal, model-based gene expression indexes (MBEI)  gene expression modelling was calculated using the PM/MM (perfect match/miss match) method. During the analysis array outliers were detected. None of the arrays exceeded 1.4% of array outliers and the probe presence call was between 36.6% and 50.3%. The experiment conditions and replicate number were electroporation (n = 13 with nEP+ = 7 and nEP- = 6), DNA injection (n = 13 with nDNA+ = 7 and nDNA- = 6) and time interval (n = 13 whereas nctrl = 3, ntime 4 hrs = 4, ntime 48 hrs = 3 and ntime 3 weeks = 3)
Raw data from the experiment are available at the gene Expression Omnibus database (GEO) in a MIAME compliant format. They are submitted under the serial code GSE6686, and under the following sample codes GSM154212, GSM154215, GSM154216, GSM154218, GSM154222, GSM154226, GSM154230, GSM154233, GSM154234, GSM154235, GSM154236, GSM154237, GSM154238.
Inferential statistics were adapted to the data by two-way ANOVA testing with a significance level of p ≤ 0.001, which have been described to give high power in multiple factor experiments . To further filter the gene lists, unpaired t-test with 50 permutations accepting a false discovery rate (FDR) less than 5% and significant level of p = 0.05 were performed. An unpaired t-test among genes differing across time resulted in 130 genes with a false discovery rate of 8.5%, therefore the significance level was changed to p ≤ 0.005, with a result of 45 differentially expressed genes and a FDR of 0%.
Descriptive statistics were applied to the significant data to obtain biologically meaningful patterns.
Unsupervised two-way hierarchical clustering was performed with a distance measure of 1 – Pearson correlation coefficient and centroid linkage with a P-value threshold of 0.05 for significant sample clusters and significant gene clusters .
A sample correlation matrix (SCM) was built on the basis of the original hierarchical unsupervised clustering, providing valuable information of how efficiently the clusters correlated. Supervised analyses were performed by applying specific issue algorithms of interest (downloaded from NetAffx™ Analysis Center). The issues of interest were cell homeostasis, ion homeostasis, ion channel activity and immune responses. To obtain significant data the data was filtered by an unpaired t-test with a significance level of p ≤ 0.05 and variation across samples between 0.50 < Standard deviation/mean < 1000.00. No FDR criterion was applied because we wanted to see a possible existing pattern. Furthermore to find a linear combination of explanatory variables and coefficient estimates within the three groups of interest (EP, Plasmid, time) the difference between the distributions of linear discriminate scores (LDA-plot) for each group was maximized.
Finally the supervised analysis tool gene ontology (GO) was applied to the data. Data from the two-way hierarchical clustering was uploaded to two software programs EASE and NetAffx™ Analysis Center, which have the abilities to calculate the GO similarity. The comparison criterions within the groups were set to the highest similarity within the genes of same gene category.
Two-step real-time reverse transcription polymerase chain reaction (RT-PCR) was performed to validate samples of differentially regulated genes in the electrotransferred muscles as determined by interferential statistic analysis. The aRNA prepared for the array hybridisation was reverse transcribed using random p(dN)6 primers (Roche, Germany) and AMV reverse transcriptase (Roche, Germany). Eleven selected genes were amplified from the cDNA using the Brilliant SYBR Green Q-PCR kit (Stratagene, AH Diagnostics, Denmark). The primers were designed using Primer3 and PCR and detection was performed by Q-PCR (MX3000P, Stratagene, AH Diagnostics, Denmark). The RNA levels were normalised to β-microglobulin and down and up-regulation in expression compared to the control muscles were determined.
Evaluation of macroscopic and histological changes and GFP expression in tibialis cranialis muscles
Respectively 4 hrs, 48 hrs and 3 weeks after treatment the mice were euthanized and whole tibialis cranialis muscles were carefully isolated and evaluated for macroscopic changes. The isolated muscles were either evaluated for GFP expression by fluorescence stereo microscopy (Nikon SMZ1500, Japan) in blinded and randomised order or fixed in 1 ml formalin buffer and imbedded in paraffin blocks following standard procedures. Sections 3 – 5 μm thick were stained with hematoxylin and eosin (Mayer-Sour) or with alkaline phosphatase-conjugated CD4 antibodies followed by counterstaining with Mayer's hematoxylin.
Test of function
The treated animals were tested for their ability to seize a grid by placing the animal on the grid and lifting the hind part by the tail. Animals with normal muscle function were given the score 0, while animals with affected function were rewarded 1 and those, who did not seize the grid at all were scored 2. For running the animals were placed alone in an empty cage. Normal running were given the score 0, while animals with a limp were given 1 and animals that did not support on the leg were rewarded 2. A combined score from the two tests were then calculated.
Measurement of force
Force was measured as previously described in detail by Clausen & Everts . In brief, isolated EDL muscles were mounted vertically with their tendons intact on a force displacement transducer (Grass FT03, W. Warwick, RI, USA) for isometric contractions in thermostatically controlled (30°C) chambers containing standard Krebs-Ringer bicarbonate buffer (pH 7.3) with (in mM): 122.1 NaCl, 25.1 NaHCO3, 2.8 KCl, 1.2 KH2PO4, 1.2 MgSO4, 1.3 CaCl2, and 5 D-glucose and gassed continuously with a mixture of 95% O2 and 5% CO2. Direct electrical stimulation was delivered via platinum electrodes on either side of the mid-portion of the muscle. Muscle length was adjusted to optimal length during repeated stimulation with single pulses. Finally force was checked using short tetanic contractions induced by supramaximal 1 ms pulses of 10 V at 90 Hz for 0.5 sec.
- EGT – DNA electrotransfer:
EP – electroporation, GO – gene ontology, HV – high voltage, LV – low voltage.
The authors would like to thank Marianne Fregil and Vibeke Uhre for excellent technical assistance and Annette Enghøj at the Department of Pathology, University of Copenhagen at Herlev Hospital, for performing the histology. The study was supported by the Danish Research Agency (22-03-0367).
- Wolff JA, Malone RW, Williams P, Chong W, Acsadi G, Jani A, Felgner PL: Direct gene transfer into mouse muscle in vivo. Science 1990, 247: 1465-1468. 10.1126/science.1690918View ArticlePubMedGoogle Scholar
- Mir LM, Bureau MF, Gehl J, Rangara R, Rouy D, Caillaud JM, Delaere P, Branellec D, Schwartz B, Scherman D: High-efficiency gene transfer into skeletal muscle mediated by electric pulses. Proc Natl Acad Sci U S A 1999, 96: 4262-4267. 10.1073/pnas.96.8.4262PubMed CentralView ArticlePubMedGoogle Scholar
- Gehl J: Electroporation: theory and methods, perspectives for drug delivery, gene therapy and research. Acta Physiol Scand 2003, 177: 437-447. 10.1046/j.1365-201X.2003.01093.xView ArticlePubMedGoogle Scholar
- Rizzuto G, Cappelletti M, Maione D, Savino R, Lazzaro D, Costa P, Mathiesen I, Cortese R, Ciliberto G, Laufer R, Monica N, Fattori E: Efficient and regulated erythropoietin production by naked DNA injection and muscle electroporation. Proc Natl Acad Sci U S A 1999, 96: 6417-6422. 10.1073/pnas.96.11.6417PubMed CentralView ArticlePubMedGoogle Scholar
- Terada Y, Tanaka H, Okado T, Inoshita S, Kuwahara M, Akiba T, Sasaki S, Marumo F: Efficient and ligand-dependent regulated erythropoietin production by naked dna injection and in vivo electroporation. Am J Kidney Dis 2001, 38: S50-S53.View ArticlePubMedGoogle Scholar
- Heller LC, Ingram SF, Lucas ML, Gilbert RA, Heller R: Effect of electrically mediated intratumor and intramuscular delivery of a plasmid encoding IFN alpha on visible B16 mouse melanomas. Technol Cancer Res Treat 2002, 1: 205-209.View ArticlePubMedGoogle Scholar
- Lucas ML, Heller R: Immunomodulation by electrically enhanced delivery of plasmid DNA encoding IL-12 to murine skeletal muscle. Mol Ther 2001, 3: 47-53. 10.1006/mthe.2000.0233View ArticlePubMedGoogle Scholar
- Bettan M, Emmanuel F, Darteil R, Caillaud JM, Soubrier F, Delaère P, Branelec D, Mahfoudi A, Duverger N, Scherman D: High-Level Protein Secretion into Blood Circulation after Electric Pulse-Mediated Gene Transfer into Skeletal Muscle. Mol Ther 2000, 2: 204-210. 10.1006/mthe.2000.0117View ArticlePubMedGoogle Scholar
- Muramatsu T, Arakawa S, Fukazawa K, Fujiwara Y, Yishida T, Sasaki R, Masuda S, Park HM: In vivo gene electroporation in skeletal muscle with special reference to the duration of gene expression. Int J Mol Med 2001, 7: 37-42.PubMedGoogle Scholar
- Bureau MF, Gehl J, Deleuze V, Mir LM, Scherman D: Importance of association between permeabilisation and electrophoretic forces for intramuscular DNA electrotransfer. Biochim Biophys Acta 2000, 1474: 353-359.View ArticlePubMedGoogle Scholar
- Satkauskas S, Bureau MF, Puc M, Mahfoudi A, Scherman D, Miklavcic D, Mir LM: Mechanisms of in vivo DNA electrotransfer: Respective contributions of cell electropermeabilisation and DNA electrophoresis. Mol Ther 2002, 5: 133-140. 10.1006/mthe.2002.0526View ArticlePubMedGoogle Scholar
- Satkauskas S, Andre F, Bureau M, Scherman D, Miklavcic D, Mir LM: Electrophoretic component of electric pulses determines the efficacy of in vivo DNA electrotransfer. Hum Gene Ther 2005, 16: 1194-1201. 10.1089/hum.2005.16.1194View ArticlePubMedGoogle Scholar
- Pavselj N, Preat V: DNA electrotransfer into the skin using a combination of one high- and one low-voltage pulse. J Control Release 2005, 106: 407-415. 10.1016/j.jconrel.2005.05.003View ArticlePubMedGoogle Scholar
- Gothelf A, Mir LM, Gehl J: Electrochemotherapy: results of cancer treatment using enhanced delivery of bleomycin by electroporation. Cancer Treat Rev 2003, 29: 371-387. 10.1016/S0305-7372(03)00073-2View ArticlePubMedGoogle Scholar
- Sersa G, Cemazar M, Rudolf Z: Electrochemotherapy: advantages and drawbacks in treatment of cancer patients. Cancer Ther 2003, 1: 133-142.Google Scholar
- Mir LM: Bases and rationale of the electrochemotherapy. Eur J Cancer Suppl 2006, 4: 38-44. 10.1016/j.ejcsup.2006.08.005View ArticleGoogle Scholar
- Marty M, Sersa G, Garbay JR, Gehl J, Collins CG, Snoj M, Billard V, Geertsen PF, Larkin JO, Miklavcic D: Electrochemotherapy - An easy, highly effective and safe treatment of cutaneous and subcutaneous metastases: Results of ESOPE (European Standard Operating Procedures of Electrochemotherapy) study. Eur J Cancer Suppl 2006, 4: 3-13. 10.1016/j.ejcsup.2006.08.002View ArticleGoogle Scholar
- Mir LM, Gehl J, Sersa G, Collins CG, Garbay J, Billard V, Geertsen PF, Rudolf Z, O'Sullivan GC, Marty M: Standard operating procedures of the electrochemotherapy: Instructions for the use of bleomycin or cisplatin administered either systemically or locally and electric pulses delivered by the CliniporatorTM by means of invasive or non-invasive electrodes. Eur J Cancer Suppl 2006, 4: 14-25. 10.1016/j.ejcsup.2006.08.003View ArticleGoogle Scholar
- Rubenstrunk A, Mahfoudi A, Scherman D: Delivery of electric pulses for DNA electrotransfer to mouse muscle does not induce the expression of stress related genes. Cell Biol Toxicol 2004, 20: 25-31. 10.1023/B:CBTO.0000021079.69089.4fView ArticlePubMedGoogle Scholar
- Schertzer JD, Plant DR, Lynch GS: Optimizing Plasmid-Based Gene Transfer for Investigating Skeletal Muscle Structure and Function. Mol Ther 2006, 13: 795-803. 10.1016/j.ymthe.2005.09.019View ArticlePubMedGoogle Scholar
- Mir LM, Moller PH, Andre F, Gehl J: Electric pulse-mediated gene delivery to various animal tissues. Advances in Genetics 54th edition. 2005, 83-114.Google Scholar
- Lefesvre P, Attema J, van Bekkum D: A comparison of efficacy and toxicity between electroporation and adenoviral gene transfer. BMC Molecular Biology 2002, 3: 3-12. 10.1186/1471-2199-3-12View ArticleGoogle Scholar
- Durieux AC, Bonnefoy R, Freyssenet D: Kinetic of transgene expression after electrotransfer into skeletal muscle: Importance of promoter origin/strength. Biochim Biophys Acta 2005,1725(3):403-409.View ArticlePubMedGoogle Scholar
- Churchill GA: Using ANOVA to analyse microarray data. Biotechniques 2004, 37: 173.PubMedGoogle Scholar
- Kakorin S, Neumann E: Kinetics of the electroporative deformation of lipid vesicles and biological cells in an electric field. Ber Bunsenges Phys Chem 1998, 102: 670-675.View ArticleGoogle Scholar
- Neumann E, Kakorin S, Toensing K: Fundamentals of electroporative delivery of drugs and gene. Bioelectrochem Bioenerg 1999, 48: 3-16. 10.1016/S0302-4598(99)00008-2View ArticlePubMedGoogle Scholar
- Teissie J, Rols MP: Manipulation of cell cytoskeleton affects the lifetime of cell membrane electropermeabilisation. Ann NY Acad Sci 1994, 720: 98-110. 10.1111/j.1749-6632.1994.tb30438.xView ArticlePubMedGoogle Scholar
- Gass GV, Chernomordik LV: Reversible large-scale deformations in the membranes of electrically-treated cells: electroinduced bleb formation. Biochim Biophys Acta 1990, 1023: 1-11. 10.1016/0005-2736(90)90002-6View ArticlePubMedGoogle Scholar
- Geiger RC, Taylor W, Glucksberg MR, Dean DA: Cyclic stretch-induced reorganization of the cytoskeleton and its role in enhanced gene transfer. Gene Ther 2006, 13: 725-731. 10.1038/sj.gt.3302693PubMed CentralView ArticlePubMedGoogle Scholar
- Chen W: Supra-physiological membrane potential induced conformational changes in K+ channel conducting system of skeletal muscle fibers. Bioelectrochemistry 2004, 62: 47-56. 10.1016/j.bioelechem.2003.10.006View ArticlePubMedGoogle Scholar
- Chen W: Evidence of electroconformational changes in membrane proteins: field-induced reductions in intra membrane nonlinear charge movement currents. Bioelectrochemistry 2004, 63: 333-335. 10.1016/j.bioelechem.2003.12.007View ArticlePubMedGoogle Scholar
- Chen W, Han Y, Chen Y, Astumian D: Electric Field-Induced Functional Reductions in the K+ Channels Mainly Resulted from Supramembrane Potential-Mediated Electroconformational Changes. Biophys J 1998, 75: 196-206.PubMed CentralView ArticlePubMedGoogle Scholar
- Hartikka J, Sukhu L, Buchner C, Hazard D, Bozoukova V, Margalith M, Nishioka WK, Wheeler CJ, Manthrope M, Sawdey M: Electroporation-facilitated delivery of plasmid DNA in skeletal muscle: Plasmid dependence of muscle damage and effect of poloxamer 188. Mol Ther 2001, 4: 407-415. 10.1006/mthe.2001.0483View ArticlePubMedGoogle Scholar
- Durieux A, Bonnefoy R, Busso T, Freyssenet D: In vivo gene electrotransfer into skeletal muscle: effects of plasmid DNA on the occurrence and extent of muscle damage. J Gene Med 2004, 6: 809-816. 10.1002/jgm.534View ArticlePubMedGoogle Scholar
- Lukacs GL, Haggie P, Seksek O, Lechardeur D, Freedman N, Verkman AS: Size-dependent DNA mobility in cytoplasm and nucleus. J Biol Chem 2000, 275: 1625-1629. 10.1074/jbc.275.3.1625View ArticlePubMedGoogle Scholar
- Lechardeur D, Sohn.K.J., Haardt M, Joshi JB, Monck M, Graham RW, Beatty B, Squire J, O'Brodowich H, Lukacs GL: Metabolic instability of plasmid DNA in the cytosol: a potential barrier to gene transfer. Gene Ther 1999, 6: 482-497. 10.1038/sj.gt.3300867View ArticlePubMedGoogle Scholar
- Pollard H, Toumaniantz G, Amos JL, Avet-Loiseau H, Guihard G, Behr JP, Escande D: Ca2+-sensitive cytosolic nucleases prevent efficient delivery to the nucleus of injected plasmids. J Gene Med 2001, 3: 153-164. 10.1002/jgm.160View ArticlePubMedGoogle Scholar
- Vaughan EE, Dean DA: Intracellular Trafficking of Plasmids during Transfection Is Mediated by Microtubules. Molecular Therapy 2005, 13: 422-428. 10.1016/j.ymthe.2005.10.004PubMed CentralView ArticlePubMedGoogle Scholar
- Ataka K, Maruyama H, Neichi T, Miyazaki J, Gejko F: Effects of erythropoietin-gene electrotransfer in rats with adenine-induced renal failure. Am J Nephrol 2003, 23: 315-323. 10.1159/000072913View ArticlePubMedGoogle Scholar
- Chenuaud P, Larchet T, Rabinowitz JE, Provost N, Cherel Y, Casadevall N, Samulski RJ, Moullier P: Autoimmune Anemia in Macaques Following Erythropoietin Gene Therapy. Blood 2004, 103: 3303-3304. 10.1182/blood-2003-11-3845View ArticlePubMedGoogle Scholar
- Maruyama H, Ataka K, Gejyo F, Higuchi N, Ito Y, Hirahara H, Imazeki I, Hirata M, Ichikawa F, Neichi T, Kikuchi H, Sugawa M, Miyazaki J: Long-term production of erythropoietin after electroporation-mediated transfer of plasmid DNA into the muscles of normal and uremic rats. Gene Ther 2001, 8: 461-468. 10.1038/sj.gt.3301412View ArticlePubMedGoogle Scholar
- Rizzuto G, Cappelletti M, Mennuni C, Wiznerowicz M, DeMartis A, Maione D, Ciliberto G, La Monica N, Fattori E: Gene electrotransfer results in a high-level transduction of rat skeletal muscle and corrects anemia of renal failure. Hum Gene Ther 2000, 11: 1891-1900. 10.1089/10430340050129503View ArticlePubMedGoogle Scholar
- Payen E, Bettan M, Rouyer-Fessard P, Beuzard Y, Scherman D: Improvement of mouse [beta]-thalassemia by electrotransfer of erythropoietin cDNA. Exp Hematol 2001, 29: 295-300. 10.1016/S0301-472X(00)00679-2View ArticlePubMedGoogle Scholar
- Samakoglu S, Fattori E, Lamartina S, Toniatti C, Stockholm D, Heard JM, Bohl D: betaMinor-globin messenger RNA accumulation in reticulocytes governs improved erythropoiesis in beta thalassemic mice after erythropoietin complementary DNA electrotransfer in muscles. Blood 2001, 97: 2213-2220. 10.1182/blood.V97.8.2213View ArticlePubMedGoogle Scholar
- Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem 1987, 162: 156-159. 10.1016/0003-2697(87)90021-2View ArticlePubMedGoogle Scholar
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: Expression index computation and outlier detection. Proc Natl Acad Sci U S A 2001, 98: 31-36. 10.1073/pnas.011404098PubMed CentralView ArticlePubMedGoogle Scholar
- Clausen T, Everts ME: K(+)-induced inhibition of contractile force in rat skeletal muscle: role of active Na(+)-K+ transport. Am J Physiol Cell Physiol 1991, 261: C799-C807.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.