Genetic variation in PARL influences mitochondrial content

Genome-wide association study of mitochondrial copy number

Mitochondrial DNA copy number (mtDNAcn) variation has been associated with increased risk of several human diseases in epidemiological studies. The quantification of mtDNAcn performed with real-time PCR is currently considered the de facto standard among several techniques. However, the heterogeneity of the laboratory methods (DNA extraction, storage, processing) used could give rise to results that are difficult to compare and reproduce across different studies. Several lines of evidence suggest that mtDNAcn is influenced by nuclear and mitochondrial genetic variability, however this relation is largely unexplored. The aim of this work was to elucidate the genetic basis of mtDNAcn variation. We performed a genome-wide association study (GWAS) of mtDNAcn in 6836 subjects from the ESTHER prospective cohort, and included, as replication set, the summary statistics of a GWAS that used 295 150 participants from the UK Biobank. We observed two novel associations with mtDNAcn variation on chromosome 19 (rs117176661), and 12 (rs7136238) that reached statistical significance at the genome-wide level. A polygenic score that we called mitoscore including all known single nucleotide polymorphisms explained 1.11% of the variation of mtDNAcn (p = 5.93 × 10-7). In conclusion, we performed a GWAS on mtDNAcn, adding to the evidence of the genetic background of this trait.

Genetic variants in nuclear DNA along with environmental factors modify mitochondrial DNA copy number: a population-based exome-wide association study

Background

Mitochondrial DNA (mtDNA) copy number has been found associated with multiple diseases, including cancers, diabetes and so on. Both environmental and genetic factors could affect the copy number of mtDNA. However, limited study was available about the relationship between genetic variants and mtDNA copy number. What’s more, most of previous studies considered only environmental or genetic factors. Therefore, it’s necessary to explore the genetic effects on mtDNA copy number with the consideration of PM_2.5 exposure and smoking.

Results

A multi-center population-based study was performed with 301 subjects from Zhuhai, Wuhan and Tianjin. Personal 24-h PM_2.5 exposure levels, smoking and mtDNA copy number were evaluated. The Illumina Human Exome BeadChip, which contained 241,305 single nucleotide variants, was used for genotyping. The association analysis was conducted in each city and meta-analysis was adopted to combine the overall effect among three cities. Seven SNPs showed significant association with mtDNA copy number with P value less than 1.00E-04 after meta-analysis. The following joint analysis of our identified SNPs showed a significant allele-dosage association between the number of variants and mtDNA copy number (P = 5.02 × 10^− 17). Further, 11 genes were identified associated with mtDNA copy number using gene-based analysis with a P value less than 0.01.

Conclusion

This study was the first attempt to evaluate the genetic effects on mtDNA copy number with the consideration of personal PM_2.5 exposure level. Our findings could provide more evidences that genetic variants played important roles in modulating the copy number of mtDNA.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5142-7) contains supplementary material, which is available to authorized users.

Single-Cell RNA-Seq of Mouse Dopaminergic Neurons Informs Candidate Gene Selection for Sporadic Parkinson Disease

Genetic variation modulating risk of sporadic Parkinson disease (PD) has been primarily explored through genome-wide association studies (GWASs). However, like many other common genetic diseases, the impacted genes remain largely unknown. Here, we used single-cell RNA-seq to characterize dopaminergic (DA) neuron populations in the mouse brain at embryonic and early postnatal time points. These data facilitated unbiased identification of DA neuron subpopulations through their unique transcriptional profiles, including a postnatal neuroblast population and substantia nigra (SN) DA neurons. We use these population-specific data to develop a scoring system to prioritize candidate genes in all 49 GWAS intervals implicated in PD risk, including genes with known PD associations and many with extensive supporting literature. As proof of principle, we confirm that the nigrostriatal pathway is compromised in Cplx1-null mice. Ultimately, this systematic approach establishes biologically pertinent candidates and testable hypotheses for sporadic PD, informing a new era of PD genetic research.

Common variants in the PARL and PINK1 genes increase the risk to leprosy in Han Chinese from South China

