Next Article in Journal
The Presence of Ganglionic Acetylcholine Receptor Antibodies in Sera from Patients with Functional Gastrointestinal Disorders: A Preliminary Study
Previous Article in Journal
A Comparative Analysis of Naïve Exosomes and Enhanced Exosomes with a Focus on the Treatment Potential in Ovarian Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Profile of Important Genes for Radiogenomics in the Amazon Indigenous Population

by
Milena Cardoso de Lima
1,
Cinthia Costa de Castro
1,
Kaio Evandro Cardoso Aguiar
1,
Natasha Monte
1,
Giovanna Gilioli da Costa Nunes
1,
Ana Caroline Alves da Costa
1,
Juliana Carla Gomes Rodrigues
1,
João Farias Guerreiro
1,2,
Ândrea Ribeiro-dos-Santos
2,
Paulo Pimentel de Assumpção
1,
Rommel Mario Rodríguez Burbano
1,
Marianne Rodrigues Fernandes
1,
Sidney Emanuel Batista dos Santos
1,2 and
Ney Pereira Carneiro dos Santos
1,*
1
Oncology Research Center, Federal University of Pará, Belém 66073-005, PA, Brazil
2
Laboratory of Human and Medical Genetics, Federal University of Pará, Belém 66075-110, PA, Brazil
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2024, 14(5), 484; https://doi.org/10.3390/jpm14050484
Submission received: 20 March 2024 / Revised: 13 April 2024 / Accepted: 18 April 2024 / Published: 30 April 2024
(This article belongs to the Section Omics/Informatics)

Abstract

:
Radiotherapy is focused on the tumor but also reaches healthy tissues, causing toxicities that are possibly related to genomic factors. In this context, radiogenomics can help reduce the toxicity, increase the effectiveness of radiotherapy, and personalize treatment. It is important to consider the genomic profiles of populations not yet studied in radiogenomics, such as the indigenous Amazonian population. Thus, our objective was to analyze important genes for radiogenomics, such as ATM, TGFB1, RAD51, AREG, XRCC4, CDK1, MEG3, PRKCE, TANC1, and KDR, in indigenous people and draw a radiogenomic profile of this population. The NextSeq 500® platform was used for sequencing reactions; for differences in the allelic frequency between populations, Fisher’s Exact Test was used. We identified 39 variants, 2 of which were high impact: 1 in KDR (rs41452948) and another in XRCC4 (rs1805377). We found four modifying variants not yet described in the literature in PRKCE. We did not find any variants in TANC1—an important gene for personalized medicine in radiotherapy—that were associated with toxicities in previous cohorts, configuring a protective factor for indigenous people. We identified four SNVs (rs664143, rs1801516, rs1870377, rs1800470) that were associated with toxicity in previous studies. Knowing the radiogenomic profile of indigenous people can help personalize their radiotherapy.

1. Introduction

Radiotherapy (RT) is one of the most important treatments for cancer. It is used in more than half of cancer patients and can have a curative or palliative goal by reducing the rate of local recurrence and improving overall survival. In RT, the tumor tissue is irradiated with ionizing radiation, which can penetrate the adjacent healthy tissue and cause toxic effects there, as the radiation damages the DNA [1,2]. In this sense, the toxic effects, adverse effects caused by the intensity of the dose, lead to losses both for the individual and due to the poor quality of life. Particular attention must be paid to childhood cancers, as these can be long-lasting, with correspondingly high impacts of radiotherapy and on the healthcare system due to the increase in public expenditure [1,3].
Radiotoxicity (RTX) may vary depending on the tissue irradiated [1]. Mucositis, xerostomia, subcutaneous fibrosis, and dysphagia [4] occur with squamous cell carcinoma; esophagitis may occur with radiation to the lung [5]. RT to the chest can cause pain, fibrosis, upper limb edema, pericarditis, and ischemic and valvular heart disease [6]. RT in the prostate can lead to dermatitis, urinary incontinence, and bladder and erectile dysfunction [7].
Radiosensitivity is a tendency to develop adverse effects in tissues under the effects of cell death due to radiation; this generates resistance to treatment and is influenced by biological or genetic factors, for example, tumor heterogeneity [8,9]. It is understood that radiosensitivity driven by genetic variation is one of the fundamental factors for radiotoxicity [10]. Given this, radiogenomics searches for biomarkers of the response to radiotherapy have been used to personalize treatments among individuals undergoing radiotherapy according to their radiogenomic profile [3,10].
RT dosing follows the guidelines of the Radiation Therapy Oncology Group, which has established the same standards for the same tumor sites. However, due to individual differences, some doses may be excessive or insufficient [11]. This scenario may be even more complicated when considering the profile of genetically diverse populations such as indigenous populations, due to their process of geographic isolation and to genetic and evolutionary factors such as genetic drift, founder effects, natural selection, and inbreeding [12].
Therefore, it is pertinent to analyze important genes related to radiogenomics in the Amazonian indigenous population, to obtain information regarding the genomic profile of this population regarding the response to radiotherapy and, in the future, to obtain pertinent results through the control of mutations. Thus, the objective of the study was to analyze the exome of the ATM, TGFB1, RAD51, AREG, XRCC4, CDK1, MEG3, PRKCE, TANC1, and KDR genes in the indigenous population and draw a genomic profile of this population, which has been poorly studied with regard to radiogenomics.
We emphasize that the Brazilian population is one of the most heterogeneous in the world, marked by indigenous, African, and European ancestry; the indigenous population contributes significantly to the formation of Latin American groups, with a share of about 30, especially in the northern region of Brazil [13,14]. In addition, studies conducted with indigenous peoples from other countries have shown high rates of cancer deaths [15]. On the other hand, there have been few genetic studies aiming to understand oncologic mechanisms and therapies that have dealt with this group. We can also point out that there have been no studies evaluating the impact of variants important for radiogenomics on the indigenous population and the population mixed with it.

2. Materials and Methods

2.1. Population Analysis for the Study

The study participants and their ethnic leaders signed the Free and Informed Consent Form. This study was approved by the National Ethics and Research Commission (CONEP) and the Research Ethics Committee of the Center for Tropical Medicine of the Federal University of Pará (UFPA) under CAAE 20654313.6.0000.5172. The population is made up of 64 indigenous people from the Amazon region of Brazil, representing 12 ethnicities: Asurini do Xingu, Arara, Araweté, Asurini do To-cantins, Awagujá, Kayapó/Xikrin, Zo’é, Wajãpi, Karipuna, Phurere, Munduruku, and Ju-runa. The population was healthy and did not undergo radiotherapy. The genetic ancestry of the population group was obtained through a panel of 61 ancestry informative markers (AIMs), which were used to identify individual ancestry and a mixture of three continents (European, African, and Amerindian). The results obtained were associated with other world populations, according to information available in version 3 of the 1000 Genomes database (available online at http://www.1000genomes.org, accessed on 25 October 2023) and the Exome Aggregation Consortium (ExAC). These populations are composed of 661 Africans (AFR), 346 Americans (AMR), 504 East Asians (EAS), 503 Europeans (EUR), and 489 South Asians (SAS).

2.2. DNA Extraction and Exome Analysis

The method used for DNA extraction was phenol–chloroform extraction [16]. Verification of the integrity of the sample was carried out using a NanoDrop 8000 spectrophotometer, and tracing was performed using the Nextera Rapid Capture Exome (Illumina®, San Diego, CA, USA) and SureSelect Human All Exon V6 kits (Agilent, Santa Clara, CA, USA). The NextSeq 500® platform (Illumina®, San Diego, CA, USA) was used for sequencing reactions using the NextSeq 500 High-output v2 300 cycle kit (Illumina®, San Diego, CA, USA).

2.3. Selection of Genes and Variants

The first criterion for variant consideration was 10 minimum coverage values, and then the impact of each variant was considered according to the SNPeff v 4.3 classification (https://pcingola.github.io/SnpEff/; accessed on 20 February 2024). In addition, we considered only variants with a significant allelic frequency in at least three continental populations. As a result of the exome analysis, 143 variants were found. After selection based on the above criteria, 39 variants remained to be followed up in this study.
Ten genes were selected for this study based on their association with radiotherapy efficacy and/or toxicity; some have been included in previous studies and are indexed in databases such as PubMed. The genes are described in Table 1.

2.4. Bioinformatics Analysis

FASTQ was used to analyze the quality of the reads (FastQCv.0.11—https://www.bioinformatics.babraham.ac.uk/projects/fastqc/; accessed on 20 February 2024), and filters were applied to the samples in order to disregard low-quality reads (fastx_tools v.0.13—http://hannonlab.cshl.edu/fastx_toolset/; accessed on 20 February 2024). The reference genome (GRCH38) was used to map and align the samples, using BWA v.0.7 (http://bio-bwa.sourceforge.net/; accessed on 20 February 2024). The file was indexed and classified according to the alignment generated by the reference genome (SAMtools v.1.2—http://sourceforge.net/projects/samtools/; accessed on 20 February 2024). From there, the alignment was processed to remove PCR duplication (Picard Tools v.1.129—http://broadinstitute.github.io/picard/; accessed on 20 February 2024) and to perform structuring quality readjustment and local realignment (GATK v. 3.2—https://www.broadinstitute.org/gatk/; accessed on 20 February 2024). Then, the results were processed to establish the reference genome variants (GATK v.3.2). ViVa1 (Viewer of Variants) software (v1)—developed by the bioinformatics team at the Federal University of Rio Grande do Norte (UFRN)—was used to analyze the variant annotations. The databases employed for variant annotations were SnpEff v.4.3., Ensembl Variant Effect Predictor (Ensembl version 99), and ClinVar (v.2018-10). For in silico predictions of pathogenicity, SIFT (v.6.2.1), PolyPhen-2 (v.2.2), LRT (November 2009), Mutation Evaluator (v.3.0), and Mutation Tester (v. 2.0) were used. FATHMM (v.2.3), PROVEAN (v.1.1.3), MetaSVM (v1.0), M-CAP (v1.4), and FATHMM-MKL (http://fathmm.biocompute.org.uk/about.html; accessed on 20 February 2024) were also used. More information about bioinformatics analysis is given in the works of Rodrigues et al. [37] and Ribeiro-dos-Santos et al. [38].

2.5. Statistical Analysis

Fisher’s Exact Test was used to determine statistically significant differences (p value ≤ 0.05) between the frequencies of the variants in the world populations (AFR, EUR, AMR, EAS, and SAS). The allele frequencies of the variants found in the study population were determined by allele enumeration. Comparisons between the allele frequencies of the variants found in the indigenous population and in the five world populations were performed using Fisher’s Exact Test. To avoid discrepancies due to different population sizes, we leveled all variant frequencies with the size of the indigenous population. The frequencies of the five populations were taken from Phase 3 of the 1000 Genomes Project (http://www.1000genomes.org; accessed on 20 February 2024) and the Exome Aggregation Consortium (ExAC). These allele frequencies were also used to calculate the dissimilarity matrix and to generate the multidimensional scaling (MDS). The analyses were performed with RStudio v. 3.5.1.

3. Results

From the analyses, we identified 39 variants distributed across the 9 genes: 1 belonging to the AREG gene, 6 to the ATM gene, 9 to the KDR gene, 4 to the MEG3 gene, 10 to the PRKCE gene, 3 to the RAD51 gene, 1 to the TGFB1 gene, 4 to the TANC1 gene, and 1 to the XRCC4 gene. Figure 1 presents the distribution of variants according to their high, modifier, or moderate impact; thus, 2 were classified as high-impact, 28 as modifiers, and 5 as moderate.
Two high-impact variants were identified, distributed in the KDR and XRCC4 genes; for the modifying variants, 3.12% (1) were in the AREG gene, 15.62% (5) in ATM, 21.87% (7) in KDR, 12.50% (4) in MEG3, 31.25% (10) in PRKCE, 9.37% (3) in RAD51, and 6.25% (2) in TANC1. Regarding the moderate variants, ATM, KDR, and TGFB1 each had one moderate variant, while two were in the TANC1 gene. The variants are described in Table 2 with information about the ID, the impact predicted by SNPeff software v 4.3, the region, the type, the nucleotide change, and the frequency presented in world populations. The low-impact variants are presented in Supplementary Table S1.
As for the variants not yet described in the literature, four were found, all in the PRKCE gene on chromosome 2; these are listed in Table 3 and are characterized by a single nucleotide change and an insertion/deletion. The variants are sequentially located at 45652078, 45652096, 45652087, and 45652092. Of these variants, three had a frequency of 0.0833 in the indigenous population and one had a frequency of 0.01677. We intend to conduct further studies to further explore the effects of these variants on protein activity to obtain additional information on radiogenomics in a poorly studied population.
Of the variants described in Table 2, rs41452948 and rs1805377 are of high impact, the first occurring in the KDR gene and the second in the XRCC4 gene, the first on chromosome 4 and the second on chromosome 5. Both are represented by a change in nucleotide from G to A. These variants showed a significant allelic frequency in the indigenous population compared to the AFR, AMR, EUR, and SAS populations and a high, modifying, and moderate impact, as shown in Table 3 Comparisons with low-impact variants can be found in Supplementary Table S2.
In addition, 43 variants showed a nonzero allelic frequency in at least three global populations. However, the differences in frequency in the five continental populations compared to the indigenous population were significant for only 18 of these alleles, as described in Table 4.
Figure 2 shows the different genotypes in the global population and the indigenous population on a multidimensional scale (MDS) determined using Fisher’s Exact Test. Looking at the genes selected for this exome, our results show that the indigenous population (INDG) has a genetic profile that differs significantly from those of the five world populations, especially those of East Asians (EAS) and Africans (AFR).

4. Discussion

The curative potential of radiotherapy is limited by the intrinsic radioresistance of tumor cells, which is related to the heterogeneity of the tumor and the surrounding microenvironment, as well as to various genetic alterations [39]. In addition to radioresistance, patients may experience clinically significant side effects, namely, acute reactions such as erythema, follicular reaction, pruritus, moist or dry desquamation, ulceration, and necrosis, or delayed reactions such as tissue fibrosis, atrophy, firm subcutaneous tissue, and subcutaneous swelling [40,41].
GWA studies are often used to identify genetic variants associated with complex and multifactorial diseases [42,43]. GWA studies linking radiogenomics to radiotherapy toxicity show that genetic variants, treatment variables, and other clinical factors are independent predictors of radiotoxicity. This contributes to the suggestion that common variants may improve traditional models of the likelihood of complications in normal tissues [44,45]. When DNA is exposed to radiation, DNA double-strand breaks occur, which favors a response to this damage, such as cellular radioresistance [46,47]. This was the first exome study to examine genes associated with the response to radiotherapy in indigenous populations, a group that is usually underrepresented in genetic studies, especially in precision medicine [48].
In our results, we found 39 variants with nonzero allele frequencies in the study population, 4 of which have never been found in other genetic databases and none of which are in the CDK1 gene. Therefore, of the 35 variants present in all five populations, we highlight two high-impact variants, rs41452948 in the KDR gene and rs1805377 in the XRCC4 gene, as presenting statically significant differences when comparing the indigenous population with African, American, European, and South Asian populations, as shown in Table 4.
We emphasize that many variants were identified in intronic regions, where they were believed not to encompass gene functions. However, intronic regions play an important role in genome regulation. For example, they influence the maturation of mRNA and transcriptional regulation through alternative splicing, which can directly influence a gene’s response to radiotoxicity [49]. The rs1805377 in the splice site acceptor + intron could have consequences such as disruption of the normal splicing process, leading to important changes in parts of the final mRNA messenger gene, which, in turn, could cause genetic disorders or disease, depending on the gene affected and the specific type of variant [50].
The variants found only in the indigenous population are all located in the PRKCE gene, which is described in the literature as a gene associated with the response to radiation [27]. Members of the PKC family phosphorylate a variety of protein targets and are involved in cell signaling pathways [26]. These variants have a modifying clinical effect; three of the variants had an allelic frequency of 0.0833 in indigenous people and one had a frequency of 0.1667. These mutations with a modifying effect could be potential markers for indigenous populations in the Amazon region. However, we emphasize the intention to conduct further studies to consolidate the knowledge of these potential new variants.
Polymorphisms in the RAD51 gene have been associated in previous studies with symptoms of heart failure and the appearance of a new primary tumor in patients with breast cancer treated with trastuzumab and RT [6], suggesting an important gene for RTX; however, such polymorphisms were not detected in the population studied. Other studies have also investigated variations in DNA repair genes—such as ATM, RAD51, and XRCC4—that cause differences in the response to radiation between individuals. These differences are characterized by increased toxicity, which normally emanates from irradiated tissue [51], but these variants were also not detected in the Indigenous population, suggesting protection from these toxicities.
We discovered four variants in the TANC1 gene, two with a moderate impact (rs34588551, rs4664277) and two with a modifying impact in the intronic region (rs34344829, rs146371641). In our analyses, we found a statistically significant difference regarding rs146371641 in all five world populations compared to the indigenous group (Table 4), which is related to the fact that TANC1 is an important gene for radiogenomics; this could show that the indigenous population has a different profile in the use of radiotherapy.
TANC1 plays a crucial role in the recruitment of fusion-capable myoblasts during myotube formation. In addition, expression of this gene has been found in various tissues such as adipose tissue and the adrenal cortex [3,31]. TANC1 was associated with late radiotherapy toxicity in a large genomic study cohort with a European population. The results suggested that this gene is involved in muscle damage processes induced by radiation. However, in our study, the variants presented in the European study were not found, which may be explained by the genomic difference between the population originating from Brazil and the European population [37]. The influencing variant rs664143, located in the ATM intron, showed significant differences in the five world populations. This variant showed a worse prognosis in patients with squamous cell carcinoma of the esophagus and non-small cell lung cancer after radiotherapy. In addition, the A allele was associated with shorter survival, while the G allele indicated an increased risk of disease progression [52].
The rs1801516 variant, which has a moderate impact on ATM, showed significant differences in the indigenous population compared to the American, European, and South Asian populations (see Table 4). This variant has been studied with regard to the effects after radiotherapy, with results suggesting that individuals with the minor allele have an increased risk of developing late fibrosis after radiotherapy [53]. In studies by McDuff et al. [54], Andreassen et al. [55], and Kerns et al. [45], the variant was associated with the risk of general toxicity, acute toxicity, late toxicity, acute skin toxicity, acute rectal toxicity, telangiectasia, and fibrosis.
The variant with a moderate impact on the KDR gene, rs1870377, showed a significant difference, especially when comparing the indigenous and European populations (Table 4). A study by Tinhofer et al. [56] showed that despite the poorer prognosis of patients with squamous cell carcinoma of the head and neck, the A > T allele of germline variant rs1870377 was associated with longer survival in risk groups of patients with stage IV squamous cell carcinoma of the oropharynx and hypopharynx who underwent chemotherapy. On the other hand, Butkiewicz et al. [57] analyzed 422 patients with squamous cell carcinoma of the head and neck undergoing radiotherapy and concluded that VEGFR2 rs1870377 TT is a significant borderline risk factor for lower local recurrence and an independent predictor of poor prognosis.
When analyzing the toxicity of the rs1870377 variant in patients undergoing 5-fluorouracil-based chemotherapy for locally advanced rectal cancer, the VEGFR2 H472Q Q/Q genotype was associated with a higher risk of grade 3 mucositis of the proximal upper gastrointestinal tract in arm 2 of the study (induction and concomitant prolonged intravenous infusion of 5-FU with radiotherapy). However, in arm 1 (bolus of 5-FU followed by prolonged intravenous infusion of 5-FU with radiotherapy), this genotype was associated with a lower risk of mucositis [22].
The variant with a moderate influence on the TGFB1 gene, rs1800470, showed statistically significant differences between the five global populations and the indigenous population (Table 4). According to Xiao Y et al. [58], the T allele of variant rs1800470 (CT/TT) is associated with the risk of developing radiation-induced pneumonia in patients with esophageal squamous cell carcinoma. In patients with oropharyngeal squamous cell carcinoma undergoing radiotherapy, the presence of the rs1800470 variant of TGFβ1 may reduce and modify the risk of death and recurrence [59].
Previous research also found mitigating differences between the genetic profile of indigenous peoples and those of the other five world populations examined here [60,61,62]. However, the physical and geographic isolation of these peoples contributes to a gap in epidemiologic and genetic information. Therefore, knowledge of their molecular profile may be of importance to indigenous populations, as well as those intermixed with them [63].
In the European world, the applications of radiogenomics concepts are more advanced, considering studies such as those by O’Sullivan, NJ; Kelly [64], which used models to predict metastases after neoadjuvant radiotherapy and showed how these concepts can be used to identify low-risk patients in order not to subject them to interventions such as tumor dissection. In addition, radiogenomics has also been used to predict genetic mutations in colorectal cancer to optimize the outcome of radiotherapy and enable targeted therapy [65].

5. Conclusions

This study investigated the presence of genetic variants involved in the process of radiotherapy toxicity in an indigenous Amazonian population. No variants already associated with late RTX were found. We found two high-impact variants and four variants that have not yet been described in the other five world populations examined and may be exclusive to this population. Future research on therapeutic approaches targeting indigenous and admixed populations is critical.
Personalized medicine research has gained much attention for its role in targeting therapies and making predictions. This can help doctors plan more appropriate treatment regimens, reduce radiation damage in high-risk patients, and optimize outcomes in low-risk patients.

Supplementary Materials

The following support information can be downloaded from: https://www.mdpi.com/article/10.3390/jpm14050484/s1. Table S1: Descriptions of low-impact AREG, ATM, KDR, MEG3, PRKCE, RAD51, TGFB1, TANC1, and XRCC4 gene variants, as well as allele frequencies in continental populations (African (AFR), American (AMR), East Asian (EAS), European (EUR), and South Asian (SAS)) described in the 1000 genome database. Table S2: Paired comparisons (p-values) with significant results of allele frequencies in the Amazonian indigenous population (INDG) and in the continental populations (African (AFR), American (AMR), East Asian (EAS), European (EUR), and South Asian (SAS)) described in the 1000 genome database.

Author Contributions

Conceptualization: M.C.d.L. and N.P.C.d.S. methodology: N.M., J.C.G.R. and N.P.C.d.S.; formal analysis: N.M. and J.C.G.R.; writing—original draft preparation: M.C.d.L., C.C.d.C., K.E.C.A., G.G.d.C.N. and A.C.A.d.C.; writing—review and editing: M.C.d.L. and N.P.C.d.S.; supervision: J.F.G., Â.R.-d.-S., P.P.d.A., R.M.R.B., M.R.F., S.E.B.d.S. and N.P.C.d.S.; project administration: N.P.C.d.S. and J.F.G.; resources: Â.R.-d.-S., P.P.d.A. and S.E.B.d.S.; validation: R.M.R.B. and M.R.F. All authors have read and agreed to the published version of the manuscript.

Funding

Pró-Reitoria de Pesquisa e Pós-Graduação (PROPESP) and the Federal University of Pará (UFPA).

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Research Ethics Committee of the Health Sciences Institute of the Federal University of Pará and by the National Research Ethics Committee (CONEP), identification number 20654313.6.0000.5172, approved on 2 March 2015.

Informed Consent Statement

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

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Acknowledgments

We acknowledge Universidade Federal do Pará (UFPA), the Oncology Research Center (NPO/UFPA), and the Human and Medical Genetics Laboratory (LGHM/UFPA).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. De Ruysscher, D.; Niedermann, G.; Burnet, N.G.; Siva, S.; Lee, A.W.; Hegi-Johnson, F. Radiotherapy toxicity. Nat. Rev. Dis. Primers 2019, 5, 13. [Google Scholar] [CrossRef] [PubMed]
  2. Early Breast Cancer Trialists’ Collaborative Group (EBCTCG); Darby, S.; McGale, P.; Correa, C.; Taylor, C.; Arriagada, R.; Clarke, M.; Cutter, D.; Davies, C.; Ewertz, M.; et al. Effect of radiotherapy after breast-conserving surgery on recurrence at 10 years and death from breast cancer at 15 years: Meta-analysis of individual patient data for 10,801 women in 17 randomized trials. Lancet 2011, 378, 1707–1716. [Google Scholar] [PubMed]
  3. Rosenstein, B.S. Radiogenomics: Identification of Genomic Predictors for Radiation Toxicity. Semin. Radiat. Oncol. 2017, 27, 300–309. [Google Scholar] [CrossRef]
  4. Thiagarajan, A.; Iyer, N.G. Genomics of radiation sensitivity in squamous cell carcinomas. Pharmacogenomics 2019, 20, 457–466. [Google Scholar] [CrossRef]
  5. Holgersson, G.; Bergström, S.; Liv, P.; Nilsson, J.; Edlund, P.; Blomberg, C.; Nyman, J.; Friesland, S.; Ekman, S.; Asklund, T.; et al. Effect of Increased Radiotoxicity on Survival of Patients with Non-small Cell Lung Cancer Treated with Curatively Intended Radiotherapy. Anticancer Res. 2015, 35, 5491–5497. [Google Scholar] [PubMed]
  6. Goričar, K.; Dugar, F.; Dolžan, V.; Marinko, T. NBN, RAD51, and XRCC3 Polymorphisms as Potential Predictive Biomarkers of Adjuvant Radiotherapy Toxicity in Early HER2-Positive Breast Cancer. Cancers 2022, 14, 4365. [Google Scholar] [CrossRef] [PubMed]
  7. Mališić, E.; Petrović, N.; Brengues, M.; Azria, D.; Matić, I.Z.; Srbljak Ćuk, I.; Kopčalić, K.; Stanojković, T.; Nikitović, M. Association of polymorphisms in TGFB1, XRCC1, XRCC3 genes and CD8 T-lymphocyte apoptosis with adverse effect of radiotherapy for prostate cancer. Sci. Rep. 2022, 12, 21306. [Google Scholar] [CrossRef] [PubMed]
  8. El-Nachef, L.; Al-Choboq, J.; Restier-Verlet, J.; Granzotto, A.; Berthel, E.; Sonzogni, L.; Ferlazzo, M.L.; Bouchet, A.; Leblond, P.; Combemale, P.; et al. Human radiosensitivity and radiosusceptibility: What are the differences? Int. J. Mol. Sci. 2021, 22, 7158. [Google Scholar] [CrossRef] [PubMed]
  9. Guhlich, M.; Hubert, L.; Mergler, C.P.N.; Rave-Fraenk, M.; Dröge, L.H.; Leu, M.; Schmidberger, H.; Rieken, S.; Hille, A.; Schirmer, M.A. Identification of Risk Loci for Radiotoxicity in Prostate Cancer by Comprehensive Genotyping of TGFB1 and TGFBR1. Cancers 2021, 13, 5585. [Google Scholar] [CrossRef]
  10. Jiang, M.; Yang, J.; Li, K.; Liu, J.; Jing, X.; Tang, M. Insights into the theranostic value of precision medicine on advanced radiotherapy to breast cancer. Int. J. Med. Sci. 2021, 18, 626–638. [Google Scholar] [CrossRef]
  11. Bai, X.; Wang, B.; Wang, S.; Wu, Z.; Gou, C.; Hou, Q. Radiotherapy dose distribution prediction for breast cancer using deformable image registration. BioMedical Engineering OnLine 2020, 19, 39. [Google Scholar] [CrossRef] [PubMed]
  12. Amorim, C.E.G.; Wang, S.; Marrero, A.R.; Salzano, F.M.; Ruiz-Linares, A.; Bortolini, M.C. X-Chromosomal Genetic Diversity and Linkage Disequilibrium Patterns in Amerindians and Non-Amerindian Populations. Am. J. Hum. Biol. 2011, 23, 299–304. [Google Scholar] [CrossRef] [PubMed]
  13. BRAZIL. Ministry of Indigenous Peoples. Data from the 2022 Census Reveals that Brazil has 1.7 Million Indigenous People. 2023. Available online: https://en.wikipedia.org/wiki/Ministry_of_Native_People_(Brazil) (accessed on 20 February 2024).
  14. Rodrigues, J.C.G.; Fernandes, M.R.; Guerreiro, J.F.; da Silva, A.L.D.C.; Ribeiro-Dos-Santos, Â.; Santos, S.; Santos, N.P.C.D. Polymorphisms of ADME-related genes and their implications for drug safety and efficacy in Amazonian Amerindians. Sci. Rep. 2019, 9, 7201. [Google Scholar] [CrossRef]
  15. Moore, S.P.; Forman, D.; Piñeros, M.; Fernández, S.M.; de Oliveira Santos, M.; Bray, F. Cancer in indigenous people in Latin America and the Caribbean: A review. Cancer Med. 2014, 3, 70–80. [Google Scholar] [CrossRef] [PubMed]
  16. Green, M.R.; Sambrook, J. Isolation of high-molecular-weight DNA using organic solvents. Cold Spring Harb. Protoc. 2017, 2017, pdb-prot093450. [Google Scholar] [CrossRef] [PubMed]
  17. Bouchet, A.; Sakakini, N.; Atifi, M.E.; Le Clec’h, C.; Bräuer-Krisch, E.; Rogalev, L.; Laissue, J.A.; Rihet, P.; Le Duc, G.; Pelletier, L. Identification of AREG and PLK1 pathway modulation as a potential key of the response of intracranial 9L tumor to microbeam radiation therapy. Int. J. Cancer 2015, 136, 2705–2716. [Google Scholar] [CrossRef] [PubMed]
  18. García, M.E.G.; Kirsch, D.G.; Reitman, Z.J. Targeting the ATM Kinase to Enhance the Efficacy of Radiotherapy and Outcomes for Cancer Patients. Semin. Radiat. Oncol. 2022, 32, 3–14. [Google Scholar] [CrossRef]
  19. Liu, J.; Tang, X.; Shi, F.; Li, C.; Zhang, K.; Liu, J.; Wang, G.; Yin, J.; Li, Z. Genetic polymorphism contributes to 131I radiotherapy-induced toxicities in patients with differentiated thyroid cancer. Pharmacogenomics 2018, 19, 1335–1344. [Google Scholar] [CrossRef]
  20. Massacci, G.; Perfetto, L.; Sacco, F. The Cyclin-dependent kinase 1: More than a cell cycle regulator. Br. J. Cancer 2023, 129, 1707–1716. [Google Scholar] [CrossRef]
  21. Huang, R.X.; Zhou, P.K. DNA damage response signaling pathways and targets for radiotherapy sensitization in cancer. Signal Transduct. Target. Ther. 2020, 5, 60. [Google Scholar] [CrossRef]
  22. Zaman, N.; Dass, S.S.; Du Parcq, P.; Macmahon, S.; Gallagher, L.; Thompson, L.; Khorashad, J.S.; LimbÄck-Stanic, C. The KDR (VEGFR-2) Genetic Polymorphism Q472H and c-KIT Polymorphism M541L Are Associated with More Aggressive Behaviour in Astrocytic Gliomas. Cancer Genom. Proteom. 2020, 17, 715–727. [Google Scholar] [CrossRef]
  23. Bohanes, P.; Rankin, C.J.; Blanke, C.D.; Winder, T.; Ulrich, C.M.; Smalley, S.R.; Rich, T.A.; Martensen, J.A.; Benson, A.B., III; Mayer, R.J.; et al. Pharmacogenetic Analysis of INT 0144 Trial: Association of Polymorphisms with Survival and Toxicity in Rectal Cancer Patients Treated with 5-FU and Radiation. Clin Cancer Res. 2015, 21, 1583–1590. [Google Scholar] [CrossRef]
  24. Ashrafizadeh, M.; Gholami, M.H.; Mirzaei, S.; Zabolian, A.; Haddadi, A.; Farahani, M.V.; Kashani, S.H.; Hushmandi, K.; Najafi, M.; Zarrabi, A.; et al. Dual relationship between long non-coding RNAs and STAT3 signaling in different cancers: New insight to proliferation and metastasis. Life Sci. 2021, 270, 119006. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, Y.; Guo, Z.; Zhao, Y.; Jin, Y.; An, L.; Wu, B.; Liu, Z.; Chen, X.; Chen, X.; Zhou, H.; et al. Genetic polymorphisms of genes in the lncRNA-p53 regulatory network are associated with concomitant toxicities and efficacy of chemoradiotherapy in patients with nasopharyngeal carcinoma. Sci. Rep. 2017, 7, 8320. [Google Scholar] [CrossRef] [PubMed]
  26. Gene. Bethesda (MD): National Library of Medicine (USA), National Center for Biotechnology Information; 2004. Available online: https://www.ncbi.nlm.nih.gov/gene/5581#summary (accessed on 27 February 2024).
  27. Pu, X.; Wang, L.; Chang, J.Y.; Hildebrandt, M.A.; Ye, Y.; Lu, C.; Skinner, H.D.; Niu, N.; Jenkins, G.D.; Komaki, R.; et al. Inflammation-related genetic variants predict toxicity following definitive radiotherapy for lung cancer. Clin. Pharmacol. Ther. 2014, 96, 609–615. [Google Scholar] [CrossRef]
  28. Greenhough, L.A.; Liang, C.C.; Belan, O.; Kunzelmann, S.; Maslen, S.; Rodrigo-Brenni, M.C.; Anand, R.; Skehel, M.; Boulton, S.J.; West, S.C. Structure and function of the RAD51B-RAD51C-RAD51D-XRCC2 tumor suppressor. Nature 2023, 619, 650–657. [Google Scholar] [CrossRef] [PubMed]
  29. Li, F.; Zhang, Y.; Shi, Y.; Liu, S. Comprehensive Analysis of Prognostic and Immune Infiltrates for RAD51 in Human Breast Cancer. Crit. Rev. Eukaryot. Gene Expr. 2021, 31, 71–79. [Google Scholar] [CrossRef]
  30. Xu, H.; Xiong, C.; Chen, Y.; Zhang, C.; Bai, D. Identification of RAD51 as a prognostic biomarker correlated with immune infiltration in hepatocellular carcinoma. Bioengineered 2021, 12, 2664–2675. [Google Scholar] [CrossRef]
  31. Data for This Paper Were Retrieved from the Alliance of Genome Resources. Version 6.0.0. Available online: https://www.alliancegenome.org (accessed on 20 February 2024).
  32. Fachal, L.; Gómez-Caamaño, A.; Barnett, G.C.; Peleteiro, P.; Carballo, A.M.; Calvo-Crespo, P.; Kerns, S.L.; Sánchez-García, M.; Lobato-Busto, R.; Dorling, L.; et al. A three-stage genome-wide association study identifies a susceptibility locus for late radiotherapy toxicity at 2q24.1. Nat. Genet. 2014, 46, 891–894. [Google Scholar] [CrossRef]
  33. Aula, H.; Skyttä, T.; Tuohinen, S.; Luukkaala, T.; Hämäläinen, M.; Virtanen, V.; Raatikainen, P.; Moilanen, E.; Kellokumpu-Lehtinen, P.L. Decreases in TGF-β1 and PDGF levels are associated with echocardiographic changes during adjuvant radiotherapy for breast cancer. Radiat. Oncol. 2018, 13, 201. [Google Scholar] [CrossRef]
  34. Delgado, B.D.; Enguix-Riego, M.V.; de Bobadilla, J.C.F.; Rivera, D.H.; Gómez, J.M.N.G.; Praena-Fernández, J.M.; Del Campo, E.R.; Gordillo, M.J.O.; Fernandez, M.D.C.F.; Guerra, J.L.L. Association of single nucleotide polymorphisms at HSPB1 rs7459185 and TGFB1 rs11466353 with radiation esophagitis in lung cancer. Radiother. Oncol. 2019, 135, 161–169. [Google Scholar] [CrossRef]
  35. Talibova, G.; Bilmez, Y.; Ozturk, S. Increased double-strand breaks in aged mouse male germ cells may result from altered expression of genes essential for homologous recombination or nonhomologous end-joining repair. Histochem. Cell Biol. 2023, 159, 127–147. [Google Scholar] [CrossRef]
  36. Wen, Y.; Dai, G.; Wang, L.; Fu, K.; Zuo, S. Silencing of XRCC4 increases the radiosensitivity of triple-negative breast cancer cells. Biosci. Rep. 2019, 39, BSR20180893. [Google Scholar] [CrossRef]
  37. Rodrigues, J.C.G.; de Souza, T.P.; Pastana, L.F.; dos Santos, A.M.R.; Fernandes, M.R.; Pinto, P.; Wanderley, A.V.; de Souza, S.J.; Kroll, J.E.; Pereira, A.L.; et al. Identification of NUDT15 gene variant in Amazonian Amerindians and admixed individuals from northern Brazil. PLoS ONE 2020, 15, e0231651. [Google Scholar] [CrossRef]
  38. Ribeiro-dos-Santos, A.M.; Vidal, A.F.; Vinasco-Sandoval, T.; Guerreiro, J.; Santos, S.; Ribeiro-dos-Santos, Â.; De Souza, S.J. Exome Sequencing of Native Populations from the Amazon Reveals Patterns on the Peopling of South America. Front. Genet. 2020, 11, 548507. [Google Scholar] [CrossRef] [PubMed]
  39. Wu, Y.; Song, Y.; Wang, R.; Wang, T. Molecular mechanisms of tumor resistance to radiotherapy. Mol. Cancer 2023, 22, 96. [Google Scholar] [CrossRef]
  40. Lee, E.; Eum, S.Y.; Slifer, S.H.; Martin, E.R.; Takita, C.; Wright, J.L.; Hines, R.B.; Hu, J.J. Association between polymorphisms in DNA damage repair genes and radiation therapy–induced early adverse skin reactions in a breast cancer population: A polygenic risk score approach. Int. J. Radiat. Oncol. Biol. Phys. 2020, 106, 948–957. [Google Scholar] [CrossRef] [PubMed]
  41. Jandu, H.K.; Veal, C.D.; Fachal, L.; Luccarini, C.; Aguado-Barrera, M.E.; Altabas, M.; Azria, D.; Baten, A.; Bourgier, C.; Bultijnck, R.; et al. Genome-wide association study of treatment-related toxicity two years following radiotherapy for breast cancer. Radiother. Oncol. 2023, 187, 109806. [Google Scholar] [CrossRef] [PubMed]
  42. Ko, D.C.; Urban, T.J. Understanding Human Variation in Susceptibility to Infectious Diseases through Clinical and Cellular GWAS. PLoS Pathog. 2013, 9, e1003424. [Google Scholar] [CrossRef]
  43. Nikpay, M.; Goel, A.; Won, H.H.; Hall, L.M.; Willenborg, C.; Kanoni, S.; Saleheen, D.; Kyriakou, T.; Nelson, C.P.; Hopewell, J.C.; et al. A comprehensive genome-wide association meta-analysis of 1000 coronary artery disease genomes. Nat. Genet. 2015, 47, 1121–1130. [Google Scholar]
  44. El Naqa, I.; Pandey, G.; Aerts, H.; Chien, J.T.; Andreassen, C.N.; Niemierko, A.; Ten Haken, R.K. Radiotherapy outcome models in the era of radiomics and radiogenomics: Uncertainties and validation. Int. J. Radiat. Oncol. Biol. Phys. 2018, 102, 1070–1073. [Google Scholar] [CrossRef]
  45. Kerns, S.L.; Fachal, L.; Dorling, L.; Barnett, G.C.; Baran, A.; Peterson, D.R.; Hollenberg, M.; Hao, K.; Narzo, A.D.; Ahsen, M.E.; et al. Radiogenomics Consortium Genome-Wide Association Study Meta-Analysis of Late Toxicity After Prostate Cancer Radiotherapy. J. Natl. Cancer Inst. 2020, 112, 179–190. [Google Scholar] [CrossRef]
  46. Ranjha, L.; Howard, S.M.; Cejka, P. Main steps in DNA double-strand break repair: An introduction to homologous recombination and related processes. Chromosoma 2018, 127, 187–214. [Google Scholar] [CrossRef] [PubMed]
  47. Zhang, J.; Wu, Q.; Zhu, L.; Xie, S.; Tu, L.; Yang, Y.; Wu, K.; Zhao, Y.; Wang, Y.; Xu, Y.; et al. SERPINE2/PN-1 regulates the DNA damage response and radioresistance by activating ATM in lung cancer. Cancer Lett. 2022, 524, 268–283. [Google Scholar] [CrossRef] [PubMed]
  48. Fernandes, M.R.; Rodrigues, J.C.G.; Maroñas, O.; Latorre-Pellicer, A.; Cruz, R.; Guerreiro, J.F.; Burbano, R.M.R.; Assumpção, P.P.D.; Ribeiro-dos-Santos, A.; Santos, S.E.B.D.; et al. Genetic Diversity of Drug-Related Genes in Native Americans of the Brazilian Amazon. Pharmacogenom. Pers. Med. 2021, 14, 117–133. [Google Scholar] [CrossRef] [PubMed]
  49. Chorev, M.; Carmel, L. The function of introns. Front. Genet. 2012, 3, 55. [Google Scholar] [CrossRef] [PubMed]
  50. Zhang, P.; Chaldebas, M.; Ogishi, M.; Al Qureshah, F.; Ponsin, K.; Feng, Y.; Rinchai, D.; Milisavljevic, B.; Han, J.E.; Moncada-Vélez, M.; et al. Genome-wide detection of human intronic AG-gain variants located between splicing branch points and canonical splice acceptor sites. Proc. Natl. Acad. Sci. USA 2023, 120, e2314225120. [Google Scholar] [CrossRef]
  51. Mumbrekar, K.D.; Goutham, H.V.; Vadhiraja, B.M.; Bola Sadashiva, S.R. Polymorphisms in double-strand break repair related genes influence radiosensitivity phenotype in lymphocytes from healthy individuals. DNA Repair 2016, 40, 27–34. [Google Scholar] [CrossRef] [PubMed]
  52. Du, Z.; Zhang, W.; Zhou, Y.; Yu, D.; Chen, X.; Chang, J.; Qiao, Y.; Zhang, M.; Huang, Y.; Wu, C.; et al. Associations of ATM Polymorphisms with Survival in Advanced Esophageal Squamous Cell Carcinoma Patients Receiving Radiation Therapy. Int. J. Radiat. Oncol. Biol. Phys. 2015, 93, 181–189. [Google Scholar] [CrossRef]
  53. Zhang, Y.; Liu, Z.; Wang, M.; Tian, H.; Su, K.; Cui, J.; Dong, L.; Han, F. Single Nucleotide Polymorphism rs1801516 in Ataxia Telangiectasia-Mutated Gene Predicts Late Fibrosis in Cancer Patients after Radiotherapy: A PRISMA-Compliant Systematic Review and Meta-Analysis. Medicine 2016, 95, e3267. [Google Scholar] [CrossRef]
  54. McDuff, S.G.R.; Bellon, J.R.; Shannon, K.M.; Gadd, M.A.; Dunn, S.; Rosenstein, B.S.; Ho, A.Y. ATM Variants in Breast Cancer: Implications for Breast Radiation Therapy Treatment Recommendations. Int. J. Radiat. Oncol. Biol. Phys. 2021, 110, 1373–1382. [Google Scholar] [CrossRef] [PubMed]
  55. Andreassen, C.N.; Rosenstein, B.S.; Kerns, S.L.; Ostrer, H.; De Ruysscher, D.; Cesaretti, J.A.; Barnett, G.C.; Dunning, A.M.; Dorling, L.; West, C.M.; et al. International Radiogenomics Consortium. Individual patient data meta-analysis shows a significant association between the ATM rs1801516 SNP and toxicity after radiotherapy in 5456 breast and prostate cancer patients. Radiother. Oncol. 2016, 121, 431–439. [Google Scholar] [CrossRef] [PubMed]
  56. Tinhofer, I.; Stenzinger, A.; Eder, T.; Konschak, R.; Niehr, F.; Endris, V.; Distel, L.; Hautmann, M.G.; Mandic, R.; Stromberger, C.; et al. Targeted next-generation sequencing identifies molecular subgroups in squamous cell carcinoma of the head and neck with distinct outcomes after concurrent chemoradiation. Ann. Oncol. 2016, 27, 2262–2268. [Google Scholar] [CrossRef] [PubMed]
  57. Butkiewicz, D.; Gdowicz-Kłosok, A.; Krześniak, M.; Rutkowski, T.; Krzywon, A.; Cortez, A.J.; Domińczyk, I.; Składowski, K. Association of Genetic Variants in ANGPT/TEK and VEGF/VEGFR with Progression and Survival in Head and Neck Squamous Cell Carcinoma Treated with Radiotherapy or Radiochemotherapy. Cancers 2020, 12, 1506. [Google Scholar] [CrossRef]
  58. Xiao, Y.; Yuan, X.; Qiu, H.; Li, Q. Single-nucleotide polymorphisms of TGFβ1 and ATM associated with radiation-induced pneumonitis: A prospective cohort study of thoracic cancer patients in China. Int. J. Clin. Exp. Med. 2015, 8, 16403–16413. [Google Scholar] [PubMed]
  59. Niu, Z.; Sun, P.; Liu, H.; Wei, P.; Wu, J.; Huang, Z.; Gross, N.D.; Shete, S.; Wei, Q.; Zafereo, M.E.; et al. Functional Genetic Variants in TGFβ1 and TGFβR1 in miRNA-Binding Sites Predict Outcomes in Patients with HPV-positive Oropharyngeal Squamous Cell Carcinoma. Clin. Cancer Res. 2023, 29, 3081–3091. [Google Scholar] [CrossRef] [PubMed]
  60. Pastana, L.F.; Silva, T.A.; Gellen, L.P.A.; Vieira, G.M.; de Assunção, L.A.; Leitão, L.P.C.; da Silva, N.M.; Coelho, R.D.C.C.; de Alcântara, A.L.; Vinagre, L.W.M.S.; et al. The Genomic Profile Associated with Risk of Severe Forms of COVID-19 in Amazonian Native American Populations. J. Pers. Med. 2022, 12, 554. [Google Scholar] [CrossRef]
  61. Aguiar, K.E.C.; Oliveira, I.D.S.; Cohen-Paes, A.D.N.; Coelho, R.D.C.C.; Vinagre, L.W.M.S.; Rodrigues, J.C.G.; Ribeiro-Dos-Santos, A.M.; De Souza, S.J.; Ribeiro-Dos-Santos, Â.; Guerreiro, J.F.; et al. Molecular Profile of Variants in CDH1, TP53, PSCA, PRKAA1 and TTN Genes Related to Gastric Cancer Susceptibility in Amazonian Indigenous Population. J. Pers. Med. 2023, 13, 1364. [Google Scholar] [CrossRef]
  62. Coelho, R.d.C.C.; Martins, C.L.e.L.P.; Pastana, L.F.; Rodrigues, J.C.G.; Aguiar, K.E.C.; Cohen-Paes, A.d.N.; Gellen, L.P.A.; Moraes, F.C.A.d.; Calderaro, M.C.L.; de Assunção, L.A.; et al. Molecular Profile of Variants Potentially Associated with Severe Forms of COVID-19 in Amazonian Indigenous Populations. Viruses 2024, 16, 359. [Google Scholar] [CrossRef]
  63. Cohen-Paes, A.; de Alcântara, A.L.; de Souza Menezes, E.; Moreira, F.C.; Fernandes, M.R.; Guerreiro, J.F.; Ribeiro-Dos-Santos, Â.; Dos Santos, S.E.B.; Santos, N.P.C.D. Characterization of DNA Polymerase Genes in Amazonian Amerindian Populations. Genes 2022, 14, 53. [Google Scholar] [CrossRef]
  64. O’Sullivan, N.J.; Kelly, M.E. Radioomics and Radiogenomics in Pelvic Oncology: Current Applications and Future Directions. Curr. Oncol. 2023, 30, 4936–4945. [Google Scholar] [CrossRef] [PubMed]
  65. Badic, B.; Tixier, F.; Cheze Le Rest, C.; Hatt, M.; Visvikis, D. Radiogenomics in Colorectal Cancer. Cancers 2021, 13, 973. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The relative contributions of variants discriminated according to high, modifier, or moderate impact in the AREG, ATM, KDR, MEG3, PRKCE, RAD51, TGFB1, TANC1, and XRCC4 genes.
Figure 1. The relative contributions of variants discriminated according to high, modifier, or moderate impact in the AREG, ATM, KDR, MEG3, PRKCE, RAD51, TGFB1, TANC1, and XRCC4 genes.
Jpm 14 00484 g001
Figure 2. Differences in allele frequencies of the variants studied in the continental populations and indigenous population, plotted on an MDS.
Figure 2. Differences in allele frequencies of the variants studied in the continental populations and indigenous population, plotted on an MDS.
Jpm 14 00484 g002
Table 1. The functions of the AREG, ATM, KDR, MEG3, PRKCE, RAD51, TGFB1, TANC1, and XRCC4 genes.
Table 1. The functions of the AREG, ATM, KDR, MEG3, PRKCE, RAD51, TGFB1, TANC1, and XRCC4 genes.
GeneDescription
AREGRelated to epidermal growth factor (EGF) and transforming growth factor alpha (TGF-alpha). The protein interacts with the EGF/TGF-alpha receptor, promotes the growth of normal epithelial cells, and inhibits the growth of carcinoma cell lines. Studies indicate that this gene is involved in tumor escape and radiation resistance and is significantly overexpressed after radiotherapy [17].
ATMFavors the interruption of the cell cycle, senescence, and apoptosis and acts on the cellular response to DNA double-strand breaks, occasionally caused by radiotherapy. Therefore, inhibiting ATM seems to increase the effectiveness of radiotherapy [18]. SNPs (rs620815 and rs11212570) in the gene are associated with an increased risk of gastrointestinal toxicity and odynophagia in individuals treated with radiotherapy [19].
CDK1Cyclin-dependent kinase (CDK1), a cell cycle regulator, is intimately involved in several cellular events vital for cell survival, such as the regulation of gene expression. CDK1 is responsible for regulating the transition between the G2 phase and mitosis [20]. CDK1 is inactivated when DNA damage occurs due to ionizing irradiation to arrest the cell cycle at the G2 checkpoint and facilitate the repair of double-strand breaks, which may be beneficial in radiotherapy treatment [21].
KDRThe kinase insert domain receptor, also known as VEGFR (vascular endothelial growth factor), acts to promote angiogenesis in normal and pathological conditions, such as cancer [22]. An SNP (rs1870377) in the gene was associated with mucositis in patients treated with chemoradiotherapy [23].
MEG3A tumor suppressor that negatively regulates proteins such as p53 and STAT3, consequently affecting cell proliferation and metastasis and promoting apoptosis [24,25]. The rs1032552 variant of MEG3 has been associated with toxicities such as dermatitis and anemia in individuals treated with chemoradiotherapy [25].
PRKCEMembers of the PKC family phosphorylate a wide variety of protein targets and are involved in cellular signaling pathways such as neuron channel activation and apoptosis [26]. Furthermore, the literature indicates that the PRKCE gene is also associated with radiation toxicity in lung cancer [27].
RAD51Acts on DNA repair and, when there is interference in the binding to ATPase and single-stranded DNA, generates genomic instability [4,28]. Gene expression is more frequent in certain tumor tissues and implies a lower survival rate, tumor progression, immunosuppression, radioresistance, and worse prognosis [29,30]. It is associated with radiotoxicity in HER2-positive breast cancer [6].
TANC1Tetratricopeptide repeat, ankyrin repeat, and coiled-coil containing 1 is a protein-coding gene, considered a scaffold component related to the regulation post-synapse [31]. It is involved in the repair of damaged muscle cells, suggesting that its biological mechanism is associated with the greater development of radiation toxicity; this is due to its potential role in the regeneration of radio-induced damage in muscle tissue [32].
TGFB1Growth factor β 1 acts to control growth, proliferation, differentiation, and apoptosis. It appears to be a marker of radiosensitivity, where its reduction during radiotherapy treatment is associated with a positive response [33]. It plays an important role in inflammation and cell proliferation, which may be related to radiation-induced fibrosis [34].
XRCC4Acts in the repair of DNA double-strand breaks [35]. Studies have indicated that reducing XRCC4 expression delays the repair of DNA damage. Furthermore, the protein is a radiosensitivity marker and a promising study object for radiogenomics [36].
Table 2. Descriptions of the variants in the AREG, ATM, KDR, MEG3, PRKCE, RAD51, TGFB1, TANC1, and XRCC4 genes according to their impact, region, variant type, and change in nucleotide.
Table 2. Descriptions of the variants in the AREG, ATM, KDR, MEG3, PRKCE, RAD51, TGFB1, TANC1, and XRCC4 genes according to their impact, region, variant type, and change in nucleotide.
ImpactGenedbSNPRegionVar TypeChange in NucleotideFrequencies
INDGAFRAMREASEURSAS
HighKDRrs41452948Protein_Structural_Interaction_
Locus
SNVG > A0.0517000.00400
HighXRCC4rs1805377Splice_Site_Acceptor+Intron SNVG > A0.66670.476 0.3330.7040.1400.172
ModifierAREGrs368667736IntronSNVC > T0.10000.2260.2420.5220.1450.242
ModifierATMrs672655IntronSNVA > G0.12500.2160.6400.4200.6170.626
ModifierATMrs58978479IntronINDELAT > A0.24590.0050.025 0.030-0.030
ModifierATMrs2066734IntronINDELTAA > T0.62300.1230.5320.3840.4510.549
ModifierATMrs664143IntronSNVA > G0.21050.7000.6830.4470.6270.679
ModifierATMrs3218681IntronINDELAA > AAA0.00000.4300.6430.4080.6250.667
ModifierMEG3rs142677044IntragenicSNVG > C0.189700.0450.01400
ModifierMEG3rs139003317IntronINDELCCCT > C0.05410.6060.173 0.0010.210 0.182
ModifierMEG3rs56363527IntronSNVC > T0.05710.0610.3210.2370.3350.276
ModifierMEG3rs7160821IntragenicSNVG > A0.19050.2160.4610.3250.3440.333
ModifierKDRrs2305946IntronSNVC > T0.0135 0.1040.1400.454 0.2380.152
ModifierKDRrs2219471IntronSNVT > C0.04690.135 0.1470.459 0.2400.121
ModifierKDRrs3816584IntronSNVA > G0.0135 0.1030.1400.454 0.2380.152
ModifierKDRrs7655964IntronSNVA > C0.64060.1080.4960.313 0.3630.306
ModifierKDRrs17085310IntronSNVG > A0.23440.001 0.1120.0740.0050.055
ModifierKDRrs3214870IntronINDELGG > GGG0.06350.1970.1460.453 0.2390.152
ModifierKDRrs7692791IntronSNVC > T0.0857 0.3900.4640.6300.5450.721
ModifierRAD51rs45455000IntronSNVT > G0.00000.1160.0810.1330.0750.126
ModifierPRKCErs60465117IntronSNVA > C0.00000.0660.0170.085 0.005 0.084
ModifierPRKCErs4953294IntronSNVG > A0.0270 0.067 0.1830.0870.2670.337
ModifierPRKCErs1987070IntronSNVC > A0.01390.152 0.1410.2190.2730.197
ModifierPRKCErs201731045IntronSNVT > C0.0806 0.0030.001000
ModifierPRKCErs2249505IntronSNVC > T0.6471 0.0140.3980.584 0.2160.195
ModifierPRKCErs10495929IntronSNVG > A0.0147 0.186 0.1330.0590.1490.215
ModifierTANC1rs34344829IntronSNVA > G0.0270 0.0230.2870.048 0.5130.186
ModifierTANC1rs146371641IntronINDELCCC > CCCC0.04050.185 0.295 0.0480.5140.186
ModifierRAD51rs45457497IntronSNVT > G0.50000.2750.3930.6440.1530.342
ModifierRAD51rs200723181IntronINDELT > TCT0.0000 0.3010.3410.569 0.1500.268
ModerateATMrs1801516Non_Synonymous_CodingSNVG > A0.00000.0080.0970.0160.1620.080
ModerateKDRrs1870377Non_Synonymous_CodingSNVT > A0.04690.0900.1310.4650.2350.149
ModerateTANC1rs34588551Non_Synonymous_CodingSNVC > T0.03230.0300.249 0.0030.3470.122
ModerateTANC1rs4664277Non_Synonymous_CodingSNVA > G0.94530.4610.610 0.6180.3730.535
ModerateTGFB1rs1800470Non_Synonymous_CodingSNVG > A0.21430.5860.4940.4450.6180.554
(-) No annotation. INDG: Indigenous Amazonian population, AFR: African population, AMR: American population, EAS: East Asian population, EUR: European population, SAS: South Asian population.
Table 3. Description of new variants found in the indigenous population from the Brazilian Amazon in genes relevant to radiotherapy.
Table 3. Description of new variants found in the indigenous population from the Brazilian Amazon in genes relevant to radiotherapy.
GeneChromosomePositionVar TypeRegion DetailedReferenceVariantImpactProtein ChangeVariant Allele
Frequency
PRKCEchr245652078SNV5UTRGCMODIFIERc.-23G > C0.0833
PRKCEchr245652096INDEL5UTRCCCCCCCAGGGTMODIFIERc.-5_-4insCCCCCAGGGT0.0833
PRKCEchr245652087SNV5UTRTCMODIFIERc.-14T > C0.0833
PRKCEchr245652092SNV5UTRGCMODIFIERc.-9G > C0.1667
Table 4. Comparison between the allelic frequencies in Indigenous Amazonian (INDG) and continental populations (AFR, AMR, EAS, EUR, and SAS) described in the 1000 Genomes and the Exome Aggregation Consortium (ExAC).
Table 4. Comparison between the allelic frequencies in Indigenous Amazonian (INDG) and continental populations (AFR, AMR, EAS, EUR, and SAS) described in the 1000 Genomes and the Exome Aggregation Consortium (ExAC).
GENEdbSPNINDG × AFR *INDG × AMR *INDG × EAS *INDG × EUR *INDG × SAS *
PRKCErs22495052.17 × 10−292.40 × 10−4-2.58 × 10−111.18 × 10−12
KDRrs23059461.06 × 10−25.43 × 10−41.23 × 10-171.56 × 10−72.88 × 10−4
TANC1rs343448296.22 × 10−96.81 × 10−20--5.17 × 10−5
KDRrs38165841.06 × 10−25.43 × 10−41.23 × 10−171.56 × 10−72.88 × 10−4
TANC1rs34588551 1.07 × 10−6-1.83 × 10−101.50 × 10−2
KDRrs22194714.08 × 10−21.72 × 10−22.51 × 10−142.82 × 10−5-
PRKCErs604651171.21 × 10−2-1.62 × 10−3-1.62 × 100−3
PRKCErs4953294-1.00 × 10−4-5.65 × 10−86.14 × 10−11
TANC1rs1463716415.63 × 10−47.86 × 10−8-6.18 × 10−185.63 × 10−4
MEG3rs1426770441.31 × 10−71.50 × 10−31.08 × 10−51.31 × 10−71.31 × 10−7
MEG3rs1390033171.90 × 10−219.17 × 10−32.28 × 10−27.49 × 10−45.70 × 10−3
KDRrs76559643.81 × 10−183.55 × 10-27.25 × 10−72.92 × 10−53.75 × 10−7
TANC1rs46642771.51 × 10−172.25 × 10−104.70 × 10−102.36 × 10−226.04 × 10−14
KDRrs170853101.31 × 10−92.04 × 10−28.51 × 10−42.07 × 10−81.49 × 10−4
PRKCErs19870702.88 × 10−45.43 × 10−46.57 × 10−74.13 × 10−95.43 × 10−6
PRKCErs2017310453.12 × 10−33.12 × 10−33.12 × 10−33.12 × 10−33.12 × 10−3
XRCC4rs18053776.73 × 10−37.42 × 10−7-4.44 × 10−176.88 × 10−15
KDRrs32148704.60 × 10−3-1.63 × 10−122.22 × 10−4-
MEG3rs56363527-1.44 × 10−71.49 × 10−43.89 × 10−87.28 × 10−6
KDRrs76927913.91 × 10−85.38 × 10−119.11 × 10−204.71 × 10−153.10 × 10−25
MEG3rs7160821-1.25 × 10−52.34 × 10−21.21 × 10−21.72 × 10−2
PRKCErs104959291.08 × 10−51.00 × 10−3-2.88 × 10−46.57 × 10−7
AREGrs3686677361.79 × 10−28.23 × 10−31.05 × 10−12-8.23 × 10−3
KDRrs1870377-4.08 × 10−21.05 × 10−145.17 × 10−51.72 × 10−2
KDRrs414529482.28 × 10−22.28 × 10−2-2.28 × 10−22.28 × 10−2
ATMrs672655-6.04 × 10−174.11 × 10−71.17 × 10−154.38 × 10−16
ATMrs1801516-8.12 × 10−4-1.21 × 10−63.12 × 10−3
ATMrs589784799.43 × 10−94.53 × 10−72.11 × 10−6-2.11 × 10−6
TGFB1rs18004705.18 × 10−91.07 × 10−52.45 × 10−42.48 × 10−108.16 × 10−8
RAD51rs454550001.02 × 10−43.12 × 10−32.39 × 10−53.12 × 10−34.99 × 10−5
ATMrs20667344.38 × 10−16-3.84 × 10−41.43 × 10−2-
ATMrs6641431.45 × 10−142.64 × 10−132.45 × 10−41.13 × 10−102.64 × 10−13
RAD51rs454574976.74 × 10−4-4.70 × 10−22.40 × 10−82.52 × 10−2
ATMrs32186814.41 × 10−195.56 × 10−327.45 × 10−186.51 × 10−01.17 × 10−33
RAD51rs2007231817.85 × 10−131.05 × 10−143.03 × 10−275.43 × 10−65.36 × 10−11
(-) No annotation. INDG: Indigenous Amazonian population, AFR: African population, AMR: American population, EAS: East Asian population, EUR: European population, SAS: South Asian population. * p-value ≤ 0.05.
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

de Lima, M.C.; de Castro, C.C.; Aguiar, K.E.C.; Monte, N.; Nunes, G.G.d.C.; Costa, A.C.A.d.; Rodrigues, J.C.G.; Guerreiro, J.F.; Ribeiro-dos-Santos, Â.; Assumpção, P.P.d.; et al. Molecular Profile of Important Genes for Radiogenomics in the Amazon Indigenous Population. J. Pers. Med. 2024, 14, 484. https://doi.org/10.3390/jpm14050484

AMA Style

de Lima MC, de Castro CC, Aguiar KEC, Monte N, Nunes GGdC, Costa ACAd, Rodrigues JCG, Guerreiro JF, Ribeiro-dos-Santos Â, Assumpção PPd, et al. Molecular Profile of Important Genes for Radiogenomics in the Amazon Indigenous Population. Journal of Personalized Medicine. 2024; 14(5):484. https://doi.org/10.3390/jpm14050484

Chicago/Turabian Style

de Lima, Milena Cardoso, Cinthia Costa de Castro, Kaio Evandro Cardoso Aguiar, Natasha Monte, Giovanna Gilioli da Costa Nunes, Ana Caroline Alves da Costa, Juliana Carla Gomes Rodrigues, João Farias Guerreiro, Ândrea Ribeiro-dos-Santos, Paulo Pimentel de Assumpção, and et al. 2024. "Molecular Profile of Important Genes for Radiogenomics in the Amazon Indigenous Population" Journal of Personalized Medicine 14, no. 5: 484. https://doi.org/10.3390/jpm14050484

APA Style

de Lima, M. C., de Castro, C. C., Aguiar, K. E. C., Monte, N., Nunes, G. G. d. C., Costa, A. C. A. d., Rodrigues, J. C. G., Guerreiro, J. F., Ribeiro-dos-Santos, Â., Assumpção, P. P. d., Burbano, R. M. R., Fernandes, M. R., dos Santos, S. E. B., & dos Santos, N. P. C. (2024). Molecular Profile of Important Genes for Radiogenomics in the Amazon Indigenous Population. Journal of Personalized Medicine, 14(5), 484. https://doi.org/10.3390/jpm14050484

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