Comparison of low and high dose ionising radiation using topological analysis of gene coexpression networks

Short- and long-term effects of radiation exposure at low dose and low dose rate in normal human VH10 fibroblasts

Introduction

Experimental studies complement epidemiological data on the biological effects of low doses and dose rates of ionizing radiation and help in determining the dose and dose rate effectiveness factor.

Methods

Human VH10 skin fibroblasts exposed to 25, 50, and 100 mGy of 137Cs gamma radiation at 1.6, 8, 12 mGy/h, and at a high dose rate of 23.4 Gy/h, were analyzed for radiation-induced short- and long-term effects. Two sample cohorts, i.e., discovery (n = 30) and validation (n = 12), were subjected to RNA sequencing. The pool of the results from those six experiments with shared conditions (1.6 mGy/h; 24 h), together with an earlier time point (0 h), constituted a third cohort (n = 12).

Results

The 100 mGy-exposed cells at all abovementioned dose rates, harvested at 0/24 h and 21 days after exposure, showed no strong gene expression changes. DMXL2, involved in the regulation of the NOTCH signaling pathway, presented a consistent upregulation among both the discovery and validation cohorts, and was validated by qPCR. Gene set enrichment analysis revealed that the NOTCH pathway was upregulated in the pooled cohort (p = 0.76, normalized enrichment score (NES) = 0.86). Apart from upregulated apical junction and downregulated DNA repair, few pathways were consistently changed across exposed cohorts. Concurringly, cell viability assays, performed 1, 3, and 6 days post irradiation, and colony forming assay, seeded just after exposure, did not reveal any statistically significant early effects on cell growth or survival patterns. Tendencies of increased viability (day 6) and reduced colony size (day 21) were observed at 12 mGy/h and 23.4 Gy/min. Furthermore, no long-term changes were observed in cell growth curves generated up to 70 days after exposure.

Discussion

In conclusion, low doses of gamma radiation given at low dose rates had no strong cytotoxic effects on radioresistant VH10 cells.

High-throughput screening strategies for space-based radiation countermeasure discovery

As humanity begins to venture further into space, approaches to better protect astronauts from the hazards found in space need to be developed. One particular hazard of concern is the complex radiation that is ever present in deep space. Currently, it is unlikely enough spacecraft shielding could be launched that would provide adequate protection to astronauts during long-duration missions such as a journey to Mars and back. In an effort to identify other means of protection, prophylactic radioprotective drugs have been proposed as a potential means to reduce the biological damage caused by this radiation. Unfortunately, few radioprotectors have been approved by the FDA for usage and for those that have been developed, they protect normal cells/tissues from acute, high levels of radiation exposure such as that from oncology radiation treatments. To date, essentially no radioprotectors have been developed that specifically counteract the effects of chronic low-dose rate space radiation. This review highlights how high-throughput screening (HTS) methodologies could be implemented to identify such a radioprotective agent. Several potential target, pathway, and phenotypic assays are discussed along with potential challenges towards screening for radioprotectors. Utilizing HTS strategies such as the ones proposed here have the potential to identify new chemical scaffolds that can be developed into efficacious radioprotectors that are specifically designed to protect astronauts during deep space journeys. The overarching goal of this review is to elicit broader interest in applying drug discovery techniques, specifically HTS towards the identification of radiation countermeasures designed to be efficacious towards the biological insults likely to be encountered by astronauts on long duration voyages.

Radiation-response in primary fibroblasts of long-term survivors of childhood cancer with and without second primary neoplasms: the KiKme study

Background

The etiology and most risk factors for a sporadic first primary neoplasm in childhood or subsequent second primary neoplasms are still unknown. One established causal factor for therapy-associated second primary neoplasms is the exposure to ionizing radiation during radiation therapy as a mainstay of cancer treatment. Second primary neoplasms occur in 8% of all cancer survivors within 30 years after the first diagnosis in Germany, but the underlying factors for intrinsic susceptibilities have not yet been clarified. Thus, the purpose of this nested case–control study was the investigation and comparison of gene expression and affected pathways in primary fibroblasts of childhood cancer survivors with a first primary neoplasm only or with at least one subsequent second primary neoplasm, and controls without neoplasms after exposure to a low and a high dose of ionizing radiation.

Methods

Primary fibroblasts were obtained from skin biopsies from 52 adult donors with a first primary neoplasm in childhood (N1), 52 with at least one additional primary neoplasm (N2+), as well as 52 without cancer (N0) from the KiKme study. Cultured fibroblasts were exposed to a high [2 Gray (Gy)] and a low dose (0.05 Gy) of X-rays. Messenger ribonucleic acid was extracted 4 h after exposure and Illumina-sequenced. Differentially expressed genes (DEGs) were computed using limma for R, selected at a false discovery rate level of 0.05, and further analyzed for pathway enrichment (right-tailed Fisher’s Exact Test) and (in-) activation (z ≥|2|) using Ingenuity Pathway Analysis.

Results

After 0.05 Gy, least DEGs were found in N0 (n = 236), compared to N1 (n = 653) and N2+ (n = 694). The top DEGs with regard to the adjusted p-value were upregulated in fibroblasts across all donor groups (SESN1, MDM2, CDKN1A, TIGAR, BTG2, BLOC1S2, PPM1D, PHLDB3, FBXO22, AEN, TRIAP1, and POLH). Here, we observed activation of p53 Signaling in N0 and to a lesser extent in N1, but not in N2+. Only in N0, DNA (excision-) repair (involved genes: CDKN1A, PPM1D, and DDB2) was predicted to be a downstream function, while molecular networks in N2+ were associated with cancer, as well as injury and abnormalities (among others, downregulation of MSH6, CCNE2, and CHUK). After 2 Gy, the number of DEGs was similar in fibroblasts of all donor groups and genes with the highest absolute log₂ fold-change were upregulated throughout (CDKN1A, TIGAR, HSPA4L, MDM2, BLOC1SD2, PPM1D, SESN1, BTG2, FBXO22, PCNA, and TRIAP1). Here, the p53 Signaling-Pathway was activated in fibroblasts of all donor groups. The Mitotic Roles of Polo Like Kinase-Pathway was inactivated in N1 and N2+. Molecular Mechanisms of Cancer were affected in fibroblasts of all donor groups. P53 was predicted to be an upstream regulator in fibroblasts of all donor groups and E2F1 in N1 and N2+. Results of the downstream analysis were senescence in N0 and N2+, transformation of cells in N0, and no significant effects in N1. Seven genes were differentially expressed in reaction to 2 Gy dependent on the donor group (LINC00601, COBLL1, SESN2, BIN3, TNFRSF10A, EEF1AKNMT, and BTG2).

Conclusion

Our results show dose-dependent differences in the radiation response between N1/N2+ and N0. While mechanisms against genotoxic stress were activated to the same extent after a high dose in all groups, the radiation response was impaired after a low dose in N1/N2+, suggesting an increased risk for adverse effects including carcinogenesis, particularly in N2+.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10020-022-00520-6.

Comparison of time and dose dependent gene expression and affected pathways in primary human fibroblasts after exposure to ionizing radiation

Background

Exposure to ionizing radiation induces complex stress responses in cells, which can lead to adverse health effects such as cancer. Although a variety of studies investigated gene expression and affected pathways in human fibroblasts after exposure to ionizing radiation, the understanding of underlying mechanisms and biological effects is still incomplete due to different experimental settings and small sample sizes. Therefore, this study aims to identify the time point with the highest number of differentially expressed genes and corresponding pathways in primary human fibroblasts after irradiation at two preselected time points.

Methods

Fibroblasts from skin biopsies of 15 cell donors were exposed to a high (2Gy) and a low (0.05Gy) dose of X-rays. RNA was extracted and sequenced 2 h and 4 h after exposure. Differentially expressed genes with an adjusted p-value < 0.05 were flagged and used for pathway analyses including prediction of upstream and downstream effects. Principal component analyses were used to examine the effect of two different sequencing runs on quality metrics and variation in expression and alignment and for explorative analysis of the radiation dose and time point of analysis.

Results

More genes were differentially expressed 4 h after exposure to low and high doses of radiation than after 2 h. In experiments with high dose irradiation and RNA sequencing after 4 h, inactivation of the FAT10 cancer signaling pathway and activation of gluconeogenesis I, glycolysis I, and prostanoid biosynthesis was observed taking p-value (< 0.05) and (in) activating z-score (≥2.00 or ≤ − 2.00) into account. Two hours after high dose irradiation, inactivation of small cell lung cancer signaling was observed. For low dose irradiation experiments, we did not detect any significant (p < 0.05 and z-score ≥ 2.00 or ≤ − 2.00) activated or inactivated pathways for both time points.

Conclusions

Compared to 2 h after irradiation, a higher number of differentially expressed genes were found 4 h after exposure to low and high dose ionizing radiation. Differences in gene expression were related to signal transduction pathways of the DNA damage response after 2 h and to metabolic pathways, that might implicate cellular senescence, after 4 h. The time point 4 h will be used to conduct further irradiation experiments in a larger sample.

Effects of High-Dose Ionizing Radiation in Human Gene Expression: A Meta-Analysis

The use of high-dose Ionizing Radiation (IR) is currently one of the most common modalities in treatment of many types of cancer. The objective of this work was to investigate the effects of high-dose ionizing radiation on healthy human tissue, utilizing quantitative analysis of gene expression. To this end, publicly available transcriptomics datasets from human samples irradiated with a high dose of radiation and non-irradiated (control) ones were selected, and gene expression was determined using RNA-Seq data analysis. Raw data from these studies were subjected to quality control and trimming. Mapping of RNA-Seq reads was performed by the partial selective alignment method, and differential gene expression analysis was conducted. Subsequently, a meta-analysis was performed to select differentially expressed genes across datasets. Based on the differentially expressed genes discovered by meta-analysis, we constructed a protein-to-protein interaction network, and we identified biological pathways and processes related to high-dose IR effects. Our findings suggest that cell cycle arrest is activated, supported by our top down-regulated genes associated with cell cycle activation. DNA repair genes are down-regulated in their majority. However, several genes implicated in the nucleotide excision repair pathway are upregulated. Nevertheless, apoptotic mechanisms seem to be activated probably due to severe high-dose-induced complex DNA damage. The significant upregulation of CDKN1A, as a downstream gene of TP53, further validates programmed cell death. Finally, down-regulation of TIMELESS, signifies a correlation between IR response and circadian rhythm. Nonetheless, high-dose IR exposure effects regarding normal tissue (radiation toxicity) and its possible long-term outcomes should be studied to a greater extend.

The Role of the Nuclear Factor κB Pathway in the Cellular Response to Low and High Linear Energy Transfer Radiation

Astronauts are exposed to considerable doses of space radiation during long-term space missions. As complete shielding of the highly energetic particles is impracticable, the cellular response to space-relevant radiation qualities has to be understood in order to develop countermeasures and to reduce radiation risk uncertainties. The transcription factor Nuclear Factor κB (NF-κB) plays a fundamental role in the immune response and in the pathogenesis of many diseases. We have previously shown that heavy ions with a linear energy transfer (LET) of 100–300 keV/µm have a nine times higher potential to activate NF-κB compared to low-LET X-rays. Here, chemical inhibitor studies using human embryonic kidney cells (HEK) showed that the DNA damage sensor Ataxia telangiectasia mutated (ATM) and the proteasome were essential for NF-κB activation in response to X-rays and heavy ions. NF-κB’s role in cellular radiation response was determined by stable knock-down of the NF-κB subunit RelA. Transfection of a RelA short-hairpin RNA plasmid resulted in higher sensitivity towards X-rays, but not towards heavy ions. Reverse Transcriptase real-time quantitative PCR (RT-qPCR) showed that after exposure to X-rays and heavy ions, NF-κB predominantly upregulates genes involved in intercellular communication processes. This process is strictly NF-κB dependent as the response is completely absent in RelA knock-down cells. NF-κB’s role in the cellular radiation response depends on the radiation quality.

Low dose of uranium induces multigenerational epigenetic effects in rat kidney

PURPOSE

A protocol of chronic exposure to low dose of uranium was established in order to distinguish the sexual differences and the developmental process that are critical windows for epigenetic effects over generations.

METHODS

Both male and female rats were contaminated through their drinking water with a non-toxic solution of uranyl nitrate for 9 months. The exposed generation (F0) and the following two generations (F1 and F2) were examined. Clinical monitoring, global DNA methylation profile and DNA methyltransferases (DNMTs) gene expression were analyzed in kidneys.

RESULTS

While the body weight of F1 males increased, a small decrease in kidney and body weight was observed in F2 males. In addition, global DNA hypermethylation profile in kidney cells was observed in F1 and F2 males. qPCR results reveal a significant increase of methyltransferase genes expression (DNMT1 and DNMT3a) for F2 females.

CONCLUSIONS

In the field of public health policy and to raise attention to generational effects for the risk assessment of the environmental exposures, low doses of uranium do not imply clinical effects on adult exposed rats. However, our results confirm the importance of the developmental windows' sensitivity in addition to the sexual dimorphisms of the offspring.

Phosphoprotein profiles of candidate markers for early cellular responses to low-dose γ-radiation in normal human fibroblast cells

Ionizing radiation causes biological damage that leads to severe health effects. However, the effects and subsequent health implications caused by exposure to low-dose radiation are unclear. The objective of this study was to determine phosphoprotein profiles in normal human fibroblast cell lines in response to low-dose and high-dose γ-radiation. We examined the cellular response in MRC-5 cells 0.5 h after exposure to 0.05 or 2 Gy. Using 1318 antibodies by antibody array, we observed ≥1.3-fold increases in a number of identified phosphoproteins in cells subjected to low-dose (0.05 Gy) and high-dose (2 Gy) radiation, suggesting that both radiation levels stimulate distinct signaling pathways. Low-dose radiation induced nucleic acid-binding transcription factor activity, developmental processes, and multicellular organismal processes. By contrast, high-dose radiation stimulated apoptotic processes, cell adhesion and regulation, and cellular organization and biogenesis. We found that phospho-BTK (Tyr550) and phospho-Gab2 (Tyr643) protein levels at 0.5 h after treatment were higher in cells subjected to low-dose radiation than in cells treated with high-dose radiation. We also determined that the phosphorylation of BTK and Gab2 in response to ionizing radiation was regulated in a dose-dependent manner in MRC-5 and NHDF cells. Our study provides new insights into the biological responses to low-dose γ-radiation and identifies potential candidate markers for monitoring exposure to low-dose ionizing radiation.

A secretome analysis reveals that PPARα is upregulated by fractionated-dose γ-irradiation in three-dimensional keratinocyte cultures

Studies have shown that γ-irradiation induces various biological responses, including oxidative stress and apoptosis, as well as cellular repair and immune system responses. However, most such studies have been performed using traditional two-dimensional cell culture systems, which are limited in their ability to faithfully represent in vivo conditions. A three-dimensional (3D) environment composed of properly interconnected and differentiated cells that allow communication and cooperation among cells via secreted molecules would be expected to more accurately reflect cellular responses. Here, we investigated γ-irradiation-induced changes in the secretome of 3D-cultured keratinocytes. An analysis of keratinocyte secretome profiles following fractionated-dose γ-irradiation revealed changes in genes involved in cell adhesion, angiogenesis, and the immune system. Notably, peroxisome proliferator-activated receptor-α (PPARα) was upregulated in response to fractionated-dose γ-irradiation. This upregulation was associated with an increase in the transcription of known PPARα target genes in secretome, including angiopoietin-like protein 4, dermokine and kallikrein-related peptide 12, which were differentially regulated by fractionated-dose γ-irradiation. Collectively, our data imply a mechanism linking γ-irradiation and secretome changes, and suggest that these changes could play a significant role in the coordinated cellular responses to harmful ionizing radiation, such as those associated with radiation therapy. This extension of our understanding of γ-irradiation-induced secretome changes has the potential to improve radiation therapy strategies.

Radiation Metabolomics: Current Status and Future Directions

Human exposure to ionizing radiation (IR) disrupts normal metabolic processes in cells and organs by inducing complex biological responses that interfere with gene and protein expression. Conventional dosimetry, monitoring of prodromal symptoms, and peripheral lymphocyte counts are of limited value as organ- and tissue-specific biomarkers for personnel exposed to radiation, particularly, weeks or months after exposure. Analysis of metabolites generated in known stress-responsive pathways by molecular profiling helps to predict the physiological status of an individual in response to environmental or genetic perturbations. Thus, a multi-metabolite profile obtained from a high-resolution mass spectrometry-based metabolomics platform offers potential for identification of robust biomarkers to predict radiation toxicity of organs and tissues resulting from exposures to therapeutic or non-therapeutic IR. Here, we review the status of radiation metabolomics and explore applications as a standalone technology, as well as its integration in systems biology, to facilitate a better understanding of the molecular basis of radiation response. Finally, we draw attention to the identification of specific pathways that can be targeted for the development of therapeutics to alleviate or mitigate harmful effects of radiation exposure.

Global Gene Expression Alterations as a Crucial Constituent of Human Cell Response to Low Doses of Ionizing Radiation Exposure

Exposure to ionizing radiation (IR) is inevitable to humans in real-life scenarios; the hazards of IR primarily stem from its mutagenic, carcinogenic, and cell killing ability. For many decades, extensive research has been conducted on the human cell responses to IR delivered at a low dose/low dose (LD) rate. These studies have shown that the molecular-, cellular-, and tissue-level responses are different after low doses of IR (LDIR) compared to those observed after a short-term high-dose IR exposure (HDIR). With the advent of high-throughput technologies in the late 1990s, such as DNA microarrays, changes in gene expression have also been found to be ubiquitous after LDIR. Very limited subset of genes has been shown to be consistently up-regulated by LDIR, including CDKN1A. Further research on the biological effects and mechanisms induced by IR in human cells demonstrated that the molecular and cellular processes, including transcriptional alterations, activated by LDIR are often related to protective responses and, sometimes, hormesis. Following LDIR, some distinct responses were observed, these included bystander effects, and adaptive responses. Changes in gene expression, not only at the level of mRNA, but also miRNA, have been found to crucially underlie these effects having implications for radiation protection purposes.

Chronic low-dose γ-irradiation of Drosophila melanogaster larvae induces gene expression changes and enhances locomotive behavior

Although radiation effects have been extensively studied, the biological effects of low-dose radiation (LDR) are controversial. This study investigates LDR-induced alterations in locomotive behavior and gene expression profiles of Drosophila melanogaster. We measured locomotive behavior using larval pupation height and the rapid iterative negative geotaxis (RING) assay after exposure to 0.1 Gy γ-radiation (dose rate of 16.7 mGy/h). We also observed chronic LDR effects on development (pupation and eclosion rates) and longevity (life span). To identify chronic LDR effects on gene expression, we performed whole-genome expression analysis using gene-expression microarrays, and confirmed the results using quantitative real-time PCR. The pupation height of the LDR-treated group at the first larval instar was significantly higher (∼2-fold increase in PHI value, P < 0.05). The locomotive behavior of LDR-treated male flies (∼3 - 5 weeks of age) was significantly increased by 7.7%, 29% and 138%, respectively (P < 0.01), but pupation and eclosion rates and life spans were not significantly altered. Genome-wide expression analysis identified 344 genes that were differentially expressed in irradiated larvae compared with in control larvae. We identified several genes belonging to larval behavior functional groups such as locomotion (1.1%), oxidation reduction (8.0%), and genes involved in conventional functional groups modulated by irradiation such as defense response (4.9%), and sensory and perception (2.5%). Four candidate genes were confirmed as differentially expressed genes in irradiated larvae using qRT-PCR (>2-fold change). These data suggest that LDR stimulates locomotion-related genes, and these genes can be used as potential markers for LDR.

Donor-specific cell-based assays in studying sensitivity to low-dose radiation: a population-based perspective

Currently, a linear no-threshold model is used to estimate health risks associated with exposure to low-dose radiation, a prevalent exposure in the general population, because the direct estimation from epidemiological studies suffers from uncertainty. This model has been criticized based on unique biology of low-dose radiation. Whether the departure from linearity is toward increased or decreased risk is intensely debated. We present an approach based on individual radiosensitivity testing and discuss how individual radiosensitivity can be assessed with the goal to develop a quantifiable measure of cellular response that can be conducted via high-throughput population testing.

Differential network analyses of Alzheimer's disease identify early events in Alzheimer's disease pathology

In late-onset Alzheimer's disease (AD), multiple brain regions are not affected simultaneously. Comparing the gene expression of the affected regions to identify the differences in the biological processes perturbed can lead to greater insight into AD pathogenesis and early characteristics. We identified differentially expressed (DE) genes from single cell microarray data of four AD affected brain regions: entorhinal cortex (EC), hippocampus (HIP), posterior cingulate cortex (PCC), and middle temporal gyrus (MTG). We organized the DE genes in the four brain regions into region-specific gene coexpression networks. Differential neighborhood analyses in the coexpression networks were performed to identify genes with low topological overlap (TO) of their direct neighbors. The low TO genes were used to characterize the biological differences between two regions. Our analyses show that increased oxidative stress, along with alterations in lipid metabolism in neurons, may be some of the very early events occurring in AD pathology. Cellular defense mechanisms try to intervene but fail, finally resulting in AD pathology as the disease progresses. Furthermore, disease annotation of the low TO genes in two independent protein interaction networks has resulted in association between cancer, diabetes, renal diseases, and cardiovascular diseases.

Intensity modulated radiotherapy induces pro-inflammatory and pro-survival responses in prostate cancer patients

Intensity modulated radiotherapy (IMRT) is one of the modern conformal radiotherapies that is widely used within the context of cancer patient treatment. It uses multiple radiation beams targeted to the tumor, however, large volumes of the body receive low doses of irradiation. Using γ-H2AX and global genome expression analysis, we studied the biological responses induced by low doses of ionizing radiation in prostate cancer patients following IMRT. By means of different bioinformatics analyses, we report that IMRT induced an inflammatory response via the induction of viral, adaptive, and innate immune signaling. In response to growth factors and immune-stimulatory signaling, positive regulation in the progression of cell cycle and DNA replication were induced. This denotes pro-inflammatory and pro-survival responses. Furthermore, double strand DNA breaks were induced in every patient 30 min after the treatment and remaining DNA repair and damage signaling continued after 18-24 h. Nine genes belonging to inflammatory responses (TLR3, SH2D1A and IL18), cell cycle progression (ORC4, SMC2 and CCDC99) and DNA damage and repair (RAD17, SMC6 and MRE11A) were confirmed by quantitative RT-PCR. This study emphasizes that the risk assessment of health effects from the out-of-field low doses during IMRT should be of concern, as these may increase the risk of secondary cancers and/or systemic inflammation.

iPSC-derived neurons as a higher-throughput readout for autism: promises and pitfalls

The elucidation of disease etiologies and establishment of robust, scalable, high-throughput screening assays for autism spectrum disorders (ASDs) have been impeded by both inaccessibility of disease-relevant neuronal tissue and the genetic heterogeneity of the disorder. Neuronal cells derived from induced pluripotent stem cells (iPSCs) from autism patients may circumvent these obstacles and serve as relevant cell models. To date, derived cells are characterized and screened by assessing their neuronal phenotypes. These characterizations are often etiology-specific or lack reproducibility and stability. In this review, we present an overview of efforts to study iPSC-derived neurons as a model for autism, and we explore the plausibility of gene expression profiling as a reproducible and stable disease marker.

Co-expression network analysis and genetic algorithms for gene prioritization in preeclampsia

BACKGROUND

In this study, we explored the gene prioritization in preeclampsia, combining co-expression network analysis and genetic algorithms optimization approaches. We analysed five public projects obtaining 1,146 significant genes after cross-platform and processing of 81 and 149 microarrays in preeclamptic and normal conditions, respectively.

METHODS

After co-expression network construction, modular and node analysis were performed using several approaches. Moreover, genetic algorithms were also applied in combination with the nearest neighbour and discriminant analysis classification methods.

RESULTS

Significant differences were found in the genes connectivity distribution, both in normal and preeclampsia conditions pointing to the need and importance of examining connectivity alongside expression for prioritization. We discuss the global as well as intra-modular connectivity for hubs detection and also the utility of genetic algorithms in combination with the network information. FLT1, LEP, INHA and ENG genes were identified according to the literature, however, we also found other genes as FLNB, INHBA, NDRG1 and LYN highly significant but underexplored during normal pregnancy or preeclampsia.

CONCLUSIONS

Weighted genes co-expression network analysis reveals a similar distribution along the modules detected both in normal and preeclampsia conditions. However, major differences were obtained by analysing the nodes connectivity. All models obtained by genetic algorithm procedures were consistent with a correct classification, higher than 90%, restricting to 30 variables in both classification methods applied.Combining the two methods we identified well known genes related to preeclampsia, but also lead us to propose new candidates poorly explored or completely unknown in the pathogenesis of preeclampsia, which may have to be validated experimentally.