Mitochondrial DNA copy number is associated with risk of head and neck squamous cell carcinoma in Chinese population

Mitochondrial DNA copy number and risk of papillary thyroid carcinoma

BACKGROUND

Mitochondrial DNA (mtDNA) copy number is associated with tumor activity and carcinogenesis. This study was undertaken to investigate mtDNA copy number in papillary thyroid cancer (PTC) tissues and to evaluate the risk of PTC development. The clinicopathological features of patients and mtDNA copy number were correlated. The value of mtDNA copy number was evaluated as a biomarker for PTC.

METHOD

DNA was extracted from 105 PTC tissues and 67 control thyroid tissues, and mtDNA copy number mtDNA oxidative damage were determined using qPCR techniques.

RESULTS

Overall, the relative mtDNA copy number was significantly higher in PTC patients (p < 0.001). The risk of developing PTC increased significantly across the tertiles of mtDNA copy number (p trend < 0.001). The higher the mtDNA copy number tertile, the greater the risk of developing PTC. Patients with follicular variants had an odds ratio of 2.09 (95% CI: 1.78-2.44) compared to those with classical variants (p < 0.001). The level of mtDNA oxidative damage in PTC was significantly elevated compared to controls (p < 0.001). The ROC analysis of mtDNA copy number indicated an area under the curve (AUC) of 77.7% (95% CI: 0.71 to 0.85, p < 0.001) for the ability of mtDNA copy number z-scores in differentiate between PTC and controls.

CONCLUSION

Our results indicated that the augmentation of mtDNA content plays a significant role during the initiation of thyroid cancer, and it might represent a potential biomarker for predicting the risk of PTC.

Different impact of vitamin D on mitochondrial activity and morphology in normal and malignant keratinocytes, the role of genomic pathway

Deregulation of mitochondria activity is one of the hallmarks of cancerogenesis and an important target for cancer therapy. Therefore, we compared the impact of an active form of vitamin D3 (1,25(OH)2D3) on mitochondrial morphology and bioenergetics in human squamous cell carcinoma (A431) and immortalized HaCaT keratinocytes. It was shown that mitochondria of cancerous A431 cells differ from that observed in HaCaT keratinocytes in terms of network, morphology, bioenergetics, glycolysis, and mitochondrial DNA copy number, while treatment of A431 with 1,25(OH)2D3 partially eliminates these differences. Furthermore, mitochondrial membrane potential, basal respiration, and mitochondrial reactive oxygen species production were decreased in A431 cells treated with 1,25(OH)2D3. Additionally, the expression and protein level of mitophagy marker PINK1 was significantly increased in A431 1,25(OH)2D3 treated cells, but not observed in treated HaCaT cells. Knockout of VDR (vitamin D receptor) or RXRA (binding partner retinoid X receptor) partially altered mitochondrial morphology and function as well as mitochondrial response to 1,25(OH)2D3. Transcriptomic analysis on A431 cells treated with 1,25(OH)2D3 revealed modulation of expression of several mitochondrial-related genes involved in mitochondrial depolarization, mitochondrial protein translation (i.e. LYRM9, MARS2), and fusion-fission (OPA1, FIS1, MFN1 and 2), however, none of the genes coded by mitochondrial DNA was affected. Interestingly, in silico analyses of nuclear-encoded mitochondrial genes revealed that they are rather activated by the secondary genomic response to 1,25(OH)2D3. Taken together, 1,25(OH)2D3 remodels mitochondrial architecture and bioenergetics through VDR-dependent and only partially RXRA-dependent activation of the genomic pathway, thus outlining a new perspective for anticancer properties of vitamin D3 in relation to mitochondria in squamous cell carcinoma.

Associations of mitochondrial DNA copy number with incident risks of gastrointestinal cancers: A prospective case-cohort study

