- Research article
- Open Access
Gene expression profiles in skeletal muscle after gene electrotransfer
BMC Molecular Biology volume 8, Article number: 56 (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
In an unsupervised cluster analysis, the EP groups clustered significantly (p = 0.0006) together as depicted from the sample correlation matrix plot in Fig. 1. No correlation within the EP or non-EP clusters, as to whether the groups had received plasmid injection or at which time following treatment the muscles had been excised, was established (additional file 1).
Performing gene ontology annotation (GO) on the cluster analysis revealed that RNA encoding proteins assembling the cytoskeleton e.g. microtubule-associated protein 7 (NM008635), cadherin-like 26 (BB748621) and sarcospan (BC021484) were down-regulated in the electroporated groups. Yet, 3 weeks after treatment the expression of these genes in the DNA+HV+LV group had returned to control levels (Table 1).
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
The plasmid injected groups clustered significantly together (p = 0.0006) in the sample correlation matrix analysis with clear separation from the non-injected groups (Fig. 2). Interestingly the LDA plot revealed that within the plasmid injected groups, two apparent sub-clusters of EP and non-EP groups were found. Likewise among the non-injected groups, the EP groups formed a distinct sub-cluster (additional file 2).
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).
Groups within each of the 4 time points clustered significantly together (Fig. 3) regardless of the treatment. The LDA plot revealed that groups from 4 and 48 hrs clustered closely together, while the control and 3 weeks groups formed individual clusters (additional file 3).
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
The gene expression levels were validated by Q-PCR for 11 selected genes (Table 2). These measurements confirmed the level of mRNA determined by microarray analysis. For example transducer of ERBB2 (AV174616) was in the array analysis found to be up-regulated at 48 hours after all treatments, indeed Q-PCR also showed that the transducer of ERBB2 was up-regulated at 48 hours e.g. in the case of DNA+HV+LV treatment by 6-fold. Similarly the expression of the solute carrier family 6 (BG069516), responsible for creatine transport, was verified by Q-PCR to be down-regulated at 4 and 48 hours.
Histological changes in tibialis cranialis muscles after DNA electrotransfer
Preparations of muscles slides were stained with HE and examined for morphological changes like integrity of fibres, central nuclei and fainting striation pattern. Furthermore, CD4 and CD8 staining were performed on parallel slides to verify infiltration of mononuclear cells. At 4 hrs, no histological changes were observed for the control, DNA injected and electroporated muscles (Fig. 4). Contrarily after DNA+HV+LV treatment, loss of fibre structure and fainting striation pattern was observed in 15% of the muscle. After 48 hrs, no changes were observed in the control and DNA injected muscles, while the electroporated muscles exhibited fainting striation pattern in a few muscle fibres. In the DNA+HV+LV treated muscles, structural changes with loss of cell integrity and striation pattern were observed in 15% of the fibres, in addition limited intra-fiber infiltration of mononuclear cells was also observed. Three weeks after treatment, the DNA+HV+LV treated muscles had around 5% of the fibres containing central nuclei, indicating regenerating myocytes. In the electroporated muscles single cells with central nuclei were also detectable, while the control and DNA injected muscles showed no signs of histological changes.
No loss in muscle performance after DNA electrotransfer
The reflex and motor function of mice were tested by assessing the running and seizing patterns. Thus, mice were allowed to run on the top of the cage, their tails were gently manually lifted, and it was observed whether they could grasp the bars with normal force. Similarly, spontaneous motoric movement was observed in the cage. Fig. 5 depicts the average score, short duration films of the test results can be found in additional file 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19. At 4 hours the mice had not fully awoken from the anaesthesia and could therefore not be tested. At 48 hours no effects were observed in the control and DNA injected mice, while slight impairments were observed in a couple of the electroporated and DNA+HV+LV treated mice. From 5 days after treatment no functional impairments in the seizing and running patterns were observed in any of the mice.
The effects of DNA electrotransfer on the force generating capacity were tested in rat EDL muscles 2 weeks after treatment and depicted in Fig. 6. Rat EDL muscles were chosen for this measurements as mice TA cannot be used in our set-up. The force generating capacity in the DNA electrotransferred and electroporated muscles did not differ significantly from the control muscles.
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.
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.1690918
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.4262
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.x
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.6417
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.
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.
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.0233
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.0117
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.
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.
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.0526
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.1194
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.003
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-2
Sersa G, Cemazar M, Rudolf Z: Electrochemotherapy: advantages and drawbacks in treatment of cancer patients. Cancer Ther 2003, 1: 133-142.
Mir LM: Bases and rationale of the electrochemotherapy. Eur J Cancer Suppl 2006, 4: 38-44. 10.1016/j.ejcsup.2006.08.005
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.002
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.003
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.4f
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.019
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.
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-12
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.
Churchill GA: Using ANOVA to analyse microarray data. Biotechniques 2004, 37: 173.
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.
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-2
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.x
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-6
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.3302693
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.006
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.007
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.
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.0483
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.534
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.1625
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.3300867
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.160
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.004
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/000072913
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-3845
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.3301412
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/10430340050129503
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-2
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.2213
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-2
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.011404098
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.
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).
PHM carried out the experimental work and drafted the manuscript. JRZ performed the microarray analysis and statistics and helped revise the manuscript. HG performed the force generation measurements. JE participated in the design of the study. JG participated in the design of the study, interpretation of results and drafting the manuscript. All authors have read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
About this article
Cite this article
Hojman, P., Zibert, J.R., Gissel, H. et al. Gene expression profiles in skeletal muscle after gene electrotransfer. BMC Molecular Biol 8, 56 (2007). https://doi.org/10.1186/1471-2199-8-56
- Green Fluorescent Protein Expression
- Force Generate Capacity
- Gene Electrotransfer
- Plasmid Injection
- Intracellular Transport Protein