Association of Mitochondrial DNA Polymerase γ Gene POLG1 Polymorphisms with Parkinsonism in Chinese Populations

Background Mitochondrial DNA polymerase gamma (POLG1) mutations were associated with levodopa-responsive Parkinsonism. POLG1 gene contains a number of common nonsynonymous SNPs and intronic regulatory SNPs which may have functional consequences. It is of great interest to discover polymorphisms variants associated with Parkinson's disease (PD), both in isolation and in combination with specific SNPs. Materials and Methods We conducted a case-control study and genotyped twenty SNPs and poly-Q polymorphisms of POLG1 gene including in 344 Chinese sporadic PD patients and 154 healthy controls. All the polymorphisms of POLG1 we found in this study were sequenced by PCR products with dye terminator methods using an ABI-3100 sequencer. Hardy-Weinberg equilibrium and linkage disequilibrium (LD) for association between twenty POLG1 SNPs and PD were calculated using the program Haploview. Principal Results We provided evidence for strong association of four intronic SNPs of the POLG1 gene (new report: c.2070-12T>A and rs2307439: c.2070-64G>A in intron 11, P = 0.00011, OR = 1.727; rs2302084: c.3105-11T>C and rs2246900: c.3105-36A>G in intron 19, P = 0.00031, OR = 1.648) with PD. However, we did not identify any significant association between ten exonic SNPs of POLG1 and PD. Linkage disequilibrium analysis indicated that c.2070-12T>A and c.2070-64G>A could be parsed into one block as Haplotype 1 as well as c.3105-11T>C and c.3105-36A>G in Haplotype 2. In addition, case and control study on association of POLG1 CAG repeat (poly-Q) alleles with PD showed a significant association (P = 0.03, OR = 2.16) of the non-10/11Q variants with PD. Although intronic SNPs associated with PD didn't influence POLG1 mRNA alternative splicing, there was a strong association of c.2070-12T>A and c.2070-64G>A with decreased POLG1 mRNA level and protein levels. Conclusions Our findings indicate that POLG1 may play a role in the pathogenesis of PD in Chinese populations.


Introduction
Parkinson's disease (PD) is one of the most frequent neurodegenerative disorders caused by loss of dopaminergic neurons in the substantia nigra [1]. 5-10% of PD cases begin at onset of the ages between 20 and 50, which is classified as young onset [2]. The majority of PD patients are sporadic cases [2]. Approximately 5-10% of patients have genetic factors, yet the etiology of PD remains unclear [3,4]. Mitochondrial dysfunction was implicated in the pathogenesis of sporadic, idiopathic PD [5,6,7]. Accumulation of somatic mitochondrial DNA (mtDNA) deletions was observed in the substantia nigra in PD and several gene products of familial Parkinsonism, e.g., mitochondrial DNA polymerase gamma gene (POLG1), had a connection to mitochondrial function [2,5,7].
Numerous evidences indicated that POLG1 mutations cosegregated with levodopa-responsive Parkinsonism in some families, suggesting that POLG1 could be responsible for some mendelian transmission of Parkinsonism [5,6,8,9]. POLG1 encodes part of the heterotrimeric mitochondrial DNA polymerase consisting of a catalytic subunit [5,10], which is essential for mtDNA replication and repair due to its 59 to 39 polymerase activity and 39 to 59 exonuclease activity. POLG1 mutations lead to uncorrect nucleotide incorporation in mtDNA, which resulted in dysfunction of the respiratory chain [11]. Homozygous knock-in mice expressing a polymerase deficiency mutant of POLG1 shown an phenotype with increased amounts of point mutations and deleted mtDNA [12]. In some cases, mitochondrial DNA primary genetic abnormalities or secondary rearrangements due to POLG1 gene mutation could directly cause Parkinsonism [5,9]. Two missense mutations in POLG1, G737R and R853W, were reported in a family with early onset Parkinsonism [13]. Another W748S missense mutation was found in a Finnish patient with a mixed disease background of Parkinsonism, chronic progressive external ophthalmoplegia, and peripheral neuropathy [14]. POLG1 contains an N-terminal polymorphic polyglutamine tract encoded by a (CAG) n tract in exon 2 which may have functional consequences [6,15]. Increased frequency of rare alleles of the POLG1 CAG-repeat (poly-Q) was found in Finnish idiopathic apparently sporadic PD patients, but conflicting reports exist. Therefore, the gene for POLG1 may represent an ideal candidate gene for Parkinsonism susceptibility. If susceptibility alleles are common in the general population and have rather small effects, association can be a more sensitive method than linkage to identify these loci. It is therefore of great interest to discover specific polymorphism variants of POLG1 associated with PD.
Here we scrutinized the genetic polymorphisms in POLG1 gene to evaluate it as potential candidate gene for PD genetic study. We performed extensive screening of POLG1 gene by direct sequencing to detect polymorphisms, and statistical analysis to examine the genetic effects on PD in Chinese populations.

Materials and Methods
Sample collection and DNA/RNA preparations A total of 498 study subjects were included in this study, comprising 344 patients with PD (175 males and 169 females) and 154 control subjects without evidence of PD (84 males and 70 females). We recruited all patients from the Department of Neurology from Sir Run Run Shaw hospital affiliated with Zhejiang University School of Medicine and received a standard neurological examination as well as a psychiatric interview. The PD patients with a mean age of 58.3869.32 years (age at onset 54.2268.52). The clinical diagnosis of PD was confirmed by the senior neurologist specializing in movement disorders according to the Consensus Statement of the Movement Disorders Society in 1998 [16]. The healthy control subjects (mean age of 59.57610.86 years old) were free of neurological or psychiatric disorders. Gender and age proportion between cases and controls was statistically matched. All subjects were with Chinese Han ethnic background. Informed consent for participation in the study was obtained from either directly or from his or her guardian in all subjects and the work received approval from the institution ethics committee and conformed to the tenets of the Declaration of Helsinki. After informed consent was obtained from all participating individuals, blood samples were collected for genomic RNA/DNA extraction using the QIAamp Blood kit (Qiagen). To avoid genomic DNA contamination in RNAs, samples were treated with DNase I (Qiagen). All samples were tested by realtime PCR, showing no amplification after 40 cycles in cDNA preparation without RT.
Another group of 110 patients with sporadic early-onset Parkinsonism was studied. All the patients were from Hubei province, Central China. There were 62 male and 48 female sporadic patients: mean age at onset, 32.766.4 years. All patients underwent a standardized neurological examination by two movement disorder specialists. 110 PD patients were selected according to the following criteria: (1) at least two of the three cardinal motor signs (resting tremor, bradykinesia, rigidity); (2) excellent response following L-Dopa therapy; (3) an age at onset #40 years; (4) Absence of extensor plantar reflexes, ophthalmoplegia, early dementia, or early autonomic failure in the family members. A control group of 130 healthy individuals from the same geographic areas was obtained. All controls were free of symptoms suggestive of PD and with a negative family history of movement disorders. Informed consent was obtained before participation into our study.

PCR and SNP screening by sequencing
To find common SNPs, we designed specific primer sets based on the POLG1 sequences (NM_002693) in the GenBank DNA database to screen the exons and the flanking regions, including the promoter regions (1.5 kb) of these genes. The first-step SNP discovery screening panel included 15 affected and 10 unaffected individuals. PCR conditions were set as follows: 2 min initial denaturation at 95uC, 30 cycles of 30 s denaturation at 94uC, 30 s annealing at 53-60uC, 30 s extension at 72uC and 10 min final extension at 72uC. Before sequencing, amplification products were incubated with shrimp alkaline phosphatase (0.5 units; Amersham) and exonuclease I (1 units; Amersham) at 37uC for 30 min, followed by heat inactivation at 80uC for 15 min. The DNA samples containing extension products were analyzed using the sequencing kit ABI PRISM dye terminator cycle and the automatic sequencing system ABI 3100 (Applied Biosystems).
Twenty SNPs within the POLG1 gene were found in twenty-five individuals, including ten exonic SNPs and ten intronic SNPs (Table 1). There is no other polymorphism with a frequency over 4% in all exons and flanking regions of POLG1. We further conducted a case-control study on these twenty SNPs of POLG1 gene from the subjects, including 344 affected PD patients and 154 unaffected individuals, by PCR-direct sequencing using the ABI-3100 sequencer. The twenty-five samples in the SNP discovery panel were included in the case-control study and calculated in the final significance analysis. Specific primers were designed based on genomic sequences obtained from the GenBank DNA database.

Sequence variation analysis
We aligned multiple sequences of POLG1 from different species using the Clustal W program.
Isolation of peripheral blood lymphocytes from patients' blood samples carefully poured onto the Ficoll solution. Harvest the ring with white blood cells without touching the Ficoll using a sterile pipette tips after centrifugation of the tubes 20 min at 1600 rpm. We diluted the peripheral blood lymphocytes with PBS and washed them twice in PBS. Peripheral blood lymphocytes were cultured in RPMI1640 medium supplemented with 10% (v/v) FBS (Invitrogen).

Quantitative mRNA levels of POLG1
The mRNA levels of POLG1 in human peripheral blood lymphocytes were measured with real-time PCR using genespecific primers as follows: 5-aagcactgtctcgaacaggg-3 (sense); 5-cactgcagctcgcaagttct-3 (antisense); For real-time PCR, 1 mL of cDNA was amplified for 45 cycles with a master mix (SYBR Green Supermix; ABI) using a thermo cycler (MJ). Melting curve analysis was done at the end of the reaction (after 45 cycles) to assess the quality of the final PCR products. The threshold cycle C(t) values were calculated by fixing the basal fluorescence at 0.05 unit. Three replicates were used for each sample, and the average C(t) value was calculated. The DC(t) values were calculated as C(t) sample -C(t) GAPDH (or b-Actin). The N-fold increase or decrease in expression was calculated by the DDCt method using the C(t) value as the reference point. MtDNA copy number quantification The measure of mtDNA copy number was performed on peripheral blood lymphocytes by quantitative RT-PCR based on SYBR-Green I fluorescence. Similar to previous analysis [17], the amount of mtDNA was compared with the amount of the nuclear gene encoding the b-globin gene, contained in the same sample. The fluorescent signal intensities were recorded and analyzed during PCR in the thermo cycler. The mtDNA/b-globin gene DNA ratio obtained in the patient's samples was expressed as a percentage of the mean value obtained in control samples, which represent the 100% value.

SDS-PAGE and Western blot analysis
Patients' blood cell lysates were mixed with equal volume of Laemmli buffer (62.5 mM Tris-HCl pH 6.8, 2% sodium dodecyl sulfate (SDS), 50 mM Dithiothreitol (DTT), 10% glycerol, 0.01% bromophenol blue), boiled for 3 min at 100uC. Electrophoresis of blood cell extracts from PD patients and health controls was done under denaturing reducing conditions using 4%-12% Nupage gels with Mops-SDS running buffer (Invitrogen). Approximately 25 mg of protein was used per lane for analysis as described in the protocol [18]. Electrophoresis was for approximately 60 minutes at 150 V. After electrophoresis, the proteins were transferred to nitrocellulose membranes (Amersham) using the Surelock X-cell transfer system (Invitrogen). Then the membrane was blocked in TBS/T buffer (16 TBS containing 0.05% Tween 20) with 5% BSA and probed with appropriate primary antibody: antibody against POLG1 (Abcam Catalog#Ab128899, 1:2000 dilution) and Rabbit secondary antibodies (Abcam, 1:2000 dilution) coupled to horseradish peroxidase. The blots were washed in TBS/T, incubated using the ECL kit (Cell Signaling Technology), and exposed to X-ray film (Kodak). GAPDH/PORIN levels were detected as loading controls by antibodies against GAPDH and PORIN (Abcam).

Statistical analysis
Genetic association and the Hardy-Weinberg equilibrium for the distribution of genotypes were tested by x 2 analysis with Yates' correction. Odds ratios (ORs) were calculated with 95% confidence intervals (CI). The pairwise linkage disequilibrium (LD) coefficient r2 was calculated using the program Haploview from the Broad Institute website (http://www.broad.mit.edu/mpg/ haploview/index.php) and German Research Center for Environmental Health's website (http://ihg2.helmholtz-muenchen.de/ ihg/snps.html). P less than 0.05 was considered statistically significant. We defined pairs to be in strong LD if one side upper 95% confidence bound on D9 is .0.98 and lower bound is .0.7. Fraction of strong LD in informative comparisons must be at least 0.95.

POLG1 exonic and intronic SNP screening
By direct DNA sequencing in twenty-five individuals, we identified twenty SNPs within the POLG1 gene, including ten exonic SNPs and ten intronic SNPs. Among twenty SNPs identified, seventeen SNPs were already registered in the dbSNP database, and three SNPs are novel (Table 1, Figure 1A). Twenty identified SNPs were selected for larger scale genotyping in 344 PD patients and 154 health controls for PD genetic association study based on frequencies and/or locations.

Association of exonic and intronic SNPs with PD
Cases and controls were genotyped for ten exonic SNPs including the five coding region mis-sense substitutions previously investigated (rs56410699: c.722C.T-P241L; rs145289229: c1174C.G-L392V; rs185645212: c.2165G.A-R722H; rs41549716: c.2492A.G-Y831C; and rs3179578: c.3689T.C-S1230F), and four alleles previously described in the UK population (rs113994095: c.1399G.A-A467T; rs113994097: c.2243G.C-W748S; rs113994099: c.2864A.G-Y955C; and rs2307441: c.3708G.T-Q1236H), and novel reported c.1830G.A-E516E. All of the exonic variant alleles were in Hardy Weinberg equilibrium. We did not identify any significant association between these ten exonic SNPs and PD (P.0.1, Table 2). To determine whether the alleles were acting in combination, we compared compound genotypes in PD cases and controls. There was no difference in the frequency of compound genotypes defined by L392V, W748S, Y831C, S1230F, and Q1236H (data not shown). However, four intronic SNPs were significantly associated with PD among ten intronic SNPs of the POLG1 gene ( Table 3) Figure 2A). Furthermore, a comparison of allelic frequencies in the haplotype 1 indicated that ht1-AA was significantly more frequent in PD individuals than in the healthy control subjects (P = 0.00011). Another pair of polymorphisms between rs2302084: c.3105-11T.C and rs2246900: c.3105-36A.G in intron 19 was detected as significant association with PD (Haplotype 2 (ht2); P = 0.00031, OR = 1.648) ( Figure 2B). To investigate the pattern of LD in the POLG1 locus, pairwise LD was measured by D9 values among nine intronic polymorphisms ( Figure 1B). Linkage disequilibrium analysis indicated that c.2070-12T.A and c.2070-64G.A in ht1 as well as c.3105-11T.C and c.3105-36A.G in ht2 could be parsed into two different blocks, which had a strong LD spine. Furthermore, we confirmed our association analysis between the two intronic haplotypes and PD in an independent case-control sample consisting of 110 individuals with PD and 130 normal controls from the Hubei province of China. The statistical Number of valid subjects who were successfully genotyped for each of SNP. 3 Analysis performed by a 262 table for each SNP using major homozygotes vs. others in cases and controls. 4 Reference group (controls) designated with an OR of 1.00.  Table 4). These results demonstrated that intronic SNPs rather than exonic SNPs of POLG1 were associated with PD development in Chinese Parkinsonism populations. Furthermore, c.2070-12T.A and c.2070-64G.A substitution in intron 11 as well as c.3105-11T.C and c.3105-36A.G in intron 19 occur highly conserved ( Figure 3A and 3B). These intronic substitutions just happened to be nearby the pyrimidinerich region of the 39 splicing acceptor. Mutations within the pyrimidine-rich intronic sequence could cause human diseases [19,20,21]. Therefore, it is of great interest and importance to examine the effect of these four polymorphisms on POLG1 RNA splicing.