Leprosy is a chronic infectious and neurological disease caused by Mycobacterium leprae, an unculturable pathogen with massive genomic decay and dependence on host metabolism. We hypothesized that mitochondrial genes PARL and PINK1 would confer risk to leprosy. Thirteen tag SNPs of PARL and PINK1 were analyzed in 3620 individuals with or without leprosy from China. We also sequenced the entire exons of PARL, PINK1 and PARK2 in 80 patients with a family history of leprosy by using the next generation sequencing technology (NGS). We found that PARL SNP rs12631031 conferred a risk to leprosy (P_adjusted = 0.019) and multibacillary leprosy (MB, P_adjusted = 0.020) at the allelic level. rs12631031 and rs7653061 in PARL were associated with leprosy and MB (dominant model, P_adjusted < 0.05) at the genotypic level. PINK1 SNP rs4704 was associated with leprosy at the genotypic level (P_adjusted = 0.004). We confirmed that common variants in PARL and PINK1 were associated with leprosy in patients underwent NGS. Furthermore, PARL and PINK1 could physically interact with each other and were involved in the highly connected network formed by reported leprosy susceptibility genes. Together, our results showed that PARL and PINK1 genetic variants are associated with leprosy.

PARL: The mitochondrial rhomboid protease

The rhomboid family comprises evolutionary conserved intramembrane proteases involved in a wide spectrum of biologically relevant activities. A mitochondrion-localized rhomboid, called PARL in mammals, and conserved in yeast and Drosophila as RBD1/PCP1 and rho-7, respectively, plays an indispensable role in cell homeostasis as illustrated by the severe phenotypes caused by its genetic ablation in the various investigated species. Although several substrates of PARL have been proposed to explain these phenotypes, there remains a lot of controversy in this important area of research. We review here the putative functions and substrates of PARL and its orthologues in different species, highlighting areas of uncertainty, and discuss its potential involvement in some prevalent diseases such as type II diabetes and Parkinson's disease.

Common variants of the PINK1 and PARL genes do not confer genetic susceptibility to schizophrenia in Han Chinese

Schizophrenia is a prevalent psychiatric disorder with a complex etiology. Mitochondrial dysfunction has been frequently reported in schizophrenia. Phosphatase and tension homologue-induced kinase 1 (PINK1) and presenilin-associated rhomboid-like protease (PARL) are mitochondrial proteins, and genetic variants of these two genes may confer genetic susceptibility to schizophrenia by influencing mitochondrial function. In this study, we conducted a two-stage genetic association study to test this hypothesis. We genotyped 4 PINK1 and 5 PARL genetic variants and evaluated the potential association of the 9 SNPs with schizophrenia in two independent case-control cohorts of 2510 Han Chinese individuals. No positive association of common genetic variants of the PINK1 and PARL genes with schizophrenia was identified in our samples after Bonferroni correction. Re-analysis of the newly updated Psychiatric Genetics Consortium (PGC) data sets confirmed our negative result. Intriguingly, one PINK1 SNP (rs10916832), which showed a marginally significant association in only Hunan samples (P = 0.032), is associated with the expression of a schizophrenia susceptible gene KIF17 according to the expression quantitative trait locus (eQTL) analysis. Our study indicated that common genetic variants of the PINK1 and PARL genes are unlikely to be involved in schizophrenia. Further studies are essential to characterize the role of the PINK1 and PARL genes in schizophrenia.

A genome-wide association study in the genetic analysis of idiopathic thrombophilia project suggests sex-specific regulation of mitochondrial DNA levels

Identifying genes that regulate mitochondrial DNA (mtDNA) levels is of interest due to an increasing number of diseases in humans that are associated with altered mtDNA levels. We searched for nuclear polymorphisms that influence mtDNA levels using a family-based genome-wide association (GWAS) method. Also, our aim was to determine if sex influences the genetic control of mtDNA levels. Two intron-polymorphisms, in the genes PARK2 and MRPL37, showed a tendency toward an association with mtDNA levels only in females and only in males, respectively. Both genes have a role in mitochondrial biogenesis and are potential candidates for the sex-specific control of mtDNA levels.

ABCC5 transporter is a novel type 2 diabetes susceptibility gene in European and African American populations

Numerous functional studies have implicated PARL in relation to type 2 diabetes (T2D). We hypothesised that conflicting human association studies may be due to neighbouring causal variants being in linkage disequilibrium (LD) with PARL. We conducted a comprehensive candidate gene study of the extended LD genomic region that includes PARL and transporter ABCC5 using three data sets (two European and one African American), in relation to healthy glycaemic variation, visceral fat accumulation and T2D disease. We observed no evidence for previously reported T2D association with Val262Leu or PARL using array and fine-map genomic and expression data. By contrast, we observed strong evidence of T2D association with ABCC5 (intron 26) for European and African American samples (P = 3E-07) and with ABCC5 adipose expression in Europeans [odds ratio (OR) = 3.8, P = 2E-04]. The genomic location estimate for the ABCC5 functional variant, associated with all phenotypes and expression data (P = 1E-11), was identical for all samples (at Chr3q 185,136 kb B36), indicating that the risk variant is an expression quantitative trait locus (eQTL) with increased expression conferring risk of disease. That the association with T2D is observed in populations of disparate ancestry suggests the variant is a ubiquitous risk factor for T2D.

Efficient mitochondrial biogenesis drives incomplete penetrance in Leber's hereditary optic neuropathy

Leber's hereditary optic neuropathy is a maternally inherited blinding disease caused as a result of homoplasmic point mutations in complex I subunit genes of mitochondrial DNA. It is characterized by incomplete penetrance, as only some mutation carriers become affected. Thus, the mitochondrial DNA mutation is necessary but not sufficient to cause optic neuropathy. Environmental triggers and genetic modifying factors have been considered to explain its variable penetrance. We measured the mitochondrial DNA copy number and mitochondrial mass indicators in blood cells from affected and carrier individuals, screening three large pedigrees and 39 independently collected smaller families with Leber's hereditary optic neuropathy, as well as muscle biopsies and cells isolated by laser capturing from post-mortem specimens of retina and optic nerves, the latter being the disease targets. We show that unaffected mutation carriers have a significantly higher mitochondrial DNA copy number and mitochondrial mass compared with their affected relatives and control individuals. Comparative studies of fibroblasts from affected, carriers and controls, under different paradigms of metabolic demand, show that carriers display the highest capacity for activating mitochondrial biogenesis. Therefore we postulate that the increased mitochondrial biogenesis in carriers may overcome some of the pathogenic effect of mitochondrial DNA mutations. Screening of a few selected genetic variants in candidate genes involved in mitochondrial biogenesis failed to reveal any significant association. Our study provides a valuable mechanism to explain variability of penetrance in Leber's hereditary optic neuropathy and clues for high throughput genetic screening to identify the nuclear modifying gene(s), opening an avenue to develop predictive genetic tests on disease risk and therapeutic strategies.

The role of PARL and HtrA2 in striatal neuronal injury after transient global cerebral ischemia

The presenilin-associated rhomboid-like (PARL) protein and high temperature requirement factor A2 (HtrA2) are key regulators of mitochondrial integrity and play pivotal roles in apoptosis. However, their roles after cerebral ischemia have not been thoroughly elucidated. To clarify these roles, mice were subjected to transient global cerebral ischemia, and striatal neuronal injury was assessed. Western blot and coimmunoprecipitation analyses revealed that PARL and processed HtrA2 localized to mitochondria, and that PARL was bound to HtrA2 in sham animals. Expression of PARL and processed HtrA2 in mitochondria significantly decreased 6 to 72 hours after ischemia, and the binding of PARL to HtrA2 disappeared after ischemia. In contrast, expression of processed HtrA2 increased 24 hours after ischemia in the cytosol, where HtrA2 was bound to X chromosome-linked inhibitor-of-apoptosis protein (XIAP). Administration of PARL small interfering RNA inhibited HtrA2 processing and worsened ischemic neuronal injury. Our results show that downregulation of PARL after ischemia is a key step in ischemic neuronal injury, and that it decreases HtrA2 processing and increases neuronal vulnerability. In addition, processed HtrA2 released into the cytosol after ischemia contributes to neuronal injury via inhibition of XIAP.

The emerging role of proteolysis in mitochondrial quality control and the etiology of Parkinson's disease

Mitochondria are highly dynamic organelles that are important for many diverse cellular processes, such as energy metabolism, calcium buffering, and apoptosis. Mitochondrial biology and dysfunction have recently been linked to different types of cancers and neurodegenerative diseases, most notably Parkinson's disease. Thus, a better understanding of the quality control systems that maintain a healthy mitochondrial network can facilitate the development of effective treatments for these diseases. In this perspective, we will discuss recent advances on two mitochondrial quality control pathways: the UPS and mitophagy, highlight how new players may be contributing to regulate these pathways. We believe the proteases involved will be key and novel regulators of mitochondrial quality control, and this knowledge will provide insights into future studies aimed to combat neurodegenerative diseases.

