Skip to main content

ORIGINAL RESEARCH article

Front. Aging Neurosci., 04 July 2022
Sec. Parkinson’s Disease and Aging-related Movement Disorders
This article is part of the Research Topic Genetic and Molecular Diversity in Parkinson's Disease View all 9 articles

Genetic Analysis of Six Transmembrane Protein Family Genes in Parkinson’s Disease in a Large Chinese Cohort

\r\nYuwen Zhao,&#x;Yuwen Zhao1,2†Kailin Zhang,&#x;Kailin Zhang1,2†Hongxu PanHongxu Pan1Yige WangYige Wang1Xiaoxia ZhouXiaoxia Zhou1Yaqin XiangYaqin Xiang1Qian XuQian Xu1Qiying SunQiying Sun3Jieqiong TanJieqiong Tan4Xinxiang YanXinxiang Yan1Jinchen Li,,Jinchen Li2,3,4Jifeng Guo,,,Jifeng Guo1,2,4,5Beisha Tang,,,Beisha Tang1,2,4,5Zhenhua Liu,,,*Zhenhua Liu1,2,4,5*
  • 1Department of Neurology, Xiangya Hospital, Central South University, Changsha, China
  • 2National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, China
  • 3Department of Geriatrics, Xiangya Hospital, Central South University, Changsha, China
  • 4Centre for Medical Genetics and Hunan Key Laboratory of Medical Genetics, School of Life Sciences, Central South University, Changsha, China
  • 5Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, China

Objectives: Parkinson’s disease (PD) is a neurodegenerative disorder with the manifestation of motor symptoms and non-motor symptoms. Previous studies have indicated the role of several transmembrane (TMEM) protein family genes in PD pathogenesis.

Materials and Methods: In order to better investigate the genetic role of PD-related TMEM protein family genes in PD, including TMEM230, TMEM59, TMEM108, TMEM163, TMEM175, and TMEM229B, 1,917 sporadic early onset PD (sEOPD) or familial PD (FPD) patients and 1,652 healthy controls were analyzed by whole-exome sequencing (WES) while 1,962 sporadic late-onset PD (sLOPD) and 1,279 healthy controls were analyzed by whole-genome sequencing (WGS). Rare and common variants for each gene were included in the analysis.

Results: One hundred rare damaging or loss of function variants of six genes were found at the threshold of MAF < 0.1%. Three rare Dmis variants of TMEM230 were specifically identified in PD. Rare missense variants of TMEM59 were statistically significantly associated with PD in the WES cohort, indicating the role of TMEM59 in FPD and sEOPD. Rare missense variants of TMEM108 were suggestively associated with PD in the WGS cohort, indicating the potential role of TMEM108 in sLOPD. The rare variant of the other three genes and common variants of six genes were not significantly associated with PD.

Conclusion: We performed a large case-control study to systematically investigate the role of several PD-related TMEM protein family genes in PD. We identified three PD-specific variants in TMEM230, the significant association of TMEM59 with FPD, and sEOPD and the suggestive association of TMEM108 with sLOPD.

Introduction

Parkinson’s disease (PD) is the second most prevalent neurodegenerative disorder characterized by different combinations of motor and non-motor symptoms, including bradykinesia, resting tremor, rigidity, postural instability, olfactory loss, sleep dysfunction, autonomic dysfunction, and cognitive impairment (Armstrong and Okun, 2020). The pathophysiology of PD is dopaminergic neuron death and α-synuclein accumulation (Bloem et al., 2021). Neuroinflammation, mitochondrial/oxidative stress pathways dysregulation, autophagy/lysosomal dysfunction, and synaptic vesicle endocytosis disruptions have been indicated in the pathogenesis of PD (Nguyen et al., 2019; Hirsch and Standaert, 2021). Genetic, epigenetic, environmental, and lifestyle risk factors contribute to PD etiology (Jankovic and Tan, 2020). Previous studies have identified more than 20 PD-causing genes through linkage analysis or homozygous mapping strategies and many PD risk loci affecting disease risk, age at onset (AAO), and progression of PD through the genome-wide association study (GWAS) (Blauwendraat et al., 2020).

Transmembrane (TMEM) proteins belong to a large protein family with the characteristic of having at least one transmembrane peptide that entirely or partially crosses the biological membrane (Marx et al., 2020). With the transmembrane property, TMEM proteins play an essential role in many cell processes, such as signal transduction, ion/small molecule transportation, and cell adhesion (Ryu et al., 2019). Especially some members of the TMEM protein family play roles in PD-related pathways, such as synaptic vesicle trafficking and lysosomal and mitochondrial function (Deng et al., 2016; Jinn et al., 2017; Wang et al., 2021a). However, many TMEM proteins are still poorly investigated (Marx et al., 2020).