Epidemiological investigations implied that mitochondrial DNA copy number (mtDNAcn) variations could trigger predisposition to multiple cancers, but evidence regarding gastrointestinal cancers (GICs) was still uncertain. We conducted a case-cohort study within the prospective Dongfeng-Tongji cohort, including incident cases of colorectal cancer (CRC, n = 278), gastric cancer (GC, n = 138), and esophageal cancer (EC, n = 72) as well as a random subcohort (n = 1173), who were followed up from baseline to the end of 2018. We determined baseline blood mtDNAcn and associations of mtDNAcn with the GICs risks were estimated by using weighted Cox proportional hazards models. Significant U-shaped associations were observed between mtDNAcn and GICs risks. Compared to subjects within the second quartile (Q2) mtDNAcn subgroup, those within the 1st (Q1), 3rd (Q3), and 4th (Q4) quartile subgroups showed increased risks of CRC (hazard ratio [HR] [95% confidence interval, CI] = 2.27 [1.47-3.52], 1.65 [1.04-2.62], and 2.81 [1.85-4.28], respectively) and total GICs (HR [95%CI] = 1.84 [1.30-2.60], 1.47 [1.03-2.10], and 2.51 [1.82-3.47], respectively], and those within Q4 subgroup presented elevated GC and EC risks (HR [95% CI] = 2.16 [1.31-3.54] and 2.38 [1.13-5.02], respectively). Similar associations of mtDNAcn with CRC and total GICs risks remained in stratified analyzes by age, gender, smoking, and drinking status. This prospective case-cohort study showed U-shaped associations between mtDNAcn and GICs risks, but further research works are needed to uncover underlying biological mechanisms.

Mitochondrial DNA copy number as a diagnostic marker and indicator of degree of severity in alopecia areata

Alopecia areata (AA) is a disorder with several etiologies. The evidence suggests that the absolute copy number of mitochondrial deoxyribonucleic acid (mtDNA), as well as proportion of mutated mtDNA copies, determines disease onset. This study aims to quantify the relative index of the mtDNA copy number in patients with AA and healthy controls and correlate the results with the existing clinical ‎information. This case-control study included 50 patients with AA and 50 age- and sex-coordinated healthy persons as controls. The severity of AA was weighed using the Severity of Alopecia Tool and Kavak's classification. The relative index of the mtDNA copy number was measured by real-time qPCR. Significant statistical difference was observed between cases and controls regarding mean mtDNA copy number, p < .001. There was significant positive correlation with SALT score (p =  0.001). A cutoff value of >1.619 N/µL could significantly diagnose AA cases (p < .001), and a cutoff value of > 1.36 N/µL could discriminate mild AA cases from those with moderate AA (p =  0.007). The relative index of mtDNA copy number is significantly elevated in AA cases and could be helpful in diagnosing and evaluating AA severity.

Nuclear and mitochondrial DNA alterations in pheochromocytomas and paragangliomas, and their potential treatment

Mitochondrial DNA (mtDNA) alterations have been reported in different types of cancers and are suggested to play important roles in cancer development and metastasis. However, there is little information about its involvement in pheochromocytomas and paragangliomas (PCCs/PGLs) formation. PCCs and PGLs are rare endocrine tumors of the chromaffin cells in the adrenal medulla and extra-adrenal paraganglia that can synthesize and secrete catecholamines. Over the last 3 decades, the genetic background of about 60% of PCCs/PGLs involving nuclear DNA alterations has been determined. Recently, a study showed that mitochondrial alterations can be found in around 17% of the remaining PCCs/PGLs. In this review, we summarize recent knowledge regarding both nuclear and mitochondrial alterations and their involvement in PCCs/PGLs. We also provide brief insights into the genetics and the molecular pathways associated with PCCs/PGLs and potential therapeutical targets.

T cell-mediated tumor killing patterns in head and neck squamous cell carcinoma identify novel molecular subtypes, with prognosis and therapeutic implications

As an important process in cancer immunotherapy, T cell-mediated tumor killing (TTK) enhances the immune response of patients. However, the role of TTK in Head and Neck Squamous Cell Carcinoma (HNSCC) patients still needs further exploration. Therefore, we comprehensively analyzed the gene expression information and clinical characteristics of 1063 HNSCC in five cohorts. Univariate regression, differential expression analysis, and gene mutation profiling were combined to identify the important genes regulating the sensitivity of tumor cells to T cell-mediated killing (GSTTK) in HNSCC. A total of 20 GSTTK were identified as important genes of HNSCC. Patients were divided into C1 and C2 subgroups (TTK patterns) and displayed significant prognostic differences. Patients with C2 subtype had dismal prognosis characteristic compared to C1 subtype in all validation cohorts. Patients with C1 subgroup exhibited robust immune profile and C1 subgroup patients were significantly enriched in metabolically relevant functions. Notably, the multi-omics analysis found that C1 subgroup have higher mutation burden and C2 subgroup patients had significantly higher copy number variation. Drug sensitivity analysis found that multiple first-line chemotherapeutic drugs were more sensitive in patients with subgroup C1. In conclusion, the establishment of GSTTK provides guidance and assistance to clinicians in the personalized management and treatment of HNSCC patients.