Mitochondrial disorders and the eye

The clinical significance of disturbed mitochondrial function in the eye has emerged since mitochondrial DNA (mtDNA) mutation was described in Leber's hereditary optic neuropathy. The spectrum of mitochondrial dysfunction has become apparent through increased understanding of the contribution of nuclear and somatic mtDNA mutations to mitochondrial dynamics and function. Common ophthalmic manifestations of mitochondrial dysfunction include optic atrophy, pigmentary retinopathy, and ophthalmoplegia. The majority of patients with ocular manifestations of mitochondrial disease also have variable central and peripheral nervous system involvement. Mitochondrial dysfunction has recently been associated with age-related retinal disease including macular degeneration and glaucoma. Therefore, therapeutic targets directed at promoting mitochondrial biogenesis and function offer a potential to both preserve retinal function and attenuate neurodegenerative processes.

Mitochondrial DNA copy number is modulated by genetic variation in the signal transducer and activator of transcription 3 (STAT3)

The regulation of mitochondrial DNA (mtDNA) copy number not only is critical for the maintenance of the normal mitochondrial function but has a strong clinical significance. A recent report revealed that the signal transducer and activator of transcription 3 (STAT3) is involved in the regulation of the mitochondrial function and is required for the optimal function of the electron transport chain. In this study, we explored whether gene variants in the STAT3 influence the leukocyte mtDNA copy number. Clinical data and blood samples were collected from 179 subjects (aged 52.8 ± 0.9 years). Mitochondrial DNA quantification using nuclear DNA (nDNA) as a reference was carried out by a real-time quantitative polymerase chain reaction method; results are presented as the mtDNA/nDNA ratio. We selected 3 tag single nucleotide polymorphisms showing a minor allele frequency greater than 10% (rs2293152 C/G, rs6503695 C/T, and rs9891119 A/C), representing 24 polymorphic sites of the STAT3 (r(2) > 0.8). We observed a significant association between mtDNA/nDNA ratio and both rs6503695 and rs9891119, adjusted by age and homeostasis model assessment index. The proportion of the total variance of the mtDNA/nDNA ratio accounted for by the rs6503695 and rs9891119 genotypes was 4.7% and 6.53%, respectively. Common variation in the STAT3 may influence mtDNA copy number.

Genome-wide linkage scan and association study of PARL to the expression of LHON families in Thailand

Leber hereditary optic neuropathy (LHON) is the most common mitochondrially inherited disease causing blindness, preferentially in young adult males. Most of the patients carry the G11778A mitochondrial DNA (mtDNA) mutation. However, the marked incomplete penetrance and the gender bias indicate some additional genetic and/or environmental factors to disease expression. Herein, we first conducted a genome-wide linkage scan with 400 microsatellite markers in 9 large Thai LHON G11778A pedigrees. Using an affecteds-only nonparametric linkage analysis, 4 regions on chromosomes 3, 12, 13 and 18 showed Zlr scores greater than 2 (P < 0.025), which is consistently significant across several linkage statistics. The most suggestive marker D3S1565 (Zlr > 2 in 10 of 16 allele sharing models tested) was then expanded to include the region 3q26.2-3q28 covering SLC7A14 (3q26.2), MFN1 (3q26.32), MRPL47 (3q26.33), MCCC1 (3q27.1), PARL (3q27.1) and OPA1 (3q28-q29). All of these candidate genes were selected from the Maestro database and had known to be localized in mitochondria. Sixty tag SNPs were genotyped in 86 cases, 211 of their relatives and 32 unrelated Thai controls, by multiplex-PCR-based Invader assay. Analyses using a powerful association testing tool that adjusts for relatedness (the M(QLS) statistic) showed the most evidence of association between two SNPs, rs3749446 and rs1402000 (located in PARL presenilins-associated rhomboid-like) and LHON expression (both P = 8.8 x 10(-5)). The mitochondrial PARL protease has been recently known to play a role with a dynamin-related OPA1 protein in preventing apoptotic events by slowing down the release of cytochrome c out of mitochondrial cristae junctions. Moreover, PARL is required to activate the intramembranous proteolyses resulting in the degradation of an accumulated pro-apoptotic protein in the outer mitochondrial membrane. Under these circumstances, variants of PARL are suggested to influence cell death by apoptosis which has long been believed to intrigue the neurodegeneration of LHON.