Next Article in Journal
Mouse Models of Musculocontractural Ehlers-Danlos Syndrome
Next Article in Special Issue
Molecular Modeling Analysis Provides Genotype–Phenotype Correlation Insights in a Patient with Ankyloblepharon-Ectodermal Dysplasia-Clefting Syndrome
Previous Article in Journal
Computational Prediction of Protein Intrinsically Disordered Region Related Interactions and Functions
Previous Article in Special Issue
Selected SNPs of FCN2 Associated with Chronic Tonsillitis in the Polish Adult Population
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mitochondria-Related TFAM and POLG Gene Variants and Associations with Tumor Characteristics and Patient Survival in Head and Neck Cancer

1
Department of Genetics and Molecular Medicine, Lithuanian University of Health Sciences, 44307 Kaunas, Lithuania
2
Institute of Oncology, Lithuanian University of Health Sciences, 44307 Kaunas, Lithuania
3
Department of Pathological Anatomy, Lithuanian University of Health Sciences, 44307 Kaunas, Lithuania
4
Department of Otorhinolaryngology, Lithuanian University of Health Sciences, 44307 Kaunas, Lithuania
*
Author to whom correspondence should be addressed.
Genes 2023, 14(2), 434; https://doi.org/10.3390/genes14020434
Submission received: 4 January 2023 / Revised: 2 February 2023 / Accepted: 7 February 2023 / Published: 8 February 2023
(This article belongs to the Special Issue Head and Neck Genetics)

Abstract

:
In 2020, 878,348 newly reported cases and 444,347 deaths related to head and neck cancer were reported. These numbers suggest that there is still a need for molecular biomarkers for the diagnosis and prognosis of the disease. In this study, we aimed to analyze mitochondria-related mitochondrial transcription factor A (TFAM) and DNA polymerase γ (POLG) single-nucleotide polymorphisms (SNPs) in the head and neck cancer patient group and evaluate associations between SNPs, disease characteristics, and patient outcomes. Genotyping was performed using TaqMan probes with Real-Time polymerase chain reaction. We found associations between TFAM gene SNPs rs11006129 and rs3900887 and patient survival status. We found that patients with the TFAM rs11006129 CC genotype and non-carriers of the T allele had longer survival times than those with the CT genotype or T-allele carriers. Additionally, patients with the TFAM rs3900887 A allele tended to have shorter survival times than non-carriers of the A allele. Our findings suggest that variants in the TFAM gene may play an important role in head and neck cancer patient survival and could be considered and further evaluated as prognostic biomarkers. However, due to the limited sample size (n = 115), further studies in larger and more diverse cohorts are needed to confirm these findings.

1. Introduction

Head and neck cancer is a cancer type that affects various locations in the upper aerodigestive epithelium, such as the lip and oral cavity, oropharynx, larynx, and pharynx. Almost all head and neck cancers form in the epithelium of airways and are squamous cell carcinoma [1]. Despite arising from one cell type, this cancer type is heterogeneous due to its occurrence in different anatomical locations. Based on Global Cancer Observatory (GLOBOCAN)-reported data, in 2020, head and neck cancers accounted for 878,348 newly reported cases and 444,347 deaths [2]. The annual incidence in the United States in 2022 was projected to be approximately 66,470 new patients and ~15,050 new deaths [3]. The most-reported risk factors include excessive smoking, alcohol use, and oncogenic virus, such as human papillomavirus (HPV), infection [4]. The risk was estimated to be ten times higher for cigarette smokers than non-smokers, and alcohol consumption was also proven to be a risk factor and to enhance the smoking effect [5,6]. In older patients, head and neck cancer is associated with smoking and alcohol use. The number of such patients has decreased due to the reduced usage of tobacco [7]. HPV-positive and HPV-negative oropharyngeal cancers are different in their molecular and clinical parameters. HPV status impacts head and neck cancer. HPV-positive cancer is usually associated with a better prognosis, as it responds to treatment (radiation or chemotherapy) better [1,8,9]. However, HPV-associated cancer is more common among younger patients reflecting HPV type 16-induced head and neck cancer after oral sex exposure [10,11]. This shows that there are still unknown molecular mechanisms underlying head and neck cancer. However, Epidermal Growth Factor Receptor (EGFR) has been studied extensively, and it is known that it is upregulated in most (up to 90%) head and neck squamous cell carcinoma patients and is associated with a poor prognosis [12]. Head and neck cancer treatment is generally complex and requires specialists from multiple fields for the best patient outcome [13]. Treatment includes different options and combinations, such as surgery, chemotherapy, radiotherapy, and EGFR inhibitors [8]. Despite all current knowledge and available treatment options, there is a great need for biomarkers that would help in the early diagnosis and prognosis of the course of the disease and the more accurate prediction of the outcome. Changes in genes at the molecular level could be used as such biomarkers for the more detailed subtyping of head and neck cancer.
Single-nucleotide polymorphisms (SNPs) are among the most studied biomarkers in various diseases and health conditions. However, nuclear genes coding for mitochondrial proteins are rarely considered critical elements in carcinogenesis. Mitochondria are essential organelles important for energy generation, apoptosis, and metabolism. Moreover, they may play a critical role in carcinogenesis [14]. Since the mitochondrial genome encodes only 13 proteins, 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs), the nuclear genome encodes other mitochondrial proteins (more than 1500) [15,16]. For example, mitochondrial transcription factor A (TFAM) is a nuclear-encoded protein essential in transcription, replication, and DNA damage sensing and repair [17,18,19,20,21,22]. In a previous study, lower expression of TFAM was found to be correlated with a decreased amount of mitochondrial DNA (mtDNA) [17]. Additionally, TFAM SNPs were reported to be associated with prostate, colorectal, breast, uterine, ovary, and cervical cancer and TFAM expression alterations in lung and breast cancer [23,24,25,26,27,28,29,30]. Another protein—DNA polymerase γ (POLG)—is the polymerase active in mitochondria and encoded by the nuclear POLG gene [31]. Mutations in this polymerase are known to cause mtDNA depletion and, therefore, decreased oxidative phosphorylation (OXPHOS) [32]. POLG mutations are usually studied in mitochondrial diseases, including ataxia, progressive external ophthalmoplegia, mitochondrial epilepsy, Alpers’ syndrome, Leigh’s syndrome, Parkinsonism, and male infertility [33,34,35,36,37,38]. They are not the usual targets in cancers, although several POLG SNPs were studied in breast, cervical, and colorectal cancer studies [23,32,39,40]. However, there is still a lack of knowledge about mitochondria-related nuclear SNPs such as TFAM and POLG in cancers. There are no previous studies of the selected TFAM and POLG polymorphisms in head and neck cancer patients. Therefore, we focused on SNPs in these genes and analyzed their associations with various disease and tumor pathomorphological characteristics, such as the stage, tumor grade, size, spread to regional lymph nodes, metastasis, and survival status in this study.

2. Materials and Methods

2.1. Study Subjects

Our study included one hundred fifteen patients diagnosed with head and neck cancer. Patients were included in this study if they were adults (older than 18 years old), were diagnosed with head and neck cancer without prior diseases, and agreed to participate in this study. After their initial diagnosis, patients were informed about the study in detail and signed the consent form. Their blood samples were collected afterward and stored until DNA extraction. Data about the condition, including the patient’s clinical data and tumor characteristics, including gender, age at diagnosis, stage, survival status, tumor localization, grade, size, regional lymph node status, and metastasis status were collected from medical records. To ensure confidentiality, all samples were depersonalized and assigned identification codes.

2.2. Genotyping

Genotyping was performed at the Institute of Oncology, Lithuanian University of Health Sciences (LUHS). DNA was purified from blood samples using the GeneJet Genomic DNA purification kit (Thermo Fisher Scientific, Waltham, MA, USA, Cat. K0721) according to the user’s manual. Selected single-nucleotide variants (rs11006132, rs11006129, rs1937, rs16912174, rs1692202, and rs3900887 in the TFAM gene and rs3087374, rs2072267, rs976072, and rs2307441 in the POLG gene) were determined with predesigned and commercially available TaqMan probe SNP Genotyping Assays (Thermo Fisher Scientific, Cat. #4351379) on a QuantStudio 3 Real-Time polymerase chain reaction (PCR) System (Thermo Fisher Scientific, Cat. #A28137). The standard genotyping PCR program was used: 95 °C for 10 min (polymerase activation), 95 °C for 15 s (denaturation), and 60 °C for 1 min (annealing and extension). The full protocol included 45 cycles of denaturation, annealing, and extension. Reactions were assembled into a 12 μL total volume and included 15 ng of DNA, 6.125 μL of TaqMan Universal Master Mix (Thermo Fisher Scientific, Waltham, MA, USA, Cat. #4304437), 0.625 μL of TaqMan SNP Genotyping Assay, and nuclease-free water. A no-template control (nuclease-free water) was used to confirm the lack of contamination in every run. Genotyping results were analyzed using an allelic discrimination plot, and genotypes were determined according to VIC and FAM fluorescence intensity after the run. The no-template controls were seen at the bottom left corner of the allelic discrimination plot for all runs, which indicated that there was no contamination.

2.3. Statistical Analysis

SPSS software version 28 (IBM Corp., Armonk, NY, USA) was used for statistical analysis. Associations between genotypes and disease or tumor parameters were evaluated with Pearson’s Chi-square or Fisher’s Exact tests and binary logistic regression with univariate and multivariate models (including various parameters). Survival analyses were performed using Kaplan-Meier analysis, and the difference between survival curves was analyzed using a log-rank test. Cox regression analysis was performed with univariate and multivariate models to check the impacts of additional factors (age at the time of diagnosis, tumor size T, status of regional lymph node involvement, and tumor grade) on survival. p < 0.05 was considered to indicate a statistically significant difference.

3. Results

3.1. Study Cohort Characteristics

This study included one hundred fifteen patients diagnosed with head and neck cancer. Most patients were males (n = 106), with only a fraction of females (n = 9). The median age in this cohort was 62 years (age range, 31–85 years). Larynx carcinoma was the most frequent head and neck cancer type among our patient group and had an 80.9% incidence rate, while the least frequent was oral cavity cancer, with only a 2.6% incidence rate (Table 1). Most patients were diagnosed with stage IV cancer (37.4%), had no metastasis (96.5%) or regional lymph node involvement (60.9%) reported, and were alive (61.7%) during the last censoring. Tumors were mostly grade G2 (75.7%), and their sizes varied from T1 to T4 in approximate distribution; the data are provided in Table 2.

3.2. TFAM and POLG Genotype and Allele Frequencies

We determined that TFAM rs11006132, rs11006129, rs1937, and rs16912174 and POLG rs3087374, rs2072267, rs976072, and rs2307441 were in Hardy-Weinberg equilibrium (HWE), and the distributions of all genotypes are presented in Table 3. Due to an amplification failure, four data points are missing from the TFAM rs3900887 genotyping data.

3.3. Association Analyses

Associations between polymorphisms and clinical or tumor characteristics data were evaluated using Pearson’s Chi-square test (Table 4 and Table 5). TFAM rs16912202 was excluded from further investigations because the TT genotype was the only genotype in our study group.
We found that TFAM rs3900887 and rs11006129 were significantly associated with the metastasis status (rs3900887 p = 0.007; rs11006129 p = 0.036) and patient survival status (rs3900887 p = 0.018; rs11006129 p = 0.019). Additionally, rs11006129 was associated with the tumor differentiation grade (p = 0.036), and rs3087374 in the POLG gene was associated with tumor size (p = 0.017). To further investigate significant associations, we grouped tumor size and differentiation grade data into subgroups; however, no statistically significant associations remained (Table 4 and Table 5). We removed metastasis status from further analysis due to the small count of metastasis (3.5%) events in our study group. No associations were found between TFAM rs11006132, rs1937, or rs16912174 or POLG rs2072267, rs976072, or rs2307441 polymorphisms and clinical data or tumor characteristics.
We continued investigating the significant associations between TFAM rs3900887 and rs11006129 and survival status; therefore, we performed logistic regression analyses. We found that patients with the TFAM rs3900887 AT genotype had an approximately 4.6 times higher death probability (odds ratio (OR) = 4.565; 95% CI, 1.554–13.414; p = 0.006) than those with the TT genotype.
Additionally, patients with the TFAM rs11006129 CT genotype had an almost 3 times higher death probability (odds ratio (OR) = 2.875; 95% CI, 1.166–7.086; p = 0.022) than those with the CC genotype. These associations also remained significant in multivariate genotype and allelic model analyses, where additional cofactors were included and are presented in Table 6.

3.4. Survival Analyses

We analyzed all genotypes and alleles in our study to determine associations with overall patient survival. Kaplan–Meier plots for all significant associations are represented in Figure 1. Our findings suggest that both the TFAM rs11006129 T allele (p = 0.003) and CT vs. CC (p = 0.003) genotype are statistically significant for patient survival (Figure 1). Non-carriers of the TFAM rs11006129 T allele and those with the CC genotype survived significantly longer than T-allele carriers. This was followed by Cox univariate and multivariate regression analyses that included additional cofactors. In the Cox univariate analysis adjusted for age at diagnosis, it was determined that TFAM rs11006129 T-allele carriers had an approximately 2.5 times higher probability of shorter survival than non-carriers (hazard ratio (HR) = 2.505; 95% CI 1.331–4.713; p = 0.004). The differences remained significant in Cox multivariate analysis, following adjustment for other cofactors (tumor size T, regional lymph node status N, tumor grade G) (Table 7).
Additionally, we discovered that the TFAM rs3900887 A allele (p = 0.004) and AT vs. TT (p = 0.001) genotype are important for patient survival (Figure 1). We found that A-allele carriers had significantly shorter survival than non-carriers. Therefore, we expected that homozygous AA genotype patients would have the lowest survival rates; however, this was not the case. We hypothesized that this may possibly be due to the small sample size of AA genotypes (n = 8) and low death rates within this genotype (n = 3). In Cox univariate analysis adjusted for age at diagnosis, it was determined that A-allele carriers had an approximately two and a half times higher probability of shorter survival than non-carriers (hazard ratio (HR) = 2.460; 95% CI 1.319–4.588; p = 0.005). The differences remained significant in the Cox multivariate analysis following adjustment for other cofactors (tumor size T, regional lymph node status N, tumor grade G) (Table 7).

4. Discussion

Head and neck cancer affects various locations, such as the oral cavity, oropharynx, larynx, and pharynx. Most head and neck cancers are squamous cell carcinoma and form in the epithelium [1]. The major risk factors are excessive tobacco and alcohol consumption; however, HPV was also recognized as one of the risk factors and was associated with a better overall prognosis [4,41]. Despite the current knowledge about this cancer type, there are still 878,348 newly reported cases and 444,347 deaths worldwide each year [2], and multidisciplinary teams are required to achieve the best outcome and ensure the quality of a patient’s life.
Single-nucleotide polymorphisms, also known as SNPs, are among the widely studied biomarkers; however, mitochondrial DNA or nuclear DNA coding for mitochondrial proteins is rarely considered in cancer research. However, mitochondria are involved in major cellular processes and are considered to be important in carcinogenesis [14]. Since the nuclear genome encodes most of the mitochondrial proteins, we chose to analyze two of them. The first gene we selected was TFAM, which encodes mitochondrial transcription factor A. It is involved in crucial cellular processes, such as transcription, translation, replication, damage sensing, and mtDNA repair [17,18,19,20,21,22]. Additionally, the levels of mtDNA transcripts depend on the copy number of mtDNA, which is important for maintaining adenosine triphosphate (ATP) production [19,42]. The second gene we chose for analysis was POLG, which encodes the α subunit of polymerase γ, the only functioning polymerase in the mitochondria [31]. It is responsible for mtDNA replication, and POLG mutations are a known cause of mtDNA depletion and decreasedOXPHOS [32]. Mutations in POLG are often found in mitochondrial diseases, such as ataxia, progressive external ophthalmoplegia, mitochondrial epilepsy, Alpers’ syndrome, Leigh’s syndrome, Parkinsonism, and male infertility [33,34,35,36,37,38].
In this study, we analyzed rs11006132, rs11006129, rs1937, rs16912174, rs1692202, and rs3900887 in the TFAM gene and rs3087374, rs2072267, rs976072, and rs2307441 in the POLG gene. We selected these SNPs to be in both coding and non-coding sequences of each gene, to be benign in the general population, and to have a minor allele frequency (MAF) lower than 0.50 in the European population. Additionally, to the best of our knowledge, they were not previously studied in head and neck cancer. We found that rs11006129 and rs3900887 in the TFAM gene were significantly associated with patient survival. The rs11006129 variant is in the intron of the TFAM gene. The total MAF of this polymorphism in the 1000 Genome Project data was 0.15. However, the pathogenicity of this variant remains unclear, as it is predicted to be benign in Varsome [43] but is not reported in ClinVar. The TFAM rs11006129 CC genotype and non-carriers of the T allele had a probability of longer survival compared with those with the CT genotype or T-allele carriers. In our previous breast cancer patient cohort study, we found that patients with the rs11006129 T allele were less likely to have estrogen-positive tumors, and carriers of the C allele were less likely to develop tumors with vascular invasion [23]. Therefore, this confirms that the T allele may be a pathogenic factor for breast and head and neck cancer patients. However, no associations between rs11006129 and clinical data were found in the cervical cancer group [27]. The rs3900887 variant is in the intron of the TFAM gene. The total MAF of this polymorphism in the 1000 Genome Project data was 0.15. It is predicted to be benign by Varsome [43], while ClinVar lacks information on this variant, and the clinical significance is still unknown. We found that TFAM rs3900887 A-allele carriers tended to have shorter survival times than non-carriers. However, in the genotype analysis, we found that heterozygous patients with the AT genotype had the shortest survival times compared to both homozygous variants (TT or AA). From the results of the allele model, we expected AA homozygous patients to have the shortest survival. However, we believe that the small sample size of AA genotypes, with only eight patients and only three death events within this group, may have impacted the association. We analyzed rs3900887 in our previous cervical and breast cancer patient cohorts and found some statistically significant associations with tumor characteristics. For example, in cervical cancer, patients with TFAM rs3900887 TT and TA genotypes had a lower risk of larger tumors than patients with the AA genotype [27]. In the breast cancer group, we found that patients with the rs3900887 TT or TA genotype were associated with an increased risk of positive regional lymph nodes. Carriers of the TT genotype had an increased risk of estrogen receptor positivity or lymphatic invasion compared with patients with the AA genotype [23]. This SNP was previously analyzed in Alzheimer’s disease; however, no associations were found [44].
No associations were observed between disease characteristics and TFAM rs1937, rs11006132, or rs16912174. These variants are a missense variant (MAF 0.10), a non-coding transcript variant (MAF 0.24), and an upstream variant (MAF 0.01), respectively. One SNP, rs1692202, which is a non-coding transcript variant, MAF 0.12, was removed from statistical analyses because it only presented with the TT genotype in our study group. We did not find any previous studies focusing on these SNPs in head and neck cancer patients and evaluating SNP associations with disease characteristics. No significant associations were found with POLG rs2072267, rs3087374, rs976072, or rs2307441. POLG rs2072267 is an intron variant, MAF 0.47; rs3087374 is a missense variant, MAF 0.08; rs976072 is a downstream variant, MAF 0.39; and rs2307441 is a missense variant, 0.04. We also analyzed POLG rs2072267 in our previous cervical cancer study and determined that patients with the AA genotype survived longer without metastasis than those with the GG genotype [40]. We also found it to be associated with progression in breast cancer patients [23]. However, POLG rs2072267 was also investigated in Parkinson’s disease, ataxia, and colorectal cancer, but no significant results were reported. Our previous cervical cancer study found that the rs3087374 CA genotype was associated with the stage and tumor size compared to the CC genotype. Still, no associations were found in breast cancer or Parkinson’s disease studies [23,40,45]. Additionally, analyses of rs976072 in our previous breast and cervical cancer studies did not show any associations [23,40]. No significant results were reported in a study on bladder cancer [46]. Despite that, it was important for pancreatic cancer [47]. Similarly, in our previous breast cancer patient study, we found that rs2307441 was associated with tumor vascular invasion and metastasis-free survival [23]. However, our previous cervical cancer group study determined no associations between rs2307441 and clinical parameters, and no associations were reported in Parkinson’s or colorectal cancer [39,40,45].
There were limiting factors in this study, including a limited number of patients, as there were only one hundred fifteen patients in our cohort, and only nine of them were females. Additionally, only ten SNPs were selected and analyzed in TFAM and POLG genes that were reported to be benign in the general population. Despite these limitations, our results suggest that mitochondrial transcription factor A, encoded by the nuclear TFAM gene, may play an important role in head and neck cancer patient survival and could potentially serve as a prognostic survival biomarker. Therefore, further studies in larger and more diverse cohorts are needed to confirm our findings. After confirming TFAM SNPs in larger cohorts, future directions could potentially include the construction of the mutation by CRISPR/Cas9 or ectopic expression analysis and the evaluation of their impact on protein expression and translational activity.

Author Contributions

Conceptualization, I.G. and R.U.; methodology, A.B.; software, I.G.; validation, A.V., E.P., V.R. and E.J.; formal analysis, I.G.; investigation, A.B.; resources, E.J., R.U., V.R., A.V., E.P. and L.P.; data curation, R.U.; writing—original draft preparation, I.G. and R.U.; writing—review and editing, E.J., V.R., A.B., A.V., E.P. and L.P.; visualization, I.G.; supervision, E.J. and R.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Kaunas Regional Biomedical Research Ethics Committee (Nos. BE-2-10 and P1-BE-2-10/2014).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We would like to acknowledge our gratitude to the patients included in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leemans, C.R.; Snijders, P.J.F.; Brakenhoff, R.H. The Molecular Landscape of Head and Neck Cancer. Nat. Rev. Cancer 2018, 18, 269–282. [Google Scholar] [CrossRef] [PubMed]
  2. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef]
  3. Siegel, R.L.; Miller, K.D.; Fuchs, H.E.; Jemal, A. Cancer Statistics, 2022. CA A Cancer J. Clin. 2022, 72, 7–33. [Google Scholar] [CrossRef] [PubMed]
  4. Shaw, R.; Beasley, N. Aetiology and Risk Factors for Head and Neck Cancer: United Kingdom National Multidisciplinary Guidelines. J. Laryngol. Otol. 2016, 130, S9–S12. [Google Scholar] [CrossRef] [PubMed]
  5. Sturgis, E.M.; Cinciripini, P.M. Trends in Head and Neck Cancer Incidence in Relation to Smoking Prevalence. Cancer 2007, 110, 1429–1435. [Google Scholar] [CrossRef]
  6. Hashibe, M.; Brennan, P.; Benhamou, S.; Castellsague, X.; Chen, C.; Curado, M.P.; Maso, L.D.; Daudt, A.W.; Fabianova, E.; Wünsch-Filho, V.; et al. Alcohol Drinking in Never Users of Tobacco, Cigarette Smoking in Never Drinkers, and the Risk of Head and Neck Cancer: Pooled Analysis in the International Head and Neck Cancer Epidemiology Consortium. JNCI J. Natl. Cancer Inst. 2007, 99, 777–789. [Google Scholar] [CrossRef]
  7. Mourad, M.; Jetmore, T.; Jategaonkar, A.A.; Moubayed, S.; Moshier, E.; Urken, M.L. Epidemiological Trends of Head and Neck Cancer in the United States: A SEER Population Study. J. Oral Maxillofac. Surg. 2017, 75, 2562–2572. [Google Scholar] [CrossRef]
  8. Argiris, A.; Karamouzis, M.V.; Raben, D.; Ferris, R.L. Head and Neck Cancer. Lancet 2008, 371, 1695–1709. [Google Scholar] [CrossRef]
  9. Licitra, L.; Perrone, F.; Bossi, P.; Suardi, S.; Mariani, L.; Artusi, R.; Oggionni, M.; Rossini, C.; Cantù, G.; Squadrelli, M.; et al. High-Risk Human Papillomavirus Affects Prognosis in Patients with Surgically Treated Oropharyngeal Squamous Cell Carcinoma. J. Clin. Oncol. 2006, 24, 5630–5636. [Google Scholar] [CrossRef]
  10. Fitzmaurice, C.; Allen, C.; Barber, R.M.; Barregard, L.; Bhutta, Z.A.; Brenner, H.; Dicker, D.J.; Chimed-Orchir, O.; Dandona, R.; Dandona, L.; et al. Global Burden of Disease Cancer Collaboration Global, Regional, and National Cancer Incidence, Mortality, Years of Life Lost, Years Lived With Disability, and Disability-Adjusted Life-Years for 32 Cancer Groups, 1990 to 2015: A Systematic Analysis for the Global Burden of Disease Study. JAMA Oncol. 2017, 3, 524–548. [Google Scholar] [CrossRef]
  11. Gillison, M.L.; Chaturvedi, A.K.; Anderson, W.F.; Fakhry, C. Epidemiology of Human Papillomavirus–Positive Head and Neck Squamous Cell Carcinoma. J. Clin. Oncol. 2015, 33, 3235–3242. [Google Scholar] [CrossRef] [PubMed]
  12. Solomon, B.; Young, R.J.; Rischin, D. Head and Neck Squamous Cell Carcinoma: Genomics and Emerging Biomarkers for Immunomodulatory Cancer Treatments. Semin. Cancer Biol. 2018, 52, 228–240. [Google Scholar] [CrossRef] [PubMed]
  13. Kusampudi, S.; Konduru, N. General Introduction to Head and Neck Cancer: Etiology, Symptoms, Diagnosis, Staging, Prevention, and Treatment. In Early Detection and Treatment of Head & Neck Cancers: Theoretical Background and Newly Emerging Research; El Assal, R., Gaudilliere, D., Connelly, S.T., Eds.; Springer International Publishing: Cham, Switzerland, 2021; pp. 1–50. ISBN 978-3-030-69852-2. [Google Scholar]
  14. Van Gisbergen, M.W.; Voets, A.M.; Starmans, M.H.W.; de Coo, I.F.M.; Yadak, R.; Hoffmann, R.F.; Boutros, P.C.; Smeets, H.J.M.; Dubois, L.; Lambin, P. How Do Changes in the MtDNA and Mitochondrial Dysfunction Influence Cancer and Cancer Therapy? Challenges, Opportunities and Models. Mutat. Res. Rev. Mutat. Res. 2015, 764, 16–30. [Google Scholar] [CrossRef] [PubMed]
  15. Garesse, R.; Vallejo, C.G. Animal Mitochondrial Biogenesis and Function: A Regulatory Cross-Talk between Two Genomes. Gene 2001, 263, 1–16. [Google Scholar] [CrossRef]
  16. Taanman, J.-W. The Mitochondrial Genome: Structure, Transcription, Translation and Replication. Biochim. Biophys. Acta-Bioenerg. 1999, 1410, 103–123. [Google Scholar] [CrossRef]
  17. Kang, D.; Kim, S.H.; Hamasaki, N. Mitochondrial Transcription Factor A (TFAM): Roles in Maintenance of MtDNA and Cellular Functions. Mitochondrion 2007, 7, 39–44. [Google Scholar] [CrossRef]
  18. Hallberg, B.M.; Larsson, N.-G. TFAM Forces MtDNA to Make a U-Turn. Nat. Struct. Mol. Biol. 2011, 18, 1179–1181. [Google Scholar] [CrossRef] [PubMed]
  19. Bonawitz, N.D.; Clayton, D.A.; Shadel, G.S. Initiation and Beyond: Multiple Functions of the Human Mitochondrial Transcription Machinery. Mol. Cell 2006, 24, 813–825. [Google Scholar] [CrossRef]
  20. Garstka, H.L.; Schmitt, W.E.; Schultz, J.; Sogl, B.; Silakowski, B.; Pérez-Martos, A.; Montoya, J.; Wiesner, R.J. Import of Mitochondrial Transcription Factor A (TFAM) into Rat Liver Mitochondria Stimulates Transcription of Mitochondrial DNA. Nucleic Acids Res. 2003, 31, 5039–5047. [Google Scholar] [CrossRef]
  21. The Mitochondrial Transcription Factor TFAM Coordinates the Assembly of Multiple DNA Molecules into Nucleoid-like Structures | Molecular Biology of the Cell. Available online: https://www.molbiolcell.org/doi/full/10.1091/mbc.e07-05-0404 (accessed on 4 October 2022).
  22. Alterations to the Expression Level of Mitochondrial Transcription Factor A, TFAM, Modify the Mode of Mitochondrial DNA Replication in Cultured Human Cells | Nucleic Acids Research | Oxford Academic. Available online: https://academic.oup.com/nar/article/34/20/5815/3100463 (accessed on 4 October 2022).
  23. Golubickaite, I.; Ugenskiene, R.; Korobeinikova, E.; Gudaitiene, J.; Vaitiekus, D.; Poskiene, L.; Juozaityte, E. The Impact of Mitochondria-Related POLG and TFAM Variants on Breast Cancer Pathomorphological Characteristics and Patient Outcomes. Biomarkers 2021, 26, 343–353. [Google Scholar] [CrossRef]
  24. Granados, J.B.; Méndez, J.P.; Feria-Bernal, G.; García-García, E.; Tejeda, M.E.; Rojano-Mejía, D.; Tapia, A.; Canto, P. Association of a TFAM Haplotype with Aggressive Prostate Cancer in Overweight or Obese Mexican Mestizo Men. Urol. Oncol. 2017, 35, e9–e111. [Google Scholar] [CrossRef] [PubMed]
  25. Guo, J.; Zheng, L.; Liu, W.; Wang, X.; Wang, Z.; Wang, Z.; French, A.J.; Kang, D.; Chen, L.; Thibodeau, S.N.; et al. Frequent Truncating Mutation of TFAM Induces Mitochondrial DNA Depletion and Apoptotic Resistance in Microsatellite-Unstable Colorectal Cancer. Cancer Res. 2011, 71, 2978–2987. [Google Scholar] [CrossRef] [PubMed]
  26. Hu, W.; Ma, S.-L.; Liu, L.-L.; Zhu, Y.-H.; Zeng, T.-T.; Li, Y.; Guan, X.-Y. Impact of Mitochondrial Transcription Factor A Expression on the Outcomes of Ovarian, Endometrial and Cervical Cancers. Am. J. Transl. Res. 2020, 12, 5343–5361. [Google Scholar] [PubMed]
  27. Golubickaite, I.; Ugenskiene, R.; Cepaite, J.; Ziliene, E.; Inciura, A.; Poskiene, L.; Juozaityte, E. Mitochondria-Related TFAM Gene Variants and Their Effects on Patients with Cervical Cancer. Biomed. Rep. 2021, 15, 1–7. [Google Scholar] [CrossRef]
  28. Xie, D.; Wu, X.; Lan, L.; Shangguan, F.; Lin, X.; Chen, F.; Xu, S.; Zhang, Y.; Chen, Z.; Huang, K.; et al. Downregulation of TFAM Inhibits the Tumorigenesis of Non-Small Cell Lung Cancer by Activating ROS-Mediated JNK/P38MAPK Signaling and Reducing Cellular Bioenergetics. Oncotarget 2016, 7, 11609–11624. [Google Scholar] [CrossRef]
  29. Gao, W.; Wu, M.-H.; Wang, N.; Ying, M.-Z.; Zhang, Y.-Y.; Hua, J.; Chuan, L.; Wang, Y.-J. Mitochondrial Transcription Factor A Contributes to Cisplatin Resistance in Patients with Estrogen Receptor-positive Breast Cancer. Mol. Med. Rep. 2016, 14, 5304–5310. [Google Scholar] [CrossRef]
  30. Peng, H.; Yang, M.; Chen, Z.; Chen, P.; Guan, C.; Xiang, X.; Cai, S.; Chen, Y.; Fang, X. Expression and Methylation of Mitochondrial Transcription Factor a in Chronic Obstructive Pulmonary Disease Patients with Lung Cancer. PLoS ONE 2013, 8, e82739. [Google Scholar] [CrossRef]
  31. Chan, S.S.L.; Copeland, W.C. DNA Polymerase γ and Mitochondrial Disease: Understanding the Consequence of POLG Mutations. Biochim. Biophys. Acta Bioenerg. 2009, 1787, 312–319. [Google Scholar] [CrossRef]
  32. Singh, K.K.; Ayyasamy, V.; Owens, K.M.; Koul, M.S.; Vujcic, M. Mutations in Mitochondrial DNA Polymerase γ Promote Breast Tumorigenesis. J. Hum. Genet. 2009, 54, 516–524. [Google Scholar] [CrossRef]
  33. Rahman, S.; Copeland, W.C. POLG-Related Disorders and Their Neurological Manifestations. Nat. Rev. Neurol. 2019, 15, 40–52. [Google Scholar] [CrossRef]
  34. Hudson, G.; Chinnery, P.F. Mitochondrial DNA Polymerase-γ and Human Disease. Hum. Mol. Genet. 2006, 15, R244–R252. [Google Scholar] [CrossRef]
  35. Davidzon, G.; Greene, P.; Mancuso, M.; Klos, K.J.; Ahlskog, J.E.; Hirano, M.; DiMauro, S. Early-Onset Familial Parkinsonism Due to POLG Mutations. Ann. Neurol. 2006, 59, 859–862. [Google Scholar] [CrossRef] [PubMed]
  36. Longley, M.J.; Graziewicz, M.A.; Bienstock, R.J.; Copeland, W.C. Consequences of Mutations in Human DNA Polymerase γ. Gene 2005, 354, 125–131. [Google Scholar] [CrossRef] [PubMed]
  37. Naviaux, R.K.; Nguyen, K.V. POLG Mutations Associated with Alpers’ Syndrome and Mitochondrial DNA Depletion. Ann. Neurol. 2004, 55, 706–712. [Google Scholar] [CrossRef] [PubMed]
  38. Van Goethem, G.; Dermaut, B.; Löfgren, A.; Martin, J.J.; Van Broeckhoven, C. Mutation of POLG Is Associated with Progressive External Ophthalmoplegia Characterized by MtDNA Deletions. Nat. Genet. 2001, 28, 211–212. [Google Scholar] [CrossRef] [PubMed]
  39. Linkowska, K.; Jawień, A.; Marszałek, A.; Malyarchuk, B.A.; Tońska, K.; Bartnik, E.; Skonieczna, K.; Grzybowski, T. Mitochondrial DNA Polymerase γ Mutations and Their Implications in MtDNA Alterations in Colorectal Cancer. Ann. Hum. Genet. 2015, 79, 320–328. [Google Scholar] [CrossRef]
  40. Golubickaite, I.; Ugenskiene, R.; Ziliene, E.; Beniusyte, J.; Inciura, A.; Poskiene, L.; Juozaityte, E. POLG Gene Variants in Cervical Cancer Patients and Their Associations with Clinical and Pathomorphological Tumor Characteristics. J. Clin. Med. 2021, 10, 1838. [Google Scholar] [CrossRef]
  41. Evidence for a Causal Association Between Human Papillomavirus and a Subset of Head and Neck Cancers | JNCI: Journal of the National Cancer Institute | Oxford Academic. Available online: https://academic.oup.com/jnci/article/92/9/709/2906131?login=true (accessed on 22 March 2022).
  42. Ekstrand, M.I.; Falkenberg, M.; Rantanen, A.; Park, C.B.; Gaspari, M.; Hultenby, K.; Rustin, P.; Gustafsson, C.M.; Larsson, N.-G. Mitochondrial Transcription Factor A Regulates MtDNA Copy Number in Mammals. Hum. Mol. Genet. 2004, 13, 935–944. [Google Scholar] [CrossRef]
  43. Kopanos, C.; Tsiolkas, V.; Kouris, A.; Chapple, C.E.; Albarca Aguilera, M.; Meyer, R.; Massouras, A. VarSome: The Human Genomic Variant Search Engine. Bioinformatics 2019, 35, 1978–1980. [Google Scholar] [CrossRef]
  44. Grupe, A.; Li, Y.; Rowland, C.; Nowotny, P.; Hinrichs, A.L.; Smemo, S.; Kauwe, J.S.K.; Maxwell, T.J.; Cherny, S.; Doil, L.; et al. A Scan of Chromosome 10 Identifies a Novel Locus Showing Strong Association with Late-Onset Alzheimer Disease. Am. J. Hum. Genet. 2006, 78, 78–88. [Google Scholar] [CrossRef] [Green Version]
  45. Gui, Y.; Xu, Z.; Lv, W.; Liu, H.; Zhao, J.-J.; Hu, X.-Y. Association of Mitochondrial DNA Polymerase γ Gene POLG1 Polymorphisms with Parkinsonism in Chinese Populations. PLoS ONE 2012, 7, e50086. [Google Scholar] [CrossRef] [PubMed]
  46. Wei, H.; Kamat, A.; Chen, M.; Ke, H.-L.; Chang, D.W.; Yin, J.; Grossman, H.B.; Dinney, C.P.; Wu, X. Association of Polymorphisms in Oxidative Stress Genes with Clinical Outcomes for Bladder Cancer Treated with Bacillus Calmette-Guérin. PLoS ONE 2012, 7, e38533. [Google Scholar] [CrossRef] [PubMed]
  47. Li, D.; Duell, E.J.; Yu, K.; Risch, H.A.; Olson, S.H.; Kooperberg, C.; Wolpin, B.M.; Jiao, L.; Dong, X.; Wheeler, B.; et al. Pathway Analysis of Genome-Wide Association Study Data Highlights Pancreatic Development Genes as Susceptibility Factors for Pancreatic Cancer. Carcinogenesis 2012, 33, 1384–1390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Overall survival of patients with different TFAM and POLG genotypes. Kaplan–Meier survival curves for overall survival of (a) TFAM rs11006129 T-allele carriers (p = 0.003); (b) patients with TFAM rs11006129 CC and CT genotypes (p = 0.003); (c) TFAM rs3900887 A-allele carriers (p = 0.011); (d) patients with TFAM rs3900887 AT and TT genotypes (p = 0.004) and TT and AT genotypes (p = 0.001). Censored cases are shown as a vertical line. The shadow with the same color represents the 95% confidence interval.
