- Research article
- Open Access
Genetically-regulated transcriptomics & copy number variation of proctitis points to altered mitochondrial and DNA repair mechanisms in individuals of European ancestry
BMC Cancer volume 20, Article number: 954 (2020)
Proctitis is an inflammation of the rectum and may be induced by radiation treatment for cancer. The genetic heritability of developing radiotoxicity and prior role of genetic variants as being associated with side-effects of radiotherapy necessitates further investigation for underlying molecular mechanisms. In this study, we investigated gene expression regulated by genetic variants, and copy number variation in prostate cancer survivors with radiotoxicity.
We investigated proctitis as a radiotoxic endpoint in prostate cancer patients who received radiotherapy (n = 222). We analyzed the copy number variation and genetically regulated gene expression profiles of whole-blood and prostate tissue associated with proctitis. The SNP and copy number data were genotyped on Affymetrix® Genome-wide Human SNP Array 6.0. Following QC measures, the genotypes were used to obtain gene expression by leveraging GTEx, a reference dataset for gene expression association based on genotype and RNA-seq information for prostate (n = 132) and whole-blood tissue (n = 369).
In prostate tissue, 62 genes were significantly associated with proctitis, and 98 genes in whole-blood tissue. Six genes - CABLES2, ATP6AP1L, IFIT5, ATRIP, TELO2, and PARD6G were common to both tissues. The copy number analysis identified seven regions associated with proctitis, one of which (ALG1L2) was also associated with proctitis based on transcriptomic profiles in the whole-blood tissue. The genes identified via transcriptomics and copy number variation association were further investigated for enriched pathways and gene ontology. Some of the enriched processes were DNA repair, mitochondrial apoptosis regulation, cell-to-cell signaling interaction processes for renal and urological system, and organismal injury.
We report gene expression changes based on genetic polymorphisms. Integrating gene-network information identified these genes to relate to canonical DNA repair genes and processes. This investigation highlights genes involved in DNA repair processes and mitochondrial malfunction possibly via inflammation. Therefore, it is suggested that larger studies will provide more power to infer the extent of underlying genetic contribution for an individual’s susceptibility to developing radiotoxicity.
Prostate cancer is one of the most prevalent diseases in older men, with 66 years being the average age at the time of cancer diagnosis . According to the SEER cancer statistics of 2019, approximately 3 million men have been previously diagnosed with prostate cancer and are still alive today. This feat can be credited to the advancement in cancer treatment which has contributed to the 5-year relative survival rate of 90% in prostate cancer survivors . Radiation therapy is one of the primary forms of treatment for prostate cancer, delivered as external beam radiotherapy (EBRT) or brachytherapy. While the dose and precision of radiation delivery to the tumor tissue has improved over the years, surrounding normal tissues get irradiated leading to clinical side effects . Proctitis is the inflammation of the rectum, which can result from receiving radiation therapy around the pelvic region such as during prostate cancer treatment . The inflammation of the rectum can either be acute or chronic. Acute proctitis appears within 3 months of receiving radiation therapy, and progression of rectal inflammation after 3 months of completing radiation therapy is identified as chronic proctitis . The development of radiation-induced chronic proctitis affects 5–20% of cancer survivors and is relatively more common  than acute proctitis which affects approximately 13% of the cancer population . As of 2016, the total population of cancer survivors in the US was estimated to be approximately 15 million, and by the year 2026 is expected to reach 20 million individuals . Given the prevalence of chronic proctitis affecting cancer individuals (5–20%), we can deduce that approximately 1–4 million cancer survivors experience proctitis from receiving radiotherapy. The goal of radiation therapy in treating cancer is to damage the DNA of cancer tissue by creating double-strand breaks (DSBs). While the cellular system is capable of repairing breaks in the DNA, strands with DSBs are difficult to restore leading to activation of apoptotic signals which ultimately kills cancer cells. Unfortunately, the normal tissue around the targeted region is also affected by DNA damage from radiation, and must rely on DNA repair mechanisms for rehabilitation of cellular functions . A recent twin-study has reported that certain SNPs and their transcriptomic influence are associated with individual radiation sensitivity and a heritability estimate of 66% . Therefore, it is vital to understand genetics underlying molecular mechanisms involved in adverse effects of radiotherapy and individual genetic variations that may induce radiation sensitivity. Genome-wide studies have been conducted to identify gene variants that may contribute towards developing radiotoxic side-effects. These studies have identified genetic variations involved in DNA repair pathways to be associated with overall radiotoxicity [3, 10]. However, the role of altered gene expression  from aggregated single nucleotide polymorphisms (SNPs) remains elusive in radiotoxicity phenotypes (e.g., proctitis). Regulatory variants are SNPs within coding regions which contribute towards tissue – specific gene expression alterations leading to variations in the phenotypic spectrum. Estimating the contribution of SNP aggregates to gene expression can be carried out using correlation weights derived from reference datasets which contain both SNP and RNA-seq information as modelled in PrediXcan . One such dataset is the GTEx project, an NIH funded initiative that stores genotype and RNA-seq data of 53 tissues from 620 donors (v7). The majority of the donors in the GTEx dataset are Caucasian, and more than 50% of the donors are over the age of 50 years. These characteristics make GTEx an excellent reference dataset to derive gene expression values from individual level SNP profiles of prostate cancer patients who have received radiation treatment.
Beyond SNPs, genetic discordance from gene dosage and structural effects can be attributed to copy number variation (CNV). CNVs are segments of duplicated DNA that are greater than 1 kb with differences in size between the two copies . In a clinical setting, testing of CNVs is relatively more common than other genetic tests  due to the majority of phenotypic changes being associated with these variations in segment size. CNVs associated with radiotoxicity phenotypes (e.g. proctitis) have not been investigated extensively and may prove to be significant contributors to phenotype risk.
We hypothesize that genotype-regulated gene expression profiles and variations in copy number will identify genetic alterations associated with a spectrum of DNA repair functions. Here, we investigate both CNV and tissue specific (prostate and whole-blood) transcriptomic profiles derived from individual-level SNPs that are associated with radiation induced proctitis in prostate cancer patients.
The overall methodology is visually summarized in Fig. S1.
Data access to study subjects
The Gene-PARE was approved by the Institutional Review Board of the Icahn School of Medicine at Mount Sinai and Florida Radiation Oncology Group (Kerns et al. [3, 14]). All patients provided informed consent under the parent study – Gene-PARE, at Icahn School of Medicine, Mount Sinai and Florida Radiation Oncology Group (Kerns et al. [3, 14]). We obtained access to anonymized individual level genotype data from Genetic Predictors of Adverse Radiotherapy Effects (Gene-PARE) (phs000772.v1.p1) via dbGaP’s authorized application, under the approval of North Texas Regional IRB protocol 2016–090. The study described here was performed under the North Texas Regional IRB (formerly the University of North Texas Health Science Center IRB), and was given “EXEMPT” status based on the criteria that our study involved data available from public repository, i.e. dbGaP and does not require approval of receiving informed consent. We analyzed prostate cancer individuals from the discovery set (N = 367), which were genotyped for 934,940 SNPs on Affymetrix® Genome-wide Human SNP Array 6.0. The dataset contains phenotypic information on prostate cancer patients who have received radiation treatment either via EBRT or brachytherapy. Out of three radiotoxicity phenotypes – erectile dysfunction, proctitis and urinary morbidity (IPSS/AUASS) – we focused our investigation on proctitis because (1) it is also prevalent in other pelvic region cancers  and (2) the data for proctitis was complete for all individuals.
We extracted SNP data from the intensity (*.CEL) files using Affymetrix® Genotyping Console using the BIRDSEED v2 algorithm and genotype call rate of 95% and default settings from the array, leaving 905,280 markers and total of 355 individuals. The files were then exported to plink  format to perform QC measures as suggested by Anderson et al. . At the individual-level filtering, we removed 5 individuals for either failing heterozygosity or having greater than 9% missing genotypes. At the IBD filter of 0.1875, we removed 29 individuals; further, 120 individuals were removed who were not of European ancestry or failed to cluster with the main patient population based on principal component (PC) analysis of PC1 and PC2 of the genetic data. After SNP-level filtering on SNP missingness, minimum allele frequency and Hardy-Weinberg equilibrium, we had 746,684 SNPs and 222 individuals. The final cohort characteristics after QC were analyzed using Fischer’s exact test for categorical variables and Student’s t-test for continuous variables (Table 1).
Genetically-regulated gene expression and GSEA
Tissue specific gene expression prediction using each individual’s genotype profile was performed using PrediXcan . Integrating expression quantitative trait loci (eQTLs) with genotype increases power to detect gene-based associations with phenotype of interest. The PrediXcan approach determines gene expression attributed to genetic variants, the aggregated effect of local (cis) genetic variants to the gene expression is tested for association with the phenotype. The PrediXcan model, derives weights of SNPs and tissue specific gene expression using lasso regression from GTEx  as reference dataset, including variance and covariance patterns based on linkage disequilibrium of cis-SNPs. We obtained model weights for prostate and whole-blood tissues for the GTEx-v7 reference panel, accessible at http://predictdb.org/; filename- GTEx-V7_HapMap-2017-11-29.tar.gz. In the GTEx (v7), there are 132 prostate tissue donors and 369 donors for whole-blood tissue. PrediXcan implements gene expression value prediction, followed by gene-based association. The z-scores identify the direction of expression  for each gene and their corresponding p-values for association testing. For prostate tissue, a total of 3113 gene and 5954 genes for whole-blood tissue were imputed for their tissue-specific gene expression based on genotype variation. Following association tests, significant genes (identified as p-value < 0.05) were investigated further by constructing tissue-specific protein-protein interaction (PPI) networks between query (significant genes) and interacting genes using the DifferentialNet  database and NetworkAnalyst3.0 . The network was filtered on betweenness centrality of 4.0 in order to reduce isolated neighboring nodes (each gene is a node). All the genes in the network were subsequently analyzed for gene set enrichment using clusterProfiler  for gene ontology and visualized in GOPlot .
Replication of transcriptomic findings with second cohort
We identified the study by Oorschot et al.  as an external cohort for replication of our transcriptomic findings. Their study recruited 200 patients who received EBRT for prostate cancer. Prior to receiving treatment, whole blood was drawn, and lymphocytes were cultured, irradiating half of the cells with 2Gy gamma rays and other half left untreated followed by genome-wide microarray analysis. The data were deposited in Gene Expression Omnibus - GSE85570 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE85570). Further details of the study have been described elsewhere . We used BART  and NCBI’s GEO2R to analyze differentially expressed genes between untreated and treated with 2 Gy of radiation groups. We tested the significance of genes overlapping our results from prostate and whole-blood tissue to those identified in the replication dataset using GeneOverlap (https://github.com/shenlab-sinai/geneoverlap) which applies Fisher’s exact test testing significant genes for whole blood tissue (n = 98), prostate tissue (n = 62) and replication cohort (n = 2973) against genome (n = 21,196) .
CNV association and GSEA
The intensity (*.CEL) files of 222 individuals from the above QC protocol were extracted for copy number analysis using Affymetrix® Genotyping Console. Copy number segments were filtered to regions (minimum genomic size of 2 kbps) with 10 marker per segment . The copy number data was exported as tab-delimited file for copy number association in CNVRuler . CNV regions were significantly associated at FDR p-value < 0.05. The significant CNV regions were visualized using Phenogram  then mapped to genes using UCSC browser for GRCh37/hg19 assembly (https://genome.ucsc.edu/). The genes within CNV regions were analyzed for functional and diseases processes and visualized using Ingenuity Pathway Analysis (QIAGEN Inc., https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis).
Genes identified in prostate tissue
In the association analysis between prostate cancer individuals who developed proctitis (cases) and who didn’t develop proctitis (controls), we found a total of 62 differentially expressed genes to be significantly associated. Based on z-score direction, 28 genes were downregulated, and 34 genes were upregulated in the prostate tissue (Table S1). We mapped the genes to tissue-specific protein-protein interaction (PPI) network to understand combined functional relationship of differentially expressed genes. Integrating PPI information with identified genes, we found an additional 22 genes which interact with 20 genes associated with proctitis based on genetically regulated expression of prostate tissue. Analyzing all the genes in the network for enriched gene ontology of biological processes, molecular functions and cellular components (Fig. 1) highlighted protein deubiquitination (ARRB2, TP53, SHMT2, BRCA1, ESR1, NEDD8, MYC), wnt signaling (MOV10, ARRB2, LRRK2, TNIK, ESR1, APP, CUL3), regulation of apoptosis signaling (ARRB2, TP53, LRRK2, BRCA1, YWHAZ, PTTG1IP), response to radiation (CIRBP, TP53, BRCA1, APP, MYC) and mitochondrial organization & apoptotic mitochondrial changes (ARRB2, TP53, LRRK2, YWHAZ) (Table S2). Specifically, gene ontology categories that included genes identified via SNP-based gene expression were ‘response to radiation’ and ‘regulation of apoptotic signaling pathway’.
Genes identified in whole blood tissue
We found a total of 98 genes to be associated with proctitis in whole blood tissue. 49 genes were upregulated, and 49 genes were downregulated (Table S3). Integrating PPI network information with the significant genes, found an additional 51 intermediate genes that interact with identified genes based on genetically regulated gene expression from whole blood. Investigating all the genes in the network for over-represented biological processes, highlighted DNA repair processes (Fig. 2) such as DNA replication (TERF2, EGFR, CDC7, BRCA1, ATRIP, RBBP8, SLX4, ORC1, ORC6, RAD50, CDK2, MCM2, DTL, RPA1, RPA2, RPA3), DNA integrity checkpoint (FBXO6, BRCA1, CDC5L, ATRIP, MDM2, ORC1, FZR1, CDK2, TP53, DTL, RPA2), nucleotide excision repair (COPS6, RBBP8, DDB1, UBC, SLX4, TP53, RPA1, RPA2, RPA3), recombinational repair (CDC7, BRCA1, RBBP8, SLX4, RAD50, RPA1, RPA2, RPA3), detection of DNA damage and response (DDB1,UBC,DTL,RPA1,RPA2,RPA3) and telomeric maintenance (TERF2, CCT5, SLX4, RAD50, TELO2, RPA1, RPA2, RPA3) (Table S4). Thus, integrating PPI information allows for the inference of gene-network as the gene expression changes attributed to genetic variants is limited in identifying all gene members of DNA damage and its related networks.
Replication dataset for transcriptomic findings
In order to replicate our genetically-regulated gene expression genes associated with radiation toxicity-proctitis, we identified the study of van Oorschot et al., whose data were deposited in gene expression omnibus (GEO) . Their study isolated and cultured lymphocytes from individuals who received radiation treatment for their prostate cancer and were assessed for radiotoxicity for a period of 2 years. In their study design they irradiated lymphocytes (collected prior to radiation treatment) with untreated and 2Gy of gamma ray, followed by microarray analysis. We analyzed differential expression of genes between the two irradiated groups. A total of 27 genes overlapped between significant genes identified in replication cohort and significant genes in whole blood tissue (overlapping p-value: 0.00028), and 5 genes overlapped with significant genes in prostate tissue. We further investigated the 27 genes that significantly overlapped with the external cohort for direction of gene expression. We found 14 genes to have concordant direction of expression between both results whereas the remaining 13 genes were discordant. The GSEA analysis of these 14 genes was enriched for Cilium Assembly, DNA repair, Organelle biogenesis and maintenance, CHEK2 PCC Network, Response to Ionizing radiation, TP53 Targets and Damaged DNA binding (Fig. 3; Supplementary file – Table S5).
CNV association and GSEA of mapped genes
We found 7 CNV regions associated with proctitis on chromosomes 1, 3, 4, 11, 12 and 15 (Fig. 4; Table S7). We identified genes within CNV regions using the UCSC browser (hg19) (Table S8). Interestingly, out of the two regions on chromosome 11 that were significant, we observed a high number of TRIM family genes (chr11:89487937–89,909,274 bps). The mapped genes from copy number regions were investigated for gene interactions using biobase knowledge of Ingenuity Pathway Analysis®. The pathway with the highest number of query genes (Fig. 5) was further analyzed for enriched disease and functional categories (Table S9). Cell-to-cell signaling interaction processes for renal and urological system, connective tissue development and function, and organismal injury were significantly associated processes, and their functional categorization included synthesis, proliferation, apoptosis and transmembrane transport. It is interesting to note, that most of these processes were dominated by TRH and TRIM-family genes. Furthermore, we also observed that the ALG1L2 gene, which was one of the significantly downregulated genes in whole-blood tissue, was also mapped to significant CNV region on chr3:129690192–129,896,364 bps which observes both gain and loss of copy, referred to as mixed regions.
Genetic susceptibility towards developing radiotoxic phenotypes is an upcoming research interest of significant clinical impact to improve the quality of life of cancer survivors . Previously, GWAS studies have been conducted to identify genetic loci associated with overall toxicity, decreased urinary stream, and erectile dysfunction [10, 29] in prostate cancer individuals who received radiotherapy . While these findings have shed much-needed light on SNP loci associated with susceptibility towards radiotoxicity, cumulative effects of exonic SNPs on gene expression and other genetic alterations such as copy number differences have not been previously studied. Here, we integrated genotyping data to identify genetic risk associated with proctitis by (1) employing genetic variant–derived gene expression of both prostate and whole-blood tissue, and (2) identifying associated genomic CNVs. The transcriptomic analyses points to several novel genes that play role in DNA-repair processes. In addition, we identified variable copy number regions that had multiple members of TRIM-family genes to be associated with proctitis. Along with novel genes identified through the analysis, the incorporation of PPI map reveals convergence of the implicated gene sets on known DNA-repair, mitochondrial, and telomeric regulation processes, highlighting their involvement with radiotoxic phenotypes (e.g. proctitis).
Six genes from both prostate and whole-blood tissue were associated with proctitis. CABLES2, ATP6AP1L and IFIT5 were under expressed, and ATRIP and TELO2 were upregulated in both tissues, however PARD6G was over expressed in prostate tissue and under expressed in whole blood tissue. CABLES2 (Cdk5 And Abl Enzyme Substrate 2), which is involved in regulation of the cell cycle, was also reported to be under expressed in lymphocytes of occupational workers who were exposed to ionizing radiation . ATP6AP1L (ATPase H+ Transporting Accessory Protein 1 Like) is critical for proton transportation and ATP synthase activity in the mitochondria; it has a paralog, ATP6AP1, which is involved in secretory granules and regulating neuroendocrine responses . IFIT5 (IFN-induced protein with tetratricopeptide repeats) has been reported to act as an enhancer in immune responses, with partial containment in mitochondria . A recent study has reported that elevated IFIT5 gene expression was correlated with interferon-γ levels in prostate cancer individuals after radiation, and demonstrated that IFN-γ stimulated epithelial-to-mesenchymal transition through the activation of JAK-STAT pathway . ATRIP [TREX1] is an ATR interacting protein that is a DNA exonuclease  that is known to initiate DNA repair pathways . ATRIP has also been reported to provide telomere protection by recruiting ATM – a key player in regulating cellular damage responses— to telomeric and DNA damage sites unaided by ATR kinase activity . Additionally, ATR responds to UV damage via the downregulation of Pin1 demonstrating anti-apoptotic activity in mitochondria . ATRIP [TREX1] has been shown to be upregulated in radiation-induced immunogenicity of tumor cells  by degrading cytosolic dsDNA and transferring cancer cells to dendritic cells under the stimulation of interferon-type1 . In addition to ATRIP, TELO2 [CLK2] also interacts with ATM to stimulate cell cycle arrest in response to radiation induced double strand breaks  via AKT activation . High TELO2 expression activity has been identified to be correlated with cell protection when exposed to high radiation doses ; conversely, TELO2 overexpression can trigger inflammation by influencing PIKKs (via mTORC1 binding ) while responding to DNA damage . PARD6G was found to have opposite direction of expression in prostate and whole-blood tissue, which could be attributed to tissue specific differences. Hypermethylation and downregulation of PARD6G was concluded to be involved in DNA repair mechanisms  in bisphenol A (BPA, a xenoestrogen) exposed human-derived breast cancer epithelial cells.
Assessment of copy number variation associated with radiation toxicity phenotypes can help identify genetic alterations that may lead to functional changes in gene expression, and thus, phenotypic variation. While CNVs have recently gathered interest in cancer diagnosis and treatment, their importance in studying radiation toxicity phenotypes remains understudied. Here, analysis of copy number data identified 7 regions to be associated with proctitis. The ALG1L2 mapped to the significant CNV region on chr3:129690192–129,896,364 (along with TRH and FAM86HP) was also found to be associated with proctitis in the transcriptomic analysis of whole-blood tissue, while the other two genes (TRH & FAM86HP) were not predicted in either of the two tissues. Based on gene ontology annotations, ALG1L2 functions as a mannosyl transferase in protein glycosylation . Another gene within chromosome 3 CNV region is TRH (Thyrotropin releasing hormone); its function includes carbohydrate and amino acid metabolism, and it has been involved in endocrine system disorders, metabolic disease, and organismal injury (Table S7). TRH has been shown to mobilize calcium from endoplasmic reticulum and mitochondria , and plays role in mitochondrial endoxidation via mitochondrial complex I and IV enzyme activity in skin samples , which is aligned with TRH’s known involvement in hair and skin development (as indicated in our IPA results, Table S7). Variants in FAM86HP (Family With Sequence Similarity 86 Member H, Pseudogene) have been reported to be associated with BMI-adjusted waist-hip ratio . In the pathway analysis (Fig. 5), we observed that FAM86HP has an indirect interaction with TGM2 (Transglutaminase 2), a stress-response gene  involved in mRNA metabolism  which has been reported to be upregulated during inflammation . It is interesting to note that TGM2 was over expressed in individuals receiving chemo & radiation therapy, suggesting its involvement in sensitivity to radiation [52, 53]. Both, TGM2 and NUPR1 (Fig. 5) have been reported to be implicated in inflammation driven primarily via the JAK/Stat and IL-17A signaling pathway . In the pathway we observe, NUPR1 has direct interaction with FAM72C/FAM27D and LOC100133315, both of which were identified from the copy number variation association. NUPR1 is known to repair double strand DNA breaks  and regulate cell cycle progression  from damage induced by gamma irradiation . The under expression of NUPR1 has been reported to result in increased ROS production thus creating a deficit in mitochondrial membrane potential. This alteration in OXPHOS activity has been associated with ER stress and triggers programmed necrosis in cancer cells . During angiogenesis in cancer cells, NUPR1 was reported to be upregulated in association with triiodothyronine thyroid hormone receptor , which is regulated with TRH - Thyrotropin releasing hormone . It is interesting to note that NUPR1 has been reported to play different roles before cancer development, during cancer progression and in response to cancer treatment.
Normal tissues receive varying amounts of radiation depending on their proximity to the tumor tissue, thus exhibiting a spectrum of toxicity effects. Genome-wide and molecular studies have shown that alterations in the DNA-damage response (DDR) from ATM (ataxia-telangiectasia mutation) mutations influence intraindividual variations to radiation toxicity . The findings in this study observe that genes such as ATRIP, NUPR1, TELO2, and TRH have multiple roles in damage detection and response mechanism. The genes and their PPI networks highlight their involvement in cell cycle arrest upon detection of DNA damage, affecting DNA replication and repair. During course of DDR, ROS from radiation seems to promote mitochondrial-induced apoptosis . Further, a constant but varying amount of inflammatory responses due to tissue injury (normal and tumor) appears to trigger ROS–induced mitochondrial oxidation , which exacerbates local inflammation in nearby tissues and propagates this DNA repair-inflammation stress cycle. Our findings report several novel genes that have been observed to be associated with known BRCA1-ATM-RAD50 damage response complex  which are activated in response to radiation, thus extending our understanding of these new players and their multifactorial roles associated with proctitis.
Our study has several limitations. The sample size is small, and our findings should be replicated and functionally validated in future studies using an animal model for targeted effects of radiotoxicity. Our study concentrated on prostate cancer survivors of Caucasian ethnicity, hence, to understand if similar or different genes are associated with proctitis developed during treatment of other cancers, it is imperative that our methods be applied to other races/ethnicities and other cancers. In the CNV regions, we mapped several noncoding regions, including lcRNA and snRNAs that would require fine mapping to further understand their involvement with proctitis. Unfortunately, there are only a handful of studies that investigate genetics underlying radiotoxicity phenotypes.
There are also many strengths to our study. Leveraging SNPs for reference transcriptomic data and copy number association, we identified several novel genes associated with proctitis – an inflammation of the rectum resulting from radiation therapy received for prostate cancer. The integration of tissue - specific PPI network aided in understanding the biological interactions between known and reported genes. We were able to replicate 27 of our genes to be significantly representative in an external cohort of prostate cancer diagnosis receiving radiation treatment. Analysis of copy number variation identified several genes in the reported regions and their pathways analysis highlighted two primary genes that have distinctive roles before cancer development and in response to cancer treatment.
In conclusion, this investigation highlights genes primarily involved in DNA repair processes and mitochondrial malfunction threaded via inflammation. The field of radiogenomics –work investigating the role of genetics in developing radiation toxicity—calls for investigation of genetic risk that can help inform dose management of radiation treatment and toxicity monitoring during treatment . We anticipate that understanding genetic data from both CNV and SNPs would contribute towards optimization of radiation treatment on an individual basis. Similar studies in the future would play strong role in early clinical interventions or periodic checkups for individuals who have high expression of DNA damage activity and alterations in copy number within specific regions.
Availability of data and materials
All results generated or analyzed during this study are included in this article (and its Supplementary Information files). Original genotype data can be accessed via dbGaP. The model weights for PrediXcan for GTEx genotype-gene expression associations can be accessed from predictdb.org, and external cohort can be accessed on Gene Expression Omnibus.
Single nucleotide polymorphism
Copy Number Variation
long coding RNA
Gene Set Enrichment Analysis
Ingenuity Pathway Analysis
External beam radiotherapy
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67:7–30. https://doi.org/10.3322/caac.21387.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69:7–34. https://doi.org/10.3322/caac.21551.
Fachal L, Gómez-Caamaño A, Barnett GC, Peleteiro P, Carballo AM, Calvo-Crespo P, 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–4. https://doi.org/10.1038/ng.3020.
Do NL, Nagle D, Poylin VY. Radiation proctitis: current strategies in management. Gastroenterol Res Pract. 2011;2011:917941. https://doi.org/10.1155/2011/917941.
Tabaja L, Sidani SM. Management of radiation proctitis. Dig Dis Sci. 2018;63:2180–8. https://doi.org/10.1007/s10620-018-5163-8.
Vanneste BGL, Van De Voorde L, de Ridder RJ, Van Limbergen EJ, Lambin P, van Lin EN. Chronic radiation proctitis: tricks to prevent and treat. Int J Color Dis. 2015;30:1293–303. https://doi.org/10.1007/s00384-015-2289-4.
Miller KD, Siegel RL, Lin CC, Mariotto AB, Kramer JL, Rowland JH, et al. Cancer treatment and survivorship statistics, 2016. CA Cancer J Clin. 2016;66:271–89. https://doi.org/10.3322/caac.21349.
West CM, Barnett GC. Genetics and genomics of radiotherapy toxicity: towards prediction. Genome Med. 2011;3:52. https://doi.org/10.1186/gm268.
Zyla J, Kabacik S, O’Brien G, Wakil S, Al-Harbi N, Kaprio J, et al. Combining CDKN1A gene expression and genome-wide SNPs in a twin cohort to gain insight into the heritability of individual radiosensitivity. Funct Integr Genomics. 2019;19:575–85. https://doi.org/10.1007/s10142-019-00658-3.
Barnett GC, Thompson D, Fachal L, Kerns S, Talbot C, Elliott RM, et al. A genome wide association study (GWAS) providing evidence of an association between common genetic variants and late radiotherapy toxicity. Radiother Oncol. 2014;111:178–85. https://doi.org/10.1016/j.radonc.2014.02.012.
Gamazon ER, Wheeler HE, Shah KP, Mozaffari SV, Aquino-Michaels K, Carroll RJ, et al. A gene-based association method for mapping traits using reference transcriptome data. Nat Genet. 2015;47:1091–8. https://doi.org/10.1038/ng.3367.
Valsesia A, Macé A, Jacquemont S, Beckmann JS, Kutalik Z. The growing importance of cnvs: new insights for detection and clinical interpretation. Front Genet. 2013;4:92. https://doi.org/10.3389/fgene.2013.00092.
Yang X, Song Z, Wu C, Wang W, Li G, Zhang W, et al. Constructing a database for the relations between CNV and human genetic diseases via systematic text mining. BMC Bioinformatics. 2018;19(Suppl 19):528. https://doi.org/10.1186/s12859-018-2526-2.
Kerns SL, Stock R, Stone N, Buckstein M, Shao Y, Campbell C, et al. A 2-stage genome-wide association study to identify single nucleotide polymorphisms associated with development of erectile dysfunction following radiation therapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2013;85:e21–8. https://doi.org/10.1016/j.ijrobp.2012.08.003.
Grodsky MB, Sidani SM. Radiation proctopathy. Clin Colon Rectal Surg. 2015;28:103–11. https://doi.org/10.1055/s-0035-1547337.
Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7. https://doi.org/10.1186/s13742-015-0047-8.
Anderson CA, Pettersson FH, Clarke GM, Cardon LR, Morris AP, Zondervan KT. Data quality control in genetic case-control association studies. Nat Protoc. 2010;5:1564–73. https://doi.org/10.1038/nprot.2010.116.
GTEx Consortium, Laboratory, Data Analysis &Coordinating Center (LDACC)—analysis working group, statistical methods groups—analysis working group, enhancing GTEx (eGTEx) groups, NIH common fund, NIH/NCI, et al. Genetic effects on gene expression across human tissues. Nature. 2017;550:204–13. https://doi.org/10.1038/nature24277.
Gaspar HA, Hübel C, Breen G. Drug Targetor: a web interface to investigate the human druggome for over 500 phenotypes. Bioinformatics. 2019;35:2515–7. https://doi.org/10.1093/bioinformatics/bty982.
Basha O, Shpringer R, Argov CM, Yeger-Lotem E. The DifferentialNet database of differential protein-protein interactions in human tissues. Nucleic Acids Res. 2018;46:D522–6. https://doi.org/10.1093/nar/gkx981.
Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019;47:W234–41. https://doi.org/10.1093/nar/gkz240.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7. https://doi.org/10.1089/omi.2011.0118.
Walter W, Sánchez-Cabo F, Ricote M. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics. 2015;31:2912–4. https://doi.org/10.1093/bioinformatics/btv300.
van Oorschot B, Uitterhoeve L, Oomen I, Ten Cate R, Medema JP, Vrieling H, et al. Prostate Cancer patients with late radiation toxicity exhibit reduced expression of genes involved in DNA double-Strand break repair and homologous recombination. Cancer Res. 2017;77:1485–91. https://doi.org/10.1158/0008-5472.CAN-16-1966.
Amaral ML, Erikson GA, Shokhirev MN. BART: bioinformatics array research tool. BMC Bioinformatics. 2018;19:296. https://doi.org/10.1186/s12859-018-2308-x.
Yu YP, Liu S, Huo Z, Martin A, Nelson JB, Tseng GC, et al. Genomic copy number variations in the genomes of leukocytes predict prostate cancer clinical outcomes PLoS One 2015;10:e0135982. doi:https://doi.org/10.1371/journal.pone.0135982.
Kim J-H, Hu H-J, Yim S-H, Bae JS, Kim S-Y, Chung Y-J. CNVRuler: a copy number variation-based case-control association analysis tool. Bioinformatics. 2012;28:1790–2. https://doi.org/10.1093/bioinformatics/bts239.
Wolfe D, Dudek S, Ritchie MD, Pendergrass SA. Visualizing genomic information across chromosomes with PhenoGram. BioData Min. 2013;6:18. https://doi.org/10.1186/1756-0381-6-18.
Kerns SL, Stock RG, Stone NN, Blacksburg SR, Rath L, Vega A, et al. Genome-wide association study identifies a region on chromosome 11q14.3 associated with late rectal bleeding following radiation therapy for prostate cancer. Radiother Oncol. 2013;107:372–6. https://doi.org/10.1016/j.radonc.2013.05.001.
Kerns SL, Dorling L, Fachal L, Bentzen S, Pharoah PDP, Barnes DR, et al. Meta-analysis of genome wide association studies identifies genetic markers of late toxicity following radiotherapy for prostate Cancer. EBioMedicine. 2016;10:150–63. https://doi.org/10.1016/j.ebiom.2016.07.022.
Fachin AL, Mello SS, Sandrin-Garcia P, Junta CM, Ghilardi-Netto T, Donadi EA, et al. Gene expression profiles in radiation workers occupationally exposed to ionizing radiation. J Radiat Res. 2009;50:61–71.
Logue MW, Smith AK, Baldwin C, Wolf EJ, Guffanti G, Ratanatharathorn A, et al. An analysis of gene expression in PTSD implicates genes involved in the glucocorticoid receptor pathway and neural responses to stress. Psychoneuroendocrinology. 2015;57:1–13. https://doi.org/10.1016/j.psyneuen.2015.03.016.
Zhang B, Liu X, Chen W, Chen L. IFIT5 potentiates anti-viral response through enhancing innate immune signaling pathways. Acta Biochim Biophys Sin Shanghai. 2013;45:867–74. https://doi.org/10.1093/abbs/gmt088.
Lo U-G, Pong R-C, Yang D, Gandee L, Hernandez E, Dang A, et al. IFNγ-induced IFIT5 promotes epithelial-to-Mesenchymal transition in prostate Cancer via miRNA processing. Cancer Res. 2019;79:1098–112. https://doi.org/10.1158/0008-5472.CAN-18-2207.
Komaki R, Ueda T, Tsuji Y, Miyawaki T, Kusuhara S, Hara S, et al. Retinal vasculopathy with cerebral leukoencephalopathy carrying TREX1 mutation diagnosed by the intracranial calcification: a case report. Rinsho Shinkeigaku. 2018;58:111–7. https://doi.org/10.5692/clinicalneurol.cn-001096.
Shigechi T, Tomida J, Sato K, Kobayashi M, Eykelenboom JK, Pessina F, et al. ATR-ATRIP kinase complex triggers activation of the Fanconi anemia DNA repair pathway. Cancer Res. 2012;72:1149–56. https://doi.org/10.1158/0008-5472.CAN-11-2904.
Subramanian L, Nakamura TM. A kinase-independent role for the Rad3(ATR)-Rad26(ATRIP) complex in recruitment of Tel1(ATM) to telomeres in fission yeast. PLoS Genet. 2010;6:e1000839. https://doi.org/10.1371/journal.pgen.1000839.
Hilton BA, Li Z, Musich PR, Wang H, Cartwright BM, Serrano M, et al. ATR plays a direct Antiapoptotic role at mitochondria, which is regulated by Prolyl Isomerase Pin1. Mol Cell. 2015;60:35–46. https://doi.org/10.1016/j.molcel.2015.08.008.
Vanpouille-Box C, Alard A, Aryankalayil MJ, Sarfraz Y, Diamond JM, Schneider RJ, et al. DNA exonuclease Trex1 regulates radiotherapy-induced tumour immunogenicity. Nat Commun. 2017;8:15618. https://doi.org/10.1038/ncomms15618.
Diamond JM, Vanpouille-Box C, Spada S, Rudqvist N-P, Chapman JR, Ueberheide BM, et al. Exosomes shuttle TREX1-sensitive IFN-stimulatory dsDNA from irradiated Cancer cells to DCs. Cancer Immunol Res. 2018;6:910–20. https://doi.org/10.1158/2326-6066.CIR-17-0581.
Garcia-Muse T, Boulton SJ. Distinct modes of ATR activation after replication stress and DNA double-strand breaks in Caenorhabditis elegans. EMBO J. 2005;24:4345–55. https://doi.org/10.1038/sj.emboj.7600896.
Nam SY, Seo HH, Park HS, An S, Kim J-Y, Yang KH, et al. Phosphorylation of CLK2 at serine 34 and threonine 127 by AKT controls cell survival after ionizing radiation. J Biol Chem. 2010;285:31157–63. https://doi.org/10.1074/jbc.M110.122044.
Brown MC, Gromeier M. MNK controls mTORC1:substrate association through regulation of TELO2 binding with mTORC1. Cell Rep. 2017;18:1444–57. https://doi.org/10.1016/j.celrep.2017.01.023.
Quek H, Lim YC, Lavin MF, Roberts TL. PIKKing a way to regulate inflammation. Immunol Cell Biol. 2018;96:8–20. https://doi.org/10.1111/imcb.1001.
Fernandez SV, Huang Y, Snider KE, Zhou Y, Pogash TJ, Russo J. Expression and DNA methylation changes in human breast epithelial cells’ ‘ after bisphenol a exposure. Int J Oncol. 2012;41:369–77. https://doi.org/10.3892/ijo.2012.1444.
Rouillard AD, Gundersen GW, Fernandez NF, Wang Z, Monteiro CD, McDermott MG, et al. The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins. Database (Oxford). 2016;2016. https://doi.org/10.1093/database/baw100.
Ronning SA, Heatley GA, Martin TF. Thyrotropin-releasing hormone mobilizes Ca2+ from endoplasmic reticulum and mitochondria of GH3 pituitary cells: characterization of cellular Ca2+ pools by a method based on digitonin permeabilization. Proc Natl Acad Sci U S A. 1982;79:6294–8. https://doi.org/10.1073/pnas.79.20.6294.
Knuever J, Poeggeler B, Gáspár E, Klinger M, Hellwig-Burgel T, Hardenbicker C, et al. Thyrotropin-releasing hormone controls mitochondrial biology in human epidermis. J Clin Endocrinol Metab. 2012;97:978–86. https://doi.org/10.1210/jc.2011-1096.
Southam L, Gilly A, Süveges D, Farmaki A-E, Schwartzentruber J, Tachmazidou I, et al. Whole genome sequencing and imputation in isolated populations identify genetic associations with medically-relevant complex traits. Nat Commun. 2017;8:15606. https://doi.org/10.1038/ncomms15606.
Agnihotri N, Kumar S, Mehta K. Tissue transglutaminase as a central mediator in inflammation-induced progression of breast cancer. Breast Cancer Res. 2013;15:202. https://doi.org/10.1186/bcr3371.
Ooko E, Kadioglu O, Greten HJ, Efferth T. Pharmacogenomic characterization and Isobologram analysis of the combination of ascorbic acid and Curcumin-two Main metabolites of Curcuma longa-in Cancer cells. Front Pharmacol. 2017;8:38. https://doi.org/10.3389/fphar.2017.00038.
Leicht DT, Kausar T, Wang Z, Ferrer-Torres D, Wang TD, Thomas DG, et al. TGM2: a cell surface marker in esophageal adenocarcinomas. J Thorac Oncol. 2014;9:872–81. https://doi.org/10.1097/JTO.0000000000000229.
Garnier D, Meehan B, Kislinger T, Daniel P, Sinha A, Abdulkarim B, et al. Divergent evolution of temozolomide resistance in glioblastoma stem cells is reflected in extracellular vesicles and coupled with radiosensitization. Neuro-Oncology. 2018;20:236–48. https://doi.org/10.1093/neuonc/nox142.
Rajamani D, Singh PK, Rottmann BG, Singh N, Bhasin MK, Kumar A. Temporal retinal transcriptome and systems biology analysis identifies key pathways and hub genes in Staphylococcus aureus endophthalmitis. Sci Rep. 2016;6:21502. https://doi.org/10.1038/srep21502.
Gironella M, Malicet C, Cano C, Sandi MJ, Hamidi T, Tauil RMN, et al. p8/nupr1 regulates DNA-repair activity after double-strand gamma irradiation-induced DNA damage. J Cell Physiol. 2009;221:594–602. https://doi.org/10.1002/jcp.21889.
Hamidi T, Cano CE, Grasso D, Garcia MN, Sandi MJ, Calvo EL, et al. Nupr1-aurora kinase a pathway provides protection against metabolic stress-mediated autophagic-associated cell death. Clin Cancer Res. 2012;18:5234–46. https://doi.org/10.1158/1078-0432.CCR-12-0026.
Santofimia-Castaño P, Lan W, Bintz J, Gayet O, Carrier A, Lomberk G, et al. Inactivation of NUPR1 promotes cell death by coupling ER-stress responses with necrosis. Sci Rep. 2018;8:16999. https://doi.org/10.1038/s41598-018-35020-3.
Chen C-Y, Wu S-M, Lin Y-H, Chi H-C, Lin S-L, Yeh C-T, et al. Induction of nuclear protein-1 by thyroid hormone enhances platelet-derived growth factor a mediated angiogenesis in liver cancer. Theranostics. 2019;9:2361–79. https://doi.org/10.7150/thno.29628.
Brent GA. Mechanisms of thyroid hormone action. J Clin Invest. 2012;122:3035–43. https://doi.org/10.1172/JCI60047.
De Ruysscher D, Niedermann G, Burnet NG, Siva S, Lee AWM, Hegi-Johnson F. Radiotherapy toxicity. Nat Rev Dis Primers. 2019;5:13. https://doi.org/10.1038/s41572-019-0064-5.
Li X, Fang F, Gao Y, Tang G, Xu W, Wang Y, et al. ROS induced by KillerRed targeting mitochondria (mtKR) enhances apoptosis caused by radiation via Cyt c/Caspase-3 pathway. Oxidative Med Cell Longev. 2019;2019:4528616. https://doi.org/10.1155/2019/4528616.
Zorov DB, Juhaszova M, Sollott SJ. Mitochondrial reactive oxygen species (ROS) and ROS-induced ROS release. Physiol Rev. 2014;94:909–50. https://doi.org/10.1152/physrev.00026.2013.
Jang ER, Lee J-S. DNA damage response mediated through BRCA1. Cancer Res Treat. 2004;36:214–21. https://doi.org/10.4143/crt.2004.36.4.214.
We would like to acknowledge Dr. Ranajit Chakraborty for his guidance and interest during the start of the project, we unfortunately lost him in September 2018. We appreciate the data received via authorized access from dbGaP- GenePARE dataset and GEO dataset – GSE85570.
We would like to acknowledge the NIH – Neurobiology of Aging T32 grant AG020494 for supporting this research. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Ethics approval and consent to participate
The Gene-PARE was approved by the Institutional Review Board of the Icahn School of Medicine at Mount Sinai and Florida Radiation Oncology Group. All patients provided informed consent under the parent study – Gene-PARE, at Icahn School of Medicine, Mount Sinai and Florida Radiation Oncology Group (see methods for references). We obtained access to anonymized individual level genotype data from Genetic Predictors of Adverse Radiotherapy Effects (Gene-PARE) (phs000772.v1.p1) via dbGaP’s authorized application, under the approval of North Texas Regional IRB protocol 2016–090. The study described here was performed under the North Texas Regional IRB (formerly the University of North Texas Health Science Center IRB), and was given “EXEMPT” status based on the criteria that our study involved data available from public repository, i.e. dbGaP and does not require approval of receiving informed consent from our institute. (dbGaP - https://www.ncbi.nlm.nih.gov/gap/).
Consent for publication
All authors on the manuscript have provided consent for publication. No individual data was used/discussed in the current manuscript.
The authors have no competing interests to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Pathak, G.A., Polimanti, R., Silzer, T.K. et al. Genetically-regulated transcriptomics & copy number variation of proctitis points to altered mitochondrial and DNA repair mechanisms in individuals of European ancestry. BMC Cancer 20, 954 (2020). https://doi.org/10.1186/s12885-020-07457-1
- DNA repair