The human PINK1 locus is regulated in vivo by a non-coding natural antisense RNA during modulation of mitochondrial function
© Scheele et al; licensee BioMed Central Ltd. 2007
Received: 13 March 2007
Accepted: 15 March 2007
Published: 15 March 2007
Mutations in the PTEN induced putative kinase 1 (PINK1) are implicated in early-onset Parkinson's disease. PINK1 is expressed abundantly in mitochondria rich tissues, such as skeletal muscle, where it plays a critical role determining mitochondrial structural integrity in Drosophila.
Herein we characterize a novel splice variant of PINK1 (svPINK1) that is homologous to the C-terminus regulatory domain of the protein kinase. Naturally occurring non-coding antisense provides sophisticated mechanisms for diversifying genomes and we describe a human specific non-coding antisense expressed at the PINK1 locus (naPINK1). We further demonstrate that PINK1 varies in vivo when human skeletal muscle mitochondrial content is enhanced, supporting the idea that PINK1 has a physiological role in mitochondrion. The observation of concordant regulation of svPINK1 and naPINK1 during in vivo mitochondrial biogenesis was confirmed using RNAi, where selective targeting of naPINK1 results in loss of the PINK1 splice variant in neuronal cell lines.
Our data presents the first direct observation that a mammalian non-coding antisense molecule can positively influence the abundance of a cis-transcribed mRNA under physiological abundance conditions. While our analysis implies a possible human specific and dsRNA-mediated mechanism for stabilizing the expression of svPINK1, it also points to a broader genomic strategy for regulating a human disease locus and increases the complexity through which alterations in the regulation of the PINK1 locus could occur.
PTEN induced putative kinase 1 gene (PINK1) is a serine-threonine kinase directly linked to a recessive form of familial parkinsonism [1–4]. A mutation at the nucleotide binding site within the kinase domain renders the protein unable to protect neuroblastoma cells from apoptosis, whereas over expression of the native peptide protects SH-SY5Y cells [1, 3]. Over expressed tagged PINK1 localizes to the mitochondria [2, 4], raising the possibility that PINK1 phosphorylates and regulates proteins involved in oxidative phosphorylation or the mitochondrial translocator pore [5, 6], processes linked to neuronal cell death [7, 8]. Indeed, PINK1 inhibits mitochondrial cytochrome c release, attenuating the general apoptosis machinery [3, 6], while mitochondrial localization is supported by the existence of an 8.5 kDa mitochondrial import tag  and in vitro mitochondrial translocation . Most recently, Gandhi et al.  produced evidence that endogenous PINK1 protein locates to the mitochondrial membranes, while three additional articles highlighted a critical role for PINK1 in Drosophila flight muscle mitochondria [10–12]. A fourth study in Drosophila show an antioxidant rescue of neurodegeneration induced by knockdown of Pink1 . It is clear, therefore, that PINK1 may play a rather central physiological role in neuronal and myocyte energy metabolism, consistent with its abundant expression in mitochondria-rich tissue  and consistent with the proposed mitochondrial basis for Parkinsonism.
The C-terminal end of the protein regulates PINK1 kinase activity . Intriguingly, a predicted novel short splice variant of PINK1 would translate into a protein sequence which represents the C-terminus portion of the protein. Furthermore, a cis-transcribed non-coding natural antisense (ncNAT) originating from the PINK1 gene locus is predicted from the human EST databases  and is now, very recently, listed in a natural antisense transcript database . According to the available well annotated databases, PINK1 ncNAT is not apparent in other species. The ncNAT demonstrates near complete sequence overlap with the short PINK1 splice variant and only partial overlap with PINK1, at the 3'end. Widespread expression of natural antisense transcripts (NAT) has, in recent times [17–19], emerged as a potential mechanism for bringing diversity and regulatory complexity to a surprisingly finite human 'protein coding' genome . While there are a handful of examples of NAT (coding or non-coding) which negatively regulate protein coding mRNA expression , few are associated with human disease  and none have been shown to directly regulate a cis-transcribed partner in a positive or concordant manner. In addition, despite the annotation of several thousands of ncNAT's in mammalian genomes, few have had their role in mammalian biology clarified [17, 20].
Greater appreciation of the mechanisms responsible for regulating the PINK1 gene locus, beyond rare mutations, is of importance for a number of reasons. For example, it is likely that alteration of PINK1 function, regardless of mechanism, would influence mitochondrial function [1, 10–12, 21]. In addition, as PINK1 appears to influence mitochondrial membrane potential  and mitochondrial dysfunction is implicated in many age related diseases , studying the regulation of the PINK1 locus during dynamic in vivo modulation of mitochondrial function is likely to be relevant to human ageing. While it is impossible to study this process in the adult human brain, it is plausible to model gains in mitochondrial function in adult human skeletal muscle  and examine the regulation of transcription from the PINK1 locus. In doing so we reveal direct evidence that there is a discordant transcriptional response between the two coding PINK1 splice variants, in vivo, during dynamic alterations in mitochondrial function, whereas there is a concordant regulation of short transcript variant of PINK1 (svPINK1) by a cis-transcribed ncNAT (naPINK1) originating from the PINK1 gene locus. Not only do we provide the first direct evidence that a ncNAT can positively regulate the abundance of a sense coding transcript in humans under physiological conditions, we also highlight that the complexity of PINK1 locus regulation is potentially greater in humans than other species.
Analysis of the PINK1 locus
An alignment of the 5' upstream regions (promoter) of PINK1 in mouse and human indicate a high degree of conservation, implying that PINK1 transcriptional regulation may be similar across these species. A search for transcription factor binding sites within the PINK1 5' region  identified PPARγ and numerous myogenic transcription factor binding sites (see Additional file 2), supporting the robust expression of PINK1 in mitochondria rich striated muscle and the recently demonstrated and devastating impact of loss of PINK1 expression in Drosophila flight muscle . We also examined the genomic sequences of DDOST and PINK1 (see Additional file 2). The intergenic region and the last exons of both PINK1 and DDOST are noticeably less conserved between human and mouse than the rest of the locus (approximately 60% conservation compared with >80%), which supports the idea that the region from where the naPINK1 originates, is specific to the human genome. Interestingly, DDOST encodes the AGE-R1 scavenger receptor, a counter regulatory mechanism for advanced glycation end-product induced oxidative stress , with a role in diabetes and ageing, processes clearly relevant to Parkinson's disease . Evidence that this 'intergenic' region between the 3' end of DDOST and the 5' end of naPINK1 is subject to post-transcriptional regulation comes from the observation that between the 3' end of DDOST and the 5' end of naPINK1, there are numerous consensus sequences for endogenous microRNA's (see Additional file 2) providing a potential mechanism through which the naPINK1 transcript could be processed from a large combined transcript (see below).
Detection of the RNA transcripts from the PINK1 locus
PINK1 qRT-PCR forward primer
PINK1 MGB probe
PINK1 qRT-PCR reverse primer
svPINK1 qRT-PCR forward primer
svPINK1 MGB probe
svPINK1 qRT-PCR reverse primer
naPINK1 qRT-PCR forward primer
naPINK1 qRT-PCR reverse primer
Northern forward primer
Northern reverse primer
5' RACE outer primer
5' RACE inner primer
RACE product internal forward primer
RACE product internal reverse primer
AS siRNA 1 sense strand
AS siRNA 1 antisense strand
AS siRNA 2 sense strand
AS siRNA 2 antisense strand
S-PINK1 siRNA (spanning two exons) sense strand
S-PINK1 siRNA (spanning two exons) antisense strand
In vivo human transcript regulation of the PINK1 genomic locus
RNAi demonstrates that naPINK1 regulates svPINK1 abundance
Discussions and conclusion
In the present study we demonstrate that a naturally expressed non coding antisense RNA and a short splice variant of PINK1 exist in human neuroblastoma cells and human skeletal muscle tissue. We further demonstrate that this NAT regulates the PINK1 locus in a manner that involves stabilizing, or promotion of, the expression of the PINK1 splice variant, svPINK1. This is the first compelling evidence that an ncNAT positively modulates a sense coding RNA at levels found in vivo in humans. Furthermore, given that svPINK1 codes for a homologue of the C-terminus of PINK1; a peptide sequence which regulates PINK1 kinase activity , then modulation of naPINK1 expression may have direct relevance in Parkinson's disease.
Teleologically, any process that regulates PINK1 activity may influence Parkinson's disease progression and thus it is important to consider the broader regulation of the PINK1 locus. The scale of the FANTOM Consortium mammalian antisense transcriptome project  enabled robust estimations, using gene ontology based classification scores, of the function of the coding genes which have cis-transcribed associated NATs. For convergent NAT species, the corresponding coding transcripts were most likely translated into catalytic proteins involved in nucleotide metabolism, while for fully overlapping cis-transcribed NAT species, many of the corresponding sense mRNAs were translated into a protein with protein binding abilities or kinase activity. These statistical associations are clearly in agreement with the naPINK1 and PINK1/svPINK1 relationships. In the present study we provide robust evidence of regulation by the naPINK1 transcript of the short variant of the Parkinson's disease gene, PINK1, under conditions where the mRNA species do not exceed their natural abundances (i.e. we do not rely on artificial over expression to study these RNA interactions).
A further point of interest is that the naPINK1 transcript is formed at a region concurrent with the 3' end of the DDOST gene. Examination reveals no obvious promoter region within the overlapping region that encodes for naPINK1 and the 3'UTR of DDOST, and so we can speculate that naPINK1 and DDOST are transcribed from a common promoter and that naPINK1 is a result of transcription read-through. Our Northern blot analysis (Fig. 1B) supports this idea as we found a large transcript (6–7 kb) in both the total RNA and isolated mRNA preparations that could represent a DDOST-naPINK1 transcript (which should be 6.4 kb). It has been postulated that neurodegenerative diseases, such as Parkinson's, are associated with oxidative stress [7, 8] and there is further association between metabolic syndrome and neurodegeneration [32–34]. DDOST encodes a receptor that mediates suppression of advanced glycation mediated cellular oxidative stress and apoptosis . While transcriptional modulation of the DDOST gene has not been studied, it's possible that metabolic syndrome induced transcriptional modulation of DDOST could subsequently result in altered naPINK1 expression.
We must also ask how naPINK1 is able to 'stabilize' the svPINK1 transcript. In general, the precise molecular mechanisms through which ncNAT leads to the regulation of a sense mRNA have not been clarified [17, 19, 12, 35]. RNA duplex formation leading to RNA editing or degradation has been suggested, while antisense induced chromosomal modifications has been demonstrated and result in gene silencing. For example, the murine ncNAT Air , which regulates the maternally expressed igf2 receptor gene (igfr), highlights a multi-gene interaction and provides a clear indication of the complexity of NAT regulated transcription. Although only over-lapping igf2r, the paternal antisense transcript Air, silences the expression of igf2r, Slc22a2 and Slc22a3 through an ill-defined epigenetic modification . There are also examples where the antisense transcript regulates transcription (N-myc), inhibits splicing (Rev-erbAα), transport (mld) or promotes degradation (bFGF). These examples all involve a reduction in mRNA expression . A short co-transcribed ncNAT from the WT1 gene loci has been associated with greater abundance of WT1 protein  yet in vitro mechanistic studies failed to clarify this effect. In our case, naPINK1 positively regulates the abundance of svPINK1, both in vivo and during siRNA knockdown experiments and thus we present a novel observation for a ncNAT. Importantly, we demonstrate that the detectable levels of naPINK1 and svPINK1 are likely to be mature polyadenylated transcripts and so the strong in vivo association may reflect a cytosolic interaction. Insight into the mechanism through which naPINK1 may regulate svPINK1 also comes from the observation that naPINK1 appears only to influence the abundance of the near full length complementary cis transcribed svPINK1 sequence (only 52 bases of the 5' end of svPINK1 are non-complementary) suggesting a potential role for dsRNA mediated stabilization. Recent data from chloroplasts, involving the mitochondrial ATP synthase gene (atpB), demonstrated that a chloroplast mutant that produced a complementary coded atpb NAT mRNA resulted in dsRNA formation, which prevented exonuclease mediated degradation of the atpB mRNA . It is also thought that the degree of RNA editing, modulates activity of the RNAi pathway, and it is plausible that this might explain the current observations.
It was recently reported that siRNA mediated transcript degradation can act in the nucleus of mammals , however this was demonstrated using a high concentration of a siRNA containing a sequence known to result in off-target effects . In our experiments, RNAi directed degradation of naPINK1 is most likely to occur in the cytosol, and we clearly demonstrate that this causes a significant reduction of svPINK1. One could argue that naPINK1 may influence the splicing of PINK1 pre-mRNA. Under this scenario, association between the 5' end of the PINK1 pre-processed mRNA with the natural antisense molecule would somehow favor the alternative splicing process to produce svPINK, starting from the intronic region between Exon 4 and Exon 5 of PINK1, implying that naPINK1 potentially masks a splicing inhibitor site. In agreement with this idea, if one examines the in vivo relationship between the ratio of PINK1:svPINK1 RNA then this is also linearly related to the expression level of the naPINK1 (see Additional file 3). Nevertheless, this correlation is driven by the svPINK1-naPINK1 relationship and it is clear that PINK1 itself does not co-vary with the naPINK1 RNA under the present experimental conditions (see Additional file 3). Given that we demonstrate that siRNA H and I do not target the pre-spliced version of svPINK1 or PINK1 (which prior to nuclear processing and export contain a common intronic region identical to the siRNA targeted region on the naPINK1 transcript) and that is possible that svPINK1-naPINK1 relationship reflects alternative splicing, then we have to consider that the synthetic siRNA may indeed direct RNAi within the nucleus (naPINK1 knock-down) with a subsequent change in PINK1 splicing, a consideration that merits more robust investigation.
The regulation of PINK1 protein expression is likely to be complex, particularly under disease conditions [9–12, 42]. Available evidence suggests that the full length protein is processed (removal of its mitochondrial localization tag [2, 9]), resulting in the detection of at least two proteins on a western blot. Neither of these proteins is consistent with the predicted size of a svPINK1 protein (they are both too large). Multiple shorter proteins have been found when GST-tagged PINK1 is over-expressed and detected using an antibody to GST  and while others have presented western blots demonstrating a single protein band , so far all commercial antibodies that we purchased reproducibly (or received as a gift) give rise to >8 discrete bands and thus are not reliable when profiling cell models. As more reliable reagents become available  it will become important to combine well designed RNAi experiments with protein detection, to study the biochemical consequence of loss of PINK1, svPINK1 and naPINK1 in human cells. From a speciation perspective, there is prior evidence that natural antisense RNA sequences can be conserved selectively in higher mammals. A cis-transcribed intronic antisense mRNA from the FAS receptor gene, named Saf, alters FAS splice variant expression when Saf is artificially over-expressed  and Saf appears only to be conserved in primates. Future experiments will also have to address the relative importance of naPINK and svPINK in relation to PINK1 function. For example, the relative abundance of naPINK approaches that of PINK1 in tumor cell lines, while naPINK1 is less abundant, although clearly present in human muscle tissue. Ongoing svPINK1 over-expression studies should resolve the function of svPINK1, while we must also decide whether svPINK1 is a truly independent 'gene' transcribed from within the PINK1 genomic loci. In addition, given the complexity of defining what merits being called a 'gene' it is reasonable to suggest that a NAT, with a proven biological function, should also be considered a gene.
In summary, we produce in vivo human data and supportive RNAi data which demonstrates that a non-coding RNA so far limited to human expression, appears to 'stabilize' a sense coding RNA. While our data is the first evidence that endogenous levels of a non-coding antisense RNA influencing a sense coded mRNA in a concordant fashion, we can not rule out the possibility that naPINK1 recruits additional cytosolic proteins  to bind and stabilize svPINK1 in a more complex scenario than direct complementary hybridization. Given the close genomic association between DDOST and naPINK specifically in the human genome, we speculate that these two disease related genes may be co-regulated, and uniquely so in humans.
Cell culture and transfection
All cell culture reagents were purchased from GIBCO Invitrogen (Sweden) unless otherwise stated. Neuroblastoma cell lines SK-N-MC, SK-N-DZ, SK-N-F1, SK-N-AS were sub-cultured in Minimal Essential Media (MEM) + Earle's supplemented with 10% Fetal Bovine serum (FBS), 2 mM L-glutamine, 1 mM sodium pyruvate, 1 × non-essential amino acids, 100 U/ml penicillin and 100 μg/ml streptomycin. Neuroblastoma cell lines SH-SY5Y and SK-N-Be(2) was sub-cultured in 50% MEM + Earle's, 50% F12 (HAM) supplemented with 10% FBS, 2 mM L-glutamine, 1 mM sodium pyruvate, 1 × non-essential amino acids, 100 U/ml penicillin and 100 μg/ml streptomycin. Neuroblastoma cell line SK-N-SH was sub-cultured in RPMI 1640 supplemented with 10% FBS, 2 mM L-glutamine, 1 × non-essential amino acids 100 U/ml penicillin and 100 μg/ml streptomycin. Cell cultures were maintained at 37°C in 5% CO2.
Prior to transfection, cells were counted using a hemacytometer and 50,000 SK-N-MC cells or 70,000 SH-SY5Y cells were seeded per well, into 6-well plates, 48 hours or 24 hours before transfection, respectively. SK-N-MC attach more slowly to the surface after splitting than SH-SY5Y and are sensitive to transfection if confluence is less than 75%. On the other hand, transfection efficiency of SH-SY5Y was improved by reducing the seeding density. Transfection was performed using 0.25% Lipofectamine 2000 (Invitrogen, Sweden) and 20 nM siRNA in 2.5 ml serum-containing cell culture media (without antibiotics) according to the manufacturer's protocol. Briefly, Lipofectamine 2000 was pre-mixed with Optimem, incubated for 5 minutes and then mixed with siRNA. Following 20 minutes of incubation, the Lipofectamine 2000 – siRNA mix was applied to cells. Media was changed 4 hours later and cells were harvested 48 hrs after transfection for RNA isolation. The siRNAs (Table S1) targeting naPINK1 were designed using an on-line design tool called BLOCK-iT™ RNAi Designer (Invitrogen) and chemically synthesized by Invitrogen. The siRNA targeting PINK1 (Table S1) was a pre-validated siRNA from Ambion (siRNA ID # 1199) and thus designed and chemically synthesized by Ambion. The control siRNA was pre-designed by Ambion to not target any gene in the genome (Silencer®Negative Control#1, Ambion)
RNA isolation and quantitative real-time PCR (qRT-PCR)
Human muscle tissue biopsies were homogenized in TRIzol (Invitrogen) using a motor-driven homogenizer (Polytron, Kinematica) and total RNA was isolated according to the manufacturer's protocol. In brief, two paired samples (before and after activity) were processed in parallel in 1 ml TRIzol at any given time. The sizes of the tissues ranged between 10 and 40 mg. Total RNA was dissolved in 30 μl RNase-free water and quantified using a Nanodrop (NanoDrop Technologies). Typically, the yield was 8–10 μg of total RNA from 20 mg of tissue. For the cell samples, total RNA was isolated and on-column DNase-treated using RNeasy Mini Kit (Qiagen) and RNase-Free DNase set (Qiagen) according to the manufacturer's protocol. Total RNA was eluted in 30 μl of RNase-free water and concentration was measured using a spectrophotometer (UltroSpec 3000, Pharmacia Biotech). The typical yield of SH-SY5Y was 700 ng and of SK-N-MC 5 μg of total RNA per well in a 6-well plate.
1 μg muscle total RNA or 200 ng neuroblastoma cell total RNA, in reaction volumes of 40 μl was reverse transcribed using RT-PCR kit (Applied Biosystems) according to the manufacturer's protocol. Typically, random hexamers were used for the first strand cDNA synthesis. The advantage of using random hexamers is the possibility to use 18S rRNA as an endogenous control in the quantitative real-time PCR (qRT-PCR) analysis. However, to ensure that mRNA and not pre-mRNA was reverse transcribed, a second cDNA synthesis was performed for a sub-set of samples, using oligo d(T)16 primer. This was used in an attempt to exclude RNA species lacking poly A tails from being reverse transcribed and only measure processed mature polyadenylated mRNA transcripts in the subsequent qRT-PCR reaction.
Detection of mRNA was performed using an ABI-PRISM® 7000 Sequence Detection system (Applied Biosystems, Foster City, CA, USA). Primers and MGB probes (Table 1) were designed using Primer Express software (Applied Biosystems). Primers were synthesized by Invitrogen and MGB probes by Applied Biosystems. A pre-optimized primer and probe assay for 18S rRNA was used as an endogenous control (Applied Biosystems). For cDNA synthesized with oligo d(T)16, β-actin was used as endogenous control. Primers and probes were pre-mixed with TaqMan Universal Master Mix (for probes) or SYBR®GREEN PCR Master Mix (for primers only) (Applied Biosystems) and applied to 96-well MicroAmp Optical plates (Applied Biosystems). cDNA aliquots of 4 μl were added in triplicates. Muscle cDNA was diluted 1:50 for detection of PINK1, 1:10 for detection of svPINK1 and naPINK1 and 1:500 for detection of 18S rRNA. Cell cDNA was diluted 1:2 for detection of all target genes and β-actin and 1:200 for detection of 18S rRNA. Samples were accessed for DNA contamination by including cDNA controls prepared without reverse transcriptase. The amplification of genomic DNA typically amounted to a maximum of <1% of the target gene when using the TRIzol protocol and even less when using RNeasy columns including DNase-treatment. Thermal cycling conditions were: 2 min at 50°C, 10 min at 95°C and 40 cycles of 15 s at 95°C and 1 min at 65°C.
Two-fold dilutions series were performed for all target genes and endogenous controls to determine the amplification efficiency (see Additional file 4). The ΔΔCt method was used to calculate relative changes in mRNA abundance. The threshold cycle (CT) for 18S was subtracted from the CT for the target gene to adjust for variations in the cDNA synthesis. For the in vivo samples, the pre-intervention values reflect baseline gene expression levels and were subtracted from the post mitochondrial biogenesis sample (induced under clinical trial conditions using an endurance training model  to calculate the increase or decrease in mRNA abundance. For the siRNA experiments, 18S rRNA or β-actin adjusted CT-values were related to the adjusted CT-value of the control-siRNA treated sample in each experiment, transformed from logarithmic to linear and then a mean from all experiments was calculated. The results are presented as percentage mRNA abundance of control siRNA treated samples (mean ± SE).
Neuroblastoma cell lines SK-N-MC, SH-SY5Y, SK-N-DZ, SK-N-F1, SK-N-AS, SK-N-Be(2) and SK-N-SH were harvested from 10 cm culture dishes with 2 ml of TRIzol (Invitrogen) and total RNA was isolated as described in the manufacturer's protocol. A pool of totalRNA from all the above cell lines was used for polyadenylated mRNA isolation using an oligo d(T) magnetic bead kit for mRNA isolation (Promega) Total RNA and mRNA was examined by Northern blotting. 10 μg of total RNA were separated on an agarose gel (1.2%) containing 20 mM MOPS (pH 7.0) and 6.7% formaldehyde. Ethidium bromide was added to the gel to examine the distribution of RNA and to determine the abundance of 18S rRNA abundance across the samples. The gel was run for 5 to 6 hours at 80 V and the RNA was then transferred by capillary blotting overnight onto a Hybond-XL membrane (Amersham Biosciences). The membrane was pre hybridized at 42°C for at least 3 h in hybridization solution containing 5 × SSC (1 × SSC is 0.15 M NaCl/0.015 M sodium citrate), 5 × Denhardt's (1 × Denhardt's is 0.02% Ficoll 400/0.02% polyvinylpyrrolidone/0.02% BSA), 0.5% (w/v) SDS, 50 mM sodium phosphate (pH 6.5), 50% formamide and 100 μg/ml sheared herring sperm DNA. Hybridization of the membrane was performed overnight at 42°C in a hybridization solution with the addition of denatured probe, [32P]dCTP-labelled cDNA, corresponding to PINK1 amplified by PCR (primers in table 1) and labeled using a random priming kit (Amersham Biosciences). After hybridization, the solution was removed and the membrane was washed twice for 30 min at room temperature in a solution of 2 × SSC and 0.2% (w/v) SDS and twice for 15 min at 50°C in a solution of 0.1 × SSC and 0.2% (w/v) SDS. The membrane was sealed in a plastic envelope and then exposed and scanned in a Fuji Film Phosphoimager. Specific mRNA signals were analyzed with the Image Gauge V3.45 program (Fuji Film).
5'-RACE and cloning of svPINK1
RNA Ligase Mediated Rapid amplification of cDNA ends (RLM-RACE kit, Ambion) was performed according to the manufacturer's protocol in order to clone svPINK1 full-length cDNA. Briefly, total RNA from SK-N-MC cells (isolated with TRIzol in the same way as for Northern blot according to above) was used as template for the reaction. A 45 base RNA adaptor oligonucleotide was ligated to the 5' end of the de-capped RNA population. The RNA was reverse transcribed using RT-PCR kit from Applied Biosystems instead of the reverse transcription reagents included in the kit, as the latter gave a very weak signal. Reverse transcription PCR was performed in a PTC-225 Peltier Thermal Cycler (MJ Research) at 25°C for 10 minutes, 48°C for 30 minutes and 95°C for 5 minutes. PCR was performed with the inner primer towards the ligated adaptor sequence (provided in RLM-RACE kit, Ambion) and the outer gene-specific primer (Table 1) designed downstream of the predicted stop codon for svPINK1. A second PCR analysis was performed using the inner primer towards the ligated adaptor sequence (provided in RLM-RACE kit, Ambion) and the inner gene-specific primer (table 1) designed downstream of the predicted stop codon for svPINK1 (see Additional file 4). The sequence was confirmed by a third PCR analysis using internal primers towards the previously amplified PCR product. The 5' internal primer targets and thus the product include the start of the svPINK1 sequence, containing the 52 bases unique to this splice variant. A clone (AK075225) containing the svPINK1 sequence was identified by database searching, purchased from National Institute of Technology and Evaluation (NITE) Department of Biotechnology Biological Resource Center (NBRC), Kazusa-kamatari, Kisarazu, Chiba, 292-0818, Japan  and used as template with the internal primers in parallel with the obtained PCR-product from the RACE experiment. The PCR products were validated by agarose gel analysis which demonstrated bands at ~854 bp for both templates, accordingly to the size of the sequence. The product obtained from 5' RACE as well as the product obtained from the AK075225 clone was cloned using standard techniques and sequence analysis confirmed the identity of these sequences.
Human in vivo skeletal muscle mitochondrial experiments
The study was approved by the ethics committee of the Karolinska Institute, Stockholm, Sweden, and informed consent was obtained from each of the volunteers. For stimulating increases in mitochondrial function twenty-four subjects trained under full clinical trial conditions on a cycle ergometer as described previously . Physiological determination of increased aerobic capacity and the muscle biopsy procedure were performed as described [23, 47]. All muscle biopsy samples were divided into several portions and stored under liquid nitrogen until analysis. RNA analysis was as described above. Mitochondrial enzyme activity was analyzed in whole muscle as described previously. Briefly, muscle biopsies were homogenized with a Potter-Elvehjem glass homogenizer in a KCl buffer (100 mM KCl, 1 mM EDTA, 50 mM Tris, 5 mM MgCl2, 1.8 mM ATP, pH 7.2) to obtain a 5% homogenate . For Complex I measurements, 20 μl of whole muscle homogenate was treated with 590 μl disruption solution 1 (5 mM KH2PO4, 5 mM MgCl2, 0.5 g × l-1 essentially fatty acid free bovine serum albumin) and then in 50 μl disruption buffer 2 (7.15 g × l-1 saponin dissolved in disruption buffer 1). Seventy-two μl of the pretreated homogenate was incubated with 75 μl of reagent (50 mM KH2PO4, 5 mM MgCl2, 5 gl-1 BSA, 0.2 mM potassium cyanide, 4.8 μg × ml-1 Antimycine A (Sigma) and 0.12 mM cytochrome c, pH7.5) for 7 min at 35°C. After the incubation, NADH was added to a final concentration of 0.15 mM, and the decrease in absorbance was monitored at 340 nm, for 1 min before and after the addition of rotenone, 2 mg × l-1. The rotenone sensitive activity was calculated. The CV for this method is 8.9%.
Statistical analysis and bioinformatics
All data are presented as mean ± se unless otherwise stated. ANOVA was utilized to analyze the human gene expression responses. When a significant F ratio was achieved, post hoc analysis was utilized to make individual comparisons. For cell based responses, unpaired t-tests were utilized. Sample size and significance level is shown for each graph. The genomic region including the PINK1 and DDOST genes was obtained from Ensembl v38 for both human (Chromosome 1: 20825000 – 20865000 bp) and mouse (Chromosome 4: 137680000 – 137710000 bp). The human sequence was aligned to the reverse complement of the mouse sequence using Kalign . The alignment was analyzed for conservation using Kalignvu. Transcription factor binding sites were predicted using Consite  with a stringent cut-off criteria of 15 bits, using all models in the Jaspar [30, 52] database. A model of the PINK1 protein structure was created using the structure of human protein kinase C theta type (KPCT_HUMAN) in the region 377–625 which is homologous to PINK1 in the region 153–501 .
non-coding natural antisense
natural antisense transcripts
short transcript variant of PINK1
non-coding natural antisense at the PINK1 locus
open reading frame
quantitative real-time PCR
RNA ligase mediated rapid amplification of cDNA ends
short interfering RNA
minimal essential media
We thank Carl Johan Sundberg for collaborations associated with the skeletal muscle studies. The present analysis was supported by the Swedish Diabetes Association (JAT) and the Swedish Research Council (CW).
- Valente EM, Abou-Sleiman PM, Caputo V, Muqit MM, Harvey K, Gispert S, Ali Z, Del Turco D, Bentivoglio AR, Healy DG, Albanese A, Nussbaum R, Gonzalez-Maldonado R, Deller T, Salvi S, Cortelli P, Gilks WP, Latchman DS, Harvey RJ, Dallapiccola B, Auburger G, Wood NW: Hereditary early-onset Parkinson's disease caused by mutations in PINK1. Science. 2004, 304 (5674): 1158-1160. 10.1126/science.1096284.PubMedView ArticleGoogle Scholar
- Silvestri L, Caputo V, Bellacchio E, Atorino L, Dallapiccola B, Valente EM, Casari G: Mitochondrial import and enzymatic activity of PINK1 mutants associated to recessive parkinsonism. Hum Mol Genet. 2005, 14 (22): 3477-3492. 10.1093/hmg/ddi377.PubMedView ArticleGoogle Scholar
- Petit A, Kawarai T, Paitel E, Sanjo N, Maj M, Scheid M, Chen F, Gu Y, Hasegawa H, Salehi-Rad S, Wang L, Rogaeva E, Fraser P, Robinson B, St George-Hyslop P, Tandon A: Wild-type PINK1 prevents basal and induced neuronal apoptosis, a protective effect abrogated by Parkinson disease-related mutations. J Biol Chem. 2005, 280 (40): 34025-34032. 10.1074/jbc.M505143200.PubMedView ArticleGoogle Scholar
- Beilina A, Van Der Brug M, Ahmad R, Kesavapany S, Miller DW, Petsko GA, Cookson MR: Mutations in PTEN-induced putative kinase 1 associated with recessive parkinsonism have differential effects on protein stability. Proc Natl Acad Sci U S A. 2005, 102 (16): 5703-5708. 10.1073/pnas.0500617102.PubMed CentralPubMedView ArticleGoogle Scholar
- Horbinski C, Chu CT: Kinase signaling cascades in the mitochondrion: a matter of life or death. Free Radic Biol Med. 2005, 38 (1): 2-11. 10.1016/j.freeradbiomed.2004.09.030.PubMedView ArticleGoogle Scholar
- Green DR, Kroemer G: The pathophysiology of mitochondrial cell death. Science. 2004, 305 (5684): 626-629. 10.1126/science.1099320.PubMedView ArticleGoogle Scholar
- Cookson MR: The biochemistry of Parkinson's disease *. Annu Rev Biochem. 2005, 74: 29-52. 10.1146/annurev.biochem.74.082803.133400.PubMedView ArticleGoogle Scholar
- Abou-Sleiman PM, Muqit MM, Wood NW: Expanding insights of mitochondrial dysfunction in Parkinson's disease. Nat Rev Neurosci. 2006, 7 (3): 207-219. 10.1038/nrn1868.PubMedView ArticleGoogle Scholar
- Gandhi S, Muqit MM, Stanyer L, Healy DG, Abou-Sleiman PM, Hargreaves I, Heales S, Ganguly M, Parsons L, Lees AJ, Latchman DS, Holton JL, Wood NW, Revesz T: PINK1 protein in normal human brain and Parkinson's disease. Brain. 2006, 129 (Pt 7): 1720-1731. 10.1093/brain/awl114.PubMedView ArticleGoogle Scholar
- Yang Y, Gehrke S, Imai Y, Huang Z, Ouyang Y, Wang JW, Yang L, Beal MF, Vogel H, Lu B: Mitochondrial pathology and muscle and dopaminergic neuron degeneration caused by inactivation of Drosophila Pink1 is rescued by Parkin. Proc Natl Acad Sci U S A. 2006, 103 (28): 10793-10798. 10.1073/pnas.0602493103.PubMed CentralPubMedView ArticleGoogle Scholar
- Clark IE, Dodson MW, Jiang C, Cao JH, Huh JR, Seol JH, Yoo SJ, Hay BA, Guo M: Drosophila pink1 is required for mitochondrial function and interacts genetically with parkin. Nature. 2006, 441 (7097): 1162-1166. 10.1038/nature04779.PubMedView ArticleGoogle Scholar
- Park J, Lee SB, Lee S, Kim Y, Song S, Kim S, Bae E, Kim J, Shong M, Kim JM, Chung J: Mitochondrial dysfunction in Drosophila PINK1 mutants is complemented by parkin. Nature. 2006, 441 (7097): 1157-1161. 10.1038/nature04788.PubMedView ArticleGoogle Scholar
- Wang D, Qian L, Xiong H, Liu J, Neckameyer WS, Oldham S, Xia K, Wang J, Bodmer R, Zhang Z: Antioxidants protect PINK1-dependent dopaminergic neurons in Drosophila. Proc Natl Acad Sci U S A. 2006, 103 (36): 13520-13525. 10.1073/pnas.0604661103.PubMed CentralPubMedView ArticleGoogle Scholar
- Unoki M, Nakamura Y: Growth-suppressive effects of BPOZ and EGR2, two genes involved in the PTEN signaling pathway. Oncogene. 2001, 20 (33): 4457-4465. 10.1038/sj.onc.1204608.PubMedView ArticleGoogle Scholar
- Ensembl Genome Browser. [http://www.ensembl.org]
- Zhang Y, Li J, Kong L, Gao G, Liu QR, Wei L: NATsDB: Natural Antisense Transcripts DataBase. Nucleic Acids Res. 2006Google Scholar
- Katayama S, Tomaru Y, Kasukawa T, Waki K, Nakanishi M, Nakamura M, Nishida H, Yap CC, Suzuki M, Kawai J, Suzuki H, Carninci P, Hayashizaki Y, Wells C, Frith M, Ravasi T, Pang KC, Hallinan J, Mattick J, Hume DA, Lipovich L, Batalov S, Engstrom PG, Mizuno Y, Faghihi MA, Sandelin A, Chalk AM, Mottagui-Tabar S, Liang Z, Lenhard B, Wahlestedt C: Antisense transcription in the mammalian transcriptome. Science. 2005, 309 (5740): 1564-1566. 10.1126/science.1112009.PubMedView ArticleGoogle Scholar
- Yelin R, Dahary D, Sorek R, Levanon EY, Goldstein O, Shoshan A, Diber A, Biton S, Tamir Y, Khosravi R, Nemzer S, Pinner E, Walach S, Bernstein J, Savitsky K, Rotman G: Widespread occurrence of antisense transcription in the human genome. Nat Biotechnol. 2003, 21 (4): 379-386. 10.1038/nbt808.PubMedView ArticleGoogle Scholar
- Lavorgna G, Dahary D, Lehner B, Sorek R, Sanderson CM, Casari G: In search of antisense. Trends Biochem Sci. 2004, 29 (2): 88-94. 10.1016/j.tibs.2003.12.002.PubMedView ArticleGoogle Scholar
- Timmons JA, Good L: Does everything now make (anti)sense?. Biochem Soc Trans. 2006, 34 (Pt 6): 1148-1150.PubMedView ArticleGoogle Scholar
- Willingham AT, Orth AP, Batalov S, Peters EC, Wen BG, Aza-Blanc P, Hogenesch JB, Schultz PG: A strategy for probing the function of noncoding RNAs finds a repressor of NFAT. Science. 2005, 309 (5740): 1570-1573. 10.1126/science.1115901.PubMedView ArticleGoogle Scholar
- Lowell BB, Shulman GI: Mitochondrial dysfunction and type 2 diabetes. Science. 2005, 307 (5708): 384-387. 10.1126/science.1104343.PubMedView ArticleGoogle Scholar
- Timmons JA, Larsson O, Jansson E, Fischer H, Gustafsson T, Greenhaff PL, Ridden J, Rachman J, Peyrard-Janvid M, Wahlestedt C, Sundberg CJ: Human muscle gene expression responses to endurance training provide a novel perspective on Duchenne muscular dystrophy. Faseb J. 2005, 19 (7): 750-760. 10.1096/fj.04-1980com.PubMedView ArticleGoogle Scholar
- Vega Genome Browser. [http://vega.sanger.ac.uk]
- Pieper U, Eswar N, Davis FP, Braberg H, Madhusudhan MS, Rossi A, Marti-Renom M, Karchin R, Webb BM, Eramian D, Shen MY, Kelly L, Melo F, Sali A: MODBASE: a database of annotated comparative protein structure models and associated resources. Nucleic Acids Res. 2006, 34 (Database issue): D291-5. 10.1093/nar/gkj059.PubMed CentralPubMedView ArticleGoogle Scholar
- Expasy Proteomics Server. [http://www.expasy.org]
- FANTOM3: Functional annotation of mouse-3. [http://fantom3.gsc.riken.go.jp]
- NATsDB: Natural Antisense transcripts Database. [http://natsdb.cbi.pku.edu.cn]
- Vlieghe D, Sandelin A, De Bleser PJ, Vleminckx K, Wasserman WW, van Roy F, Lenhard B: A new generation of JASPAR, the open-access repository for transcription factor binding site profiles. Nucleic Acids Res. 2006, 34 (Database issue): D95-7. 10.1093/nar/gkj115.PubMed CentralPubMedView ArticleGoogle Scholar
- Lu C, He JC, Cai W, Liu H, Zhu L, Vlassara H: Advanced glycation endproduct (AGE) receptor 1 is a negative regulator of the inflammatory response to AGE in mesangial cells. Proc Natl Acad Sci U S A. 2004, 101 (32): 11767-11772. 10.1073/pnas.0401588101.PubMed CentralPubMedView ArticleGoogle Scholar
- Hutvagner G: Small RNA asymmetry in RNAi: function in RISC assembly and gene regulation. FEBS Lett. 2005, 579 (26): 5850-5857. 10.1016/j.febslet.2005.08.071.PubMedView ArticleGoogle Scholar
- Fishel MA, Watson GS, Montine TJ, Wang Q, Green PS, Kulstad JJ, Cook DG, Peskind ER, Baker LD, Goldgaber D, Nie W, Asthana S, Plymate SR, Schwartz MW, Craft S: Hyperinsulinemia provokes synchronous increases in central inflammation and beta-amyloid in normal adults. Arch Neurol. 2005, 62 (10): 1539-1544. 10.1001/archneur.62.10.noc50112.PubMedView ArticleGoogle Scholar
- Craft S, Watson GS: Insulin and neurodegenerative disease: shared and specific mechanisms. Lancet Neurol. 2004, 3 (3): 169-178. 10.1016/S1474-4422(04)00681-7.PubMedView ArticleGoogle Scholar
- Ramasamy R, Vannucci SJ, Yan SS, Herold K, Yan SF, Schmidt AM: Advanced glycation end products and RAGE: a common thread in aging, diabetes, neurodegeneration, and inflammation. Glycobiology. 2005, 15 (7): 16R-28R. 10.1093/glycob/cwi053.PubMedView ArticleGoogle Scholar
- Knee R, Murphy PR: Regulation of gene expression by natural antisense RNA transcripts. Neurochem Int. 1997, 31 (3): 379-392. 10.1016/S0197-0186(96)00108-8.PubMedView ArticleGoogle Scholar
- Sleutels F, Zwart R, Barlow DP: The non-coding Air RNA is required for silencing autosomal imprinted genes. Nature. 2002, 415 (6873): 810-813.PubMedView ArticleGoogle Scholar
- O'Neill MJ: The influence of non-coding RNAs on allele-specific gene expression in mammals. Hum Mol Genet. 2005, 14 Spec No 1: R113-20. 10.1093/hmg/ddi108.PubMedView ArticleGoogle Scholar
- Moorwood K, Charles AK, Salpekar A, Wallace JI, Brown KW, Malik K: Antisense WT1 transcription parallels sense mRNA and protein expression in fetal kidney and can elevate protein levels in vitro. J Pathol. 1998, 185 (4): 352-359. 10.1002/(SICI)1096-9896(199808)185:4<352::AID-PATH119>3.0.CO;2-#.PubMedView ArticleGoogle Scholar
- Nishimura Y, Kikis EA, Zimmer SL, Komine Y, Stern DB: Antisense transcript and RNA processing alterations suppress instability of polyadenylated mRNA in chlamydomonas chloroplasts. Plant Cell. 2004, 16 (11): 2849-2869. 10.1105/tpc.104.026203.PubMed CentralPubMedView ArticleGoogle Scholar
- Robb GB, Brown KM, Khurana J, Rana TM: Specific and potent RNAi in the nucleus of human cells. Nat Struct Mol Biol. 2005, 12 (2): 133-137. 10.1038/nsmb886.PubMedView ArticleGoogle Scholar
- Judge AD, Sood V, Shaw JR, Fang D, McClintock K, MacLachlan I: Sequence-dependent stimulation of the mammalian innate immune response by synthetic siRNA. Nat Biotechnol. 2005, 23 (4): 457-462. 10.1038/nbt1081.PubMedView ArticleGoogle Scholar
- Muqit MM, Abou-Sleiman PM, Saurin AT, Harvey K, Gandhi S, Deas E, Eaton S, Payne Smith MD, Venner K, Matilla A, Healy DG, Gilks WP, Lees AJ, Holton J, Revesz T, Parker PJ, Harvey RJ, Wood NW, Latchman DS: Altered cleavage and localization of PINK1 to aggresomes in the presence of proteasomal stress. J Neurochem. 2006, 98 (1): 156-169. 10.1111/j.1471-4159.2006.03845.x.PubMedView ArticleGoogle Scholar
- Nakajima A, Kataoka K, Hong M, Sakaguchi M, Huh NH: BRPK, a novel protein kinase showing increased expression in mouse cancer cell lines with higher metastatic potential. Cancer Lett. 2003, 201 (2): 195-201. 10.1016/S0304-3835(03)00443-9.PubMedView ArticleGoogle Scholar
- Deng H, Jankovic J, Guo Y, Xie W, Le W: Small interfering RNA targeting the PINK1 induces apoptosis in dopaminergic cells SH-SY5Y. Biochem Biophys Res Commun. 2005, 337 (4): 1133-1138.PubMedView ArticleGoogle Scholar
- Akimoto T, Pohnert SC, Li P, Zhang M, Gumbs C, Rosenberg PB, Williams RS, Yan Z: Exercise stimulates Pgc-1alpha transcription in skeletal muscle through activation of the p38 MAPK pathway. J Biol Chem. 2005, 280 (20): 19587-19593. 10.1074/jbc.M408862200.PubMedView ArticleGoogle Scholar
- National Institute of Technology and Evaluation. [http://www.nbrc.nite.go.jp]
- Timmons JA, Gustafsson T, Sundberg CJ, Jansson E, Hultman E, Kaijser L, Chwalbinska-Moneta J, Constantin-Teodosiu D, Macdonald IA, Greenhaff PL: Substrate availability limits human skeletal muscle oxidative ATP regeneration at the onset of ischemic exercise. J Clin Invest. 1998, 101 (1): 79-85.PubMed CentralPubMedView ArticleGoogle Scholar
- Fredriksson K, Radell P, Eriksson LI, Hultenby K, Rooyackers O: Effect of prolonged mechanical ventilation on diaphragm muscle mitochondria in piglets. Acta Anaesthesiol Scand. 2005, 49 (8): 1101-1107. 10.1111/j.1399-6576.2005.00718.x.PubMedView ArticleGoogle Scholar
- Lassmann T, Sonnhammer EL: Kalign--an accurate and fast multiple sequence alignment algorithm. BMC Bioinformatics. 2005, 6: 298-10.1186/1471-2105-6-298.PubMed CentralPubMedView ArticleGoogle Scholar
- Sandelin A, Wasserman WW, Lenhard B: ConSite: web-based prediction of regulatory elements using cross-species comparison. Nucleic Acids Res. 2004, 32 (Web Server issue): W249-52. 10.1093/nar/gkh372.PubMed CentralPubMedView ArticleGoogle Scholar
- Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004, 32 (Database issue): D91-4. 10.1093/nar/gkh012.PubMed CentralPubMedView ArticleGoogle Scholar