ORIGINAL RESEARCH article

Front. Genet., 02 August 2021

Sec. Predictive Toxicology and Exposomics

Volume 12 - 2021 | https://doi.org/10.3389/fgene.2021.710143

Low-Dose Radiation Can Cause Epigenetic Alterations Associated With Impairments in Both Male and Female Reproductive Cells

  • 1. Hong Kong Branch of the Southern Marine Science and Engineering Guangdong Laboratory, Hong Kong, China

  • 2. Department of Chemistry, City University of Hong Kong, Kowloon Tong, Hong Kong

  • 3. School of Biological Sciences, The University of Hong Kong, Pok Fu Lam, Hong Kong

  • 4. State Key Laboratory of Marine Pollution, City University of Hong Kong, Kowloon Tong, Hong Kong

  • 5. Department of Physics, City University of Hong Kong, Kowloon Tong, Hong Kong

  • 6. School of Life Sciences, Hong Kong Bioinformatics Centre, The Chinese University of Hong Kong, Shatin, Hong Kong

  • 7. Department of Biology, Hong Kong Baptist University, Kowloon Tsai, Hong Kong

  • 8. Laboratory of Environmental Pollution and Integrative Omics, Guilin Medical University, Guilin, China

  • 9. Department of Science and Environmental Studies, The Education University of Hong Kong, Tai Po, Hong Kong

Article metrics

View details

22

Citations

7,8k

Views

2,3k

Downloads

Abstract

Humans are regularly and continuously exposed to ionizing radiation from both natural and artificial sources. Cumulating evidence shows adverse effects of ionizing radiation on both male and female reproductive systems, including reduction of testis weight and sperm count and reduction of female germ cells and premature ovarian failure. While most of the observed effects were caused by DNA damage and disturbance of DNA repairment, ionizing radiation may also alter DNA methylation, histone, and chromatin modification, leading to epigenetic changes and transgenerational effects. However, the molecular mechanisms underlying the epigenetic changes and transgenerational reproductive impairment induced by low-dose radiation remain largely unknown. In this study, two different types of human ovarian cells and two different types of testicular cells were exposed to low dose of ionizing radiation, followed by bioinformatics analysis (including gene ontology functional analysis and Ingenuity Pathway Analysis), to unravel and compare epigenetic effects and pathway changes in male and female reproductive cells induced by ionizing radiation. Our findings showed that the radiation could alter the expression of gene cluster related to DNA damage responses through the control of MYC. Furthermore, ionizing radiation could lead to gender-specific reproductive impairment through deregulation of different gene networks. More importantly, the observed epigenetic modifications induced by ionizing radiation are mediated through the alteration of chromatin remodeling and telomere function. This study, for the first time, demonstrated that ionizing radiation may alter the epigenome of germ cells, leading to transgenerational reproductive impairments, and correspondingly call for research in this new emerging area which remains almost unknown.

Introduction

Humans are regularly and continuously exposed to ionizing radiation from both natural and artificial sources. Natural background radiation (NBR) may come from radon and its daughter products, crust, cosmic rays, soil, water, food etc. According to the United States report in 2008, the average world population exposure to natural background radiation was about 2.4 millisievert (mSv) per annum, whereas additional radiation contributed from medical diagnosis is estimated at about 0.6 mSv annually (UNSCEAR, 2010). Radiation injuries, such as epithelial and stromal lesion, vascular lesions, fibrosis, and neoplasia may occur upon irradiation (Fajardo, 2005), and level of injuries mainly depends on the radiation dose, duration, and cell cycle stage. In general, male and female gametes are more susceptible to radiation which may affect the composition and biological integrity of their proteins, lipids, and nucleic acids, and hence compromising their capacity to produce normal embryos (García-Rodríguez et al., 2018). Indeed, the adverse effects of ionizing radiation on the reproductive system of mammals (including human being) have been clearly demonstrated by numerous studies. Radiation, for example, has been shown to reduce testis weight and sperm count in male rodent (Haines et al., 2002). In human, a radiation dose as low as 1 Gy can disrupt gonadotropin levels, reduce the number of spermatocytes and spermatids, while a higher radiation dose of 10 Gy can kill all spermatogonial stem cells, leading to permanent azoospermia (Rowley et al., 1974; Clifton and Bremner, 1983). In vivo studies in female rats demonstrated that radiation can reduce a number of germ cells and increase synaptonemal complex (SC) fragmentation (Pujol et al., 1996, 1997). Female cancer patients after irradiation treatments often experience premature ovarian failure, infertility, uterine damage, premature deliveries, and miscarriage (Critchley et al., 1992; Wallace et al., 2003).

At cellular level, ionizing radiation can induce DNA damage: (i) directly, via energy deposition to DNA molecules; and (ii) indirectly, via the attack from reactive oxygen species generated by other ionized molecules, such as free radicals (Ward, 1988). DNA breaks, if failed to repair, can increase the risk of mutagenesis and carcinogenesis (Goodhead, 1994; Ward, 1995; Little, 2000). Haploid germ cells are more susceptible to DNA mutation induced by radiation due to the absence of sister chromatid. DNA repair is thus performed by non-homologous end joining, which is generally accepted as an error-prone pathway (Ahmed et al., 2015; Wdowiak et al., 2019).

Apart from DNA damage, emerging evidence showed that ionizing radiation may further alter DNA methylation pattern (Miousse et al., 2017), modify histone and chromatin structure (Rogakou et al., 1998; Pogribny et al., 2005; Kumar R. et al., 2013), and miRNA expression levels (Halimi et al., 2012; Metheetrairut and Slack, 2013). These induced changes in paternal epigenome can be passed on and affect fertilization and embryogenesis in subsequent generations (Jenkins and Carrell, 2012), thereby causing transgenerational effects. Non-coding RNA is suggested to play a role in epigenetic inheritance by the trafficking of microRNA-containing vesicles to sperm via blood stream (Paris et al., 2015; Grandjean et al., 2016; Sharma et al., 2018; Szatmári et al., 2019). In zebrafish, it has been shown that differential methylation of specific genes in F0 caused by parental radiation exposure can be transmitted to their F3 offspring, which had never been exposed to radiation before throughout their life cycle (Kamstra et al., 2018). Nevertheless, existing knowledge on radiation-induced epigenetics alterations on reproductive system, especially female, are limited. Irradiated (2.5 Gy) mice showed an upregulation of miR-29 family in male germline, resulting in reduced expression of de novo DNA methyltransferase 3a and hypomethylation of interspersed nuclear elements associated to chromatin modification (Filkowski et al., 2010). Increased phosphorylation of histone H2AX (γ-H2AX) was observed in immature spermatozoa of cranially exposed (20 Gy) rats via bystander effect. Significant reduction in global DNA methylation was also found in testes and mature sperm cells (Tamminga et al., 2008). On the contrary, significant increase in global methylation was found in spermatozoa of a human exposed to occupational radiation (Kumar D. et al., 2013). Nevertheless, the epigenetic and transgenerational effects induced by radiation on the reproductive system are still not fully understood, and the underlying mechanisms remained obscure and need further studies (Skrzypek et al., 2019; Dubrova and Sarapultseva, 2020).

Using two different types of human ovarian cells and two different types of testicular cells as models, we carried out a series of in vitro experiments followed by bioinformatics analysis (including gene ontology functional analysis and Ingenuity Pathway Analysis), to unravel and compare epigenetic effects and pathway changes in male and female reproductive cells caused by environmentally relevant dose of ionizing radiation. Specifically, we hypothesize that: (a) ionizing radiation could alter the expression of gene cluster related to DNA damage response, and (b) reproductive impairment caused by ionizing radiation is gender specific and controlled by different gene networks.

Materials and Methods

Ovarian and Testicular Cell Culture and Ionizing Radiation Exposure

Two human ovarian cancer cells (SKOV3 and COV434) and two mouse testicular germ cells (GC-1 and TM4) were cultured under the conditions described in Supplementary Table 1. Only human and mouse cell lines (not primary cells) were used, and these were purchased from an international company. In accordance with the national legislation and the institutional requirements, the Human Research Ethics Committee of The University of Hong Kong waived the requirement for ethical approval and written informed consent for participants in this study. The cells were cultured at 37°C under 95% air and 5% carbon dioxide. For the ionizing radiation exposure, the cells were seeded onto six-well plate 1 day before exposure to 10 cGy of X-ray (320 kV, 2 mA) for 1 min (X-RAD 320 X-ray system).

Cell Viability Test

Cells were seeded in a 96-well plate (24 replicate wells for each treatment). After ionizing radiation exposure, cell viability was measured by the MTT assay (Sigma). Colorimetric reaction was measured at 570 nm.

RNA Isolation

After radiation exposure, total RNA of the cells was extracted using TRIzol Reagent (Invitrogen) according to the manufacturer’s instructions. Briefly, the cells were lysed in 1 ml of TRIzol. Then, 200 μl of chloroform was added, and the sample was centrifuged at 12,000 × g for 15 min. Next, 500 μl of the aqueous phase was mixed with 500 μl of isopropanol and stored in −20°C. After overnight precipitation, the mixture was centrifuged at 12,000 × g for 20 min. The RNA pellet was then washed twice using ice-cold 70% ethanol, followed by resuspension in RNAse-free distilled water. RNA quality was assessed using an Agilent 2100 Bioanalyzer system, and samples with an RNA integrity number (RIN) greater than eight was used for RNA library construction.

RNA Sequencing and Bioinformatics Analysis

RNA (cDNA) libraries (three biological replicates from each treatment) of cells were constructed as previously described (Li et al., 2021) and sequenced by the Beijing Genomics Institute (Wuhan, China). Single-end 50-bp read-length reads were sequenced on a BGISEQ-500RS sequencer. Sequence reads were dynamically trimmed according to the q algorithm of BWA (Li et al., 2021). At least 50 million quality-trimmed reads were obtained in each sample. Quality-trimmed sequence reads were mapped to human genome reference (GRCh38/hg38) for SKOV3 and COV434 cell lines, and mouse genome reference (GRCm39/mm39) for GC-1 and TM4 cell lines. Read-count data was then subjected to differential expression analysis using the edgeR package (Robinson et al., 2010). Genes with FDR < 0.05 were considered differentially expressed genes (DEGs). Furthermore, Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, and Ingenuity Pathway Analysis (IPA®, QIAGEN1) were used to decipher the biological effects and possible epigenetic effect of ionizing radiation on the female and male reproductive systems.

Data Availability

Sequencing data of transcriptome sequencing that support the findings of this study have been deposited in the NCBI BioProject database2 with the BioProject accession codes PRJNA730377.

Statistical Analysis

For bioinformatics analysis, all pathways, diseases, or biofunctions with p < 0.05 was considered statistically significant. In MTT assay, results were compared using paired Student’s t-test. All statistical analyses were performed using GraphPad Prism 3.02 (GraphPad Software Inc.), and results with p < 0.05 was considered statistically significant.

Results

Environmental Relevant Dose of Ionizing Radiation Has no Effect on the Viability of the Female and Male Reproductive Cells

MTT assay was employed to determine the cytotoxicity of the tested ionizing radiation (10 cGy). No significant change in cell viability was found after ionizing radiation exposure (Figure 1), suggesting that the level of radiation used in this study had no cytotoxic effect on both female and male reproductive cells.

FIGURE 1

Ionizing Radiation Induced Differential Gene Expression in Both Female and Male Reproductive Cells

In an attempt to understand the biological functions altered by the ionizing radiation, a comparative transcriptomic analysis was conducted. Deep sequencing of RNA libraries derived from control and ionizing radiation treatment groups of each cell lines generated at least 42 million quality-trimmed clean reads (Supplementary Table 2). A total of 51.2 Gb quality-trimmed bases were obtained from the transcriptome sequencing (Supplementary Table 2). Over 95% of sequencing reads could be mapped to the reference genome (Supplementary Table 2). Comparative transcriptomic analysis of ovarian cancer SKOV3 cells found a total 1,144 differentially expressed genes (DEGs), including 574 upregulated and 570 downregulated genes in the ionizing radiation group as compared with the control group (Figure 2A and Supplementary Table 3). Four thousand three hundred ninety-nine DEGs, including 1,575 upregulated and 2,824 downregulated genes were found in ovarian cancer COV434 cells after ionizing radiation exposure (Figure 2B and Supplementary Table 4). Upregulation of 24 genes and downregulation of 35 genes were common in SKOV3 and COV434 after ionizing radiation (Figure 2C and Table 1). Comparative transcriptomic analysis further revealed a smaller number of DEGs in male reproductive cell lines compared with female ovarian cancer cell lines after ionizing radiation. In testicular TM4 germ cell, a total of 783 DEGs, including 278 upregulated and 505 downregulated genes, were found in ionizing radiation group as compared with control group (Figure 3A and Supplementary Table 5). In another testicular GC-1 germ cell, we found 248 DEGs including 150 upregulated and 98 downregulated genes under ionizing radiation (Figure 3B and Supplementary Table 6). Notably, TM4 and GC-1 cells shared 26 upregulated and 32 downregulated genes (Figure 3C and Table 2).

TABLE 1

Common
ABHD14A-ACY16.321928095
AGAP6−0.26995382
AHNAK−0.036692267
ANKRD52−0.125579585
BHLHE400.159297821
CCNB10.118534794
CCT80.112181258
CELSR1−0.169768603
CEP170B−0.203081668
CLK1−0.581517732
CYR610.50372215
DAG1−0.137058218
DNAJB10.421062267
FLNA−0.137776264
FLNB−0.101035897
HCFC1−0.256265359
HMGA2−0.279167524
HSPA1A0.674635288
HSPA1B0.802124194
HSPA80.522838677
HSPE1-MOB41.696165063
HSPG2−0.131589484
ILF3−0.117865207
KMT2D−0.128633535
LARP1−0.130260559
LDHA0.070238521
LENG8−0.167621054
MAP3K14−0.246671519
MSH5−1.036220188
MTRNR2L10.086443006
MTRNR2L20.055259141
MTRNR2L60.366614351
MTRNR2L80.07506049
MYC0.465690304
NACA0.066625449
NBPF11−2.273922722
NDST1−0.324568416
NOL9−0.227755207
NOMO30.261103909
PKD1−0.212287115
PLEC−0.211283151
POLE−0.149276115
PPP1R15A0.341577734
PRRC2B−0.102852978
RHOB0.332811579
RPS178.515699838
SH3BP4−0.306432533
SMAD3−0.169052451
SNRNP200−0.082025905
SRCAP−0.18936661
SRRM2−0.08339425
TCP10.131439798
TJP1−0.138449491
TMEM189-UBE2V1−4.906890596
TNRC18−0.247618398
TPT10.082415032
TPX20.109746859
TRIO−0.102527386
UBR4−0.100565996

Common deregulated genes in SKOV3 and COV434 after radiation exposure.

TABLE 2

CommonGenes
6820431F20Rik−0.148053365
Actb0.015429673
Actg10.035451843
Amotl20.117501459
Atad5−0.144339873
Atrx−0.155713413
Bclaf1−0.091555479
Cd209c−0.538806202
Cd2ap−0.108090558
Cited20.074274825
Cmah−0.055197597
Col1a10.115874807
Ctgf0.213118333
Cxcr20.239034728
Cyr610.184896033
Ddx17−0.121821168
Dst−0.118465395
Dusp10.227403063
Egr10.294032038
Eif10.081156225
Eif4a2−0.147533967
Flnb0.044000894
Fn10.03142685
Fubp1−0.054877449
Hspa80.027844383
Ier20.169029369
Ier30.226920302
Luc7l3−0.306875903
Mat2a−0.098084338
Mki67−0.187693922
Myc0.09462868
Ncl−0.055239947
Npm1−0.02494927
Pcm1−0.107611547
Phlda10.066301904
Plk20.129986464
Psip1−0.107652929
Ptgs20.085702
Ranbp2−0.106473682
Rhob0.100308303
Rplp00.014678559
Scaf11−0.053057537
Serpine10.060359909
Setd2−0.085843794
Smarca5−0.063291942
Smc2−0.18366615
Smc3−0.112639541
Smc4−0.145773918
Sparc0.081296273
Srsf10−0.063547185
Thbs10.048740214
Tmem254c−7.996389141
Top2a−0.085094243
Trpm7−0.156997671
Ttc14−0.263145409
U2surp−0.094367453
Ubc0.050980317
Zfc3h1−0.170572408

Common deregulated genes in TM4 and GC-1 after radiation exposure.

FIGURE 2

FIGURE 3

Ionizing Radiation Altered Biological Processes and Signaling Pathways in the Reproductive Cells

The common DEGs were then subjected to GO and KEGG enrichment analysis, to elucidate the alteration of biological functions and signaling pathways in female and male reproductive cells caused by the ionizing radiation. The result of GO analysis in the female reproductive cells showed that the ionizing radiation would trigger DNA damage response, leading to the alteration of many biological processes related to cell apoptosis, cell growth, cell cycle arrest, protein stabilization and folding, cell-cell adhesion, gene expression, and in utero embryonic development (Figure 4A and Table 3). More importantly, our results showed that ionizing radiation can alter the biological processes closely related to epigenetic regulation such as establishment of protein localization to telomere, telomerase RNA localization to Cajal body, and heterochromatin assembly (Figure 4A and Table 3). Results of KEGG pathway analysis further highlighted that radiation causes alteration of MAPK signaling pathway, protein processing in endoplasmic reticulum, spliceosome, antigen processing and presentation, estrogen signaling pathway, and cell cycle (Figure 4B and Table 4). In the male reproductive cells, the ionizing radiation could similarly cause a DNA damage response (Figure 4C and Table 5), leading to alterations of different biological processes such as cell apoptosis, cell proliferation, cell cycle, and cell death (Figure 4C and Table 5). Although the responses were similar to those observed in female reproductive cells, the gene clusters involved in the processes were largely different (Table 6) and only shared the induction of MYC and CYR61 (Table 6). More importantly, the ionizing radiation resulted in the modulation of chromatin remodeling such as chromosome organization, meiotic chromosome segregation, and meiotic chromosome condensation, which was considered a major event leading to epigenetic modification in male reproductive cells (Figure 4C and Table 5). The chromatin remodeling is controlled by a family of structural maintenance of chromosome protein including SMC2, SMC3, and SMC4 and SJQ/SMB-related matrix-associated actin-dependent regulator of chromatin subfamily A member 5 (SMARCA5) (Table 5). The result of KEGG pathway analysis highlighted the alteration of focal adhesion, Hippo signaling pathway, ECM-receptor interaction, and MAPK signaling pathway (Figure 4D and Table 7). Taken together, our data suggested that ionizing radiation caused a similar biological alteration in both the female and male reproductive cells through the regulation of different gene clusters.

TABLE 3

TermCountP-valueGenes
GO:2001022∼positive regulation of response to DNA damage stimulus20.026451MYC, HMGA2
GO:0043066∼negative regulation of apoptotic process50.048141SMAD3, FLNA, HMGA2, CYR61, HSPA1B
GO:0030308∼negative regulation of cell growth30.037132SMAD3, SH3BP4, HSPA1A
GO:0007050∼cell cycle arrest30.016601SMAD3, MYC, PKD1
GO:0045787∼positive regulation of cell cycle20.092951CCNB1, MYC
GO:0050821∼protein stabilization52.61E-04SMAD3, TCP1, FLNA, CCT8, HCFC1
GO:0006457∼protein folding40.003713DNAJB1, HSPA8, TCP1, CCT8
GO:0098609∼cell–cell adhesion83.37E-07TJP1, DNAJB1, LDHA, LARP1, FLNB, HCFC1, PLEC, HSPA1A
GO:0010628∼positive regulation of gene expression50.015824PPP1R15A, HSPA8, SMAD3, HMGA2, HCFC1
GO:0042993∼positive regulation of transcription factor import into nucleus20.038247SMAD3, FLNA
GO:0001701∼in utero embryonic development50.005915KMT2D, CCNB1, SMAD3, MYC, PKD1
GO:1904851∼positive regulation of establishment of protein localization to telomere20.021693TCP1, CCT8
GO:1904874∼positive regulation of telomerase RNA localization to Cajal body20.035899TCP1, CCT8
GO:0032212∼positive regulation of telomere maintenance via telomerase20.075063TCP1, CCT8
GO:0031507∼heterochromatin assembly20.009699HMGA2, TNRC18
GO:0007339∼binding of sperm to zona pellucida30.003459TCP1, CCT8, HSPA1B
GO:0060428∼lung epithelium development20.019306HMGA2, PKD1
GO:1904871(positive regulation of protein localization to Cajal body20.019306TCP1, CCT8
GO:0060236(regulation of mitotic spindle organization20.019306TPX2, PKD1
GO:0001649(osteoblast differentiation30.030878SMAD3, SNRNP200, CYR61
GO:0051085(chaperone mediated protein folding requiring cofactor20.031186DNAJB1, HSPA8
GO:0044319(wound healing, spreading of cells20.031186FLNA, CYR61
GO:0006468(protein phosphorylation50.050757CLK1, ILF3, CCNB1, TRIO, MAP3K14
GO:0060441(epithelial tube branching involved in lung morphogenesis20.05683DAG1, HMGA2
GO:0045216(cell—cell junction organization20.05683SMAD3, FLNA
GO:0043484∼regulation of RNA splicing20.068266CLK1, AHNAK
GO:0001837∼epithelial to mesenchymal transition20.072803FLNA, HMGA2
GO:0042177∼negative regulation of protein catabolic process20.081811SMAD3, FLNA
GO:1901998∼toxin transport20.081811TCP1, CCT8
GO:0090307∼mitotic spindle assembly20.081811TPX2, FLNA
GO:0043085∼positive regulation of catalytic activity20.090733HSPA8, MYC
GO:0045944∼positive regulation of transcription from RNA polymerase II promoter60.092546KMT2D, SMAD3, MYC, HMGA2, PKD1, CYR61
GO:0048565∼digestive tract development20.097369CCNB1, PKD1

Gene ontology enrichment analysis on SKOV3 and COV434 common deregulated genes.

TABLE 4

TermCountP-valueGenes
mmu04010:MAPK signaling pathway77.20E-05HSPA8, MYC, FLNA, FLNB, MAP3K14, HSPA1B, HSPA1A
mmu04141:protein processing in endoplasmic reticulum50.001406PPP1R15A, DNAJB1, HSPA8, HSPA1B, HSPA1A
mmu03040:spliceosome40.006952HSPA8, SNRNP200, HSPA1B, HSPA1A
mmu04612:antigen processing and presentation30.024577HSPA8, HSPA1B, HSPA1A
mmu04915:estrogen signaling pathway30.034176HSPA8, HSPA1B, HSPA1A
mmu04110:cell cycle30.052334CCNB1, SMAD3, MYC
mmu05134:legionellosis30.012357HSPA8, HSPA1B, HSPA1A
mmu05164:influenza A40.013793DNAJB1, HSPA8, HSPA1B, HSPA1A
mmu05205:proteoglycans in cancer40.021748MYC, FLNA, FLNB, HSPG2
mmu05132:Salmonella infection30.022385TJP1, FLNA, FLNB
mmu05145:toxoplasmosis30.038769HSPA8, HSPA1B, HSPA1A
mmu04144:endocytosis40.04145HSPA8, SMAD3, HSPA1B, HSPA1A
mmu05166:HTLV-I infection40.047645SMAD3, MYC, MAP3K14, POLE
mmu05162:measles30.06166HSPA8, HSPA1B, HSPA1A

KEGG enrichment analysis on SKOV3 and COV434 common deregulated genes.

TABLE 5

TermCountp-valueGenes
GO:2001022∼positive regulation of response to DNA damage stimulus20.031782BCLAF1, MYC
GO:0006974∼cellular response to DNA damage stimulus50.03441TOP2A, ATAD5, MYC, ATRX, SMC3
GO:0043066∼negative regulation of apoptotic process103.51E-05NPM1, CITED2, DUSP1, NCL, PLK2, CXCR2, FN1, THBS1, CYR61, IER3
GO:0043065∼positive regulation of apoptotic process85.22E-05TOP2A, BCLAF1, DUSP1, TRPM7, PTGS2, PHLDA1, CYR61, RHOB
GO:0043280∼positive regulation of cysteine-type endopeptidase activity involved in apoptotic process30.008723MYC, CYR61, CTGF
GO:0008284∼positive regulation of cell proliferation70.004907NPM1, MYC, CXCR2, FN1, PTGS2, THBS1, CTGF
GO:0071364∼cellular response to epidermal growth factor stimulus30.004982COL1A1, NCL, MYC
GO:0007049∼cell cycle60.033428DUSP1, MKI67, SMC3, SMC4, SMC2, CD2AP
GO:0051383∼kinetochore organization20.008768SMC4, SMC2
GO:0010942∼positive regulation of cell death30.008723EGR1, PTGS2, CTGF
GO:0010628∼positive regulation of gene expression70.001063EGR1, HSPA8, CITED2, FUBP1, SERPINE1, FN1, CTGF
GO:0045893∼positive regulation of transcription, DNA-templated70.006565DDX17, COL1A1, EGR1, NPM1, CITED2, MYC, SMARCA5
GO:0010629∼negative regulation of gene expression50.007523NPM1, CITED2, MYC, SERPINE1, CTGF
GO:0045944∼positive regulation of transcription from RNA polymerase II promoter90.007997DDX17, TOP2A, EGR1, CITED2, NCL, MYC, ATRX, PSIP1, CYR61
GO:0051276∼chromosome organization41.87E-04MYC, SMC3, SMC4, SMC2
GO:0045132∼meiotic chromosome segregation20.014571SMC4, SMC2
GO:0010032∼meiotic chromosome condensation20.014571SMC4, SMC2
GO:0006338∼chromatin remodeling30.022492MYC, ATRX, SMARCA5
GO:0001525∼angiogenesis76.84E-05SETD2, NCL, SERPINE1, FN1, PTGS2, CTGF, RHOB
GO:0032355∼response to estradiol52.18E-04COL1A1, DUSP1, MYC, PTGS2, CTGF
GO:0051591∼response to cAMP44.48E-04COL1A1, SPARC, MAT2A, DUSP1
GO:0030261∼chromosome condensation39.85E-04TOP2A, SMC4, SMC2
GO:0048661∼positive regulation of smooth muscle cell proliferation40.001666EGR1, MYC, PTGS2, THBS1
GO:0042060∼wound healing40.002641COL1A1, SPARC, SERPINE1, FN1
GO:0045766∼positive regulation of angiogenesis40.005367SERPINE1, CXCR2, THBS1, RHOB
GO:0043085∼positive regulation of catalytic activity30.005828HSPA8, NPM1, MYC
GO:0010595∼positive regulation of endothelial cell migration30.006425SPARC, THBS1, RHOB
GO:0035914∼skeletal muscle cell differentiation30.011335EGR1, CITED2, MYC
GO:0007155∼cell adhesion60.013561DST, FN1, THBS1, CYR61, CTGF, RHOB
GO:0010757∼negative regulation of plasminogen activation20.014571SERPINE1, THBS1
GO:0048146∼positive regulation of fibroblast proliferation30.015584MYC, SERPINE1, FN1
GO:0051592∼response to calcium ion30.016973SPARC, DUSP1, THBS1
GO:0071636∼positive regulation of transforming growth factor beta production20.017461PTGS2, THBS1
GO:0043434∼response to peptide hormone30.0194COL1A1, SPARC, CTGF
GO:0051918∼negative regulation of fibrinolysis20.020341SERPINE1, THBS1
GO:0046599∼regulation of centriole replication20.020341NPM1, PLK2
GO:0051384∼response to glucocorticoid30.020409SPARC, DUSP1, PTGS2
GO:0030335∼positive regulation of cell migration40.021653COL1A1, FN1, THBS1, CYR61
GO:0009749∼response to glucose30.021963EGR1, THBS1, CTGF
GO:0071347∼cellular response to interleukin-130.023026MYC, SERPINE1, FN1
GO:1904628∼cellular response to phorbol 13-acetate 12-myristate20.023214MYC, RPLP0
GO:0007613∼memory30.024659PLK2, TRPM7, PTGS2
GO:0000375∼RNA splicing, via transesterification reactions20.026078SCAF11, SRSF10
GO:0030728∼ovulation20.031782MYC, PTGS2
GO:0008380∼RNA splicing40.033622HSPA8, SCAF11, LUC7L3, SRSF10
GO:0006334(nucleosome assembly30.037334NPM1, ATRX, SMARCA5
GO:0003197(endocardial cushion development20.037453CITED2, THBS1
GO:0007076(mitotic chromosome condensation20.037453SMC4, SMC2
GO:0030336(negative regulation of cell migration30.038648CITED2, SERPINE1, RHOB
GO:0030194(positive regulation of blood coagulation20.040277SERPINE1, THBS1
GO:0050921(positive regulation of chemotaxis20.043092FN1, THBS1
GO:0006376(mRNA splice site selection20.043092LUC7L3, SRSF10
GO:0007067(mitotic nuclear division40.047549SMC3, SMC4, SMC2, CD2AP
GO:0046907(intracellular transport20.048698RANBP2, DST
GO:0070542(response to fatty acid20.054272PTGS2, CTGF
GO:0006259(DNA metabolic process20.057047TOP2A, MKI67
GO:0046697(decidualization20.068067CITED2, PTGS2
GO:0007126(meiotic nuclear division20.068067MKI67, SMC3
GO:0030511(positive regulation of transforming growth factor beta receptor signaling pathway20.068067CITED2, THBS1
GO:0043044(ATP-dependent chromatin remodeling20.070802SMARCA5, ACTB
GO:0001895(retina homeostasis20.070802ACTB, ACTG1
GO:0048025(negative regulation of mRNA splicing, via spliceosome20.070802NPM1, SRSF10
GO:0010667∼negative regulation of cardiac muscle cell apoptotic process20.07353HSPA8, NPM1
GO:0016925∼protein sumoylation20.076249TOP2A, RANBP2
GO:0042493∼response to drug40.077109COL1A1, MAT2A, PTGS2, THBS1
GO:0007052∼mitotic spindle organization20.078961PLK2, SMC3
GO:0007569∼cell aging20.084361NPM1, CITED2
GO:0001937∼negative regulation of endothelial cell proliferation20.084361SPARC, THBS1
GO:0044344∼cellular response to fibroblast growth factor stimulus20.084361COL1A1, MYC
GO:0070372∼regulation of ERK1 and ERK2 cascade20.087049FN1, CYR61
GO:0033574∼response to testosterone20.089729DUSP1, THBS1
GO:0006446∼regulation of translational initiation20.092402EIF4A2, EIF1
GO:0006351∼transcription, DNA-templated100.095659DDX17, EGR1, HSPA8, BCLAF1, SETD2, CITED2, MYC, FUBP1, ATRX, PSIP1
GO:0051301∼cell division40.09664SMC3, SMC4, SMC2, CD2AP

Gene ontology enrichment analysis on GC-1 and TM4 common deregulated genes.

TABLE 6

SKOV3/COV434Biological processesTM4/GC-1
MYC, HMGA2DNA damage responseBCLAF1, TOP2A, ATAD5, MYC, ATRX, SMC3
SMAD3, FLNA, HMGA2, CYR61, HSPA1BCell apoptosisNPM1, CITED2, DUSP1, NCL, PLK2, CXCR2, FN1, THBS1, CYR61, IER3, MYC, CTGF, TOP2A, BCLAF1, TRPM7, PTGS2, PHLDA1, RHOB
SMAD3, MYC, PKD1, CCNB1Cell cycleDUSP1, MKI67, SMC3, SMC4, SMC2, CD2AP, MYC
SMAD3, SH3BP4, HSPA1ACell proliferationNPM1, MYC, CXCR2, FN1, PTGS2, THBS1, CTGF, COL1A1, NCL,
PPP1R15A, HSPA8, SMAD3, HMGA2, HCFC1, FLNAGene expressionEGR1, HSPA8, CITED2, FUBP1, SERPINE1, FN1, CTGF, NPM1, MYC, DDX17, COL1A1, SMARCA5

Different gene clusters involved in similar biological processes.

TABLE 7

TermCountp-valueGenes
mmu04510:focal adhesion67.98E-04COL1A1, FN1, FLNB, THBS1, ACTB, ACTG1
mmu04390:hippo signaling pathway50.002025MYC, SERPINE1, ACTB, ACTG1, CTGF
mmu04512:ECM-receptor interaction30.040363COL1A1, FN1, THBS1
mmu04010:MAPK signaling pathway40.061762HSPA8, DUSP1, MYC, FLNB
mmu05205:proteoglycans in cancer77.28E-05COL1A1, MYC, FN1, FLNB, THBS1, ACTB, ACTG1
mmu05100:bacterial invasion of epithelial cells40.002739FN1, ACTB, CD2AP, ACTG1
mmu04145:phagosome40.023537CD209C, THBS1, ACTB, ACTG1
mmu05132:Salmonella infection30.032375FLNB, ACTB, ACTG1
mmu05222:small cell lung cancer30.037085MYC, FN1, PTGS2
mmu04919:thyroid hormone signaling pathway30.064117MYC, ACTB, ACTG1
mmu04611:platelet activation30.081647COL1A1, ACTB, ACTG1
mmu03040:spliceosome30.0838HSPA8, SRSF10, U2SURP

KEGG enrichment analysis on GC-1 and TM4 common deregulated genes.

FIGURE 4

Ionizing Radiation Caused Deregulation of Gene Network and Reproductive Disorder

In order to understand the reproductive diseases and to delineate the gene networks involved in ionizing radiation-altered biological processes, IPA was conducted. The results of disease analysis showed that the radiation exposure could potentially lead to many female reproductive disorders and diseases such as ovary growth impairment, genital tract cancer, and ovarian cancer (Table 8). A similar finding was observed in the male reproductive system that radiation could cause tumorigenesis of reproductive tract and develop malignant neoplasm and endometriosis of the male genital organ (Table 9). Canonical pathway analysis of IPA further highlighted the involvement of many components from different cellular levels in radiation-mediated reproductive impairment. In the female reproductive cells, tight junction protein 1 (TJP1) was found to be downregulated and associated with the deregulation of filamin A, alpha (FLNA) and heat shock 70 kDa protein 1 (HSPA1A) (Figure 5A). Also, a group of enzymes such as lactate dehydrogenase A (LDHA), DNA polymerase Epsilon (POLE), and small nuclear ribonucleoprotein U5 Subunit 200 (SNRP200) (Figure 5A) were highlighted. More importantly, a cluster of transcription factors including DnaJ heat shock protein family (Hsp40) member B1 (DNAJB1), mothers against decapentaplegic homolog 3 (SMAD3) and MYC proto-oncogene, and BHLH transcription factor (MYC) were involved in the ionizing radiation-mediated reproductive impairment (Figure 5A). Our result showed that MYC is a key mediator that directly controls different kinases and enzymes (Figure 5A). In the male reproductive cells, more candidates were found to be involved in reproductive impairment caused by the ionizing radiation (Figure 5B). The gene networking showed the induction of thrombospondin-1 (THBS1), fibronectin 1 (FN1), collagen, type I, alpha 1 (COL1A1), and serpin family E member 1 (SERPINE1) in extracellular matrix (Figure 5B). Also, radiation exposure could lead to downregulation of membrane protein transient receptor potential cation channel subfamily M member 7 (TRPM7). Similar to the result in the female reproductive cells, MYC was the major mediator connecting different enzymes such as methionine adenosyltransferase 2A (MAT2), Ras homolog family member B (RHOB), ubiquitin C (UBC), prostaglandin-endoperoxide synthase 2 (PTGS2), and heat shock protein family A member 8 (HSPA8), as well as a large number of transcription factors including nucleophosmin (NPM1), Cbp/p300-interacting transactivator 2 (CITED2), and ATP-dependent helicase (ATRX) (Figure 5B). Taken together, our results suggested a gender-specific gene network involved in the ionizing radiation-mediated reproductive impairment.

TABLE 8

Diseases or functions annotationp-value# moleculesMolecules
Genital tract cancer2.64E-0846ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, BHLHE40, CCN1, CCNB1, CELSR1, CEP170B, CLK1, FLNA, FLNB, HCFC1, HMGA2, HSPA1A/HSPA1B, HSPA8, HSPE1-MOB4, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NACA, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, PPP1R15A, PRRC2B, RHOB, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TJP1, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Genital tumor4.45E-0846ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, BHLHE40, CCN1, CCNB1, CELSR1, CEP170B, CLK1, FLNA, FLNB, HCFC1, HMGA2, HSPA1A/HSPA1B, HSPA8, HSPE1-MOB4, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NACA, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, PPP1R15A, PRRC2B, RHOB, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TJP1, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Female genital tract cancer1.02E-0740ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, CCN1, CCNB1, CELSR1, CLK1, FLNA, FLNB, HCFC1, HMGA2, HSPA1A/HSPA1B, HSPE1-MOB4, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, RHOB, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TJP1, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Tumorigenesis of reproductive tract1.79E-0740ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, CCN1, CCNB1, CELSR1, CLK1, FLNA, FLNB, HCFC1, HMGA2, HSPA1A/HSPA1B, HSPE1-MOB4, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, RHOB, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TJP1, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Female genital tract adenocarcinoma2.10E-0737ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, CCNB1, CELSR1, CLK1, FLNA, FLNB, HCFC1, HMGA2, HSPA1A/HSPA1B, HSPE1-MOB4, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Endometrioid endometrial adenocarcinoma2.98E-0733ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, CELSR1, FLNA, FLNB, HCFC1, HSPA1A/HSPA1B, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Endometrial cancer3.27E-0735ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, CCNB1, CELSR1, FLNA, FLNB, HCFC1, HSPA1A/HSPA1B, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TJP1, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Development of genital tumor3.67E-0739ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, CCN1, CCNB1, CELSR1, CLK1, FLNA, FLNB, HCFC1, HMGA2, HSPA1A/HSPA1B, HSPE1-MOB4, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TJP1, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Endometrial carcinoma5.26E-0734ABHD14A-ACY1, AGAP6 (includes others), AHNAK, ANKRD52, CELSR1, FLNA, FLNB, HCFC1, HSPA1A/HSPA1B, HSPG2, ILF3, KMT2D, LARP1, LENG8, MAP3K14, MSH5, MYC, NBPF10 (includes others), NDST1, NOMO1 (includes others), PKD1, PLEC, POLE, SH3BP4, SMAD3, SNRNP200, SRCAP, SRRM2, TJP1, TMEM189-UBE2V1, TNRC18, TPX2, TRIO, UBR4
Ovarian tumor6.44E-0518AHNAK, ANKRD52, CCN1, CCNB1, CLK1, HCFC1, HMGA2, HSPE1-MOB4, KMT2D, MYC, PKD1, POLE, RHOB, SMAD3, SNRNP200, SRCAP, TNRC18, TPX2
Ovarian cancer1.63E-0417AHNAK, ANKRD52, CCN1, CCNB1, CLK1, HCFC1, HMGA2, HSPE1-MOB4, KMT2D, MYC, PKD1, POLE, RHOB, SNRNP200, SRCAP, TNRC18, TPX2
Growth of genital organ4.23E-044KMT2D, MYC, NDST1, SMAD3
Growth of ovary6.45E-043KMT2D, MYC, SMAD3
Ovarian carcinoma1.76E-0314AHNAK, ANKRD52, CLK1, HCFC1, HMGA2, HSPE1-MOB4, KMT2D, MYC, PKD1, POLE, SNRNP200, SRCAP, TNRC18, TPX2

Female reproductive diseases caused by radiation exposure.

TABLE 9

Diseases or functions annotationp-value# moleculesMolecules
Malignant neoplasm of male genital organ9.69E-0730ATAD5, ATRX, BCLAF1, CD2AP, COL1A1, DDX17, DUSP1, EGR1, FLNB, FN1, HSPA8, IER3, LUC7L3, MKI67, MYC, NPM1, PHLDA1, PLK2, PSIP1, PTGS2, RGPD4 (includes others), SETD2, SMC3, SMC4, SPARC, SRSF10, TOP2A, TRPM7, UBC, ZFC3H1
Endometriosis3.09E-069ACTB, CITED2, CXCR2, DUSP1, EGR1, FN1, PLK2, PTGS2, TOP2A
Prostate cancer9.27E-0628ATAD5, ATRX, CD2AP, COL1A1, DDX17, DUSP1, EGR1, FLNB, FN1, HSPA8, IER3, LUC7L3, MKI67, MYC, NPM1, PHLDA1, PLK2, PSIP1, PTGS2, RGPD4 (includes others), SETD2, SMC3, SMC4, SPARC, SRSF10, TOP2A, UBC, ZFC3H1
Tumorigenesis of reproductive tract9.69E-0636ACTB, ACTG1, ATAD5, ATRX, BCLAF1, CD2AP, COL1A1, CXCR2, DUSP1, EIF4A2, FLNB, FN1, FUBP1, IER3, LUC7L3, MKI67, MYC, NPM1, PHLDA1, PLK2, PSIP1, PTGS2, RGPD4 (includes others), RHOB, SCAF11, SERPINE1, SETD2, SMARCA5, SMC3, SMC4, THBS1, TOP2A, TRPM7, U2SURP, UBC, ZFC3H1
Genital tumor1.16E-0541ACTB, ACTG1, ATAD5, ATRX, BCLAF1, CD2AP, COL1A1, CXCR2, DDX17, DUSP1, EGR1, EIF4A2, FLNB, FN1, FUBP1, HSPA8, IER3, LUC7L3, MKI67, MYC, NPM1, PHLDA1, PLK2, PSIP1, PTGS2, RGPD4 (includes others), RHOB, SCAF11, SERPINE1, SETD2, SMARCA5, SMC3, SMC4, SPARC, SRSF10, THBS1, TOP2A, TRPM7, U2SURP, UBC, ZFC3H1
Morphology of reproductive system1.68E-0511ATRX, CXCR2, DUSP1, EGR1, FUBP1, MYC, PLK2, PTGS2, SERPINE1, SETD2, THBS1
Prostatic carcinoma9.71E-0524ATAD5, ATRX, CD2AP, COL1A1, DDX17, DUSP1, FLNB, FN1, HSPA8, LUC7L3, MKI67, MYC, NPM1, PHLDA1, PLK2, PTGS2, RGPD4 (includes others), SETD2, SMC3, SMC4, SRSF10, TOP2A, UBC, ZFC3H1
Development of genital tumor1.75E-0433ACTB, ACTG1, ATAD5, ATRX, BCLAF1, CD2AP, COL1A1, CXCR2, EIF4A2, FLNB, FN1, FUBP1, IER3, LUC7L3, MKI67, MYC, NPM1, PHLDA1, PLK2, PSIP1, PTGS2, RGPD4 (includes others), SCAF11, SETD2, SMARCA5, SMC3, SMC4, THBS1, TOP2A, TRPM7, U2SURP, UBC, ZFC3H1

Male reproductive diseases caused by radiation exposure.

FIGURE 5

Discussion

The present study unraveled that epigenetic changes associated with reproductive impairment can be elicited by a low dose of ionizing radiation (i.e., 100 mGy). The dose we used in this study is considered a low-dose radiation because fetal absorbed about 4.9 cGy per procedure during abdominal CT scan, and for pelvis it was as high as 7.9 cGy (Groen et al., 2012). A cohort study of diagnostic medical radiation workers in South Korea showed that the workers absorbed a mean cumulative badge dose of 7.20 mSv (Lee et al., 2021). While the environmental level of ionizing radiation is about 3 mGy/year, high levels of background radiation was commonly reported in some areas around the world (Hendry et al., 2009).

In this study, we determined and compared the effects of radiation on different types of male and female reproductive cells. GC-1 spg is an immortalized type B spermatogonia line and displays characteristics of a stage between type B spermatogonia and primary spermatocytes. TM4 is an immortalized Sertoli cell line and is a somatic cell-type essential for nurturing germ cell development. Due to the low availability of freshly isolated primary ovarian reproductive cells, two ovarian cancer lines derived from ovarian surface epithelium and granulosa cells, respectively, are used. Our comparative transcriptomic analysis revealed that environmental relevant level of ionizing radiation is sufficient to induced DNA damage, leading to cell cycle arrest and cell apoptosis in both female and male reproductive cells. It is not surprising that the genes involved were largely different according to cell types, and only MYC and CYR61 were shared in both female and male reproductive cells. MYC is a transcription factor and contributes to the development of many cancer types (Godwin et al., 2021). In addition, it plays an important role in DNA damage response repair. For instance, MYC formed foci with γ-H2AX to phosphorylate ATM and to mediate DNA-PKcs activity (Cui et al., 2015), leading to enhanced chromosomal and chromatid breaks in response to γ-ray ionizing radiation (Li et al., 2012). In addition, MYC was essential for DNA damage-induced apoptosis through the control of the p53 tumor suppressor protein (Phesse et al., 2014). The other deregulated gene, cysteine-rich angiogenic inducer 61 (CYR61), is an extracellular matrix-associated signaling protein of the CCN intercellular signaling protein family (Lau, 2011). CYR61 is a cell apoptosis and senescence inducer involved in DNA damage response through the regulation of p53 upon genotoxic stress (Feng et al., 2008; Morrison et al., 2009). In addition, CYR61 was a downstream effector of the Hippo signaling pathway that is one the important cell signaling pathways of DNA damage response (Shome et al., 2020). Other than the common DNA damage responses, our result also suggested that the ionizing radiation could cause gene deregulation related to impairment of utero embryonic development in the female reproductive system. Although it has been reported that abdominopelvic ionizing irradiation (>5 Gy) increased the risk of unfavorable neonatal outcomes such as fetal malformation and disturbances of growth or development (Kumar and DeJesus, 2021), our result showed that even a much lower level of ionizing radiation (10 cGy) could also pose adverse effect on embryonic development. In addition, our data highlighted that ionizing radiation could alter estrogen signaling pathway in female reproductive cells. Estrogen signaling pathway is one of the most important pathways involved in steroid hormones and reproductive regulation in mammals (Saito and Cui, 2018), and estrogens are also associated with tumor development, particularly breast and ovary cancers (Ferreira et al., 2009). Studies also demonstrated that estrogen is also related to the DNA repair pathways and DNA integrity (Jiménez-Salazar et al., 2021; Pescatori et al., 2021).

Our results further highlighted that ionized radiation could alter a spectrum of male reproductive function-related cell signaling such as ECM-receptor interaction, Hippo signaling pathway, and MAPK signaling pathway. Extracellular matrix (ECM) receptor interaction pathway has been reported to be associated with spermatogenesis (Gong et al., 2018). In addition, ECM remodeling is required for the testicular development and maturation (Slongo et al., 2002). Hippo signaling pathway is responsible for the controls of organ size through the regulation of cell proliferation and apoptosis (Zhang et al., 2009). Yes-associated protein, the downstream effector of Hippo signaling pathway was reported to modulate the decline of germline stem cells and niche cells (Francis et al., 2019). Also, Hippo signaling cascade was found to be associated with the pubertal development of male reproductive tract and spermatogenesis in sheep (Zhang et al., 2019). MAPK signaling pathway is the key mediator that controls the phosphorylation of many downstream effectors, leading to modulate different cellular functions, including cell proliferation, differentiation, and migration (Joerger-Messerli et al., 2021). Oxidative stress-mediated p38 MAPK signaling pathway was associated with the blood-testis barrier-related junction protein and promoting apoptosis in mice testes (Huang et al., 2021). Moreover, activation of the MAPK signaling pathway was involved in the molecular mechanism of apoptosis in spermatogonia cells (Park et al., 2020). As such, deregulation of these pathways found in the present study suggested different aspects of reproductive impairment caused by the ionizing radiation.

Beside the cell signaling pathways related to reproductive functions, our result also demonstrated that ionizing radiation may cause epigenetic modification in both female and male reproductive cells. In female reproductive cells, telomerase, telomere-binding proteins, and heterochromatin assembly were disturbed by ionizing radiation. In male reproductive cells, chromosome organization such as meiotic chromosome segregation and condensation was interfered. Telomere, a repetitive DNA sequencing at the end of chromosome, protects chromosome from progressive degradation. The length of telomere is controlled by a group of telomere-binding proteins and the enzymatic activity of telomerase. Recent research showed that dysregulated telomerase activation was associated with epigenetic, transcriptional, and posttranscriptional modifications (Dogan and Forsyth, 2021). Lister-Shimauchi et al. (2021) demonstrated that deficiency for telomerase in Caenorhabditis elegans resulted in transgenerational shortening of telomeres. In addition, histone modifications including activation markers (H3K4me1 and H3K4me3) and silencing marks (H3K27me3) at distal promoters were telomere length dependent, suggesting that the epigenetic state of telomere-distal promoters could be influenced by telomere length (Mukherjee et al., 2018). Moreover, meiotic chromosome segregation is essential for the maintenance of genomic integrity of gametes and requires functional centromeres that are required for a precise epigenetic inheritance (Mahlke and Nechemia-Arbely, 2020; Tan et al., 2020).

In conclusion, we show that environmental relevant dose of ionizing radiation can alter the expression of gene cluster related to DNA damage response through the control of MYC, which agreed with other studies reporting amplification of MYC gene in somatic cells and cancer cells (Sawey et al., 1987; Mothersill et al., 1991; Käcker et al., 2013; Ginter et al., 2014). The ionizing radiation-mediated reproductive impairment is not only gender specific but also connected with different gene networks and pathways. Due to the limited number of germ cells and their extended period of developmental processes, research on the radiation-induced effects of female germline has been hindered (Skrzypek et al., 2019). In this study, we have unraveled possible pathways in ovarian cells altered by radiation, thus providing insights on the mechanisms underlying the perturbed reproductive functions. More importantly, our findings suggested that ionizing radiation can interfere telomere and chromatin remodeling, leading to possible epigenetic changes. Further investigation is warranted to elucidate whether these induced modifications can lead to transgenerational reproductive impairments.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

Sequencing data of transcriptome sequencing that support the findings of this study have been deposited in the NCBI BioProject database2 with the BioProject accession codes PRJNA730377.

Ethics statement

In accordance with the national legislation and the institutional requirements, the Human Research Ethics Committee of The University of Hong Kong waived the requirement for ethical approval and written informed consent for participants in this study.

Author contributions

RW, AW, WL, KWY, RK, and JC contributed to conception and design of the study. CL and YY performed the experiments. KL organized the database. TC, XL, and NT performed the bioinformatics and statistical analysis. KL, RW, and KNY wrote the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding

This project was supported by the Hong Kong Branch of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (SMSEGL20SC02). KL was supported by the Hong Kong SAR, Macao SAR, and Taiwan Province Talented Young Scientist Program of Guangxi.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.710143/full#supplementary-material

References

  • 1

    AhmedE.ScherthanH.deRooijD. (2015). DNA double strand break response and limited repair capacity in mouse elongated spermatids.Int. J. Mol. Sci.162992329935. 10.3390/ijms161226214

  • 2

    CliftonD. K.BremnerW. J. (1983). The effect of testicular X-irradiation on spermatogenesis in man.J. Androl.4387392. 10.1002/j.1939-4640.1983.tb00765.x

  • 3

    CritchleyH. O. D.WallaceW. H.ShaletS. M.MamtoraH.HigginsonJ.AndersonD. C. (1992). Abdominal irradiation in childhood; the potential for pregnancy.Br J Obstet Gynaecol.99392394.

  • 4

    CuiF.FanR.ChenQ.HeY.SongM.ShangZ.et al (2015). The involvement of c-Myc in the DNA double-strand break repair via regulating radiation-induced phosphorylation of ATM and DNA-PKcs activity.Mol. Cell. Biochem.4064351. 10.1007/s11010-015-2422-2

  • 5

    DoganF.ForsythN. R. (2021). Telomerase regulation: a role for epigenetics.Cancers (Basel)13:1213. 10.3390/cancers13061213

  • 6

    DubrovaY. E.SarapultsevaE. I. (2020). Radiation-induced transgenerational effects in animals.Int. J. Radiat. Biol.10.1080/09553002.2020.1793027[Epub ahead of print].

  • 7

    FajardoL. F. (2005). The pathology of ionizing radiation as defined by morphologic patterns.Acta Oncol. (Madr.)441322. 10.1080/02841860510007440

  • 8

    FengP.WangB.RenE. C. (2008). Cyr61/CCN1 is a tumor suppressor in human hepatocellular carcinoma and involved in DNA damage response.Int. J. Biochem. Cell Biol.4098109. 10.1016/j.biocel.2007.06.020

  • 9

    FerreiraA. M.WestersH.AlbergariaA.SerucaR.HofstraR. M. W. (2009). Estrogens, MSI and lynch syndrome-associated tumors.Biochim. Biophys. Acta1796194200. 10.1016/j.bbcan.2009.06.004

  • 10

    FilkowskiJ. N.IlnytskyyY.TammingaJ.KoturbashI.GolubovA.BagnyukovaT.et al (2010). Hypomethylation and genome instability in the germline of exposed parents and their progeny is associated with altered miRNA expression.Carcinogenesis3111101115. 10.1093/carcin/bgp300

  • 11

    FrancisD.ChananaB.FernandezB.GordonB.MakT.PalaciosI. M. (2019). YAP/Yorkie in the germline modulates the age-related decline of germline stem cells and niche cells.PLoS One14:e0213327. 10.1371/journal.pone.0213327

  • 12

    García-RodríguezA.GosálvezJ.AgarwalA.RoyR.JohnstonS. (2018). DNA damage and repair in human reproductive cells.Int. J. Mol. Sci.20:31. 10.3390/ijms20010031

  • 13

    GinterP. S.MosqueraJ. M.MacDonaldT. Y.D’AlfonsoT. M.RubinM. A.ShinS. J. (2014). Diagnostic utility of MYC amplification and anti-MYC immunohistochemistry in atypical vascular lesions, primary or radiation-induced mammary angiosarcomas, and primary angiosarcomas of other sites.Hum. Pathol.45709716. 10.1016/j.humpath.2013.11.002

  • 14

    GodwinI.AntoN. P.BavaS. V.BabuM. S.JineshG. G. (2021). Targeting K-Ras and apoptosis-driven cellular transformation in cancer.Cell Death Discov.7:80.

  • 15

    GongJ.ZhangQ.WangQ.MaY.DuJ.ZhangY.et al (2018). Identification and verification of potential piRNAs from domesticated yak testis.Reproduction155117127.

  • 16

    GoodheadD. T. (1994). Initial events in the cellular effects of ionizing radiations: clustered damage in DNA.Int. J. Radiat. Biol.65717. 10.1080/09553009414550021

  • 17

    GrandjeanV.FourréS.De AbreuD. A. F.DerieppeM.-A.RemyJ.-J.RassoulzadeganM. (2016). RNA-mediated paternal heredity of diet-induced obesity and metabolic disorders.Sci. Rep.5:18193.

  • 18

    GroenR. S.BaeJ. Y.LimK. J. (2012). Fear of the unknown: ionizing radiation exposure during pregnancy.Am. J. Obstet. Gynecol.206456462. 10.1016/j.ajog.2011.12.001

  • 19

    HainesG. A.HendryJ. H.DanielC. P.MorrisI. D. (2002). Germ cell and dose-dependent DNA damage measured by the comet assay in murine Spermatozoaa after testicular X-irradiation1.Biol. Reprod.67854861. 10.1095/biolreprod.102.004382

  • 20

    HalimiM.AsghariS. M.SaririR.MoslemiD.ParsianH. (2012). Cellular response to ionizing radiation: a microRNA story.Int. J. Mol. Cell. Med.1178184.

  • 21

    HendryJ. H.SimonS. L.WojcikA.SohrabiM.BurkartW.CardisE.et al (2009). Human exposure to high natural background radiation: what can it teach us about radiation risks?J. Radiol. Prot.29A29A42.

  • 22

    HuangW.LiuM.XiaoB.ZhangJ.SongM.LiY.et al (2021). Aflatoxin B1 disrupts blood-testis barrier integrity by reducing junction protein and promoting apoptosis in mice testes.Food Chem. Toxicol.148:111972. 10.1016/j.fct.2021.111972

  • 23

    JenkinsT. G.CarrellD. T. (2012). The sperm epigenome and potential implications for the developing embryo.Reproduction143727734. 10.1530/rep-11-0450

  • 24

    Jiménez-SalazarJ. E.Damian-FerraraR.ArteagaM.BatinaN.Damián-MatsumuraP. (2021). Non-genomic actions of estrogens on the DNA repair pathways are associated with chemotherapy resistance in breast cancer.Front. Oncol.11:631007. 10.3389/fonc.2021.631007

  • 25

    Joerger-MesserliM. S.ThomiG.HaeslerV.KellerI.RenzP.SurbekD. V.et al (2021). Human Wharton’s Jelly mesenchymal stromal cell-derived small extracellular vesicles drive oligodendroglial maturation by restraining MAPK/ERK and notch signaling pathways.Front. Cell Dev. Biol.9:622539. 10.3389/fcell.2021.622539

  • 26

    KäckerC.MarxA.MössingerK.SvehlaF.SchneiderU.HogendoornP. C.et al (2013). High frequency of MYC gene amplification is a common feature of radiation-induced sarcomas. Further results from EORTC STBSG TL 01/01.Genes Chromosomes Cancer529398. 10.1002/gcc.22009

  • 27

    KamstraJ. H.HuremS.MartinL. M.LindemanL. C.LeglerJ.OughtonD.et al (2018). Ionizing radiation induces transgenerational effects of DNA methylation in zebrafish.Sci. Rep.8:15373.

  • 28

    KumarD.SalianS. R.KalthurG.UppangalaS.KumariS.ChallapalliS.et al (2013). Semen abnormalities, sperm DNA damage and global hypermethylation in health workers occupationally exposed to ionizing radiation.PLoS One8:e69927. 10.1371/journal.pone.0069927

  • 29

    KumarR.DeJesusO. (2021). “Radiation Effects on the Fetus,” in StatPearls [Internet], (Treasure Island, FL: StatPearls Publishing).

  • 30

    KumarR.HorikoshiN.SinghM.GuptaA.MisraH. S.AlbuquerqueK.et al (2013). Chromatin modifications and the DNA damage response to ionizing radiation.Front. Oncol.2:214. 10.3389/fonc.2012.00214

  • 31

    LauL. F. (2011). CCN1/CYR61: the very model of a modern matricellular protein.Cell. Mol. Life Sci.6831493163. 10.1007/s00018-011-0778-3

  • 32

    LeeW. J.KoS.BangY. J.ChoeS. A.ChoiY.PrestonD. L. (2021). Occupational radiation exposure and cancer incidence in a cohort of diagnostic medical radiation workers in South Korea.Occup. Environ. Med.10.1136/oemed-2021-107452[Epub ahead of print].

  • 33

    LiR.HuangC.HoC. H.LeungC. T.KongR. Y. C.LiangX.et al (2021). The use of glutathione to reduce oxidative stress status and its potential for modifying the extracellular matrix organization in cleft lip.Free Radic. Biol. Med.164130138. 10.1016/j.freeradbiomed.2020.12.455

  • 34

    LiZ.OwonikokoT. K.SunS. Y.RamalingamS. S.DoetschP. W.XiaoZ. Q.et al (2012). c-Myc suppression of DNA double-strand break repair.Neoplasia1411901202.

  • 35

    Lister-ShimauchiE. H.DinhM.MaddoxP.AhmedS. (2021). Gametes deficient for Pot1 telomere binding proteins alter levels of telomeric foci for multiple generations.Commun. Biol.4:158.

  • 36

    LittleJ. B. (2000). Radiation carcinogenesis.Carcinogenesis21397404. 10.1093/carcin/21.3.397

  • 37

    MahlkeM. A.Nechemia-ArbelyY. (2020). Guarding the genome: CENP-A-chromatin in health and cancer.Genes (Basel)11:810. 10.3390/genes11070810

  • 38

    MetheetrairutC.SlackF. J. (2013). MicroRNAs in the ionizing radiation response and in radiotherapy.Curr. Opin. Genet. Dev.231219. 10.1016/j.gde.2013.01.002

  • 39

    MiousseI. R.KutanziK. R.KoturbashI. (2017). Effects of ionizing radiation on DNA methylation: from experimental biology to clinical applications.Int. J. Radiat. Biol.93457469. 10.1080/09553002.2017.1287454

  • 40

    MorrisonJ. P.TonT. V.CollinsJ. B.SwitzerR. C.LittleP. B.MorganD. L.et al (2009). Gene expression studies reveal that DNA damage, vascular perturbation, and inflammation contribute to the pathogenesis of carbonyl sulfide neurotoxicity.Toxicol. Pathol.37502511. 10.1177/0192623309335631

  • 41

    MothersillC.SeymourC. B.O’BrienA. (1991). Induction of c-myc oncoprotein and of cellular proliferation by radiation in normal human urothelial cultures.Anticancer Res.1116091612.

  • 42

    MukherjeeA. K.SharmaS.SenguptaS.SahaD.KumarP.HussainT.et al (2018). Telomere length-dependent transcription and epigenetic modifications in promoters remote from telomere ends.PLoS Genet.14:e1007782. 10.1371/journal.pgen.1007782

  • 43

    ParisL.GiardulloP.LeonardiS.TannoB.MeschiniR.CordelliE.et al (2015). Transgenerational inheritance of enhanced susceptibility to radiation-induced medulloblastoma in newborn Ptch1+/– mice after paternal irradiation.Oncotarget63609836112. 10.18632/oncotarget.5553

  • 44

    ParkH.-J.LeeR.YooH.HongK.SongH. (2020). Nonylphenol induces apoptosis through ROS/JNK signaling in a spermatogonia cell line.Int. J. Mol. Sci.22:307. 10.3390/ijms22010307

  • 45

    PescatoriS.BerardinelliF.AlbanesiJ.AscenziP.MarinoM.AntocciaA.et al (2021). A tale of ice and fire: the dual role for 17β-estradiol in balancing DNA damage and genome integrity.Cancers (Basel).13:1583. 10.3390/cancers13071583

  • 46

    PhesseT. J.MyantK. B.ColeA. M.RidgwayR. A.PearsonH.MuncanV.et al (2014). Endogenous c-Myc is essential for p53-induced apoptosis in response to DNA damage in vivo.Cell Death Differ.21956966. 10.1038/cdd.2014.15

  • 47

    PogribnyI.KoturbashI.TryndyakV.HudsonD.StevensonS. M.SedelnikovaO.et al (2005). Fractionated low-dose radiation exposure leads to accumulation of DNA damage and profound alterations in DNA and histone methylation in the murine thymus.Mol. Cancer Res.3553561. 10.1158/1541-7786.mcr-05-0074

  • 48

    PujolR.CusidóL.RubioA.EgozcueJ.GarciaM. (1996). Effect of X-rays on germ cells in female fetuses of Rattus norvegicus irradiated at three different times of gestation.Mutat. Res. Mol. Mech. Mutagen.356247253. 10.1016/0027-5107(96)00067-x

  • 49

    PujolR.CusidóL.RubioA.EgozcueJ.GarciaM. (1997). X-ray-induced synaptonemal complex damage during meiotic prophase in female fetuses of Rattus norvegicus.Mutat. Res. Mol. Mech. Mutagen.379127134. 10.1016/s0027-5107(97)00115-2

  • 50

    RobinsonM. D.McCarthyD. J.SmythG. K. (2010). EdgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data.Bioinformatics.2613940. 10.1093/bioinformatics/btp616

  • 51

    RogakouE. P.PilchD. R.OrrA. H.IvanovaV. S.BonnerW. M. (1998). DNA double-stranded breaks induce histone H2AX phosphorylation on serine 139.J. Biol. Chem.27358585868. 10.1074/jbc.273.10.5858

  • 52

    RowleyM. J.LeachD. R.WarnerG. A.HellerC. G. (1974). Effect of graded doses of ionizing radiation on the human testis.Radiat. Res.59665678. 10.2307/3574084

  • 53

    SaitoK.CuiH. (2018). Emerging roles of estrogen-related receptors in the brain: potential interactions with estrogen signaling.Int. J. Mol. Sci.19:1091. 10.3390/ijms19041091

  • 54

    SaweyM. J.HoodA. T.BurnsF. J.GarteS. J. (1987). Activation of c-myc and c-K-ras oncogenes in primary rat tumors induced by ionizing radiation.Mol. Cell. Biol.7932935. 10.1128/mcb.7.2.932

  • 55

    SharmaU.SunF.ConineC. C.ReichholfB.KukrejaS.HerzogV. A.et al (2018). Small RNAs are trafficked from the epididymis to developing mammalian sperm.Dev. Cell46481.e494.e.

  • 56

    ShomeD.vonWoedtkeT.RiedelK.MasurK. (2020). The HIPPO transducer YAP and its targets CTGF and Cyr61 drive a paracrine signalling in cold atmospheric plasma-mediated wound healing.Oxid. Med. Cell. Longev.2020:4910280.

  • 57

    SkrzypekM.WdowiakA.PanasiukL.StecM.SzczygiełK.ZybałaM.et al (2019). Effect of ionizing radiation on the female reproductive system.Ann. Agric. Environ. Med.26606616.

  • 58

    SlongoM. L.ZampieriM.OnistoM. (2002). Expression of matrix metalloproteases (MMP-2, MT1 -MMP) and their tissue inhibitor (TIMP-2) by rat sertoli cells in culture: implications for spermatogenesis.Biol. Chem.383235239.

  • 59

    SzatmáriT.PersaE.KisE.BenedekA.HargitaiR.SáfrányG.et al (2019). Extracellular vesicles mediate low dose ionizing radiation-induced immune and inflammatory responses in the blood.Int. J. Radiat. Biol.951222. 10.1080/09553002.2018.1450533

  • 60

    TammingaJ.KoturbashI.BakerM.KutanziK.KathiriaP.PogribnyI. P.et al (2008). Paternal cranial irradiation induces distant bystander DNA damage in the germline and leads to epigenetic alterations in the offspring.Cell Cycle712381245. 10.4161/cc.7.9.5806

  • 61

    TanH. L.ZengY. B.ChenE. S. (2020). N-terminus does not govern protein turnover of Schizosaccharomyces pombe CENP-A.Int. J. Mol. Sci.21:6175. 10.3390/ijms21176175

  • 62

    UNSCEAR (2010). Sources and Effects of Ionizing Radiation. Available online at: http://www.unscear.org/unscear/en/publications/2008_1.html(accessed January 2021).

  • 63

    WallaceW. H. B.ThomsonA. B.KelseyT. W. (2003). The radiosensitivity of the human oocyte.Hum. Reprod.18117121. 10.1093/humrep/deg016

  • 64

    WardJ. F. (1988). DNA damage produced by ionizing radiation in mammalian cells: identities, mechanisms of formation, and reparability.Prog. Nucleic Acid Res. Mol. Biol.3595125. 10.1016/S0079-6603(08)60611-X

  • 65

    WardJ. F. (1995). Radiation mutagenesis: the initial DNA lesions responsible.Radiat. Res.142:362. 10.2307/3579145

  • 66

    WdowiakA.SkrzypekM.StecM.PanasiukL. (2019). Effect of ionizing radiation on the male reproductive system.Ann. Agric. Environ. Med.26210216. 10.26444/aaem/106085

  • 67

    ZhangG.-M.ZhangT. T.AnS. Y.El-SamahyM. A.YangH.WanY.et al (2019). Expression of Hippo signaling pathway components in Hu sheep male reproductive tract and spermatozoa.Theriogenology126239248. 10.1016/j.theriogenology.2018.12.029

  • 68

    ZhangL.YueT.JiangJ. (2009). Hippo signaling pathway and organ size control.Fly (Austin)36873. 10.4161/fly.3.1.7788

Summary

Keywords

environmental radiation, epigenetic, reproductive impairments, testicular, ovarian

Citation

Leung CT, Yang Y, Yu KN, Tam N, Chan TF, Lin X, Kong RYC, Chiu JMY, Wong AST, Lui WY, Yuen KWY, Lai KP and Wu RSS (2021) Low-Dose Radiation Can Cause Epigenetic Alterations Associated With Impairments in Both Male and Female Reproductive Cells. Front. Genet. 12:710143. doi: 10.3389/fgene.2021.710143

Received

31 May 2021

Accepted

07 July 2021

Published

02 August 2021

Volume

12 - 2021

Edited by

Pradyumna Kumar Mishra, ICMR-National Institute for Research in Environmental Health, India

Reviewed by

Kailash Manda, Institute of Nuclear Medicine & Allied Sciences (DRDO), India; Ravindra M. Samartha, Bhopal Memorial Hospital & Research Centre, India

Updates

Copyright

*Correspondence: Rudolf Shiu Sun Wu, Keng Po Lai,

This article was submitted to Toxicogenomics, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics