ionizing radiation is a well - known carcinogen in humans , and breast is one of the most sensitive organs to radiogenic cancer .
the rate of breast cancer in postwar japan was among the lowest in the world , but breast cancer contributed a disproportionately large fraction of the radiation - related cancer burden in atomic bomb survivors [ 2 , 3 ] .
the data from the hiroshima and nagasaki survivors provides strong evidence for increased breast cancer following single acute doses of 20  cgy and linearity with increasing dose [ 36 ] .
also , an increase in the incidence of breast cancer has been observed in areas affected by the chernobyl accident , which resulted in radioactive contamination of large areas of belarus and ukraine . a twofold increase in risk
was observed when comparing the most ( > 40  msv cumulative dose ) and least contaminated regions .
interestingly , the increase appeared 10 years after exposure and was most prominent in women exposed at younger age .
more than 50,000 women in the united states have been treated with chest radiation ( 20  gy ) for a pediatric or young adult cancer .
children treated from cancer with radiotherapy have a 2.9 relative risk of subsequent malignancy compared to those who were not [ 8 , 9 ] .
a systematic review of 14 studies concluded that risk of breast cancer increased as early as 8 years following chest radiation and did not plateau with increasing length of follow - up .
studies estimating low - dose radiation - induced cancer risk from diagnostic x - rays and ct scans have found a small but significant increased lifetime risk [ 11 , 12 ] . while the benefits of diagnostic x - ray and ct scans outweigh potential individual lifetime risk , their use should be justified and alternatives considered .
we know remarkably little of molecular mechanisms that may be protective or risky for breast cancer after exposure to low - dose ionizing radiation ( ldir ) .
identification of transcriptomic changes induced by ldir in mammary tissue will be valuable to elucidate the molecular mechanisms associated with radiation - induced breast cancer .
long noncoding rnas ( lncrnas ) , which initially were thought of as transcriptional noise , are emerging as key regulators of a multitude of cellular processes by taking part in epigenetic , transcriptional , and posttranscriptional regulation of gene expression [ 13 , 14 ] .
the lncrnas have a weaker evolutionary constraint and lower levels of expression compared to the protein - coding transcripts [ 15 , 16 ] but exhibit more tissue specific expression than the protein - coding genes .
recently , a number of studies have shown that lncrna expression can be deregulated in human cancers [ 17 , 18 ] . as the functions of individual lncrnas in cancer are beginning to be elucidated ,
they are being categorized and referred to as either tumor suppressor or oncogenic lncrnas , in the same way as traditional protein - coding cancer genes .
however , the relevance of lncrnas in the cell and tissue response to ionizing radiation has not yet been characterized . in this study , we used agilent sureprint g3 microarrays to profile lncrna and mrna from mammary glands of balb / c mice 2 , 4 , and 8 weeks after irradiation and of spret / eij mice 4 weeks after irradiation with 10  cgy of x - radiation .
we identified lncrna and mrna expression signatures for each time point after irradiation in comparison to sham .
of the total 1338 lncrnas identified to be differentially expressed after ldir in either balb / c or spret / eij , 1337 had a significantly correlated expression pattern with at least one mrna that was also differentially expressed after ldir .
our results indicate lncrnas may exert a partial or key role in the regulation of coding rna expression induced by radiation .
balb / c and spret / eij mice were purchased from jackson laboratory , housed four per cage under a 12  hr light and 12  hr dark cycle , and fed with lab diet 5008 chow and water ad libitum .
the mice were irradiated whole body at 8 - 9 weeks of age to a single dose of 10  cgy using a precision x - ray inc rad320 320  kvp x - ray machine , operated at 300  kv , 2  ma .
mammary tissues were collected for gene expression profile at 2 , 4 , and 8 weeks after irradiation .
all animal experiments were performed at lawrence berkeley national laboratory and the study was carried out in strict accordance with the guide for the care and use of laboratory animals of the national institutes of health .
the animal use protocol was approved by the animal welfare and research committee of the lawrence berkeley national laboratory .
agilent sureprint g3 mouse ge 8x60k microarrays were used according to the manufacturer 's protocol ( arrays contained 39,430 entrez gene rnas and 16,251 lncrnas ) .
genes that were differentially expressed between sham and irradiated were identified by the unpaired student 's t - test with a p value cut - off of 0.05 ( p value < 0.001 for baseline strain comparison ) and a fold - change criteria of more than 1.5 . for each of the 8 experimental treatment groups , the average expression values of the 3 biological replicates
significantly correlated pairs of mrna and lncrna were calculated using a standard permutation test . in brief , for each potential mrna and lncrna pair , the 8 mrna values were randomly rearranged and a correlation coefficient was calculated between the 8 mrnas and 8 lncrnas values .
the p value reported in this work represents the percentage of permutations leading to a higher correlation than the original correlation between the 8 mrnas and 8 lncrnas values . in other words ,
the lower the p value is , the more likely the mrna and lncrna pair is not randomly associated .
gene lists were annotated with biological functions using ingenuity pathway analysis ( ipa ) , kegg pathway analysis ( http://bioinfo.vanderbilt.edu/webgestalt/ ) and david go gene ontology ( http://david.abcc.ncifcrf.gov/ ; p  0.05 ) .
annotations for the various shapes used in the ipa networks in figures 13 are shown in figure s3 in supplementary material available online at http://dx.doi.org/10.1155/2015/461038 .
to identify potential lncrnas and mrnas that may determine susceptibility to radiation - induced breast cancer , we profiled two inbred strains of mice with differing genetic susceptibilities : balb / c mice as more sensitive and spret / eij as more resistant .
balb / c mice carry two dna - pkcs polymorphisms with reduced catalytic subunit activity and defective nonhomologous - end - joining of double strand breaks .
we identified 195 lncrnas as upregulated and 95 lncrnas as downregulated in balb / c in comparison to spret / eij ( fold - change 1.5 ; p value < 0.001 ) ( figure 1(a ) ; table s1 ) .
additionally , 582 mrnas were upregulated and 402 mrnas were downregulated in balb / c in comparison to spret / eij ( fold - change 1.5 ; p value < 0.001 ) ( figure 1(b ) ; table s1 ) .
gene ontology analyses of differentially expressed genes between balb / c and spret / eij showed significant enrichment for metabolic processes ( p = 5.9e  06 ) , ion binding ( p = 3.00e  08 ) , and chemokine signaling ( p = 0.02 ) ( figure 1(c ) ; table s2 ) .
our analyses identified significant strain differences in gene expression between balb / c and spret / eij mammary tissues .
we found significant differences in the expression of a number of chemokines including cxcl10 , ccl6 , and ccl25 , which were expressed at higher levels in mammary tissues of the more sensitive and susceptible balb / c mice and which have previously been associated with breast cancer progression when over expressed . to identify lncrna and mrna expression changes induced by low - dose ionizing radiation ( ldir ) , we profiled lncrna and mrna expression from mammary glands of balb / c mice 2 , 4 , and 8 weeks after irradiation and of spret / eij mice 4 weeks after irradiation with 10  cgy of x - radiation . for balb / c mice ,
a total of 357 , 480 , and 335 lncrnas and 550 , 911 , and 389 coding rnas were identified to be differentially expressed at weeks 2 , 4 , and 8 after ir in comparison to sham ( fold - change 1.5 ; p value < 0.05 ) , respectively ( figure 2(a ) ; table s1 ) . for spret
/ eij , a total of 327 lncrnas and 424 mrnas were identified as differentially expressed at week 4 after irradiation in comparison to sham ( fold - change 1.5 ; p value < 0.05 ) ( figure 3(a ) ; table s1 ) .
few coding - rnas and lncrnas were found to be differentially expressed at different time points ( figure s1(a ) and s1(b ) ) and between balb / c and spret / eij ( figure s1(c ) ) . to determine the biological functions associated with the ldir response , we excluded genes whose levels fluctuate in the mouse mammary gland across the estrous cycle .
we recently mapped transcript - level changes across the estrous cycle in the murine mammary gland using rna sequencing and defined a comprehensive estrous variable gene signature of 3893 genes whose levels fluctuate in mammary glands of balb / c mice . comparison of our mapped ldir genes in mammary glands of balb / c and spret / eij mice with the estrous signature revealed an approximate 20% overlap in balb / c ( figure 2(b ) ) and 9% overlap in spret / eij mice ( figure 3(b ) ) .
nonoverlapping and differentially expressed ldir genes for each of the time points were then computationally mapped to biological functions , pathways , upstream regulators , and networks .
these analyses suggested that the ldir response signatures in mammary glands of balb / c mice transitions between time points and is distinct from the ldir response in spret / eij mice .
two weeks after ldir exposure pathways and biological functions significantly enriched in mammary glands of balb / c mice compared to sham irradiated mice included chemokine signaling ( p = 0.01 ) , ccr3 signaling in eosinophils ( p = 0.05 ) , cellular movement ( 3.92e  04
< p < 3.41e  02 ) , and cell death and survival ( 1.29e  03
gene interaction networks were enriched for tissue and endocrine system injury ( figure 2(c ) top panel ) and significant predicted upstream regulators ( table s3 ) include gli2 ( p = 5.29e  04 ) and satb1 ( p = 1.98e  03 ) . similar to the two - week ldir response , gli2 was predicted to be an upstream regulator ( p = 4.52e  03 ; table s3 ) .
gata3 and stat6 were among other significant upstream regulators associated with the four - week low - dose response ( table s3 ) .
we furthermore observed that the mammary gland of balb / c mice four weeks after ldir was enriched for inflammatory response genes ( 1.84e  05
< p < 1.45e  02 ) , cell - cell signaling ( 2.79e  05
< p < 1.60e  02 ) , while gene interaction networks were enriched for lipid metabolism ( figure 2(c ) middle panel ) .
interestingly , similar responses were observed in mammary glands of spret mice at 4  weeks after ldir including inflammatory response functions ( 1.69e  03
< p < 4.49e  02 ) , cell - cell signaling ( 5.54e  04
< p < 4.49e  02 ) , suggesting that the functional response is similar across strains and is independent of the gene transcript response . at 8 weeks after ldir we observed downregulation of genes involved in mammary gland development including progesterone receptor , prolactin , amphiregulin , and wnt4 ( figure 2(c ) bottom panel ; table s3 ) . to identify lncrnas potentially regulating the expression of coding rnas in response to radiation ,
correlation coefficients on expression data were calculated for each coding rna and lncrna that were identified as differentially expressed after ldir .
a permutation - based algorithm was then used to determine which correlations were statistically significant ( p < 0.05 ; table s4 ) .
we observed that nearly all ldir modulated lncrnas were correlated with at least one of the ldir modulated coding mrnas ( figure 4(a ) ) . to determine whether these correlations were driven by estrous variations , we only considered genes whose expression levels were not overlapping with our previously determined estrous signature ( figure 2(b ) ) . again , we observed that nearly all differentially expressed lncrnas were correlated with at least one differentially expressed mrna suggesting that estrous cycling does not affect the strong correlation between lncrna and mrna expression after ldir . to test the robustness of these correlations , we compared the number of lncrnas associated with at least one mrna at three different p values ( p < 5e  02 , p < 5e  03 , and p < 5e  04 ) . at p
< 5e  02 or p < 5e  03 , nearly all ( 97100% ) differentially expressed lncrnas were found to be correlated with at least one mrna ( table s5 ) . at p
< 5e  04 , corresponding to a correlation coefficient > 0.9 , we still observed a significant fraction ( 6381% ; table s5 ) of lncrnas correlated with mrnas .
representative correlation networks of lncrnas are shown in figure 4(b ) ( p < 5e  03 ) for each of the timepoints .
we furthermore observed that the same lncrna correlates with different gene sets across different time points ( figure s2 ) .
taken together these data show that ldir induces coordinated changes in lncrna and mrnas and suggests a critical role for lncrnas in mediating the low - dose radiation response .
in this study , we demonstrate that genetic background strongly influences the expression of lncrnas and their response to low - dose radiation by transcriptomic analysis of mouse mammary glands using microarrays that contain both lncrnas and coding rnas .
we have identified a number of lncrnas that are significantly changed after exposure to ldir at three different timepoints after radiation exposure .
moreover , the changes in the expression of lncrnas are significantly correlated with the expression of coding rnas , suggesting that lncrnas may coordinate the tissue response to radiation via regulation of coding mrnas .
however , the specific regulatory mechanism of this control requires further investigation , and knock - out and overexpression of the lncrna genes in mice and other model systems should be performed to increase our understanding of the regulatory mechanisms in response to ldir .
 long noncoding rnas ( lncrnas ) are emerging as key regulators of diverse cell functions and processes . 
 however , the relevance of lncrnas in the cell and tissue response to ionizing radiation has not yet been characterized . 
 here we used microarray profiling to determine lncrna and mrna expression in mammary glands of balb / c and spret / eij mice after low - dose ionizing radiation ( ldir ) exposure . 
 we found that unirradiated mammary tissues of these strains differed significantly in baseline expressions of 290 lncrnas . 
 ldir exposure ( 10  cgy ) induced a significant change in the expression of many lncrnas . 
 the vast majority of lncrnas identified to be differentially expressed after ldir in either balb / c or spret / eij had a significantly correlated expression pattern with at least one ldir responsive mrna . 
 functional analysis revealed that the response to ldir in balb / c mice is highly dynamic with enrichment for genes involved in tissue injury , inflammatory responses , and mammary gland development at 2 , 4 , and 8 weeks after ldir , respectively . 
 our study demonstrates that genetic background strongly influences the expression of lncrnas and their response to radiation and that lncrnas may coordinate the tissue response to ldir exposure via regulation of coding mrnas . 
1. Introduction
2. Materials and Methods
3. Results and Discussion
4. Conclusions