Roles of Mitochondria in Oral Squamous Cell Carcinoma Therapy: Friend or Foe?

Oral squamous cell carcinoma (OSCC) therapy is unsatisfactory, and the prevalence of the disease is increasing. The role of mitochondria in OSCC therapy has recently attracted increasing attention, however, many mechanisms remain unclear. Therefore, we elaborate upon relative studies in this review to achieve a better therapeutic effect of OSCC treatment in the future. Interestingly, we found that mitochondria not only contribute to OSCC therapy but also promote resistance, and targeting the mitochondria of OSCC via nanoparticles is a promising way to treat OSCC.

Insights regarding mitochondrial DNA copy number alterations in human cancer (Review)

Mitochondria are the critical organelles involved in various cellular functions. Mitochondrial biogenesis is activated by multiple cellular mechanisms which require a synchronous regulation between mitochondrial DNA (mtDNA) and nuclear DNA (nDNA). The mitochondrial DNA copy number (mtDNA-CN) is a proxy indicator for mitochondrial activity, and its alteration reflects mitochondrial biogenesis and function. Despite the precise mechanisms that modulate the amount and composition of mtDNA, which have not been fully elucidated, mtDNA-CN is known to influence numerous cellular pathways that are associated with cancer and as well as multiple other diseases. In addition, the utility of current technology in measuring mtDNA-CN contributes to its extensive assessment of diverse traits and tumorigenesis. The present review provides an overview of mtDNA-CN variations across human cancers and an extensive summary of the existing knowledge on the regulation and machinery of mtDNA-CN. The current information on the advanced methods used for mtDNA-CN assessment is also presented.

Genetic Alterations in Mitochondrial DNA Are Complementary to Nuclear DNA Mutations in Pheochromocytomas

Simple Summary

Mitochondrial DNA (mtDNA) alterations have been reported to play important roles in cancer development and metastasis. However, there is scarce information about pheochromocytomas and paragangliomas (PCCs/PGLs) formation. To determine the potential roles of mtDNA alterations in PCCs/PGLs, we analyzed a panel of 26 nuclear susceptibility genes and the entire mtDNA sequence of 77 human tumors, using NGS. We also performed an analysis of copy-number alterations, large mtDNA deletion, and gene/protein expression. Our results revealed that 53.2% of the tumors harbor a mutation in the susceptibility genes and 16.9% harbor complementary mitochondrial mutations. Large deletions and depletion of mtDNA were found in 26% and 87% of tumors, respectively, accompanied by a reduced expression of the mitochondrial biogenesis markers (PCG1α, NRF1, and TFAM). Furthermore, P62 and LC3a gene expression suggested increased mitophagy, which is linked to mitochondrial dysfunction. These finding suggest a complementarity and a potential contributing role in PCCs/PGLs tumorigenesis.

Abstract

Background: Somatic mutations, copy-number variations, and genome instability of mitochondrial DNA (mtDNA) have been reported in different types of cancers and are suggested to play important roles in cancer development and metastasis. However, there is scarce information about pheochromocytomas and paragangliomas (PCCs/PGLs) formation. Material: To determine the potential roles of mtDNA alterations in sporadic PCCs/PGLs, we analyzed a panel of 26 nuclear susceptibility genes and the entire mtDNA sequence of seventy-seven human tumors, using next-generation sequencing, and compared the results with normal adrenal medulla tissues. We also performed an analysis of copy-number alterations, large mtDNA deletion, and gene and protein expression. Results: Our results revealed that 53.2% of the tumors harbor a mutation in at least one of the targeted susceptibility genes, and 16.9% harbor complementary mitochondrial mutations. More than 50% of the mitochondrial mutations were novel and predicted pathogenic, affecting mitochondrial oxidative phosphorylation. Large deletions were found in 26% of tumors, and depletion of mtDNA occurred in more than 87% of PCCs/PGLs. The reduction of the mitochondrial number was accompanied by a reduced expression of the regulators that promote mitochondrial biogenesis (PCG1α, NRF1, and TFAM). Further, P62 and LC3a gene expression suggested increased mitophagy, which is linked to mitochondrial dysfunction. Conclusion: The pathogenic role of these finding remains to be shown, but we suggest a complementarity and a potential contributing role in PCCs/PGLs tumorigenesis.

Spectrum of mitochondrial genomic variation in parathyroid neoplasms

PURPOSE

Mitochondrial DNA (mtDNA) variations have been implicated in various cancer types. Several attempts have been made in benign parathyroid adenoma (PA); however, mtDNA variants and copy number alterations have not been investigated in parathyroid carcinoma (PC). The present study aimed to explore the variations in the mitochondrial genome in parathyroid neoplasms.

METHODS

In this study, ultradeep targeted sequencing of mtDNA was performed in 12 cases with PC and 25 cases with PA. Germline and somatic variants in the mitochondrial genome were analysed according to clinical features. The mtDNA copy number relative to nuclear DNA was measured.

RESULTS

A total of 821 germline variants and 23 somatic mutations were identified. All 160 germline nonsynonymous variants in the coding sequence (CDS) were predicted to be likely benign, and germline variants frequently occurred (26.3%) in the displacement loop region. Tumour somatic mutation in a CDS with heteroplasmic factor (HF) >0.8 showed a higher mtDNA copy number (p = 0.039). Using Spearman's rank test, tumour mtDNA copy number revealed a positive correlation with serum levels of intact parathyroid hormone and calcium.

CONCLUSION

Our results demonstrate that null recurrent mtDNA mutational hotspots were PC specific. An increased mtDNA copy number in parathyroid neoplasms was associated with somatic CDS mutations with a high HF and major clinical features. Larger cohort investigations should be performed in future analyses.

Identification of a Nuclear Mitochondrial-Related Multi-Genes Signature to Predict the Prognosis of Bladder Cancer

INTRODUCTION

Bladder cancer (BC) is one of the most prevalent urinary cancers, and its management is still a problem causing recurrence and progression, elevating mortality.

MATERIALS AND METHODS

We aimed at the nuclear mitochondria-related genes (MTRGs), collected from the MITOMAP: A Human Mitochondrial Genome Database. Meanwhile, the expression profiles and clinical information of BC were downloaded from the Cancer Genome Atlas (TCGA) as a training group. The univariate, multivariate, and the least absolute shrinkage and selection operator (LASSO) Cox regression analyses were used to construct a nuclear mitochondrial-related multi-genes signature and the prognostic nomogram.

RESULTS

A total of 17 nuclear MTRGs were identified to be correlated with the overall survival (OS) of BC patients, and a nuclear MTRGs signature based on 16 genes expression was further determined by the LASSO Cox regression analysis. Based on a nuclear MTRGs scoring system, BC patients from the TCGA cohort were divided into high- and low- nuclear MTRGs score groups. Patients with a high nuclear MTRGs score exhibited a significantly poorer outcome (median OS: 92.90 vs 20.20 months, p<0.0001). The nuclear MTRGs signature was further verified in three independent datasets, namely, GSE13507, GSE31684, and GSE32548, from the Gene Expression Omnibus (GEO). The BC patients with a high nuclear MTRGs score had significantly worse survival (median OS in GSE13507: 31.52 vs 98.00 months, p<0.05; GSE31684: 32.85 months vs unreached, p<0.05; GSE32548: unreached vs unreached, p<0.05). Furthermore, muscle-invasive bladder cancer (MIBC) patients had a significantly higher nuclear MTRGs score (p<0.05) than non-muscle-invasive bladder cancer (NMIBC) patients. The integrated signature outperformed each involved MTRG. In addition, a nuclear MTRGs-based nomogram was constructed as a novel prediction prognosis model, whose AUC values for OS at 1, 3, 5 years were 0.76, 0.75, and 0.75, respectively, showing the prognostic nomogram had good and stable predicting ability. Enrichment analyses of the hallmark gene set and KEGG pathway revealed that the E2F targets, G2M checkpoint pathways, and cell cycle had influences on the survival of BC patients. Furthermore, the analysis of tumor microenvironment indicated more CD8+ T cells and higher immune score in patients with high nuclear MTRGs score, which might confer sensitivity to immune checkpoint inhibitors.

CONCLUSIONS

Not only could the signature and prognostic nomogram predict the prognosis of BC, but it also had potential therapeutic guidance.

Mitochondrial Heteroplasmy Shifting as a Potential Biomarker of Cancer Progression

Cancer is a serious health problem with a high mortality rate worldwide. Given the relevance of mitochondria in numerous physiological and pathological mechanisms, such as adenosine triphosphate (ATP) synthesis, apoptosis, metabolism, cancer progression and drug resistance, mitochondrial genome (mtDNA) analysis has become of great interest in the study of human diseases, including cancer. To date, a high number of variants and mutations have been identified in different types of tumors, which coexist with normal alleles, a phenomenon named heteroplasmy. This mechanism is considered an intermediate state between the fixation or elimination of the acquired mutations. It is suggested that mutations, which confer adaptive advantages to tumor growth and invasion, are enriched in malignant cells. Notably, many recent studies have reported a heteroplasmy-shifting phenomenon as a potential shaper in tumor progression and treatment response, and we suggest that each cancer type also has a unique mitochondrial heteroplasmy-shifting profile. So far, a plethora of data evidencing correlations among heteroplasmy and cancer-related phenotypes are available, but still, not authentic demonstrations, and whether the heteroplasmy or the variation in mtDNA copy number (mtCNV) in cancer are cause or consequence remained unknown. Further studies are needed to support these findings and decipher their clinical implications and impact in the field of drug discovery aimed at treating human cancer.

Association Between Mitochondrial DNA Copy Number and Head and Neck Squamous Cell Carcinoma: A Systematic Review and Dose-Response Meta-Analysis

Background

The association between mitochondrial DNA (mtDNA) copy number and head and neck squamous cell carcinoma (HNSCC) risk remains unclear. Therefore, we aimed to evaluate the relationship between mtDNA copy number and HNSCC risk.

Material/Methods

We searched PubMed, Web of Science, and EMBASE until August 2020. Studies that assessed the association between mtDNA copy number and HNSCC as the outcome of interest were included. We performed a 2-class and dose-response meta-analysis to assess the association between cancer risk and mtDNA.

Results

Eight articles (2 cohort studies and 6 case-control studies) with a total of 3913 patients were included in our meta-analysis. The overall results showed that mean mtDNA copy number level from 9 studies was 0.71 higher in patients with cancer than in non-cancer controls (the standardized mean differences (SMD) 0.71, 95% CI: 0.28–1.15, P<0.001). However, when 4 studies were pooled by dichotomizing mtDNA copy number at the median value into high- and low-content groups, no significant association between mtDNA content and overall cancer risk was found (odds ratio (OR)=0.87, 95% CI: 0.52–1.44, P=0.584). Furthermore, we observed a non-linear association from 3 studies between increased mtDNA copy number levels (P for nonlinearity <0.001).

Conclusions

The elevated mtDNA copy number could predict the risk of HNSCC as a biomarker. Moreover, there was non-linear relationship of risk between HNSCC and mtDNA copy number.

Epigenetic alteration of mitochondrial biogenesis regulatory genes in arsenic exposed individuals (with and without skin lesions) and in skin cancer tissues: A case control study

Chronic arsenic toxicity has become a global concern due to its adverse pathophysiological outcome and carcinogenic potential. It is already established that arsenic induced reactive oxygen species alters mitochondrial functionality. Major regulatory genes for mitochondrial biogenesis, i.e., PGC1α, Tfam, NRF1and NRF2 are located in the nucleus. As a result, mitochondria-nucleus crosstalk is crucial for proper mitochondrial function. This previous hypothesis led us to investigateinvolvement of epigenetic alteration behindenhanced mitochondrial biogenesis in chronic arsenic exposure. An extensive case-control study was conducted with 390 study participants (unexposed, exposed without skin lesion, exposed with skin lesion and exposed skin tumour) from highly arsenic exposed areas ofWest Bengal, India. Methylation specific PCRrevealed significant promoter hypomethylation oftwo key biogenesis regulatory genes, PGC1αandTfam in arsenic exposed individuals and also in skin tumour tissues. Linear regression analysis indicated significant negative correlation between urinary arsenic concentration and promoter methylation status. Increased expression of biogenesis regulatory genes wasobtained by quantitative real-time PCR analysis. Moreover, altered mitochondrial fusion-fission regulatory gene expression was also observed in skin tumour tissues. miR663, having tumour suppressor gene like function was known to be epigenetically regulated through mitochondrial retrograde signal. Promoter hypermethylation with significantly decreased expression of miR663 was found in skin cancer tissues compared to non-cancerous control tissue. In conclusion, results indicated crucial role of epigenetic alteration in arsenic induced mitochondrial biogenesis and arsenical skin carcinogenesis for the first time. However, further mechanistic studies are necessary for detailed understanding of mitochondria-nucleus crosstalk in arsenic perturbation.