Orthologs of human circulating miRNAs associated with hepatocellular carcinoma are elevated in mouse plasma months before tumour detection.

Research examining the potential for circulating miRNA to serve as markers for preneoplastic lesions or early-stage hepatocellular carcinoma (HCC) is hindered by the difficulties of obtaining samples from asymptomatic individuals. As a surrogate for human samples, we identified hub miRNAs in gene co-expression networks using HCC-bearing C3H mice. We confirmed 38 hub miRNAs as associated with HCC in F2 hybrid mice derived from radiogenic HCC susceptible and resistant founders. When compared to a panel of 12 circulating miRNAs associated with human HCC, two had no mouse ortholog and 7 of the remaining 10 miRNAs overlapped with the 38 mouse HCC hub miRNAs. Using small RNA sequencing data generated from serially collected plasma samples in F2 mice, we examined the temporal levels of these 7 circulating miRNAs and found that the levels of 4 human circulating markers, miR-122-5p, miR-100-5p, miR-34a-5p and miR-365-3p increased linearly as the time approaching HCC detection neared, suggesting a correlation of miRNA levels with oncogenic progression. Estimation of change points in the kinetics of the 4 circulating miRNAs suggested the changes started 17.5 to 6.8 months prior to HCC detection. These data establish these 4 circulating miRNAs as potential sentinels for preneoplastic lesions or early-stage HCC.

Transcriptomic analysis links hepatocellular carcinoma (HCC) in HZE ion irradiated mice to a human HCC subtype with favorable outcomes.

High-charge, high-energy ion particle (HZE) radiations are extraterrestrial in origin and characterized by high linear energy transfer (high-LET), which causes more severe cell damage than low-LET radiations like γ-rays or photons. High-LET radiation poses potential cancer risks for astronauts on deep space missions, but the studies of its carcinogenic effects have relied heavily on animal models. It remains uncertain whether such data are applicable to human disease. Here, we used genomics approaches to directly compare high-LET radiation-induced, low-LET radiation-induced and spontaneous hepatocellular carcinoma (HCC) in mice with a human HCC cohort from The Cancer Genome Atlas (TCGA). We identified common molecular pathways between mouse and human HCC and discovered a subset of orthologous genes (mR-HCC) that associated high-LET radiation-induced mouse HCC with a subgroup (mrHCC2) of the TCGA cohort. The mrHCC2 TCGA cohort was more enriched with tumor-suppressing immune cells and showed a better prognostic outcome than other patient subgroups.

Associations between lipids in selected brain regions, plasma miRNA, and behavioral and cognitive measures following 28Si ion irradiation.

The space radiation environment consists of multiple species of charged particles, including 28Si ions, that may impact brain function during and following missions. To develop biomarkers of the space radiation response, BALB/c and C3H female and male mice and their F2 hybrid progeny were irradiated with 28Si ions (350 MeV/n, 0.2 Gy) and tested for behavioral and cognitive performance 1, 6, and 12 months following irradiation. The plasma of the mice was collected for analysis of miRNA levels. Select pertinent brain regions were dissected for lipidomic analyses and analyses of levels of select biomarkers shown to be sensitive to effects of space radiation in previous studies. There were associations between lipids in select brain regions, plasma miRNA, and cognitive measures and behavioral following 28Si ion irradiation. Different but overlapping sets of miRNAs in plasma were found to be associated with cognitive measures and behavioral in sham and irradiated mice at the three time points. The radiation condition revealed pathways involved in neurodegenerative conditions and cancers. Levels of the dendritic marker MAP2 in the cortex were higher in irradiated than sham-irradiated mice at middle age, which might be part of a compensatory response. Relationships were also revealed with CD68 in miRNAs in an anatomical distinct fashion, suggesting that distinct miRNAs modulate neuroinflammation in different brain regions. The associations between lipids in selected brain regions, plasma miRNA, and behavioral and cognitive measures following 28Si ion irradiation could be used for the development of biomarker of the space radiation response.

EGFR Mutations Compromise Hypoxia-Associated Radiation Resistance through Impaired Replication Fork-Associated DNA Damage Repair.

EGFR signaling has been implicated in hypoxia-associated resistance to radiation or chemotherapy. Non-small cell lung carcinomas (NSCLC) with activating L858R or ΔE746-E750 EGFR mutations exhibit elevated EGFR activity and downstream signaling. Here, relative to wild-type (WT) EGFR, mutant (MT) EGFR expression significantly increases radiosensitivity in hypoxic cells. Gene expression profiling in human bronchial epithelial cells (HBEC) revealed that MT-EGFR expression elevated transcripts related to cell cycle and replication in aerobic and hypoxic conditions and downregulated RAD50, a critical component of nonhomologous end joining and homologous recombination DNA repair pathways. NSCLCs and HBEC with MT-EGFR revealed elevated basal and hypoxia-induced γ-H2AX-associated DNA lesions that were coincident with replication protein A in the S-phase nuclei. DNA fiber analysis showed that, relative to WT-EGFR, MT-EGFR NSCLCs harbored significantly higher levels of stalled replication forks and decreased fork velocities in aerobic and hypoxic conditions. EGFR blockade by cetuximab significantly increased radiosensitivity in hypoxic cells, recapitulating MT-EGFR expression and closely resembling synthetic lethality of PARP inhibition.Implications: This study demonstrates that within an altered DNA damage response of hypoxic NSCLC cells, mutant EGFR expression, or EGFR blockade by cetuximab exerts a synthetic lethality effect and significantly compromises radiation resistance in hypoxic tumor cells. Mol Cancer Res; 15(11); 1503-16. ©2017 AACR.

Original Research: A case-control genome-wide association study identifies genetic modifiers of fetal hemoglobin in sickle cell disease.

Sickle cell disease (SCD) is a group of inherited blood disorders that have in common a mutation in the sixth codon of the ß-globin (HBB) gene on chromosome 11. However, people with the same genetic mutation display a wide range of clinical phenotypes. Fetal hemoglobin (HbF) expression is an important genetic modifier of SCD complications leading to milder symptoms and improved long-term survival. Therefore, we performed a genome-wide association study (GWAS) using a case-control experimental design in 244 African Americans with SCD to discover genetic factors associated with HbF expression. The case group consisted of subjects with HbF≥8.6% (133 samples) and control group subjects with HbF≤£3.1% (111 samples). Our GWAS results replicated SNPs previously identified in an erythroid-specific enhancer region located in the second intron of the BCL11A gene associated with HbF expression. In addition, we identified SNPs in the SPARC, GJC1, EFTUD2 and JAZF1 genes as novel candidates associated with HbF levels. To gain insights into mechanisms of globin gene regulation in the HBB locus, linkage disequilibrium (LD) and haplotype analyses were conducted. We observed strong LD in the low HbF group in contrast to a loss of LD and greater number of haplotypes in the high HbF group. A search of known HBB locus regulatory elements identified SNPs 5' of δ-globin located in an HbF silencing region. In particular, SNP rs4910736 created a binding site for a known transcription repressor GFi1 which is a candidate protein for further investigation. Another HbF-associated SNP, rs2855122 in the cAMP response element upstream of Gγ-globin, was analyzed for functional relevance. Studies performed with siRNA-mediated CREB binding protein (CBP) knockdown in primary erythroid cells demonstrated γ-globin activation and HbF induction, supporting a repressor role for CBP. This study identifies possible molecular determinants of HbF production.

Elucidation of changes in molecular signalling leading to increased cellular transformation in oncogenically progressed human bronchial epithelial cells exposed to radiations of increasing LET.

The early transcriptional response and subsequent induction of anchorage-independent growth after exposure to particles of high Z and energy (HZE) as well as γ-rays were examined in human bronchial epithelial cells (HBEC3KT) immortalised without viral oncogenes and an isogenic variant cell line whose p53 expression was suppressed but that expressed an active mutant K-RAS(V12) (HBEC3KT-P53KRAS). Cell survival following irradiation showed that HBEC3KT-P53KRAS cells were more radioresistant than HBEC3KT cells irrespective of the radiation species. In addition, radiation enhanced the ability of the surviving HBEC3KT-P53RAS cells but not the surviving HBEC3KT cells to grow in anchorage-independent fashion (soft agar colony formation). HZE particle irradiation was far more efficient than γ-rays at rendering HBEC3KT-P53RAS cells permissive for soft agar growth. Gene expression profiles after radiation showed that the molecular response to radiation for HBEC3KT-P53RAS, similar to that for HBEC3KT cells, varies with radiation quality. Several pathways associated with anchorage independent growth, including the HIF-1α, mTOR, IGF-1, RhoA and ERK/MAPK pathways, were over-represented in the irradiated HBEC3KT-P53RAS cells compared to parental HBEC3KT cells. These results suggest that oncogenically progressed human lung epithelial cells are at greater risk for cellular transformation and carcinogenic risk after ionising radiation, but particularly so after HZE radiations. These results have implication for: (i) terrestrial radiation and suggests the possibility of enhanced carcinogenic risk from diagnostic CT screens used for early lung cancer detection; (ii) enhanced carcinogenic risk from heavy particles used in radiotherapy; and (iii) for space radiation, raising the possibility that astronauts harbouring epithelial regions of dysplasia or hyperplasia within the lung that contain oncogenic changes, may have a greater risk for lung cancers based upon their exposure to heavy particles present in the deep space environment.

Distinct transcriptome profiles identified in normal human bronchial epithelial cells after exposure to γ-rays and different elemental particles of high Z and energy.

BACKGROUND: Ionizing radiation composed of accelerated ions of high atomic number (Z) and energy (HZE) deposits energy and creates damage in cells in a discrete manner as compared to the random deposition of energy and damage seen with low energy radiations such as γ- or x-rays. Such radiations can be highly effective at cell killing, transformation, and oncogenesis, all of which are concerns for the manned space program and for the burgeoning field of HZE particle radiotherapy for cancer. Furthermore, there are differences in the extent to which cells or tissues respond to such exposures that may be unrelated to absorbed dose. Therefore, we asked whether the energy deposition patterns produced by different radiation types would cause different molecular responses. We performed transcriptome profiling using human bronchial epithelial cells (HBECs) after exposure to γ-rays and to two different HZE particles (28Si and 56Fe) with different energy transfer properties to characterize the molecular response to HZE particles and γ-rays as a function of dose, energy deposition pattern, and time post-irradiation. RESULTS: Clonogenic assay indicated that the relative biological effectiveness (RBE) for 56Fe was 3.91 and for 28Si was 1.38 at 34% cell survival. Unsupervised clustering analysis of gene expression segregated samples according to the radiation species followed by the time after irradiation, whereas dose was not a significant parameter for segregation of radiation response. While a subset of genes associated with p53-signaling, such as CDKN1A, TRIM22 and BTG2 showed very similar responses to all radiation qualities, distinct expression changes were associated with the different radiation species. Gene enrichment analysis categorized the differentially expressed genes into functional groups related to cell death and cell cycle regulation for all radiation types, while gene pathway analysis revealed that the pro-inflammatory Acute Phase Response Signaling was specifically induced after HZE particle irradiation. A 73 gene signature capable of predicting with 96% accuracy the radiation species to which cells were exposed, was developed. CONCLUSIONS: These data suggest that the molecular response to the radiation species used here is a function of the energy deposition characteristics of the radiation species. This novel molecular response to HZE particles may have implications for radiotherapy including particle selection for therapy and risk for second cancers, risk for cancers from diagnostic radiation exposures, as well as NASA's efforts to develop more accurate lung cancer risk estimates for astronaut safety. Lastly, irrespective of the source of radiation, the gene expression changes observed set the stage for functional studies of initiation or progression of radiation-induced lung carcinogenesis.

Progenitor cell line (hPheo1) derived from a human pheochromocytoma tumor.

BACKGROUND: Pheochromocytomas are rare tumors generally arising in the medullary region of the adrenal gland. These tumors release excessive epinephrine and norepinephrine resulting in hypertension and cardiovascular crises for which surgery is the only definitive treatment. Molecular mechanisms that control tumor development and hormone production are poorly understood, and progress has been hampered by the lack of human cellular model systems. To study pheochromocytomas, we developed a stable progenitor pheochromocytoma cell line derived from a primary human tumor. METHODS: After IRB approval and written informed consent, human pheochromocytoma tissue was excised, minced, dispersed enzymatically, and cultured in vitro. Primary pheochromocytoma cells were infected with a lentivirus vector carrying the catalytic subunit of human telomerase reverse transcriptase (hTERT). The hTERT immortalized cells (hPheo1) have been passaged >300 population doublings. The resulting cell line was characterized morphologically, biochemically and for expression of neuroendocrine properties. The expression of marker enzymes and proteins was assessed by immunofluorescence staining and immunoblotting. Telomerase activity was determined by using the telomeric repeat amplification protocol (TRAP) assay. RESULTS: We have established a human pheochromocytoma precursor cell line that expresses the neuroendocrine marker, chromogranin A, when differentiated in the presence of bone morphogenic protein 4 (BMP4), nerve growth factor (NGF), and dexamethasone. Phenylethanolamine N-methyltransferase (PNMT) expression is also detected with this differentiation regimen. CD-56 (also known as NCAM, neural cell adhesion molecule) is expressed in these cells, but CD31 (also known as PECAM-1, a marker of endothelial cells) is negative. CONCLUSIONS: We have maintained hTERT-immortalized progenitor cells derived from a pheochromocytoma (hPheo1) in culture for over 300 population doublings. This progenitor human cell line is normal diploid except for a deletion in the p16 region and has inducible neuroendocrine biomarkers. These cells should be a valuable reagent for studying mechanisms of tumor development and for testing novel therapeutic approaches.

Molecular characterisation of murine acute myeloid leukaemia induced by 56Fe ion and 137Cs gamma ray irradiation.

Exposure to sparsely ionising gamma- or X-ray irradiation is known to increase the risk of leukaemia in humans. However, heavy ion radiotherapy and extended space exploration will expose humans to densely ionising high linear energy transfer (LET) radiation for which there is currently no understanding of leukaemia risk. Murine models have implicated chromosomal deletion that includes the hematopoietic transcription factor gene, PU.1 (Sfpi1), and point mutation of the second PU.1 allele as the primary cause of low-LET radiation-induced murine acute myeloid leukaemia (rAML). Using array comparative genomic hybridisation, fluorescence in situ hybridisation and high resolution melt analysis, we have confirmed that biallelic PU.1 mutations are common in low-LET rAML, occurring in 88% of samples. Biallelic PU.1 mutations were also detected in the majority of high-LET rAML samples. Microsatellite instability was identified in 42% of all rAML samples, and 89% of samples carried increased microsatellite mutant frequencies at the single-cell level, indicative of ongoing instability. Instability was also observed cytogenetically as a 2-fold increase in chromatid-type aberrations. These data highlight the similarities in molecular characteristics of high-LET and low-LET rAML and confirm the presence of ongoing chromosomal and microsatellite instability in murine rAML.

Defining molecular and cellular responses after low and high linear energy transfer radiations to develop biomarkers of carcinogenic risk or therapeutic outcome.

The variability in radiosensitivity across the human population is governed in part by genetic factors. The ability to predict therapeutic response, identify individuals at greatest risk for adverse clinical responses after therapeutic radiation doses, or identify individuals at high risk for carcinogenesis from environmental or medical radiation exposures has a medical and economic impact on both the individual and society at large. As radiotherapy incorporates particles, particularly particles larger than protons, into therapy, the need for such discriminators, (i.e., biomarkers) will become ever more important. Cellular assays for survival, DNA repair, or chromatid/chromosomal analysis have been used to identify at-risk individuals, but they are not clinically applicable. Newer approaches, such as genome-wide analysis of gene expression or single nucleotide polymorphisms and small copy number variations within chromosomes, are examples of technologies being applied to the discovery process. Gene expression analysis of primary or immortalized human cells suggests that there are distinct gene expression patterns associated with radiation exposure to both low and high linear energy transfer radiations and that those most radiosensitive are discernible by their basal gene expression patterns. However, because the genetic alterations that drive radio response may be subtle and cumulative, the need for large sample sizes of specific cell or tissue types is required. A systems biology approach will ultimately be necessary. Potential biomarkers from cell lines or animal models will require validation in a human setting where possible and before being considered as a credible biomarker some understanding of the molecular mechanism is necessary.

Irreparable complex DNA double-strand breaks induce chromosome breakage in organotypic three-dimensional human lung epithelial cell culture.

DNA damage and consequent mutations initiate the multistep carcinogenic process. Differentiated cells have a reduced capacity to repair DNA lesions, but the biological impact of unrepaired DNA lesions in differentiated lung epithelial cells is unclear. Here, we used a novel organotypic human lung three-dimensional (3D) model to investigate the biological significance of unrepaired DNA lesions in differentiated lung epithelial cells. We showed, consistent with existing notions that the kinetics of loss of simple double-strand breaks (DSBs) were significantly reduced in organotypic 3D culture compared to kinetics of repair in two-dimensional (2D) culture. Strikingly, we found that, unlike simple DSBs, a majority of complex DNA lesions were irreparable in organotypic 3D culture. Levels of expression of multiple DNA damage repair pathway genes were significantly reduced in the organotypic 3D culture compared with those in 2D culture providing molecular evidence for the defective DNA damage repair in organotypic culture. Further, when differentiated cells with unrepaired DNA lesions re-entered the cell cycle, they manifested a spectrum of gross-chromosomal aberrations in mitosis. Our data suggest that downregulation of multiple DNA repair pathway genes in differentiated cells renders them vulnerable to DSBs, promoting genome instability that may lead to carcinogenesis.

Loss of p15/Ink4b accompanies tumorigenesis triggered by complex DNA double-strand breaks.

DNA double-strand breaks (DSBs) are the most deleterious lesion inflicted by ionizing radiation. Although DSBs are potentially carcinogenic, it is not clear whether complex DSBs that are refractory to repair are more potently tumorigenic compared with simple breaks that can be rapidly repaired, correctly or incorrectly, by mammalian cells. We previously demonstrated that complex DSBs induced by high-linear energy transfer (LET) Fe ions are repaired slowly and incompletely, whereas those induced by low-LET gamma rays are repaired efficiently by mammalian cells. To determine whether Fe-induced DSBs are more potently tumorigenic than gamma ray-induced breaks, we irradiated 'sensitized' murine astrocytes that were deficient in Ink4a and Arf tumor suppressors and injected the surviving cells subcutaneously into nude mice. Using this model system, we find that Fe ions are potently tumorigenic, generating tumors with significantly higher frequency and shorter latency compared with tumors generated by gamma rays. Tumor formation by Fe-irradiated cells is accompanied by rampant genomic instability and multiple genomic changes, the most interesting of which is loss of the p15/Ink4b tumor suppressor due to deletion of a chromosomal region harboring the CDKN2A and CDKN2B loci. The additional loss of p15/Ink4b in tumors derived from cells that are already deficient in p16/Ink4a bolsters the hypothesis that p15 plays an important role in tumor suppression, especially in the absence of p16. Indeed, we find that reexpression of p15 in tumor-derived cells significantly attenuates the tumorigenic potential of these cells, indicating that p15 loss may be a critical event in tumorigenesis triggered by complex DSBs.

SNS-032 prevents hypoxia-mediated glioblastoma cell invasion by inhibiting hypoxia inducible factor-1alpha expression.

Hypoxia and hypoxia inducible factor-1alpha (HIF-1alpha) play a critical role in glioblastoma (GBM) which is characterized by highly aggressive and widespread cell invasion into adjacent normal brain tissue. The purpose of this study was to investigate the effect of the novel aminothiazole com-pound SNS-032 in glioblastoma cell invasion under hypoxic condition. SNS-032 is a potent and selective inhibitor of cyclin-dependent kinases 2, 7 and 9 and inhibits both cell cycle and transcription. We analyzed the effect of SNS-032 (0.5 microM) on HIF-1alpha expression and its major trans-regulating factors including COX-2, VEGF, MMP-2 and uPAR that are involved in cellular invasion in tumor hypoxia. Our observations demonstrate SNS-032: i) inhibited hypoxia-induced U87MG cell invasion and among all the other inhibitors tested, SNS-032 is the most effective, ii) blocked HIF-1alpha mediated transcription of COX-2, MMP-2, VEGF and uPAR expression in U87MG cells in response to hypoxia, iii) blocked HIF-1alpha expression by a proteasome independent pathway. The effects were similar to those observed with HIF-1alpha siRNA which prevented cellular invasion by blocking HIF-1alpha expression and its downstream effectors. Taken together, our data suggest that SNS-032 prevents hypoxia-mediated U87MG cell invasion by blocking the expression of HIF-1alpha and its trans-regulating factors. Our results present an opportunity in controlling highly invasive tumors such as glioblastoma using this novel class of compounds.

Enhanced identification and biological validation of differential gene expression via Illumina whole-genome expression arrays through the use of the model-based background correction methodology.

Despite the tremendous growth of microarray usage in scientific studies, there is a lack of standards for background correction methodologies, especially in single-color microarray platforms. Traditional background subtraction methods often generate negative signals and thus cause large amounts of data loss. Hence, some researchers prefer to avoid background corrections, which typically result in the underestimation of differential expression. Here, by utilizing nonspecific negative control features integrated into Illumina whole genome expression arrays, we have developed a method of model-based background correction for BeadArrays (MBCB). We compared the MBCB with a method adapted from the Affymetrix robust multi-array analysis algorithm and with no background subtraction, using a mouse acute myeloid leukemia (AML) dataset. We demonstrated that differential expression ratios obtained by using the MBCB had the best correlation with quantitative RT-PCR. MBCB also achieved better sensitivity in detecting differentially expressed genes with biological significance. For example, we demonstrated that the differential regulation of Tnfr2, Ikk and NF-kappaB, the death receptor pathway, in the AML samples, could only be detected by using data after MBCB implementation. We conclude that MBCB is a robust background correction method that will lead to more precise determination of gene expression and better biological interpretation of Illumina BeadArray data.

Genome wide expression analysis of the CCR4-NOT complex indicates that it consists of three modules with the NOT module controlling SAGA-responsive genes.

Of the nine known members of the CCR4-NOT complex, CCR4/CAF1 are most important in mRNA deadenylation whereas the NOT1-5 proteins are most critical for transcriptional repression. Whole genome microarray analysis using deletions in seven of the CCR4-NOT genes was used to determine the overall mRNA expression patterns that are affected by members of the yeast CCR4-NOT complex. Under glucose conditions, ccr4 and caf1 displayed a high degree of similarity in the manner that they affected gene expression. In contrast, the not deletions were similar in the way they affected genes, but showed no correlation with that of ccr4/caf1. A number of groups of functionally related proteins were specifically controlled by the CCR4/CAF1 or NOT modules. Importantly, the NOT proteins preferentially affected SAGA-controlled gene expression. Also, both the CCR4/CAF1 and NOT group of proteins shared much greater similarities in their effects on gene expression during the stress of glucose deprivation. BTT1, a member of the nascent polypeptide association complex that binds the ribosome, was shown to be a tenth member of the CCR4-NOT complex, interacting through CAF130. Microarray analysis indicated that BTT1 and CAF130 correlate very highly in their control of gene expression and preferentially repress genes involved in ribosome biogenesis. These results indicate that distinct portions of the CCR4-NOT complex control a number of different cellular processes.

Gene expression changes in normal human skin fibroblasts induced by HZE-particle radiation.

Studies have shown that radiation exposure affects global gene expression in mammalian cells. However, little is known about the effects of HZE particles on gene expression. To study these effects, human skin fibroblasts were irradiated with HZE particles of different energies and LETs. The data obtained from these experiments indicate that changes in gene expression are dependent on the energy of the radiation source. Particles with the highest energy, i.e. iron, induced the biggest expression changes in terms of numbers of genes and magnitudes of changes. Many genes were found to undergo significant expression changes after HZE-particle irradiation, including CDKN1A/p21, MDM2, TNFRSF6/fas, PCNA and RAD52. Unlike X rays, HZE particles expose cells to two types of radiation: primary ions and delta rays. We hypothesized that the biological effects of delta rays, which are secondary electron emissions, should resemble the effects of X rays. To explore this idea, gene expression changes between cells that had been irradiated with HZE particles and X rays were compared. The results support our hypothesis since the number of genes that commonly changed after exposure to both radiations increased as a function of particle energy.

Gene expression profiles of normal human fibroblasts after exposure to ionizing radiation: a comparative study of low and high doses.

Several types of cellular responses to ionizing radiation, such as the adaptive response or the bystander effect, suggest that low-dose radiation may possess characteristics that distinguish it from its high-dose counterpart. Accumulated evidence also implies that the biological effects of low-dose and high-dose ionizing radiation are not linearly distributed. We have investigated, for the first time, global gene expression changes induced by ionizing radiation at doses as low as 2 cGy and have compared this to expression changes at 4 Gy. We applied cDNA microarray analyses to G1-arrested normal human skin fibroblasts subjected to X irradiation. Our data suggest that both qualitative and quantitative differences exist between gene expression profiles induced by 2 cGy and 4 Gy. The predominant functional groups responding to low-dose radiation are those involved in cell-cell signaling, signal transduction, development and DNA damage responses. At high dose, the responding genes are involved in apoptosis and cell proliferation. Interestingly, several genes, such as cytoskeleton components ANLN and KRT15 and cell-cell signaling genes GRAP2 and GPR51, were found to respond to low-dose radiation but not to high-dose radiation. Pathways that are specifically activated by low-dose radiation were also evident. These quantitative and qualitative differences in gene expression changes may help explain the non-linear correlation of biological effects of ionizing radiation from low dose to high dose.