Dose-specific transcriptional responses in thyroid tissue in mice after (131)I administration

Development of a machine learning framework for radiation biomarker discovery and absorbed dose prediction

Background

Molecular radiation biomarkers are an emerging tool in radiation research with applications for cancer radiotherapy, radiation risk assessment, and even human space travel. However, biomarker screening in genome-wide expression datasets using conventional tools is time-consuming and underlies analyst (human) bias. Machine Learning (ML) methods can improve the sensitivity and specificity of biomarker identification, increase analytical speed, and avoid multicollinearity and human bias.

Aim

To develop a resource-efficient ML framework for radiation biomarker discovery using gene expression data from irradiated normal tissues. Further, to identify biomarker panels predicting radiation dose with tissue specificity.

Methods

A strategic search in the Gene Expression Omnibus database identified a transcriptomic dataset (GSE44762) for normal tissues radiation responses (murine kidney cortex and medulla) suited for biomarker discovery using an ML approach. The dataset was pre-processed in R and separated into train and test data subsets. High computational cost of Genetic Algorithm/k-Nearest Neighbor (GA/KNN) mandated optimization and 13 ML models were tested using the caret package in R. Biomarker performance was evaluated and visualized via Principal Component Analysis (PCA) and dose regression. The novelty of ML-identified biomarker panels was evaluated by literature search.

Results

Caret-based feature selection and ML methods vastly improved processing time over the GA approach. The KNN method yielded overall best performance values on train and test data and was implemented into the framework. The top-ranking genes were Cdkn1a, Gria3, Mdm2 and Plk2 in cortex, and Brf2, Ccng1, Cdkn1a, Ddit4l, and Gria3 in medulla. These candidates successfully categorized dose groups and tissues in PCA. Regression analysis showed that correlation between predicted and true dose was high with R² of 0.97 and 0.99 for cortex and medulla, respectively.

Conclusion

The caret framework is a powerful tool for radiation biomarker discovery optimizing performance with resource-efficiency for broad implementation in the field. The KNN-based approach identified Brf2, Ddit4l, and Gria3 mRNA as novel candidates that have been uncharacterized as radiation biomarkers to date. The biomarker panel showed good performance in dose and tissue separation and dose regression. Further training with larger cohorts is warranted to improve accuracy, especially for lower doses.

Age and sex effects across the blood proteome after ionizing radiation exposure can bias biomarker screening and risk assessment

Molecular biomarkers of ionizing radiation (IR) exposure are a promising new tool in various disciplines: they can give necessary information for adaptive treatment planning in cancer radiotherapy, enable risk projection for radiation-induced survivorship diseases, or facilitate triage and intervention in radiation hazard events. However, radiation biomarker discovery has not yet resolved the most basic features of personalized medicine: age and sex. To overcome this critical bias in biomarker identification, we quantitated age and sex effects and assessed their relevance in the radiation response across the blood proteome. We used high-throughput mass spectrometry on blood plasma collected 24 h after 0.5 Gy total body irradiation (15 MV nominal photon energy) from male and female C57BL/6 N mice at juvenile (7-weeks-old) or adult (18-weeks-old) age. We also assessed sex and strain effects using juvenile male and female BALB/c nude mice. We showed that age and sex created significant effects in the proteomic response regarding both extent and functional quality of IR-induced responses. Furthermore, we found that age and sex effects appeared non-linear and were often end-point specific. Overall, age contributed more to differences in the proteomic response than sex, most notably in immune responses, oxidative stress, and apoptotic cell death. Interestingly, sex effects were pronounced for DNA damage and repair pathways and associated cellular outcome (pro-survival vs. pro-apoptotic). Only one protein (AHSP) was identified as a potential general biomarker candidate across age and sex, while GMNN, REG3B, and SNCA indicated some response similarity across age. This low yield advocated that unisex or uniage biomarker screening approaches are not feasible. In conclusion, age- and sex-specific screening approaches should be implemented as standard protocol to ensure robustness and diagnostic power of biomarker candidates. Bias-free molecular biomarkers are a necessary progression towards personalized medicine and integral for advanced adaptive cancer radiotherapy and risk assessment.

Age-related long-term response in rat thyroid tissue and plasma after internal low dose exposure to 131I

131I is used clinically for therapy, and may be released during nuclear accidents. After the Chernobyl accident papillary thyroid carcinoma incidence increased in children, but not adults. The aims of this study were to compare 131I irradiation-dependent differences in RNA and protein expression in the thyroid and plasma of young and adult rats, and identify potential age-dependent biomarkers for 131I exposure. Twelve young (5 weeks) and twelve adult Sprague Dawley rats (17 weeks) were i.v. injected with 50 kBq 131I (absorbed dose to thyroid = 0.1 Gy), and sixteen unexposed age-matched rats were used as controls. The rats were killed 3-9 months after administration. Microarray analysis was performed using RNA from thyroid samples, while LC-MS/MS analysis was performed on proteins extracted from thyroid tissue and plasma. Canonical pathways, biological functions and upstream regulators were analysed for the identified transcripts and proteins. Distinct age-dependent differences in gene and protein expression were observed. Novel biomarkers for thyroid 131I exposure were identified: (PTH), age-dependent dose response (CA1, FTL1, PVALB (youngsters) and HSPB6 (adults)), thyroid function (Vegfb (adults)). Further validation using clinical samples are needed to explore the role of the identified biomarkers.

Long-term transcriptomic and proteomic effects in Sprague Dawley rat thyroid and plasma after internal low dose 131I exposure

Background

Radioiodide (¹³¹I) is commonly used to treat thyroid cancer and hyperthyroidis.¹³¹I released during nuclear accidents, have resulted in increased incidence of thyroid cancer in children. Therefore, a better understanding of underlying cellular mechanisms behind ¹³¹I exposure is of great clinical and radiation protection interest. The aim of this work was to study the long-term dose-related effects of ¹³¹I exposure in thyroid tissue and plasma in young rats and identify potential biomarkers.

Materials and methods

Male Sprague Dawley rats (5-week-old) were i.v. injected with 0.5, 5.0, 50 or 500 kBq ¹³¹I (D_thyroid ca 1–1000 mGy), and killed after nine months at which time the thyroid and blood samples were collected. Gene expression microarray analysis (thyroid samples) and LC-MS/MS analysis (thyroid and plasma samples) were performed to assess differential gene and protein expression profiles in treated and corresponding untreated control samples. Bioinformatics analyses were performed using the DAVID functional annotation tool and Ingenuity Pathway Analysis (IPA). The gene expression microarray data and LC-MS/MS data were validated using qRT-PCR and ELISA, respectively.

Results

Nine ¹³¹I exposure-related candidate biomarkers (transcripts: Afp and RT1-Bb, and proteins: ARF3, DLD, IKBKB, NONO, RAB6A, RPN2, and SLC25A5) were identified in thyroid tissue. Two dose-related protein candidate biomarkers were identified in thyroid (APRT and LDHA) and two in plasma (DSG4 and TGM3). Candidate biomarkers for thyroid function included the ACADL and SORBS2 (all activities), TPO and TG proteins (low activities). ¹³¹I exposure was shown to have a profound effect on metabolism, immune system, apoptosis and cell death. Furthermore, several signalling pathways essential for normal cellular function (actin cytoskeleton signalling, HGF signalling, NRF2-mediated oxidative stress, integrin signalling, calcium signalling) were also significantly regulated.

Conclusion

Exposure-related and dose-related effects on gene and protein expression generated few expression patterns useful as biomarkers for thyroid function and cancer.

The IRI-DICE hypothesis: ionizing radiation-induced DSBs may have a functional role for non-deterministic responses at low doses