Association of POLG1 poly-Q alleles with PD
In order to study whether the possible effect of length alleles of the POLG1 poly-Q tract in 498 Chinese PD patients vs. controls, we amplified POLG1 exon 2 and performed sequencing on an ABI 3100 automatic sequencer. We found variations of the poly-Q tract ranging from 6Q up to 14Q in our recruited materials ( Table 5). The genotype distributions for the two most common alleles (10Q and 11Q) did not deviate from Hardy-Weinberg equilibrium in PD patients or healthy controls. For statistical analysis the alleles were divided into two main groups; 10/11Q variants (10 and 11 repeats) and non-10/11Q variants (6, 7, 8, 9, 12, 13 and 14 repeats). We found a significant association (P = 0.03, OR = 2.16) of the non-10/11Q variants with PD (Table 6). These results indicated that non-10Q/11Q viriants of POLG1 gene other than the conserved 10/11Q allele might increase susceptibility to PD in Chinese populations.
Intronic SNPs associated with PD don't influence POLG1 mRNA alternative splicing Substantial evidences indicated that misregulation of alternative splicing was associated with human diseases, including PD [22,23]. To explore the molecular mechanism of association between intronic SNP variants of POLG1 and PD, the ability of the intron to influence pre-mRNA splicing was tested.
C.2070-12T.A and c.2070-64G.A in the intron 11 were located in the 64 bp upstream of exon 12, and haplotype 2 (c.3105-11T.C and c.3105-36A.G in intron 19) was located in the 36 bp upstream of exon 20. Both haplotype SNPs are within several serine/arginine-rich protein-binding motifs ( Figure 4A and 4B). Among these, a predicted SF2/ASF-binding site implicated in splicing scored highest (ESE finder), which was deleted by the T allele instead of the C allele (CAGCGTA to CAGTGTA) of haplotype 1 in intron 11 of POLG1. However, RT-PCR amplification of transcripts from exon 11 to exon 12 as well as exon 19 to exon 20 of POLG1 did not reveal any splice variants, regardless of C or T alleles ( Figure 4C and 4D). This result argues against altered splicing as primary mechanism, because if intron 11 and intron 19 were to affect pre-mRNA splicing, we would expect to see different allelic RNA ratios in the blood samples of these patients carrying haplotype 1 of c.2070-12T.A and c.2070-64G.A as well as haplotype 2 of c.3105-11T.C and c.3105-36A.G. Therefore, intronic polymorphism variants associated with PD couldn't influence POLG1 mRNA splicing. However, because of lack of primary brain tissues from the patients, we could not further examine the effect of intronic SNPs associated with PD on RNA splicing.

Intronic SNPs associated with PD affect POLG1 mRNA and protein levels
In addition to alternative splicing, mutations within the pyrimidine-rich intronic sequence could also affect mRNA transcription [22,23]. Although POLG1 expression varies widely among individuals, the contribution of genetic factors remains uncertain. Here we isolated 74 peripheral blood lymphocytes from blood samples of PD patients (six PD patients or healthy controls carrying three different genotypes of haplotype 1 or haplotype 2) to measure POLG1 mRNA levels by quantitative RT-PCR. POLG1 mRNA levels from six PD patients' peripheral blood lymphocytes with the main G/G+T/T genotype of haplotype1 in intron 11 had 1.7 fold (95%CI: 1.1-2.8) higher levels than A/A+A/A carriers combined (P = 0.003 in PD patients and P = 0.004 in healthy controls), with no interactions between genotypes and sex ( Figure 5A). To further normalize the POLG1 mRNA levels we observed, another endogenous reference control, b-Actin, was used to detect similar decreasing effect on POLG1 mRNA expression in haplotype 1-A/A+A/A genotype carriers ( Figure S1A and S1B). After adjusting for sex ratio in each genotype, intronic c.2070-12T.A and c.2070-64G.A polymorphisms in intron 11 remained significantly associated with decreased POLG1 expression (1.67-fold G/G+T/T over G/A+T/A and A/A+A/A, 95%CI: 1.11-2.46, P,0.014) (data not shown). Decreased POLG1 levels in whole blood protein levels was also detected using 4 G/ G+T/T carriers and 4 PD patients carrying A/A+A/A genotype in haplotype 1 ( Figure 5C). We also detected the effect of haplotype 2 (c.3105-11T.C and c.3105-36A.G) in intron 19 on POLG1 mRNA transcription in peripheral blood lymphocytes and protein levels in whole blood cells, and found there was no significant association between POLG1 mRNA levels, protein levels and different genotypes of haplotype 2 ( Figure 5B and 5D). Furthermore, isolated peripheral blood lymphocytes were cultured for three days in RPMI1640 medium supplemented with 10% FBS. We found a similar decrease of POLG1 mRNA levels in haplotype 1 A/A+A/A in intron 11 of the POLG1 gene (P = 0.003 in PD patients; and P = 0.002 in healthy controls, Figure 5E and F). There were no other mutations of POLG1 gene found and no intronic region methylation in the same region found by sequencing and bisulfate-sequencing PCR (unpublished data). We concluded that intronic c.2070-12T.A and c.2070-64G.A in the intron 11 were associated with decreased POLG1 mRNA level and protein levels in human blood lymphocytes. However, POLG1 expression maybe differently regulated in different human tissues such as blood and brain tissue, differential POLG1 levels due to intronic polyporphisms still need to be confirmed by primary brain tissues. Given the significant reduction in POLG1 expression in ht1-A/A+A/A carriers, we further investigated the potential effect on mtDNA copy number. The results of qRT-PCR showed a Table 4. Analysis of association of POLG1 intronic haplotypes with PD in an independent case-control sample. reduction of the mtDNA/b-globin gene DNA ratio of 69%68% (P = 0.009 in ten PD patients carrying ht1-A/A+A/A and ht1-G/ G+T/T) and 73%68% (P = 0.011 in ten healthy controls carrying ht1-A/A+A/A and ht1-G/G+T/T) in ht1-A/A+A/A carriers compared with ht1-G/G+T/T carriers in both PD patients and age-matched control individuals respectively ( Figure 5G). These data indicated that c.2070-12T.A and c.2070-64G.A in the intron 11 changed mtDNA replication.

Discussion
Since POLG1 is essential for mtDNA stability, dysfunction of which has been demonstrated in PD patients, it is a good  candidate for PD susceptibility gene. Numerous previous association studies between POLG1 polymorphisms and PD were reported in several ethnic groups, but the results were inconsistent. Parkinsonism-associated POLG1 variant Y831C was found to be associated with PD in Finnish populations [14]. Recently two genome-wide association studies (GWAS) performed in a North American, UK [24] and German material of Caucasian origin with sporadic PD and North American material of Caucasian origin with familial PD [25], screened for a missense polymorphism rs3087374 (Q1236H) in POLG1, but no significant association with PD was found. In this study, we identified two POLG1 intronic haplotypes containing two neighboring SNPs significantly associated with Parkinsonism in Chinese populations. To our knowledge, this is the first report regarding association between intronic haplotype SNPs of POLG1 and PD. We examined twenty POLG1 SNPs on the association with PD and found that four polymorphisms (c.2070-12T.A and c.2070-64G.A in intron 11; c.3105-11T.C and c.3105-36A.G in intron 19) were significantly associated with PD. The association between haplotype 1,2 and PD was further confirmed in an independent case-control samples consisting of 110 individuals with PD and 130 normal controls from the Hubei province of China. However, there was no significant association between PD and ten exonic SNPs. These results were consistant with previous report [26]. Therefore, no additional SNP was detected in other introns of the gene that could account for the observed associations in the population, indicating that disease risk was not associated with SNPs in other parts of POLG1 gene, such as the 59 promoter region, the 39UTR, or the other exons. Furthermore, we found a significant association (P = 0.03, OR = 2.16) of the non-10/11Q variants with PD from a further case and control study on association of POLG1 poly-Q alleles. Our results supported previous finding that POLG1 poly-Q alleles other than the conserved 10/11Q allele may increase susceptibility to PD in Chinese populations [6,15]. These results indicated that POLG1 polymorphisms may exert a substantial influence on the inheritance of Parkinsonism in Chinese populations.
Although positive associations of exonic polymorphisms or mutation screening of POLG1 with PD were reported, we failed to observe any association with susceptibility to PD for exonic polymorphisms of POLG1. Several factors could provide plausible explanations, e.g., (1) possible genetic differences between different  (2) Given the limited size of our study cohort, especially normal controls (n = 154), we are not able to exclude the possibility that rare POLG1 alleles are associated with a small increased risk of the disorder; (3) different study designs and possible biases in recruiting case and controls among studies; (4) and the statistical approach in the previous study did not account for multiple significance testings, raising the possibility of a false positive result. Taking all together, our observations indicated that POLG1 amino acid substitutions were of course recessive and did not cause PD, but cause other severe neurological phenotypes. Exonic splicing enhancers (ESEs), intronic splicing enhancers (ISEs), exonic splicing silencers (ESSs), and intronic splicing silencers (ISSs) were the known cis-acting sequence elements for promoting splice-site recognition [22,27,28]. These sequences tend to be short (typically ,5-10 nt in length) and consist of relatively degenerate consensus sequences that are recognized by transacting splicing factors, such as serine/arginine-rich domain proteins (SR proteins) and nuclear ribonucleoprotein proteins (hnRNPs), which can regulate the use of a splice site [23,29]. Intronic SNPs in the intron-exon boundaries were reported to have consequences for transcription factor binding sites and native splicing sites [27,28,30]. We initially studied possible alternative splicing of POLG1 mRNA. Semiquantitative RT-PCR indicated that endogenous E12 and E19 were not skipped by alternative splicing in blood from patients carrying c.2070-12T.A and c.2070-64G.A substitution in intron 11, or c.3105-11T.C and c.3105-36A.G in intron 19. Further study indicated that intronic c.2070-12T.A and c.2070-64G.A of POLG1 associated with decreased its mRNA level and protein levels. However, mechanistic studies with minigenes of POLG1 in transfected cells will be needed to discover multifactorial regulation of alternative splicing, including identifying positive and negative regulatory motifs in exon 12, exon 19, and in flanking introns. There are many reasons for intronic variants regulating mRNA transcription. The SF2/ ASF protein was reported to have additional functions beyond splicing, such as regulating translation and stabilizing mRNA [31]. To test whether SF2/ASF was involved in intronic SNPs regulated POLG1 mRNA, the expression of SF2/ASF mRNA in patients' blood was measured using real-time PCR, to account for tissue-specific effects of intronic SNPs on POLG1 expression. However, similar SF2/ASF expression levels were detected in PD patients and health controls carrying c.2070-12T.A and c.2070-64G.A, or c.3105-11T.C and c.3105-36A.G of POLG1, respectively (data not shown). Further study will be needed for the intronic regions containing regulatory elements, like enhancer/attenuator elements that regulate transcription.
In summary, we examined the genetic association of POLG1 with the risk of PD, and found the association of intronic haplotype variants of POLG1 significantly linked to PD in Chinese populations. The results of this study could be helpful in understanding the important function of POLG1 in PD development and in developing drugs to treat PD. Figure S1 Another endogenous reference control, b-Actin, is used to detect changes in POLG1 expression. (A) Quantitative RT-PCR of POLG1 mRNA from patients' peripheral blood lymphocytes with the main G/G+T/T, A/ A+A/A genotypes and heterogeneous G/A+T/A genotype carriers of haplotype 1. (B) Quantitative RT-PCR of POLG1 gene shown there is no effect of haplotype 2 genotypes on POLG1 mRNA transcription. (TIF) Figure 5. Intronic SNPs of POLG1 gene associate with decreased POLG1 mRNA and protein levels. (A) Quantitative RT-PCR of POLG1 mRNA from patients' peripheral blood lymphocytes with the main G/G+T/T, A/A+A/A genotypes and heterogeneous G/A+T/A genotype carriers of haplotype 1. (B) Quantitative RT-PCR of POLG1 mRNA shown there was no effect of haplotype 2 genotypes on POLG1 mRNA transcription. Panel C and D shows western blot of POLG1 protein levels from eight age and gender-matched patients carrying haplotype 1 genotypes (G/G+T/T, A/A+A/A) as well as six patients carrying haplotype 2 genotypes (A/A+T/T, G/G+C/C). PORIN protein was detected as a mitochondrial loading control. The signal ratio was obtained through comparing POLG1 protein band intensity with GAPDH protein intensity. (E and F) POLG1 mRNA was detected by quantitative RT-PCR from isolated and cultured patients' peripheral blood lymphocytes from different individuals' carrying with different genotypes of haplotype 1 or haplotype 2. (G) Intronic polymorphisms effect of POLG1 haplotype 1 and haplotype 2 on mtDNA copy number. doi:10.1371/journal.pone.0050086.g005