Several TMEM proteins have been reported to be associated with PD, for example, TMEM230, TMEM175, TMEM163, TMEM229B, TMEM108, and TMEM59. TMEM230 mutations were reported to cause autosomal dominant PD (Deng et al., 2016). Moreover, mutant TMEM230 could impair synaptic vesicle trafficking, disrupt mitochondria transport, and induce apoptotic cell death (Deng et al., 2016; Wang et al., 2020, 2021a). A genome-wide significant locus (rs34311866) located around the TMEM175 gene was reported in a large-scale GWAS meta-analysis study (Nalls et al., 2014), and TMEM175 deficiency resulted in lysosomal and mitochondrial dysfunction and α-synuclein aggregation (Jinn et al., 2017). Another two loci located around TMEM163 or TMEM229B were also identified to be associated with PD in a large-scale meta-analysis of the GWAS study (Nalls et al., 2014). In contrast, a large longitudinal study found the rs138073281 located around TMEM108 linked to cognitive progression in PD (Zhang et al., 2021). Specifically, the association of TMEM163 with PD was further confirmed in other cohorts (Kia et al., 2021; Rudakou et al., 2021). TMEM59 was reported within the PARK10 locus indicating a possible role in PD (Beecham et al., 2015). Moreover, over-expression of TMEM59 in a Drosophila PD model ameliorated the phenotype of shortened lifespan, impaired locomotor activity, and dopaminergic neuron loss (Zhang et al., 2018). Although the above evidence indicated that those TMEM proteins family genes might be associated with PD, few studies systematically analyze whether these genes are associated with early onset or late-onset PD populations, respectively, considering that early onset and late-onset cases may have certain genetic heterogeneity.

This study analyzed rare and common variants of TMEM230, TMEM175, TMEM163, TMEM229B, TMEM108, and TMEM59 with gene-based and allele-based, respectively, through next-generation sequencing in a large Chinese PD cohort to better explore the genetic contribution of these six TMEM protein family genes in Chinese PD patients.

Materials and Methods

Participants

Three thousand eight hundred seventy-nine unrelated PD patients and 2,931 healthy controls were enrolled in this study. Blood samples for DNA extraction and clinical data of each participant were collected. All the participants have been divided into two cohorts: 1,917 sporadic early onset PD (sEOPD) or familial PD (FPD) probands and 1,652 healthy controls were in one cohort, and 1,962 sLOPD probands and 1,279 healthy controls were in another cohort (Supplementary Table 1). PD diagnosis was made according to the Movement Disorders Society (MDS) clinical diagnostic criteria by neurologists with expertise in movement disorders at the Department of Neurology, Xiangya Hospital, and other centers from Parkinson’s Disease & Movement Disorders Multicenter Database and Collaborative Network in China (PD-MDCNC).1 This study conforms with the World Medical Association Declaration of Helsinki. Each participant signed the written consent. The Ethics Committee of Xiangya Hospital, Central South University, approved the study. It is worth noting that the PD patients harboring the pathogenic or likely pathogenic variants identified in sEOPD or FPD patients were ruled out in this study, according to the principle of our previous study (Zhao et al., 2020).

Gene Selection

We included six TMEM protein family genes to perform the genetic analysis in this study (Supplementary Table 2). TMEM230 was reported to cause autosomal dominant PD; TMEM175, TMEM163, and TMEM229B were found to be risk loci of PD in a large meta-analysis GWAS study; TMEM108 was associated with cognitive progression in PD; TMEM59 was functionally related to PD phenotype. To investigate the role of these six TMEM protein family genes in PD, we conducted the genetic analysis by next-generation sequencing in FPD, sEOPD, sLOPD patients, and healthy controls from a large Chinese cohort.

Sequencing and Bioinformatic Analysis

We sequenced the first cohort with 1,917 sEOPD or FPD patients and 1,652 healthy controls by whole-exome sequencing (WES) with an average sequencing depth of 123-fold using the HiSeq X10 platform; the second cohort with 1,962 sLOPD and 1,279 healthy controls were sequenced by whole-genome sequencing (WGS) with an average sequencing depth of 12-fold using the HiSeq Nova platform. Variants calling for those six TMEM protein family genes were done by the following methods or algorithms: Burrows-Wheeler Aligner-MEM algorithm (Li, 2014), Picard,2 and Genome Analysis Toolkit (Van der Auwera et al., 2013). Annotation for the variants was accomplished by ANNOVAR and VarCards. The pathogenicity of missense variants was predicted by the Combined Annotation Dependent Depletion (CADD) algorithm, similar to a previous publication. Variants from TMEM230, TMEM175, TMEM163, TMEM229B, TMEM108, and TMEM59 were then filtered by the same strategy used in our previous study (Pan et al., 2020). After the variants filtering process, rare missense variants, rare damaging missense (Dmis) variants, and loss-of-function (LoF) variants were used for the gene-based burden analysis. The single-variant logistic association analysis included all common variants in the coding regions and upstream/downstream 2000 bp of each gene. We conducted Sanger sequencing on several random rare damaging variants, and the results showed that these variants are actual (Supplementary Figure 1). Moreover, we performed Sanger sequencing for the available family member of corresponding families with rare damaging variants predicted by two algorithms (CADD score over 12.37 and ReVe score over 0.7). ReVe is another algorithm used for predicting the deleteriousness of the missense variants (Li et al., 2018).

Statistical Analysis

We gathered the rare missense variants, rare Dmis variants, or LoF variants, respectively, and performed gene-based burden analysis for each TMEM gene enrolled. We used the sequence kernel association test (SKAT) R package to run the optimized sequence kernel association test (SKAT-O) and calculate the association of different categories of variants from each gene with PD. We included sex and the first five principal components of ancestry as covariates to account for the potential confounding effects in the analysis for the WES cohort, and we also included an extra covariate, age, in the analysis for the WGS cohort. The gene-based burden analysis was done for each gene in each cohort at the threshold of minor allele frequency (MAF) < 1% or MAF < 0.1%. Bonferroni correction was used to adjust the multiple testing, and P < 0.0083 was regarded as statistically significant as six independent tests were performed in this study. The result with P-value between 0.0083 and 0.05 was considered “suggestively” significant in analyzing of the rare variants.

For the analysis of the common variants, PLINK was used to perform the logistic regression analysis in the WES cohort and WGS cohort, respectively, with the similar covariates mentioned above. Similarly, P < 0.0014 was considered statistically significant in the WES cohort (36 independent tests), whereas P < 0.000035 in the WGS cohort (1,444 independent tests).

Results

This study mined the WES or WGS data of 3,879 PD patients and 2,931 healthy controls. There were 1,917 sEOPD or FPD patients and 1,652 healthy controls in the first cohort with WES. The sEOPD or FPD patients were at the mean age of 52.22 ± 9.03 years with a mean AAO of 46.31 ± 8.40 years, and the healthy controls were at the mean age of 62.03 ± 12.59 at the time of enrollment. There were 1,962 sLOPD and 1,279 healthy controls for WGS in the second cohort. The mean age for the sLOPD patients and their healthy matching controls was 66.76 ± 7.08 years and 62.32 ± 7.11 years, respectively. In addition, the mean AAO for sLOPD at the enrollment was 61.88 ± 6.93 years (Supplementary Table 1).

A total of 118 rare missense and eight rare LoF variants were found in the WES cohort, whereas 111 rare missense and seven rare LoF variants were detected in the WGS cohort, at the threshold of MAF < 1%. Furthermore, 102 and 94 rare missense variants, and seven and six rare LoF variants were found in the WES and WGS data, respectively, at the threshold of MAF < 0.1% (Supplementary Table 3). On the other hand, for the WES cohort, there were 36 common variants of these six genes included in the single-variant logistic association analysis. In contrast, 1,444 common variants of these six genes included in the WGS cohort. However, they were not significantly associated with PD in our cohorts.

We detected nine missense variants, including six Dmis variants for TMEM230 with MAF < 1% in the WES, whereas 10 missense variants and six Dmis variants were identified in the WGS cohort. At the threshold of MAF < 0.1%, six or seven missense variants were found in the WES or WGS cohort, with four of them being Dmis variants in each cohort (Table 1; Figure 1A). Gene-based burden analysis showed no significant association of these variants with PD. However, three rare Dmis variants, including p.S174F, p.Y165C, and p.V9L were identified in four sporadic PD patients from the WGS cohort but not in control individuals of our two cohorts (Table 2).

TABLE 1
www.frontiersin.org

Table 1. Analysis of rare variant burden within targeted genes in Parkinson’s disease.

FIGURE 1
www.frontiersin.org

Figure 1. Schematic diagram of six PD-related transmembrane (TMEM) protein family genes. Schematic diagram showing the transmembrane regions of six PD-related TMEM protein family genes and the rare damaging missense and loss-of-function variants of TMEM230 (A), TMEM59 (B), TMEM108 (C), TMEM163 (D), TMEM175 (E), and TMEM229B (F) in both cohorts.

TABLE 2
www.frontiersin.org

Table 2. Parkinson’s disease (PD) patient specific variants of TMEM230.

We identified 14 rare missense variants, including 12 Dmis variants and one LoF variant for TMEM59 in the WES cohort with MAF < 1%. Gene-based burden analysis showed that the rare missense variants of TMEM59 in the WES cohort were suggestively associated with PD (P = 0.035). When the MAF was limited to less than 0.1%, the association of 13 rare missense variants of TMEM59 with PD became statistically significant (P = 0.007). At the same time, the 11 rare Dmis variants (P = 0.016) and the 12 Dmis plus LoF variants (P = 0.030) of TMEM59 both became suggestively associated with PD in the WES cohort (Table 1; Figure 1B). Furthermore, 24 PD cases carried rare Dmis or LoF variants of TMEM59, so we listed their corresponding clinical manifestations in Supplementary Table 4. All those patients showed good responses to L-dopa treatment, and 10 of 24 cases showed dyskinesia. Besides, 13 and 12 rare non-synonymous variants were identified in the WGS cohort at the threshold of MAF < 1% and 0.1%, respectively. However, the association of these variants with PD was not statistically significant (Table 1; Figure 1B).

Twenty-seven rare missense variants were identified in both WES and WGS cohorts at MAF < 1%, and 24 rare missense variants were identified in both cohorts at MAF < 0.1% for TMEM108. 17 Dmis variants were found in the WES cohort with MAF < 1% or MAF < 0.1%, while 12 Dmis and one LoF variant were found in the WGS cohort with two MAF thresholds (Table 1; Figure 1C). Besides, gene-based burden analysis showed that the rare missense variants of TMEM108 from the WGS cohort with MAF < 0.1% were suggestively associated with PD (P = 0.014).

TMEM175 had the most missense variants at two different thresholds of MAF. Fifty-one in the WES cohort and 49 in the WGS cohort with MAF < 1%, and 42 in the WES cohort and 40 in the WGS cohort with MAF < 0.1%. However, the association of these variants with PD was not significantly different (Table 1; Figure 1E). Around or less than 10 rare missense variants were identified in both WES and WGS cohorts with two MAF thresholds for TMEM163 and TMEM229B (Table 1; Figures 1D,F). Nevertheless, gene-based burden analysis did not find the association of these genes with PD, either.

Moreover, our study identified 20 PD families’ probands carrying rare predicted Dmis (MAF < 0.01 and CADD score over 12.37) or LoF variants. There are seven families with damaging variants predicted by two algorithms (CADD score over 12.37 and ReVe score over 0.7) discovered in FPD, and their corresponding pedigrees are shown in Supplementary Figure 2. Then, we conducted Sanger sequencing for the available family member of the corresponding family (AR-035); however, this variant did not co-segregate within the family.

Discussion

To systematically investigate the genetic burden of six TMEM protein family genes, including TMEM230, TMEM175, TMEM163, TMEM229B, TMEM108, and TMEM59, in PD, we conducted a large case-control study in PD from the Chinese population, including FPD patients, sporadic early onset, and late-onset PD patients. Moreover, we mined the WES and WGS data, including rare and common variants covering coding, intron, and 2000 bp upstream and downstream of each gene, to explore these six genes. Then we found PD patient-specific rare variants for TMEM230, the significant association of rare missense variants of TMEM59 with PD, and a suggestive association of rare missense variants of TMEM108 with PD. However, no common variant of these genes was found to be significantly associated with PD.

We identified several PD-specific rare non-synonymous variants for TMEM230, but no significant association was found between this gene with PD. The p.Y165C variant was previously identified in a Taiwanese sporadic PD patient, while p.S174F and p.V9L variants have not previously been reported in PD patients or healthy controls (Quadri et al., 2017). According to the American College of Medical Genetics and Genomics (ACMG) guidelines, these variants were predicted to be uncertain significance due to a lack of supporting evidence. Hence, verifications of these variants in more studies and related functional investigations are warranted to confirm their pathogenicity in the future. TMEM230 mutations were reported to cause autosomal dominant PD with the potential mechanism of impairing synaptic vesicle trafficking, disrupting mitochondria transport, and inducing apoptotic cell death (Deng et al., 2016; Wang et al., 2020, 2021a). Of note, TMEM230 protein is present in SNCA-positive Lewy bodies and Lewy neurites, indicating its potentially vital role in PD pathogenesis (Deng et al., 2016). There have been several subsequent verifications of the pathogenicity of the TMEM230 gene in PD. However, many of them have not been verified (Yan et al., 2017; Wang et al., 2021b), which may be due to the relatively small contribution of the TMEM230 gene in PD patients. In our recent large-cohort study, we identified one pathogenic variant (c.429delT/p.P144Qfs*2) of TMEM230, supporting the pathogenic role of TMEM230 in PD (Zhao et al., 2020).

Then, we explored the role of common and rare variants in the other five TMEM family genes located within known PD risk loci. This study is the first to report that rare variants of TMEM59 were significantly enriched in FPD and sEOPD patients, indicating TMEM59 might play an essential role in PD. TMEM59 has differentially expressed in substantia nigra pars compacta astrocytes in a PD mouse model (Khan et al., 2022). Besides, TMEM59 interacted with TREM2, which was also related to PD, and regulated microglia function (Liu et al., 2020). TMEM59 was found to mediate PD-related pathways, such as autophagy and dopamine system dysfunction (Boada-Romero et al., 2013; Liu et al., 2017); moreover, overexpression of TMEM59 in a Drosophila PD model ameliorated the phenotype of shortened lifespan, impaired locomotor activity, and loss of dopaminergic neurons by degrading a-synuclein (Zhang et al., 2018). Combining all the evidence, we supported that the rare variants of TMEM59 might predictably increase the risk of PD. Since this is the first study to report the association of TMEM59 with PD, genetic analysis studies for TMEM59 in PD from larger cohorts and different populations are warranted to verify our findings.

We also found the suggestive association of TMEM108 with PD. TMEM108 was reported to be enriched in the dentate gyrus and involved in regulating spine development and glutamatergic activity in the granule neurons of the dentate gyrus (Jiao et al., 2017). Moreover, TMEM108 was found to play a vital role in adult neurogenesis in the dentate gyrus (Yu et al., 2019). A large longitudinal study showed that TMEM108 variant rs138073281 was associated with cognitive progression in PD (Zhang et al., 2021). However, no association of the TMEM108 gene with PD was found in previous studies, and no functional study has been conducted to investigate the role of TMEM108 in PD pathogenesis. Further research for TMEM108 in PD should focus on these genetic analysis studies and functional investigations.

Our study identified several rare non-synonymous TMEM175, TMEM163, and TMEM229B variants, but no association of these genes with PD was found. TMEM175 is a component of a lysosomal K+ channel, and its deficiency results in impaired lysosomal and mitochondrial function (Jinn et al., 2017; Pergel et al., 2021). Variants of TMEM175 have been identified to be associated with PD in a GWAS study, indicating its role in PD (Nalls et al., 2014). Further studies found that both rare and common variants of TMEM175 were correlated to PD (Hopfner et al., 2020; Rudakou et al., 2021), and different common variants of TMEM175 may have opposite influences on stress-induced damage and α-synuclein accumulation in neurons, which complicated the specific role of TMEM175 in PD (Wie et al., 2021). It may explain why we could not find the association of TMEM175 with PD at the gene level. TMEM163 is a zinc efflux transporter and belongs to the cation diffusion facilitator (CDF) protein family (Sanchez et al., 2019; Styrpejko and Cuajungco, 2021). The mutations of CDF family member, SLC30A10, were reported to cause parkinsonism (Quadri et al., 2012). TMEM163 is highly expressed in the cortex and cerebellum, and its expression was positively associated with PD risk (Cuajungco et al., 2014; Kia et al., 2021). The loci around TMEM163 was reported to associate with PD in GWAS studies (Nalls et al., 2014; Bandrés-Ciga et al., 2016), but our study did not verify the association between TMEM163 with PD. Hence further studies with larger sample sizes and different ethnic populations are needed. TMEM229B is a novel protein with unclear biological function. We did not find the association of TMEM229B with PD, and the association of TMEM229B locus with PD was controversial in previous studies (Nalls et al., 2014; Chang et al., 2017). So, TMEM229B may not play a vital role in PD, and functional studies for TMEM229B and its mutants are warranted to illustrate its role in PD better.

In addition, several studies investigated the gene expression profile in the human postmortem substantia nigra of PD patients, but the expression level of the six genes investigated in our study was not significantly different between PD and controls (Kelly et al., 2019; Nido et al., 2020). We also searched the six genes in BrainEXP-NPD, a website showing expression profiling in human brains for six neuropsychiatric disorders, including PD.3 No expression level for TMEM230 was found in any brain region on the website. No significant difference between PD and controls for other five genes regarding the expression level in different brain regions was found (Supplementary Figure 3; Supplementary Table 5).

Our study has several limitations. Firstly, the statistical power of the analysis in our study was limited, even though we had a large case-control cohort. Secondly, non-coding variants cannot be detected by WES, so these variants were not analyzed in our sEOPD and FPD patients. Finally, our findings for the association of TMEM family genes with PD were only from genetic analysis and have not been further verified by functional studies.

In conclusion, we performed a large case-control study to systematically investigate the role of PD-related TMEM protein family genes in PD using WES and WGS with gene-based and allele-based burden analysis. We found three PD patient-specific variants of TMEM230, which need further verification for pathogenicity. Moreover, we identified the significant association of TMEM59 with FPD and sEOPD, and a suggestively significant association of TMEM108 with sLOPD. However, we did not find any significant association of common variants of TMEM family genes with PD. Our study shed light on the gene burden of TMEM protein family genes in PD, and furthermore extensive genetic analysis studies and functional experiments for these genes are warranted.

Data Availability Statement

According to national legislation/guidelines, specifically the Administrative Regulations of the China on Human Genetic Resources (http://www.gov.cn/zhengce/content/2019-06/10/content_5398829.htm, http://english.www.gov.cn/policies/latest_releases/2019/06/10/content_281476708945462.htm), no additional raw data is available at this time. Data of this project can be accessed after an approval application to the China National GeneBank (CNGB) (https://db.cngb.org/cnsa/). Please refer to https://db.cngb.org/, or email: CNGBdb@cngb.org for detailed application guidance. The accession code CNP0003134 should be included in the application.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of Xiangya Hospital of Central South University in China. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

YZ: design, methodology, and writing – original draft. KZ: design, execution, and writing – original draft. HP: methodology and writing. YW: blood samples collection and DNA extraction. XZ and QX: patients enrollment. YX: blood samples collection. QS: clinical data collection. JT: DNA extraction and writing – review. XY and JG: patients enrollment and clinical data collection. JL: methodology and writing – review. BT: conception and patients enrollment. ZL: design, conception, writing – review and editing, and supervision. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Key Research and Development Program of China (Grant Nos. 2016YFC306000, 2021YFC2501204, and 2021YFC2502100), the National Natural Science Foundation of China (Grant Nos. 82001359 and U20A20355), Hunan Innovative Province Construction Project (Grant Nos. 2019SK2335 and 2021SK1010), and the Scientific Research Project of Hunan Provincial Health Commission (Grant No. 202203074637).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We are grateful to all the individuals that participant in this study. We are also grateful for resources from the Bioinformatics Center, Xiangya Hospital, Central South University and the High Performance Computing Center of Central South University.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.889057/full#supplementary-material

Footnotes

  1. ^ http://pd-mdcnc.com/
  2. ^ http://broadinstitute.github.io/picard/
  3. ^ http://brainexpnpd.org:8088/BrainEXPNPD/index.html

References

Armstrong, M. J., and Okun, M. S. (2020). Diagnosis and treatment of Parkinson disease: a review. JAMA 323, 548–560. doi: 10.1001/jama.2019.22360

PubMed Abstract | CrossRef Full Text | Google Scholar

Bandrés-Ciga, S., Price, T. R., Barrero, F. J., Escamilla-Sevilla, F., Pelegrina, J., Arepalli, S., et al. (2016). Genome-wide assessment of Parkinson’s disease in a Southern Spanish population. Neurobiol. Aging 45, .e213.e3–213.e9. doi: 10.1016/j.neurobiolaging.2016.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Beecham, G. W., Dickson, D. W., Scott, W. K., Martin, E. R., Schellenberg, G., Nuytemans, K., et al. (2015). PARK10 is a major locus for sporadic neuropathologically confirmed Parkinson disease. Neurology 84, 972–980. doi: 10.1212/wnl.0000000000001332

PubMed Abstract | CrossRef Full Text | Google Scholar

Blauwendraat, C., Nalls, M. A., and Singleton, A. B. (2020). The genetic architecture of Parkinson’s disease. Lancet Neurol. 19, 170–178. doi: 10.1016/s1474-4422(19)30287-x

CrossRef Full Text | Google Scholar

Bloem, B. R., Okun, M. S., and Klein, C. (2021). Parkinson’s disease. Lancet 397, 2284–2303. doi: 10.1016/s0140-6736(21)00218-x

CrossRef Full Text | Google Scholar

Boada-Romero, E., Letek, M., Fleischer, A., Pallauf, K., Ramón-Barros, C., and Pimentel-Muiños, F. X. (2013). TMEM59 defines a novel ATG16L1-binding motif that promotes local activation of LC3. EMBO J. 32, 566–582. doi: 10.1038/emboj.2013.8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, D., Nalls, M. A., Hallgrímsdóttir, I. B., Hunkapiller, J., van der Brug, M., Cai, F., et al. (2017). A meta-analysis of genome-wide association studies identifies 17 new Parkinson’s disease risk loci. Nat. Genet. 49, 1511–1516. doi: 10.1038/ng.3955

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuajungco, M. P., Basilio, L. C., Silva, J., Hart, T., Tringali, J., Chen, C. C., et al. (2014). Cellular zinc levels are modulated by TRPML1-TMEM163 interaction. Traffic 15, 1247–1265. doi: 10.1111/tra.12205

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, H. X., Shi, Y., Yang, Y., Ahmeti, K. B., Miller, N., Huang, C., et al. (2016). Identification of TMEM230 mutations in familial Parkinson’s disease. Nat. Genet. 48, 733–739. doi: 10.1038/ng.3589

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirsch, E. C., and Standaert, D. G. (2021). Ten unsolved questions about neuroinflammation in Parkinson’s disease. Mov. Disord. 36, 16–24. doi: 10.1002/mds.28075

PubMed Abstract | CrossRef Full Text | Google Scholar

Hopfner, F., Mueller, S. H., Szymczak, S., Junge, O., Tittmann, L., May, S., et al. (2020). Rare variants in specific lysosomal genes are associated with Parkinson’s disease. Mov. Disord. 35, 1245–1248. doi: 10.1002/mds.28037

PubMed Abstract | CrossRef Full Text | Google Scholar

Jankovic, J., and Tan, E. K. (2020). Parkinson’s disease: etiopathogenesis and treatment. J. Neurol. Neurosurg. Psychiatry 91, 795–808. doi: 10.1136/jnnp-2019-322338

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiao, H. F., Sun, X. D., Bates, R., Xiong, L., Zhang, L., Liu, F., et al. (2017). Transmembrane protein 108 is required for glutamatergic transmission in dentate gyrus. Proc. Natl. Acad. Sci. U. S. A. 114, 1177–1182. doi: 10.1073/pnas.1618213114

PubMed Abstract | CrossRef Full Text | Google Scholar

Jinn, S., Drolet, R. E., Cramer, P. E., Wong, A. H., Toolan, D. M., Gretzula, C. A., et al. (2017). TMEM175 deficiency impairs lysosomal and mitochondrial function and increases α-synuclein aggregation. Proc. Natl. Acad. Sci. U. S. A. 114, 2389–2394. doi: 10.1073/pnas.1616332114

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelly, J., Moyeed, R., Carroll, C., Albani, D., and Li, X. (2019). Gene expression meta-analysis of Parkinson’s disease and its relationship with Alzheimer’s disease. Mol. Brain 12:16. doi: 10.1186/s13041-019-0436-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, A. H., Lee, L. K., and Smith, D. J. (2022). Single cell analysis of gene expression in the substantia nigra pars compacta of a pesticide-induced mouse model of Parkinson’s disease. bioRxiv [Preprint]. doi: 10.1101/2022.02.18.481079

CrossRef Full Text | Google Scholar

Kia, D. A., Zhang, D., Guelfi, S., Manzoni, C., Hubbard, L., Reynolds, R. H., et al. (2021). Identification of candidate Parkinson disease genes by integrating genome-wide association study, expression, and epigenetic data sets. JAMA Neurol. 78, 464–472. doi: 10.1001/jamaneurol.2020.5257

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H. (2014). Toward better understanding of artifacts in variant calling from high-coverage samples. Bioinformatics 30, 2843–2851. doi: 10.1093/bioinformatics/btu356

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Shi, L., Zhang, K., Zhang, Y., Hu, S., Zhao, T., et al. (2018). VarCards: an integrated genetic and clinical database for coding variants in the human genome. Nucleic Acids Res. 46, D1039–D1048. doi: 10.1093/nar/gkx1039

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Shi, J., Lin, R., and Wen, T. (2017). Dopamine and dopamine receptor D1 associated with decreased social interaction. Behav Brain Res. 324, 51–57. doi: 10.1016/j.bbr.2017.01.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Ning, J., Zheng, X., Meng, J., Han, L., Zheng, H., et al. (2020). TMEM59 interacts with TREM2 and modulates TREM2-dependent microglial activities. Cell Death Dis. 11:678. doi: 10.1038/s41419-020-02874-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Marx, S., Dal Maso, T., Chen, J. W., Bury, M., Wouters, J., Michiels, C., et al. (2020). Transmembrane (TMEM) protein family members: poorly characterized even if essential for the metastatic process. Semin Cancer Biol. 60, 96–106. doi: 10.1016/j.semcancer.2019.08.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Nalls, M. A., Pankratz, N., Lill, C. M., Do, C. B., Hernandez, D. G., Saad, M., et al. (2014). Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson’s disease. Nat. Genet. 46, 989–993. doi: 10.1038/ng.3043

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, M., Wong, Y. C., Ysselstein, D., Severino, A., and Krainc, D. (2019). Synaptic, mitochondrial, and lysosomal dysfunction in Parkinson’s disease. Trends Neurosci. 42, 140–149. doi: 10.1016/j.tins.2018.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nido, G. S., Dick, F., Toker, L., Petersen, K., Alves, G., Tysnes, O. B., et al. (2020). Common gene expression signatures in Parkinson’s disease are driven by changes in cell composition. Acta Neuropathol. Commun. 8:55. doi: 10.1186/s40478-020-00932-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, H. X., Zhao, Y. W., Mei, J. P., Fang, Z. H., Wang, Y., Zhou, X., et al. (2020). GCH1 variants contribute to the risk and earlier age-at-onset of Parkinson’s disease: a two-cohort case-control study. Transl. Neurodegener. 9:31. doi: 10.1186/s40035-020-00212-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Pergel, E., Veres, I., Csigi, G. I., and Czirják, G. (2021). Translocation of TMEM175 lysosomal potassium channel to the plasma membrane by dynasore compounds. Int. J. Mol. Sci. 22:10515. doi: 10.3390/ijms221910515

PubMed Abstract | CrossRef Full Text | Google Scholar

Quadri, M., Breedveld, G. J., Chang, H. C., Yeh, T. H., Guedes, L. C., Toni, V., et al. (2017). Mutations in TMEM230 are not a common cause of Parkinson’s disease. Mov. Disord. 32, 302–304. doi: 10.1002/mds.26900

PubMed Abstract | CrossRef Full Text | Google Scholar

Quadri, M., Federico, A., Zhao, T., Breedveld, G. J., Battisti, C., Delnooz, C., et al. (2012). Mutations in SLC30A10 cause parkinsonism and dystonia with hypermanganesemia, polycythemia, and chronic liver disease. Am. J. Hum. Genet. 90, 467–477. doi: 10.1016/j.ajhg.2012.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Rudakou, U., Yu, E., Krohn, L., Ruskey, J. A., Asayesh, F., Dauvilliers, Y., et al. (2021). Targeted sequencing of Parkinson’s disease loci genes highlights SYT11, FGF20 and other associations. Brain 144, 462–472. doi: 10.1093/brain/awaa401

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryu, H., Fuwad, A., Yoon, S., Jang, H., Lee, J. C., Kim, S. M., et al. (2019). Biomimetic membranes with transmembrane proteins: state-of-the-art in transmembrane protein applications. Int. J. Mol. Sci. 20:1437. doi: 10.3390/ijms20061437

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanchez, V. B., Ali, S., Escobar, A., and Cuajungco, M. P. (2019). Transmembrane 163 (TMEM163) protein effluxes zinc. Arch. Biochem. Biophys. 677:108166. doi: 10.1016/j.abb.2019.108166

PubMed Abstract | CrossRef Full Text | Google Scholar

Styrpejko, D. J., and Cuajungco, M. P. (2021). Transmembrane 163 (TMEM163) protein: a new member of the zinc efflux transporter family. Biomedicines 9:220. doi: 10.3390/biomedicines9020220

PubMed Abstract | CrossRef Full Text | Google Scholar

Van der Auwera, G. A., Carneiro, M. O., Hartl, C., Poplin, R., Del Angel, G., Levy-Moonshine, A., et al. (2013). From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Curr. Protoc. Bioinformatics 43, 11.10.1–11.10.33. doi: 10.1002/0471250953.bi1110s43

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Guo, G., Zhang, J., Aebez, N., Liu, Z., Liu, C. F., et al. (2021a). Mutant-TMEM230-induced neurodegeneration and impaired axonal mitochondrial transport. Hum. Mol. Genet. 30, 1535–1542. doi: 10.1093/hmg/ddab128

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Whelan, E., Liu, Z., Liu, C. F., and Smith, W. W. (2021b). Controversy of TMEM230 associated with Parkinson’s disease. Neuroscience 453, 280–286. doi: 10.1016/j.neuroscience.2020.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Wu, T., Zhang, J., Guo, G., He, X., Pei, Z., et al. (2020). Transmembrane protein 230 mediates a poly(ADP-ribose) polymerase-1-linked apoptosis. Front. Aging Neurosci. 12:235. doi: 10.3389/fnagi.2020.00235

PubMed Abstract | CrossRef Full Text | Google Scholar

Wie, J., Liu, Z., Song, H., Tropea, T. F., Yang, L., Wang, H., et al. (2021). A growth-factor-activated lysosomal K(+) channel regulates Parkinson’s pathology. Nature 591, 431–437. doi: 10.1038/s41586-021-03185-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, W., Tang, B., Zhou, X., Lei, L., Li, K., Sun, Q., et al. (2017). TMEM230 mutation analysis in Parkinson’s disease in a Chinese population. Neurobiol. Aging 49, 219.e1–219.e3. doi: 10.1016/j.neurobiolaging.2016.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Z., Lin, D., Zhong, Y., Luo, B., Liu, S., Fei, E., et al. (2019). Transmembrane protein 108 involves in adult neurogenesis in the hippocampal dentate gyrus. Cell Biosci. 9:9. doi: 10.1186/s13578-019-0272-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Feng, R., Li, Y., Gan, L., Zhou, F., Meng, S., et al. (2018). Degradation of alpha-synuclein by dendritic cell factor 1 delays neurodegeneration and extends lifespan in Drosophila. Neurobiol. Aging 67, 67–74. doi: 10.1016/j.neurobiolaging.2018.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W., Wang, R., Fang, H., Ma, X., Li, D., Liu, T., et al. (2021). Influence of low tumor content on tumor mutational burden estimation by whole-exome sequencing and targeted panel sequencing. Clin. Transl. Med. 11:e415. doi: 10.1002/ctm2.415

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Qin, L., Pan, H., Liu, Z., Jiang, L., He, Y., et al. (2020). The role of genetics in Parkinson’s disease: a large cohort study in Chinese mainland population. Brain 143, 2220–2234. doi: 10.1093/brain/awaa167

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: TMEM protein family genes, genetic analysis, Parkinson’s disease, burden analysis, rare variants, common variants

Citation: Zhao Y, Zhang K, Pan H, Wang Y, Zhou X, Xiang Y, Xu Q, Sun Q, Tan J, Yan X, Li J, Guo J, Tang B and Liu Z (2022) Genetic Analysis of Six Transmembrane Protein Family Genes in Parkinson’s Disease in a Large Chinese Cohort. Front. Aging Neurosci. 14:889057. doi: 10.3389/fnagi.2022.889057

Received: 03 March 2022; Accepted: 15 June 2022;
Published: 04 July 2022.

Edited by:

Siti Aishah Sulaiman, National University of Malaysia, Malaysia

Reviewed by:

Kuo-Hsuan Chang, Linkou Chang Gung Memorial Hospital, Taiwan
Jolanta Dorszewska, Poznań University of Medical Sciences, Poland

Copyright © 2022 Zhao, Zhang, Pan, Wang, Zhou, Xiang, Xu, Sun, Tan, Yan, Li, Guo, Tang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhenhua Liu, liuzhenhua@csu.edu.cn

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.