Low-dose ionizing radiation (IR) responses remain an unresolved issue in radiation biology and risk assessment. Accurate knowledge of low-dose responses is important for estimation of normal tissue risk in cancer radiotherapy or health risks from occupational or hazard exposure. Cellular responses to low-dose IR appear diverse and stochastic in nature and to date no model has been proposed to explain the underlying mechanisms. Here, we propose a hypothesis on IR-induced double-strand break (DSB)-induced cis effects (IRI-DICE) and introduce DNA sequence functionality as a submicron-scale target site with functional outcome on gene expression: DSB induction in a certain genetic target site such as promotor, regulatory element, or gene core would lead to changes in transcript expression, which may range from suppression to overexpression depending on which functional element was damaged. The DNA damage recognition and repair machinery depicts threshold behavior requiring a certain number of DSBs for induction. Stochastically distributed persistent disruption of gene expression may explain-in part-the diverse nature of low-dose responses until the repair machinery is initiated at increased absorbed dose. Radiation quality and complexity of DSB lesions are also discussed. Currently, there are no technologies available to irradiate specific genetic sites to test the IRI-DICE hypothesis directly. However, supportive evidence may be achieved by developing a computational model that combines radiation transport codes with a genomic DNA model that includes sequence functionality and transcription to simulate expression changes in an irradiated cell population. To the best of our knowledge, IRI-DICE is the first hypothesis that includes sequence functionality of different genetic elements in the radiation response and provides a model for the diversity of radiation responses in the (very) low dose regimen.

Dissecting the Role of Thyrotropin in the DNA Damage Response in Human Thyrocytes after 131I, γ Radiation and H2O2

BACKGROUND

The early molecular events in human thyrocytes after 131I exposure have not yet been unravelled. Therefore, we investigated the role of TSH in the 131I-induced DNA damage response and gene expression in primary cultured human thyrocytes.

METHODS

Following exposure of thyrocytes, in the presence or absence of TSH, to 131I (β radiation), γ radiation (3 Gy), and hydrogen peroxide (H2O2), we assessed DNA damage, proliferation, and cell-cycle status. We conducted RNA sequencing to profile gene expression after each type of exposure and evaluated the influence of TSH on each transcriptomic response.

RESULTS

Overall, the thyrocyte responses following exposure to β or γ radiation and to H2O2 were similar. However, TSH increased 131I-induced DNA damage, an effect partially diminished after iodide uptake inhibition. Specifically, TSH increased the number of DNA double-strand breaks in nonexposed thyrocytes and thus predisposed them to greater damage following 131I exposure. This effect most likely occurred via Gα q cascade and a rise in intracellular reactive oxygen species (ROS) levels. β and γ radiation prolonged thyroid cell-cycle arrest to a similar extent without sign of apoptosis. The gene expression profiles of thyrocytes exposed to β/γ radiation or H2O2 were overlapping. Modulations in genes involved in inflammatory response, apoptosis, and proliferation were observed. TSH increased the number and intensity of modulation of differentially expressed genes after 131I exposure.

CONCLUSIONS

TSH specifically increased 131I-induced DNA damage probably via a rise in ROS levels and produced a more prominent transcriptomic response after exposure to 131I.

The CD73/Ado System-A New Player in RT Induced Adverse Late Effects

Radiotherapy (RT) is a central component of standard treatment for many cancer patients. RT alone or in multimodal treatment strategies has a documented contribution to enhanced local control and overall survival of cancer patients, and cancer cure. Clinical RT aims at maximizing tumor control, while minimizing the risk for RT-induced adverse late effects. However, acute and late toxicities of IR in normal tissues are still important biological barriers to successful RT: While curative RT may not be tolerable, sub-optimal tolerable RT doses will lead to fatal outcomes by local recurrence or metastatic disease, even when accepting adverse normal tissue effects that decrease the quality of life of irradiated cancer patients. Technical improvements in treatment planning and the increasing use of particle therapy have allowed for a more accurate delivery of IR to the tumor volume and have thereby helped to improve the safety profile of RT for many solid tumors. With these technical and physical strategies reaching their natural limits, current research for improving the therapeutic gain of RT focuses on innovative biological concepts that either selectively limit the adverse effects of RT in normal tissues without protecting the tumor or specifically increase the radiosensitivity of the tumor tissue without enhancing the risk of normal tissue complications. The biology-based optimization of RT requires the identification of biological factors that are linked to differential radiosensitivity of normal or tumor tissues, and are amenable to therapeutic targeting. Extracellular adenosine is an endogenous mediator critical to the maintenance of homeostasis in various tissues. Adenosine is either released from stressed or injured cells or generated from extracellular adenine nucleotides by the concerted action of the ectoenzymes ectoapyrase (CD39) and 5′ ectonucleotidase (NT5E, CD73) that catabolize ATP to adenosine. Recent work revealed a role of the immunoregulatory CD73/adenosine system in radiation-induced fibrotic disease in normal tissues suggesting a potential use as novel therapeutic target for normal tissue protection. The present review summarizes relevant findings on the pathologic roles of CD73 and adenosine in radiation-induced fibrosis in different organs (lung, skin, gut, and kidney) that have been obtained in preclinical models and proposes a refined model of radiation-induced normal tissue toxicity including the disease-promoting effects of radiation-induced activation of CD73/adenosine signaling in the irradiated tissue environment. However, expression and activity of the CD73/adenosine system in the tumor environment has also been linked to increased tumor growth and tumor immune escape, at least in preclinical models. Therefore, we will discuss the use of pharmacologic inhibition of CD73/adenosine-signaling as a promising strategy for improving the therapeutic gain of RT by targeting both, malignant tumor growth and adverse late effects of RT with a focus on fibrotic disease. The consideration of the therapeutic window is particularly important in view of the increasing use of RT in combination with various molecularly targeted agents and immunotherapy to enhance the tumor radiation response, as such combinations may result in increased or novel toxicities, as well as the increasing number of cancer survivors.

Deconvolution of expression microarray data reveals 131I-induced responses otherwise undetected in thyroid tissue

High-throughput gene expression analysis is increasingly used in radiation research for discovery of damage-related or absorbed dose-dependent biomarkers. In tissue samples, cell type-specific responses can be masked in expression data due to mixed cell populations which can preclude biomarker discovery. In this study, we deconvolved microarray data from thyroid tissue in order to assess possible bias from mixed cell type data. Transcript expression data [GSE66303] from mouse thyroid that received 5.9 Gy from ¹³¹I over 24 h (or 0 Gy from mock treatment) were deconvolved by cell frequency of follicular cells and C-cells using csSAM and R and processed with Nexus Expression. Literature-based signature genes were used to assess the relative impact from ionizing radiation (IR) or thyroid hormones (TH). Regulation of cellular functions was inferred by enriched biological processes according to Gene Ontology terms. We found that deconvolution increased the detection rate of significantly regulated transcripts including the biomarker candidate family of kallikrein transcripts. Detection of IR-associated and TH-responding signature genes was also increased in deconvolved data, while the dominating trend of TH-responding genes was reproduced. Importantly, responses in biological processes for DNA integrity, gene expression integrity, and cellular stress were not detected in convoluted data–which was in disagreement with expected dose-response relationships–but upon deconvolution in follicular cells and C-cells. In conclusion, previously reported trends of ¹³¹I-induced transcriptional responses in thyroid were reproduced with deconvolved data and usually with a higher detection rate. Deconvolution also resolved an issue with detecting damage and stress responses in enriched data, and may reduce false negatives in other contexts as well. These findings indicate that deconvolution can optimize microarray data analysis of heterogeneous sample material for biomarker screening or other clinical applications.

Gene expression profiling of rat livers after continuous whole-body exposure to low-dose rate of gamma rays

PURPOSE

To study gene expression modulation in response to continuous whole-body exposure to low-dose-rate gamma radiation and improve our understanding of the mechanism of this impact at the molecular basis.

MATERIALS AND METHODS

cDNA microarray method with complete pooling of samples was used to study expression changes in the transcriptome profile of livers from rats treated with prolonged low-dose-rate ionizing radiation (IR) relative to that of sham-irradiated rats.

RESULTS

Of the 209 genes that were two-fold-up or down-regulated, 143 were known genes of which 27 were found in previous literatures to be modulated by IR. Remarkably, there were a significant number of differentially expressed genes involved in hepatic lipid metabolism.

CONCLUSION

This study showed changes in transcriptome profile of livers from low-dose irradiated rats when compared with that of sham-irradiated ones. This study will be useful for studying the metabolic changes of human exposed for long term to cosmic ray such as in space and in polar regions.

Transcriptional response to 131I exposure of rat thyroid gland

Humans are exposed to ¹³¹I in medical diagnostics and treatment but also from nuclear accidents, and better knowledge of the molecular response in thyroid is needed. The aim of the study was to examine the transcriptional response in thyroid tissue 24 h after ¹³¹I administration in rats. The exposure levels were chosen to simulate both the clinical situation and the case of nuclear fallout. Thirty-six male rats were i.v. injected with 0–4700 kBq ¹³¹I, and killed at 24 h after injection (D_thyroid = 0.0058–3.0 Gy). Total RNA was extracted from individual thyroid tissue samples and mRNA levels were determined using oligonucleotide microarray technique. Differentially expressed transcripts were determined using Nexus Expression 3.0. Hierarchical clustering was performed in the R statistical computing environment. Pathway analysis was performed using the Ingenuity Pathway Analysis tool and the Gene Ontology database. T4 and TSH plasma concentrations were measured using ELISA. Totally, 429 differentially regulated transcripts were identified. Downregulation of thyroid hormone biosynthesis associated genes (e.g. thyroglobulin, thyroid peroxidase, the sodium-iodine symporter) was identified in some groups, and an impact on thyroid function was supported by the pathway analysis. Recurring downregulation of Dbp and Slc47a2 was found. Dbp exhibited a pattern with monotonous reduction of downregulation with absorbed dose at 0.0058–0.22 Gy. T4 plasma levels were increased and decreased in rats whose thyroids were exposed to 0.057 and 0.22 Gy, respectively. Different amounts of injected ¹³¹I gave distinct transcriptional responses in the rat thyroid. Transcriptional response related to thyroid function and changes in T4 plasma levels were found already at very low absorbed doses to thyroid.

Gene signature of the post-Chernobyl papillary thyroid cancer

Purpose

Following the nuclear accidents in Chernobyl and later in Fukushima, the nuclear community has been faced with important issues concerning how to search for and diagnose biological consequences of low-dose internal radiation contamination. Although after the Chernobyl accident an increase in childhood papillary thyroid cancer (PTC) was observed, it is still not clear whether the molecular biology of PTCs associated with low-dose radiation exposure differs from that of sporadic PTC.

Methods

We investigated tissue samples from 65 children/young adults with PTC using DNA microarray (Affymetrix, Human Genome U133 2.0 Plus) with the aim of identifying molecular differences between radiation-induced (exposed to Chernobyl radiation, ECR) and sporadic PTC. All participants were resident in the same region so that confounding factors related to genetics or environment were minimized.

Results

There were small but significant differences in the gene expression profiles between ECR and non-ECR PTC (global test, p < 0.01), with 300 differently expressed probe sets (p < 0.001) corresponding to 239 genes. Multifactorial analysis of variance showed that besides radiation exposure history, the BRAF mutation exhibited independent effects on the PTC expression profile; the histological subset and patient age at diagnosis had negligible effects. Ten genes (PPME1, HDAC11, SOCS7, CIC, THRA, ERBB2, PPP1R9A, HDGF, RAD51AP1, and CDK1) from the 19 investigated with quantitative RT-PCR were confirmed as being associated with radiation exposure in an independent, validation set of samples.

Conclusion

Significant, but subtle, differences in gene expression in the post-Chernobyl PTC are associated with previous low-dose radiation exposure.

Electronic supplementary material

The online version of this article (doi:10.1007/s00259-015-3303-3) contains supplementary material, which is available to authorized users.

Circadian rhythm influences genome-wide transcriptional responses to (131)I in a tissue-specific manner in mice

Background

Circadian variation of gene expression is often neglected when ionizing radiation-induced effects are studied, whether in animal models or in cell culture. This study characterized diurnal variation of genome-wide transcriptional regulation and responses of potential biomarkers and signature genes in normal mouse tissues at 24 h after i.v. administration of ¹³¹I.

Methods

Female BALB/c nude mice were i.v. injected with 90 kBq ¹³¹I at 9:00 a.m., 12:00 p.m., or 3:00 p.m. and killed after 24 h (n = 4/group). Paired control groups were mock-treated (n = 3–4/group). The kidneys, liver, lungs, spleen, and thyroid were excised, snap-frozen, and stored at −80 °C until extraction of total RNA. RNA microarray technology was used for genome-wide expression analysis. Enriched biological processes were categorized after cellular function. Signature genes for ionizing radiation and thyroid hormone-induced responses were taken from the literature. Absorbed dose was estimated using the Medical Internal Radiation Dose (MIRD) formalism.

Results

The thyroid received an absorbed dose of 5.9 Gy and non-thyroid tissues received 0.75–2.2 mGy over 24 h. A distinct peak in the total number of significantly regulated transcripts was observed at 9:00 a.m. in the thyroid, but 3 h later in the kidney cortex, kidney medulla, and liver. Transcriptional regulation in the lungs and spleen was marginal. Associated cellular functions generally varied in quality and response strength between morning, noon, and afternoon. In the thyroid, 25 genes were significantly regulated at all investigated times of day, and 24 thereof showed a distinct pattern of pronounced down-regulation at 9:00 a.m. and comparatively weak up-regulation at later times. Eleven of these genes belonged to the species-specific kallikrein subfamily Klk1b. Responses in signature genes for thyroid hormone-induced responses were more frequent than for ionizing radiation, and trends persisted irrespective of time of day.

Conclusion

Diurnal variation of genome-wide transcriptional responses to 90 kBq ¹³¹I was demonstrated for the thyroid, kidney cortex and medulla, and liver, whereas variation was only marginal in the lungs and spleen. Overall, significant detection of potential biomarkers and signature genes was validated at each time of day, although direction of regulation and fold-change differed between morning, noon, and afternoon. These findings suggest that circadian rhythm should be considered in radiation research and that biological and analytical endpoints should be validated for circadian robustness.

Electronic supplementary material

The online version of this article (doi:10.1186/s13550-015-0150-y) contains supplementary material, which is available to authorized users.

Gene expression signature in mouse thyroid tissue after (131)I and (211)At exposure

BACKGROUND

(131)I and (211)At are used in nuclear medicine and accumulate in the thyroid gland and may impact normal thyroid function. The aim of this study was to determine transcriptional profile variations, assess the impact on cellular activity, and identify genes with biomarker properties in thyroid tissue after (131)I and (211)At administration in mice.

METHODS

To further investigate thyroid tissue transcriptional responses to (131)I and (211)At administration, we generated a new transcriptional dataset that includes re-evaluated raw intensity values from our previous (131)I and (211)At studies. Differential transcriptional profiles were identified by comparing treated and mock-treated samples using Nexus Expression 3.0 software. Further data analysis was performed using R/Bioconductor and IPA.

RESULTS

A total of 1144 genes were regulated. Hierarchical clustering subdivided the groups into two clusters containing the lowest and highest absorbed dose levels, respectively, and revealed similar transcriptional regulation patterns for many kallikrein-related genes. Twenty-seven of the 1144 genes were recurrently regulated after (131)I and (211)At exposure and divided into six clusters. Several signalling pathways were affected, including calcium, integrin-linked kinase, and thyroid cancer signalling, and the peroxisomal proliferator-activated receptor network.

CONCLUSIONS

Substantial changes in transcriptional regulation were shown in (131)I and (211)At-treated samples, and 27 genes were identified as potential biomarkers for (131)I and (211)At exposure. Clustering revealed distinct differences between transcriptional profiles of both similar and different exposures, demonstrating the necessity for better understanding of radiation-induced effects on cellular activity. Additionally, ionizing radiation-induced changes in kallikrein gene expression and identified canonical pathways should be further assessed.

Potential Biomarkers for Radiation-Induced Renal Toxicity following 177Lu-Octreotate Administration in Mice

The kidneys are one of the main dose-limiting organs in peptide receptor radionuclide therapy and due to large inter-individual variations in renal toxicity, biomarkers are urgently needed in order to optimize therapy and reduce renal tissue damage. The aim of this study was to investigate the transcriptional, functional, and morphological effects on renal tissue after ¹⁷⁷Lu-octreotate administration in normal mice, and to identify biomarkers for radiation induced renal toxicity.

Transcriptional Response in Mouse Thyroid Tissue after 211At Administration: Effects of Absorbed Dose, Initial Dose-Rate and Time after Administration

BACKGROUND

211At-labeled radiopharmaceuticals are potentially useful for tumor therapy. However, a limitation has been the preferential accumulation of released 211At in the thyroid gland, which is a critical organ for such therapy. The aim of this study was to determine the effect of absorbed dose, dose-rate, and time after 211At exposure on genome-wide transcriptional expression in mouse thyroid gland.

METHODS

BALB/c mice were i.v. injected with 1.7, 7.5 or 100 kBq 211At. Animals injected with 1.7 kBq were killed after 1, 6, or 168 h with mean thyroid absorbed doses of 0.023, 0.32, and 1.8 Gy, respectively. Animals injected with 7.5 and 100 kBq were killed after 6 and 1 h, respectively; mean thyroid absorbed dose was 1.4 Gy. Total RNA was extracted from pooled thyroids and the Illumina RNA microarray platform was used to determine mRNA levels. Differentially expressed transcripts and enriched GO terms were determined with adjusted p-value <0.01 and fold change >1.5, and p-value <0.05, respectively.

RESULTS

In total, 1232 differentially expressed transcripts were detected after 211At administration, demonstrating a profound effect on gene regulation. The number of regulated transcripts increased with higher initial dose-rate/absorbed dose at 1 or 6 h. However, the number of regulated transcripts decreased with mean absorbed dose/time after 1.7 kBq 211At administration. Furthermore, similar regulation profiles were seen for groups administered 1.7 kBq. Interestingly, few previously proposed radiation responsive genes were detected in the present study. Regulation of immunological processes were prevalent at 1, 6, and 168 h after 1.7 kBq administration (0.023, 0.32, 1.8 Gy).

Transcriptional response in normal mouse tissues after i.v. (211)At administration - response related to absorbed dose, dose rate, and time

Background

In cancer radiotherapy, knowledge of normal tissue responses and toxicity risks is essential in order to deliver the highest possible absorbed dose to the tumor while maintaining normal tissue exposure at non-critical levels. However, few studies have investigated normal tissue responses in vivo after ²¹¹At administration. In order to identify molecular biomarkers of ionizing radiation exposure, we investigated genome-wide transcriptional responses to (very) low mean absorbed doses from ²¹¹At in normal mouse tissues.

Methods

Female BALB/c nude mice were intravenously injected with 1.7 kBq ²¹¹At and killed after 1 h, 6 h, or 7 days or injected with 105 or 7.5 kBq and killed after 1 and 6 h, respectively. Controls were mock-treated. Total RNA was extracted from tissue samples of kidney cortex and medulla, liver, lungs, and spleen and subjected to microarray analysis. Enriched biological processes were categorized after cellular function based on Gene Ontology terms.

Results

Responses were tissue-specific with regard to the number of significantly regulated transcripts and associated cellular function. Dose rate effects on transcript regulation were observed with both direct and inverse trends. In several tissues, Angptl4, Per1 and Per2, and Tsc22d3 showed consistent transcript regulation at all exposure conditions.

Conclusions

This study demonstrated tissue-specific transcriptional responses and distinct dose rate effects after ²¹¹At administration. Transcript regulation of individual genes, as well as cellular responses inferred from enriched transcript data, may serve as biomarkers in vivo. These findings expand the knowledge base on normal tissue responses and may help to evaluate and limit side effects of radionuclide therapy.

Electronic supplementary material

The online version of this article (doi:10.1186/s13550-014-0078-7) contains supplementary material, which is available to authorized users.