Figure 1. Overall survival of patients with different TFAM and POLG genotypes. Kaplan–Meier survival curves for overall survival of (a) TFAM rs11006129 T-allele carriers (p = 0.003); (b) patients with TFAM rs11006129 CC and CT genotypes (p = 0.003); (c) TFAM rs3900887 A-allele carriers (p = 0.011); (d) patients with TFAM rs3900887 AT and TT genotypes (p = 0.004) and TT and AT genotypes (p = 0.001). Censored cases are shown as a vertical line. The shadow with the same color represents the 95% confidence interval.
Genes 14 00434 g001
Table 1. Frequency of head and neck cancer localization in our cohort.
Table 1. Frequency of head and neck cancer localization in our cohort.
LocalizationFrequency
Larynx93 (80.9%)
Pharynx14 (12.2%)
Paranasal and nasal5 (4.3%)
Oral cavity3 (2.6%)
Table 2. The distribution of tumor and clinical characteristics.
Table 2. The distribution of tumor and clinical characteristics.
Characteristics Frequency
Grade114 (12.2%)
287 (75.7%)
313 (11.3%)
41 (0.9%)
SizeT136 (31.3%)
T222 (19.1%)
T327 (23.5%)
T430 (26.1%)
Cancerous regional lymph nodesN070 (60.9%)
N145 (39.1%)
MetastasisAbsent111 (96.5%)
Present4 (3.5%)
StageI32 (27.8%)
II23 (20.0%)
III17 (14.8%)
IV43 (37.4%)
Death factAlive71 (61.7%)
Deceased44 (38.3%)
Table 3. The distribution of TFAM and POLG genotype and allele frequencies.
Table 3. The distribution of TFAM and POLG genotype and allele frequencies.
Gene and SNPCountGenotype FrequencyAllele Frequency
TFAM
rs11006132
56AA—0.49
AG—0.47
GG—0.04
A—0.72
G—0.28
54
5
TFAM
rs11006129
89CC—0.77
CT—0.23
C—0.89
T—0.11
26
TFAM
rs1937
83GG—0.72
CG—0.27
CC—0.01
G—0.86
C—0.14
31
1
TFAM
rs16912174
107TT—0.93
GT—0.07
T—0.97
G—0.03
8
TFAM
rs1692202 *
115TT—1.00T—1.00
TFAM
rs3900887 *
84TT—0.76
AT—0.17
AA—0.07
T—0.84
A—0.16
19
8
POLG
rs3087374
9520CC—0.83
AC—0.17
A—0.09
C—0.91
POLG
rs2072267
12AA—0.10
AG—0.54
GG—0.36
A—0.37
G—0.63
62
41
POLG
rs976072
20AA—0.17
AG—0.56
GG—0.27
A—0.45
G—0.55
64
31
POLG
rs2307441
97TT—0.84
CT—0.16
T—0.92
C—0.08
18
The table presents the genotype distribution within our study cohort. Gene and single-nucleotide polymorphism (SNP) rs numbers are provided. Patient count together with the genotype frequency is provided for each genotype within the cohort, together with the allele frequency for every polymorphism. * Genotype distribution is not in Hardy-Weinberg equilibrium.
Table 4. TFAM gene genotype associations with disease and tumor characteristics.
Table 4. TFAM gene genotype associations with disease and tumor characteristics.
Variablep-Value
rs3900887rs11006132rs11006129rs1937rs16912174
Stage (I, II, III, IV)0.8150.1450.5770.5340.455
T (T1, T2, T3, T4)0.4430.1490.8240.6340.480
T (T1, T2 and T3, T4 groups)0.1740.8770.2360.5960.349
N (negative vs. positive)0.7510.9980.4370.6270.596
M (negative vs. positive)0.0070.5500.0360.9770.746
Differentiation grade G (G1, G2, G3, G4)0.0660.2840.0360.9480.147
Differentiation grade G (G1, G2 vs. G3, G4)0.4860.1320.3410.6940.658
Survival status0.0180.6910.0190.3210.638
T—tumor size; N—regional lymph node status; M—metastasis status; G—tumor grade.
Table 5. POLG gene genotype associations with disease and tumor characteristics.
Table 5. POLG gene genotype associations with disease and tumor characteristics.
Variablep-Value
rs3087374rs2307441rs2072267rs976072
Stage (I, II, III, IV)0.0540.7560.3140.102
T (T1, T2, T3, T4)0.0170.7880.4940.934
T (T1, T2 and T3, T4 groups)0.4200.4150.2090.958
N (negative vs. positive)0.1990.4010.3360.996
M (negative vs. positive)0.1390.4990.4220.437
Differentiation grade G (G1, G2, G3, G4)0.8970.8020.9190.833
Differentiation grade G (G1, G2 vs. G3, G4)0.4550.3140.8730.585
Survival status0.0760.1020.4490.646
T—tumor size; N—regional lymph node status; M—metastasis status; G—tumor grade.
Table 6. Univariate and multivariate logistic regression analyses between TFAM rs3900887 and rs11006129 and fact of death.
Table 6. Univariate and multivariate logistic regression analyses between TFAM rs3900887 and rs11006129 and fact of death.
DependentSNPCovariatesModel No. 1Model No. 2
OR95% CIpOR95% CIp
Survival statusrs3900887AT vs. TT4.5651.55413.4140.0065.8911.69820.4470.005
Age *0.9890.9810.9960.0020.9420.9150.9710.001
T 1.3440.8572.1080.198
N 3.0411.0318.9650.044
G (G1, G2 vs. G3, G4) 4.6461.26417.0790.021
Survival statusrs3900887A-allele carriers vs. non-carriers2.9401.1987.2140.0193.6131.27410.2500.016
Age *0.9890.9820.9960.0030.9860.9421.0320.551
T 1.6511.0502.5960.030
N 2.8871.0218.1650.046
G (G1, G2 vs. G3, G4) 6.5891.59327.2620.009
Survival statusrs11006129CT vs. CC2.8751.1667.0860.0223.6791.28810.5110.015
Age *0.9890.9820.9960.0010.9860.9431.0320.546
T 1.6651.0572.6220.028
N 2.9281.0458.2050.041
G (G1, G2 vs. G3, G4) 6.9381.67828.6860.007
Survival statusrs11006129T-allele carriers vs. non-carriers2.8751.1667.0860.0223.9301.36211.3410.011
Age *0.9890.9820.9960.0010.9890.9451.0350.624
T 1.6191.0372.5280.034
N 1.9851.1453.4420.015
G (G1, G2 vs. G3, G4) 6.9971.67929.1620.008
* Age at the time of diagnosis; OR—odds ratio; CI—confidence interval; T—tumor size; N—regional lymph node status; G—tumor grade. Model No. 1—logistic regression analysis adjusted for age at diagnosis. Model No. 2—logistic regression analysis adjusted for age at diagnosis, tumor size, regional lymph node status, and tumor grade.
Table 7. Cox multivariate model with additional cofactors for patients’ overall survival.
Table 7. Cox multivariate model with additional cofactors for patients’ overall survival.
VariablesHR95% CIp-Value
rs11006129
T-allele carriers2.9821.543–5.7650.001
Age *0.9990.965–1.0350.958
T1.7041.239–2.3420.001
N1.6040.794–3.2390.188
G (G1, G2 vs. G3, G4)2.7431.291–5.8310.009
rs3900887
A-allele carriers2.7511.437–5.2660.002
Age *1.0010.967–1.0380.934
T1.6721.221–2.2890.001
N1.5370.756–3.1270.235
G (G1, G2 vs. G3, G4)2.6341.236–5.6120.012
* Age at the time of diagnosis; HR—hazard ratio; CI—confidence interval; T—tumor size; N—regional lymph node status; G—tumor grade.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Golubickaite, I.; Ugenskiene, R.; Bartnykaite, A.; Poskiene, L.; Vegiene, A.; Padervinskis, E.; Rudzianskas, V.; Juozaityte, E. Mitochondria-Related TFAM and POLG Gene Variants and Associations with Tumor Characteristics and Patient Survival in Head and Neck Cancer. Genes 2023, 14, 434. https://doi.org/10.3390/genes14020434

AMA Style

Golubickaite I, Ugenskiene R, Bartnykaite A, Poskiene L, Vegiene A, Padervinskis E, Rudzianskas V, Juozaityte E. Mitochondria-Related TFAM and POLG Gene Variants and Associations with Tumor Characteristics and Patient Survival in Head and Neck Cancer. Genes. 2023; 14(2):434. https://doi.org/10.3390/genes14020434

Chicago/Turabian Style

Golubickaite, Ieva, Rasa Ugenskiene, Agne Bartnykaite, Lina Poskiene, Aurelija Vegiene, Evaldas Padervinskis, Viktoras Rudzianskas, and Elona Juozaityte. 2023. "Mitochondria-Related TFAM and POLG Gene Variants and Associations with Tumor Characteristics and Patient Survival in Head and Neck Cancer" Genes 14, no. 2: 434. https://doi.org/10.3390/genes14020434

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop