Because of imperfect discrimination against ribonucleoside triphosphates by the replicative DNA polymerases, large numbers of ribonucleotides are incorporated into the eukaryotic nuclear genome during S-phase. Ribonucleotides, by far the most common DNA lesion in replicating cells, destabilize the DNA, and an evolutionarily conserved DNA repair machinery, ribonucleotide excision repair (RER), ensures ribonucleotide removal. Whereas complete lack of RER is embryonically lethal, partial loss-of-function mutations in the genes encoding subunits of RNase H2, the enzyme essential for initiation of RER, cause the SLE-related type I interferonopathy Aicardi-Goutières syndrome. Here, we demonstrate that selective inactivation of RER in mouse epidermis results in spontaneous DNA damage and epidermal hyperproliferation associated with loss of hair follicle stem cells and hair follicle function. The animals developed keratinocyte intraepithelial neoplasia and invasive squamous cell carcinoma with complete penetrance, despite potent type I interferon production and skin inflammation. These results suggest that compromises to RER-mediated genome maintenance might represent an important tumor-promoting principle in human cancer.

Significance: Selective inactivation of ribonucleotide excision repair by loss of RNase H2 in the murine epidermis results in spontaneous DNA damage, type I interferon response, skin inflammation, and development of squamous cell carcinoma. Cancer Res; 78(20); 5917–26. ©2018 AACR.

Genome integrity is continuously challenged by a multitude of different hazards that cause a broad spectrum of different lesions in DNA (1). Spontaneous hydrolysis, ionizing radiation, exogenous chemicals, and endogenous metabolites, including reactive oxygen species, inflict cumulative damage to DNA (1). Fortunately, various DNA repair pathways revert most of this damage, thereby dramatically slowing down age-related deterioration of DNA integrity and suppressing the development of cancer (2). In addition to continuous, time-dependent harm to DNA, other forms of damage are introduced already during replication by errors the replicative polymerases make, despite their remarkable fidelity. Such replication errors include incorporation of deoxyribonucleotides (dNTP) that are not complementary to the template, which occurs every 104–106 bases, depending on the DNA polymerase. Most mismatched nucleotides are removed by the proofreading activity that some DNA polymerases possess, whereas the remaining mismatches are detected and corrected by DNA mismatch repair (MMR), which is directly coupled to the replication machinery (3–5). Postreplicative MMR is an important tumor-suppressive principle as genetic defects in MMR proteins cause cancer predisposition syndromes, including Lynch syndrome (6) and biallelic mismatch repair deficiency (BMMRD) syndrome (7).

Another form of replication error is the incorporation of nucleotides containing the correct base, but the wrong sugar. For example, discrimination between incoming dNTPs and ribonucleotides (rNTP) is imperfect for the three DNA polymerases that replicate the vast majority of the undamaged eukaryotic nuclear genome (8), resulting in frequent incorporation of rNMPs into genomic DNA (9, 10). The probability for this mistake is increased by the far greater cellular abundance of rNTPs compared with dNTPs (8, 11) such that on average, more than 106 rNMPs are incorporated during one round of replication of a mammalian genome (12). Newly incorporated rNMPs destabilize DNA (10) and pose a major threat to genome integrity due to their reactive 2′OH group. A highly conserved repair pathway, ribonucleotide excision repair (RER), ensures their removal (10, 13, 14). Like MMR, RER may be directly coupled to replication and results in rapid postreplicative repair of rNMPs. The first step in RER is detection of rNMPs embedded in the DNA double helix by ribonuclease H2 (RNase H2). RNase H2 then incises the rNMP-containing strand immediately 5′ of the rNMP. In vitro, this is followed by strand displacement synthesis by polδ or polϵ, flap cleavage by Fen1 or Exo1, and ligation of the nick (14). Interestingly, during their brief transient presence in DNA, rNMPs seem to serve important physiologic functions in DNA metabolism (8, 13, 15). They may provide a strand discrimination signal in MMR for the nascent leading strand (16, 17), wherein RNase H2-mediated incisions, occurring only in the newly synthesized rNMP-containing strand, may serve as entry points for MMR enzymes. Moreover, nicking by RNase H2 could be crucial for relaxing leading strand torsional stress (18). rNMPs might also be incorporated during repair synthesis occurring in noncycling cells that contain very high rNTP:dNTP ratios.

RNase H2-deficient cells are incapable of executing normal RER, carry high numbers of rNMPs in their DNA and spontaneously mount a DNA damage response (DDR), reflecting the DNA-destabilizing effects of rNMPs (13, 15, 19). rNMPs alter DNA structure and thereby potentially affect multiple key DNA transactions. rNMPs increase the rate of spontaneous hydrolysis and were shown to trigger a mutagenic repair pathway in yeast depending on DNA topoisomerase 1, which, like RNase H2, can incise at single rNMPs, potentially generating unligatable ends in repetitive DNA that result in short deletion mutations (20). Topoisomerase 1 cleavage at rNTPs embedded in the DNA can also result in double-strand breaks and replication stress ensuing from abnormally long-lived topoisomerase 1-DNA complexes (21). Moreover, bypassing rNMPs in DNA may pose a problem for replicative polymerases, which may result in replication fork stalling, strand breaks, and replicative stress. In addition to genome protection by RNase H2-mediated removal of single rNMPs, RNase H2 contributes to genome stability also by resolution of R-loops (22), preventing damage resulting from collisions of RNA polymerases or replication forks with R-loops.

Loss of RNase H2 causes genome instability in yeast with high rates of gross chromosomal rearrangements and copy number variations (13). Mammalian cells lacking the enzyme display increased numbers of strand breaks, activation of DNA damage responses (DDR), and chromosomal aberrations (12, 23). Abnormal mitotic chromosome segregation associated with this damage was shown to occur in RNase H2-deficient cells, leading to formation of micronuclei (24, 25). Similarly, DNA damage can lead to herniation of chromatin into the cytosol (26). Micronuclei and herniated chromatin expose large amounts of chromosomal DNA to the cytosolic DNA sensor cGAS that activates STING and thereby triggers an inflammatory response, including upregulation of IFN-stimulated genes (ISG) and inflammatory cytokine production (24). Unbalanced type I IFN responses are likely responsible for chronic inflammation in patients suffering from the rare monogenic autoimmune disorder Aicardi-Goutières syndrome (AGS), which can be caused by partial loss-of-function of RNase H2 or by mutations in other enzymes involved in nucleic acid metabolism (27, 28).

In mouse models, the effects of targeted inactivation of RNase H2 in vivo were found to critically depend on the degree to which RNase H2 activity was reduced. The RNase H2 complex is composed of three proteins, RNASEH2A, B, and C, all of which are essential for RNase H2 activity (29, 30). Biallelic null or hypomorphic mutations in RNase H2 genes that cause massive reduction of activity resulted in embryonic or perinatal lethality associated with spontaneous DNA damage and activation of DDRs (12, 23, 31). In contrast, a homozygous mutation reducing RNase H2 activity to about 30% of wild-type levels triggered a spontaneous low level IFN response, but was not associated with gross pathology (32).

To elucidate the consequences of severe RER deficiency for the homeostasis of a tissue with rapid cell turnover in adult animals, we bypassed embryonic lethality of ubiquitous RNase H2 deficiency by conditional gene inactivation. We demonstrate that complete loss of RNase H2 in skin epithelium results in spontaneous DNA damage and epidermal hyperproliferation that progresses to skin cancer with 100% penetrance, despite potent spontaneous type I IFN responses, demonstrating that RNase H2-dependent RER is an important tumor-suppressive principle.

Mice

Rnaseh2bFL/FL (23), K14-Cre (33), Ifnar1−/− (generated from Ifnar1FL mice (34) by ubiquitous deletion of the loxP-flanked fragment), and Trp53−/− (34) were either on the C57BL/6NJ background or backcrossed for at least five generations to C57BL/6NJ. Mice were housed at the Experimental Center, Medical Faculty Carl Gustav Carus, TU Dresden, under specific pathogen-free conditions. All procedures were in accordance with institutional guidelines on animal welfare and were approved by the Landesdirektion Dresden (permit number 24-1/2013-12, 88/2017).

Quantification of mRNAs by quantitative RT-PCR

Total RNA was isolated using the NucleoSpin Kit (Macherey-Nagel) according to the manufacturer's instructions. cDNA was generated using the RevertAid H Minus First Strand cDNA Synthesis Kit (Thermo Fisher Scientific). Quantification of transcripts was performed on a Mx3005P QPCR System (Agilent Technologies) using the Maxima SYBR Green/ROX qPCR Master Mix (Thermo Fisher Scientific). The following oligos were used: Oasl1 up, 5′-GCAATCCACAGCGATATCC-3′; Oasl1 down, 5′-CAACTGCTCACTGTCCACGG-3′; Ifi44 up, 5′-GAGTCACTCATTCTCGGACTCCGC-3′; Ifi44 down, 5′-GAGCGGGCATTGAAGTAAGGGC-3′; Viperin up, 5′-TTCGCCCGCATCAAAGCCGT-3′; Viperin down, 5′-AGGGGGCAGCGGAAGTCGAT-3′; Rnaseh2b up, 5′-AGGTTTCCAGGGACAAGGAAGAGGA-3′; Rnaseh2b down, 5′-GTCAATGAAGCTGGAGGTTCTGGAAG-3′; Tbp1 up, 5′-TGACCCAGATCATGTTTGAGACCTTCA-3′; and Tbp1 down, 5′-GGAGTCCATCACAATGCCTGTGG-3′. All samples were run in triplicates. qRT-PCR data are displayed as fold change compared with means of the respective control groups ± SD.

Assay for RNase H2 activity

Cell lysates were prepared and assayed for specific cleavage of an 18-bp double-stranded DNA substrate containing a single ribonucleotide in one strand as described previously (27, 29). RNase H2–specific activity was determined by subtracting the cellular activity against a sequence-matched DNA duplex without rNTPs. Cell lysate protein concentration was determined and lysates were added to the reaction mix at a final protein concentration of 100 ng/μL.

Isolation and flow cytometric analysis of epidermal cells

Back skin was excised and adipose tissue was removed. The skin was placed on a layer of trypsin solution (0.25%, Life Technologies) with the epidermal side facing the trypsin solution. After incubation for 2 hours, the epidermis was separated from the dermis using blunt forceps. The epidermis was minced with scalpels and digested in 10 mL trypsin solution (0.25%, Life Technologies) for another hour. The digest was stopped by adding an equal volume of medium [DMEM (2/3) + HAM's F12 (1/3)] supplemented with 10% calcium-free FBS and the cells were strained through a 40-μm nylon mesh. Cells were then centrifuged for 8 minutes at 500× g at room temperature and the pellet was resuspended in 2 mL FACS buffer. Epidermal stem cell populations were stained according to Jensen and colleagues (35). For flow cytometry, cells were stained using the following antibodies: anti–CD45-FITC (1:400, eBioscience), anti–CD49f-PE (1:200, eBioscience), anti–CD34-eF660 (1:50, eBioscience), anti–Sca1-PerCP-Cy5.5 (1:200, eBioscience), and anti–CD117-APC-Cy7 (1:800, BioLegend). Cell sorting and analysis was performed on an ARIA III cell sorter (BD Biosciences) and data were recorded using DIVA software (BD Biosciences) and analyzed using FlowJo (FlowJo).

Transcriptome analysis

Keratinocytes (CD49f+) and hematopoietic cells (CD45+) were isolated from epidermal skin cell suspensions by FACS (see above). Total RNA was isolated using the RNeasy Mini Kit+ (Qiagen). mRNA libraries were prepared and subjected to deep sequencing on an IlluminaHighSeq. Reads were mapped to the reference genome mm10 (Ensembl Version 75) using gsnap (v.2014-12-06 for analysis of CD45CD49f+ cells, v.2014-12-17 for analysis of CD45+ cells) and genes were counted with featurecount v.1.4.6. Differentially expressed genes (DEG) were identified using DESeq2 v1.8.1 (36). Clustering of DEGs into specific pathways was investigated using Kyoto Encyclopedia of Genes and Genomes (http://www.ge-nome.jp/kegg/pathway.html), Reactome (https://reactome.org/) and Interferome V2.0 (www.interferome.org/) databases as well as the Ingenuity Pathway Analysis (IPA) software (Qiagen). RNA-seq data reported in this article are accessible as a super series from the Gene Expression Omnibus database under accession number GSE115005.

Histologic analysis

Formalin-fixed and paraffin-embedded skin sections were deparaffinized and rehydrated, subjected to antigen retrieval (20 minutes 98°C in sodium citrate pH 6.0) and washed three times in TBST (1xTBS, 0.1% Triton X-100).

For quantification of phosphorylated histone H2A.X (γH2A.X), sections were blocked for 1 hour (10% goat serum in TBST) and incubated at 4°C overnight with a phospho-histone H2A.X (pSer139) antibody [Cell Signaling Technology, 1:50 in 1% goat serum/tris-buffered saline with 0.05% Tween20 (TBST)]. Next, sections were washed (TBST) and incubated with a goat anti–rabbit-AF488 antibody (1:500, Thermo Fisher Scientific) for 4 hours at room temperature in the dark. After washing, nuclei were counterstained with DAPI in the mounting solution. Images were recorded on a Zeiss Axiovert ApoTome II and analyzed using Zen software (Zeiss). γH2A.X foci in at least 220 keratinocyte nuclei of the interfollicular epidermis per section were counted.

For quantification of Ki67-positive cells, endogenous peroxidase activity was blocked (3% H2O2, 20 minutes in the dark), unspecific epitopes were blocked (PBS 10% BSA, 0.5% Triton X-100 for 1 hour at room temperature) and sections were washed and incubated at 4°C overnight with anti-Ki67-Biotin (eBioscience, 1:100 in PBS 1% BSA, 0.05% Triton X-100). After washing, sections were incubated for 1 hour at room temperature with horseradish peroxidase–conjugated Streptavidin (Dako, 1:300, PBS, 1% BSA, 0.05% Triton X-100). After addition of DAB mix (Vector Labs) and counterstaining with hematoxylin, Ki67-positive keratinocytes were counted in interfollicular epidermis.

For quantification of apoptotic cells positive for activated caspase-3, antigen retrieval was performed in PBS with 1 mmol/L EDTA, 0.05% Tween 20, pH 8.0 for 20 minutes in a pressure cooker), sections were blocked (PBS, 0.1% Tween 20, 5% normal goat serum), and incubated with rabbit anti–active-caspase-3 AF835 (R&D, 1:600) at 4°C overnight, followed by washing and incubation for 1 hour at room temperature with goat anti-rabbit AF488 (ab150077, Abcam, 1:500). Nuclei were stained with DAPI. Images were recorded on a Keyence BZ-X710 microscope and positive cells were counted in interfollicular epidermis and follicular infundibula.

Histologic evaluation was performed by a professional dermato-histopathologist (J. Wenzel, Bonn, Germany) in a blinded fashion.

Statistical analysis

Unless stated otherwise, significance was calculated by unpaired, two-sided Student t test, ***, P < 0.001; **, P < 0.01; *, P < 0.05.

Complete loss of RNase H2 in the epidermis results in epithelial hyperproliferation and loss of stem cells

We inactivated the Rnaseh2b gene selectively in the epidermis by crossing Rnaseh2bFLOX mice (23) to the K14-Cre line that deletes loxP-flanked DNA in all basal cells of skin epithelium including hair follicles (23, 33). Rnaseh2bFL/FLK14-Cre mice (“Rnaseh2EKO” mice) showed strong reduction of Rnaseh2b transcripts and RNase H2 activity in total epidermis or FACS-purified epidermal cells, as expected (Supplementary Fig. S1A–S1C). The animals were conspicuous already few days after birth, showing significant hyperpigmentation, most prominently of ears, snout, tails, and paws (Fig. 1A), a known sign of ongoing DNA damage responses in mouse epidermis (37, 38). Although whiskers were hypomorphic and reduced in number already in young mice (Fig. 1A), the fur initially appeared macroscopically normal. However, Rnaseh2EKO mice began to lose hair at about 12 weeks and were almost completely nude by 20 weeks (Fig. 1A).

Figure 1.

Epithelial hyperproliferation and stem cell exhaustion in mice lacking RNase H2 in the epidermis (Rnaseh2EKO). A, Macroscopic phenotype of a 15-week-old Rnaseh2EKO mouse and a control (Rnaseh2bFL/FLCre-negative) littermate. B, Representative histology of back skin from a 12-week-old Rnaseh2EKO and a control (Rnaseh2bFL/FLCre-negative) littermate. Left, hematoxylin-eosin (H&E) staining showing focal epithelial thickening and hyperkeratosis (thickening of cornified layer) of mutant epidermis (see Supplementary Material for detailed histopathologic findings); right, Ki67 immunostaining showing increased proliferation of mutant epidermis. Scale bars, 40 μm. Graph shows numbers of Ki67+ interfollicular epidermal cells per high power field (40×); mean ± SD; *, P = 0.014 (5 weeks); **, P = 0.0039 (10–12 weeks); ***, P = 0.0007 (18–21 weeks); *, P = 0.0101 (30 weeks). C, Detection of apoptotic cells in interfollicular epidermis by immunostaining for active caspase-3 in back skin sections of Rnaseh2EKO mice and control (Rnaseh2bFL/FLCre-negative) littermates. Representative images are shown. Scale bar, 50 μm. D, Quantification of active caspase-3 per mm epidermal length (n = 6–10/group and timepoint); means ± SD, **, P = 0.0073 (10–13 weeks); ***, P < 0.0001 (18–21 weeks). E, Quantification of “active” hair follicles (i.e., follicles extending beyond the dermal collagen into the dermal adipose layer, thus including follicles in anagen II–VI and early catagen stages) per low power (10×) field in hematoxylin and eosin–stained back skin sections of Rnaseh2EKO mice and control (Rnaseh2bFL/FLCre-negative) littermates. Means ± SD; *, P = 0.0391 (10–13 weeks); *, P = 0.0335 (18–21 weeks). F, Flow cytometric quantification of epithelial stem cell populations in back skin of Rnaseh2EKO mice and control (Rnaseh2bFL/FLCre-negative) littermates. See Supplementary Fig. S1 for gating. SCs, stem cells; HFSC, hair follicle stem cells; IFE, interfollicular epidermis. Mean and SD; **, P = 0.0011 (suprabasal bulge SCs, 20 weeks); *, P = 0.0205 (bulge SCs, 10–12 weeks); **, P = 0.0033 (bulge SCs, 20 weeks); *, P = 0.0239 (IFE+ infundibulum SCs, 20 weeks).

Figure 1.

Epithelial hyperproliferation and stem cell exhaustion in mice lacking RNase H2 in the epidermis (Rnaseh2EKO). A, Macroscopic phenotype of a 15-week-old Rnaseh2EKO mouse and a control (Rnaseh2bFL/FLCre-negative) littermate. B, Representative histology of back skin from a 12-week-old Rnaseh2EKO and a control (Rnaseh2bFL/FLCre-negative) littermate. Left, hematoxylin-eosin (H&E) staining showing focal epithelial thickening and hyperkeratosis (thickening of cornified layer) of mutant epidermis (see Supplementary Material for detailed histopathologic findings); right, Ki67 immunostaining showing increased proliferation of mutant epidermis. Scale bars, 40 μm. Graph shows numbers of Ki67+ interfollicular epidermal cells per high power field (40×); mean ± SD; *, P = 0.014 (5 weeks); **, P = 0.0039 (10–12 weeks); ***, P = 0.0007 (18–21 weeks); *, P = 0.0101 (30 weeks). C, Detection of apoptotic cells in interfollicular epidermis by immunostaining for active caspase-3 in back skin sections of Rnaseh2EKO mice and control (Rnaseh2bFL/FLCre-negative) littermates. Representative images are shown. Scale bar, 50 μm. D, Quantification of active caspase-3 per mm epidermal length (n = 6–10/group and timepoint); means ± SD, **, P = 0.0073 (10–13 weeks); ***, P < 0.0001 (18–21 weeks). E, Quantification of “active” hair follicles (i.e., follicles extending beyond the dermal collagen into the dermal adipose layer, thus including follicles in anagen II–VI and early catagen stages) per low power (10×) field in hematoxylin and eosin–stained back skin sections of Rnaseh2EKO mice and control (Rnaseh2bFL/FLCre-negative) littermates. Means ± SD; *, P = 0.0391 (10–13 weeks); *, P = 0.0335 (18–21 weeks). F, Flow cytometric quantification of epithelial stem cell populations in back skin of Rnaseh2EKO mice and control (Rnaseh2bFL/FLCre-negative) littermates. See Supplementary Fig. S1 for gating. SCs, stem cells; HFSC, hair follicle stem cells; IFE, interfollicular epidermis. Mean and SD; **, P = 0.0011 (suprabasal bulge SCs, 20 weeks); *, P = 0.0205 (bulge SCs, 10–12 weeks); **, P = 0.0033 (bulge SCs, 20 weeks); *, P = 0.0239 (IFE+ infundibulum SCs, 20 weeks).

Close modal

Histologically, Rnaseh2EKO epidermis was normal at the age of 5 weeks except for discrete hyperkeratosis (thickening of the cornified layer). In older animals (10–13 weeks and 18–21 weeks), focal thickening of the epithelium and moderate hyperkeratosis suggested increased epithelial proliferation (Fig. 1B). Immunostaining for the cell-cycle marker Ki67 (Fig. 1B) revealed increased numbers of proliferating epithelial cells in Rnaseh2EKO compared with control skin in all age groups tested. More rapid proliferation was paralleled by a strong increase in the frequency of apoptotic keratinocytes as detected by staining for active caspase-3 (Fig. 1C and D). Loss of hair follicles did not account for the almost complete boldness of Rnaseh2EKO mice as substantial numbers of follicles were still present even at 30 weeks of age (Supplementary Fig. S1D). In their active (anagen) phase, hair follicles grow to extend well into the dermal adipose tissue paralleled by thickening of the adipose tissue, whereas resting (telogen) follicles are usually confined to the dermal collagen layer. Numbers of anagen follicles were not different from control numbers at 5 weeks of age. At 10–13 weeks, we found no follicles extending into the adipose layer in control skin, indicating a synchronous resting phase of the hair cycle in this part of the body. In contrast, all Rnaseh2EKO mice featured numerous “active” follicles at this timepoint, demonstrating disturbed hair cycle regulation (Fig. 1E). This was also reflected by increased thickness of Rnaseh2EKO skin at 10 to 13 weeks as determined by skin fold measurements (Supplementary Fig. S1E). Loss of hair follicle function in older Rnaseh2EKO mice suggested exhaustion of epithelial regenerative capacity. In accordance with this notion, we found that the number of hair follicle stem cells, in particular bulge stem cells (CD49f+CD34+Sca1lo) was strongly reduced in older mutants compared with controls, as determined by flow cytometric analysis of epidermal single cell suspensions (Fig. 1F; Supplementary Fig. S1F). In contrast, other undifferentiated epidermal cell populations were not reduced in mutant skin (Fig. 1F; Supplementary Fig. S1F).

Collectively, loss of RNase H2 in the epidermis results in more rapid epithelial cell turnover due to enhanced epithelial cell proliferation and apoptosis, associated with reduction of hair follicle stem cell numbers and loss of hair follicle function.

Spontaneous type I IFN response and inflammation in Rnaseh2EKO skin

Histology and flow cytometry revealed enhanced leukocyte infiltration of Rnaseh2EKO skin (Fig. 2A and B), which was accentuated around hair follicles (Fig. 2A). Inflammatory cell invasion into the epithelium was associated with degeneration of the outermost epithelial cells of the hair follicle (Fig. 2A), an inflammation pattern (“interface dermatitis”) typical of skin lesions in SLE (39).

Figure 2.

Loss of RNase H2 in the epidermis results in spontaneous type I IFN response and inflammation. A, Representative Rnaseh2EKO back skin section (age 10–12 weeks) immunostained for CD45 (purple) for identification of hematopoietic cells (see Supplementary Materials and Methods for detailed histopathologic findings). Original magnification, ×200. B, Quantification of CD45+ leukocytes in single cell suspensions prepared from the epidermal compartment of back skin from 10–12 week-old Rnaseh2EKO (EKO) mice (n = 10) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 9). Mean and SD; ***, P < 0.0001. C, Comparison of ISG transcript levels determined by qRT-PCR in RNA extracted from total epidermis of 10–12 and 20-week-old Rnaseh2EKO (EKO) mice to the mean obtained for control (Rnaseh2bFL/FLCre-negative) littermates (10–12 wk, n = 6; 20 wk, n = 3), which was set to 1 (dotted line).10–12 weeks: **, P = 0.0018 (Oasl1); *, P = 0.0105 (Ifi44); *, P = 0.0162 (Viperin); 20 weeks: *, P = 0.0242 (Oasl1); **, P = 0.0061 (Ifi44); *, P = 0.0246 (Viperin). D, Comparison of gene expression profiles of FACS-purified CD45+ leukocytes from epidermis of 12-week-old Rnaseh2EKO (EKO) mice (n = 7) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 6) based on RNA sequencing. Black, significantly (P < 0.05) deregulated genes; red, significantly deregulated ISGs. E, Macroscopic phenotype of an 11-week-old Ifnar1−/−Rnaseh2EKO mouse (note hyperpigmentation of ears, tail, and paws, and almost complete hair loss) and a control (Ifnar1−/−Rnaseh2bFL/FLCre-negative) littermate. F, Quantification of CD45+ leukocytes in the epidermis of 10–12 week-old Ifnar1−/−Rnaseh2EKO (EKO) mice (n = 7) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 7). Mean and SD; ***, P = 0.0009. G, Flow cytometric quantification of hair follicle stem cell populations in back skin of 10–12 week-old Ifnar1−/−Rnaseh2EKO mice (n = 7) and control (Ifnar1−/−Rnaseh2bFL/FLCre-negative) littermates (n = 7). See Supplementary Fig. S1F for gating. Mean and SD; *, P = 0.0209; ***, P < 0.0001.

Figure 2.

Loss of RNase H2 in the epidermis results in spontaneous type I IFN response and inflammation. A, Representative Rnaseh2EKO back skin section (age 10–12 weeks) immunostained for CD45 (purple) for identification of hematopoietic cells (see Supplementary Materials and Methods for detailed histopathologic findings). Original magnification, ×200. B, Quantification of CD45+ leukocytes in single cell suspensions prepared from the epidermal compartment of back skin from 10–12 week-old Rnaseh2EKO (EKO) mice (n = 10) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 9). Mean and SD; ***, P < 0.0001. C, Comparison of ISG transcript levels determined by qRT-PCR in RNA extracted from total epidermis of 10–12 and 20-week-old Rnaseh2EKO (EKO) mice to the mean obtained for control (Rnaseh2bFL/FLCre-negative) littermates (10–12 wk, n = 6; 20 wk, n = 3), which was set to 1 (dotted line).10–12 weeks: **, P = 0.0018 (Oasl1); *, P = 0.0105 (Ifi44); *, P = 0.0162 (Viperin); 20 weeks: *, P = 0.0242 (Oasl1); **, P = 0.0061 (Ifi44); *, P = 0.0246 (Viperin). D, Comparison of gene expression profiles of FACS-purified CD45+ leukocytes from epidermis of 12-week-old Rnaseh2EKO (EKO) mice (n = 7) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 6) based on RNA sequencing. Black, significantly (P < 0.05) deregulated genes; red, significantly deregulated ISGs. E, Macroscopic phenotype of an 11-week-old Ifnar1−/−Rnaseh2EKO mouse (note hyperpigmentation of ears, tail, and paws, and almost complete hair loss) and a control (Ifnar1−/−Rnaseh2bFL/FLCre-negative) littermate. F, Quantification of CD45+ leukocytes in the epidermis of 10–12 week-old Ifnar1−/−Rnaseh2EKO (EKO) mice (n = 7) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 7). Mean and SD; ***, P = 0.0009. G, Flow cytometric quantification of hair follicle stem cell populations in back skin of 10–12 week-old Ifnar1−/−Rnaseh2EKO mice (n = 7) and control (Ifnar1−/−Rnaseh2bFL/FLCre-negative) littermates (n = 7). See Supplementary Fig. S1F for gating. Mean and SD; *, P = 0.0209; ***, P < 0.0001.

Close modal

Because RNase H2 deficiency was shown to be associated with spontaneous activation of type I IFN responses (31, 32), we quantified transcript levels of type I IFN-inducible genes in Rnaseh2EKO and control skin by qRT-PCR analysis. ISG mRNAs were strongly increased in total epidermis and in purified keratinocytes from young and older mutant mice (Fig. 2C; Supplementary Fig. S2A). Comparison of transcriptomes of leukocytes from control versus Rnaseh2EKO skin showed upregulation of multiple ISGs in the cells from the mutants, indicating exposure to type I IFN (Fig. 2D; Supplementary Table S1).

Type I IFNs were shown to regulate stem cell function in the hematopoietic system (40). We therefore addressed what role the chronic IFN response in Rnaseh2EKO epidermis played and crossed Rnaseh2EKO mice to a type I IFN receptor knock out (Ifnar1−/−) line (41; Supplementary Fig. S2B). Hyperpigmentation, hair loss (Fig. 2E), hyperkeratosis and focal epidermal thickening (Supplementary Fig. S2C), and inflammatory leukocyte infiltration (Fig. 2F) of Rnaseh2EKO skin were not mitigated by additional inactivation of type I IFN signaling. Rnaseh2EKOIfnar1−/− mice showed a reduction of hair follicle stem cell numbers compared with controls (Fig. 2G), similar to IFNAR-competent Rnaseh2EKO mice (Fig. 1E).

Collectively, the skin inflammation of Rnaseh2EKO mice is associated with type I IFN production; however, leukocyte infiltration, epidermal hyperproliferation, and loss of hair follicle function occur independent of type I IFN.

Loss of RNase H2 in the epidermis results in spontaneous DNA damage and skin cancer

Spontaneous DNA damage and activation of DNA damage responses was shown in RNase H2-deficient embryos and cells (12, 23). Moreover, the hyperpigmented skin of Rnaseh2EKO mice is suggestive of chronic DNA damage responses in the epidermis (37, 38). Immunostaining of skin sections for γH2AX and 53BP1 repair foci indeed revealed increased numbers of γH2AX and 53BP1 foci in Rnaseh2EKO compared with control epidermal cells (Fig. 3A; Supplementary Fig. S3A), consistent with increased frequencies of double-strand breaks in RNase H2-deficient yeast (21). In line with this finding, transcriptome analysis of flow cytometrically purified Rnaseh2EKO keratinocytes showed upregulation of several p53-inducible genes compared with control cells (Fig. 3B; Supplementary Table S2), indicating ongoing DNA damage responses.

Figure 3.

Loss of RNase H2 in the epidermis results in spontaneous DNA damage and skin cancer. A, Immunostaining of γH2AX repair foci (green) in back skin of 10–13 week-old Rnaseh2EKO (EKO) mice (n = 7) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 6). Left, representative result for EKO and control (DAPI counterstaining). Scale bars, 5 μm. Right, quantification of γH2AX foci per nucleus. At least 220 nuclei (interfollicular epidermis) were counted per animal. Mean and SD. ***, P < 0.0001 (one-way ANOVA). B, Comparison of gene expression profiles of FACS-purified keratinocytes from epidermis of 9–12 week-old Rnaseh2EKO (EKO) mice (n = 3) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 3) based on RNA sequencing. Black, significantly (P < 0.05) deregulated genes; red, significantly deregulated p53-inducible genes. C, Incidence of macroscopic lesions (ulcerations and tumors) in Rnaseh2EKO mice. Graph indicates the timepoints at which euthanasia was required for these lesions. D, Representative histology of the lesions represented in C, Top, hematoxylin and eosin–stained section of Rnaseh2EKO skin adjacent to ulceration showing characteristic features of KIN, that is, cancer that has not yet broken through the basement membrane. Scale bars, 200 μm, left; 100 μm, right. Ki67 immunostaining (bottom right) demonstrates proliferating Ki67+ cells in all strata of the epithelium as is characteristic of KIN. Scale bar, 40 μm. Macroscopic tumor in the genital area of a 36-week-old Rnaseh2EKO mouse histologically (bottom left) proved to be invasive SCC (Supplementary Table S3, animal 62166). Scale bar, 200 μm. See Supplementary Materials and Methods for detailed histopathologic findings. E, Incidence of invasive SCC and incidence of cancer of all stages (KIN and/or SCC) in Rnaseh2EKO mice. *, P = 0.0131 (log-rank test).

Figure 3.

Loss of RNase H2 in the epidermis results in spontaneous DNA damage and skin cancer. A, Immunostaining of γH2AX repair foci (green) in back skin of 10–13 week-old Rnaseh2EKO (EKO) mice (n = 7) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 6). Left, representative result for EKO and control (DAPI counterstaining). Scale bars, 5 μm. Right, quantification of γH2AX foci per nucleus. At least 220 nuclei (interfollicular epidermis) were counted per animal. Mean and SD. ***, P < 0.0001 (one-way ANOVA). B, Comparison of gene expression profiles of FACS-purified keratinocytes from epidermis of 9–12 week-old Rnaseh2EKO (EKO) mice (n = 3) and control (Rnaseh2bFL/FLCre-negative) littermates (n = 3) based on RNA sequencing. Black, significantly (P < 0.05) deregulated genes; red, significantly deregulated p53-inducible genes. C, Incidence of macroscopic lesions (ulcerations and tumors) in Rnaseh2EKO mice. Graph indicates the timepoints at which euthanasia was required for these lesions. D, Representative histology of the lesions represented in C, Top, hematoxylin and eosin–stained section of Rnaseh2EKO skin adjacent to ulceration showing characteristic features of KIN, that is, cancer that has not yet broken through the basement membrane. Scale bars, 200 μm, left; 100 μm, right. Ki67 immunostaining (bottom right) demonstrates proliferating Ki67+ cells in all strata of the epithelium as is characteristic of KIN. Scale bar, 40 μm. Macroscopic tumor in the genital area of a 36-week-old Rnaseh2EKO mouse histologically (bottom left) proved to be invasive SCC (Supplementary Table S3, animal 62166). Scale bar, 200 μm. See Supplementary Materials and Methods for detailed histopathologic findings. E, Incidence of invasive SCC and incidence of cancer of all stages (KIN and/or SCC) in Rnaseh2EKO mice. *, P = 0.0131 (log-rank test).

Close modal

Starting from week 12, Rnaseh2EKO mice develop chronic ulcerations of the skin, most frequently in the dorsal neck area, affecting most animals by the age of 1 year (Fig. 3C). All Rnaseh2EKO mice had to be sacrificed between 23 and 55 weeks of age because of these ulcerations, and/or because of large tumors in various locations, often in the mandibular region close to the ear (Fig. 3C). Histology revealed (Fig. 3D; Supplementary Table S3) that all ulcerations occurred in neoplastic skin that was classified keratinocyte intraepithelial neoplasia (KIN; ref. 42), that is, malignant epidermal growth that has not yet breached the epithelial basement membrane, also called “carcinoma in situ”. In three of these cases, the tumor had focally broken through the basement membrane and was thus classified invasive squamous cell carcinoma (SCC). All macroscopic tumors histologically proved to be invasive SCC, except for two, which were KIN producing large masses of cornified material (Fig. 3D; Supplementary Table S3). Figure 3E shows that 100% of Rnaseh2EKO mice developed cancer of at least the KIN stage by about 50 weeks of age. By this time, 60% of the animals had already progressed to invasive SCC (Fig. 3E). Additional inactivation of type I IFN signaling resulted in a slight (albeit insignificant) acceleration of neoplastic skin ulceration (Supplementary Fig. S3B), most likely reflecting loss of stimulation of antitumor immunity by type I IFN.

Collectively, we show that RNase H2-deficient hyperproliferative epidermis features spontaneous DNA damage as demonstrated by increased numbers of repair foci and increased transcript levels of p53-inducible genes, and progresses to cancer in 100% of the animals within the first year of life.

Additional loss of p53 in Rnaseh2EKO mice enhances the epidermal IFN response and accelerates hyperproliferation and carcinogenesis

To determine the effect of p53-dependent DNA damage responses on the Rnaseh2EKO phenotype, we crossed Rnaseh2EKO to Trp53−/− mice (34) and observed potent effects of Trp53 gene dose. Rnaseh2EKO mice heterozygous for the Trp53null allele showed ameliorated hyperpigmentation compared with p53-competent controls (Supplementary Fig. S4A) and no sign of alopecia until the age of 40 weeks. All Rnaseh2EKO Trp53+/− mice that were allowed to age had to be sacrificed because of macroscopic ulcerations or skin tumors before the age of 40 weeks. Histology showed that these lesions were invasive SCC of various stages and grades (Supplementary Table S4; Fig. 4A).

Figure 4.

Additional loss of p53 in Rnaseh2EKO mice enhances the epidermal IFN response and accelerates hyperproliferation and cancerogenesis. A, Incidence of invasive SCC (left) and of cancer of all stages (KIN and/or SCC; right) in Rnaseh2EKO mice that were Trp53WT/− or Trp53−/−. Cancer incidence of Trp53WT/WT Rnaseh2EKO mice from Fig. 3E is shown in light gray for comparison. ***, P < 0.0001, both (log-rank test). B, Macroscopic phenotype of a 10-week-old Trp53−/−Rnaseh2EKO mouse and an 11-week-old control (Trp53−/−Rnaseh2bFL/FLCre-negative). Note normal pigmentation (tail and paws) of the mutant, in contrast to Trp53wtwtRnaseh2EKO mice (Fig. 1A). Except for an inflamed, hyperkeratotic lesion in the neck with circumscribed hair loss, the fur is largely normal. C, Representative hematoxylin and eosin–stained neck skin section of a Trp53−/−Rnaseh2EKO mouse sacrificed for severe itching and erosions. Note typical features of KIN, including massive thickening of epidermis and severe dysplasia of the epithelium (irregular nuclear shape; numerous mitotic figures also in upper strata). See Supplementary Materials and Methods for detailed histopathologic findings. Scale bar, left, 50 μm; inset, 20 μm. D, Immunostaining of neck skin sections of 10–12 week-old Trp53−/−Rnaseh2EKO mice (n = 4) and control (Trp53−/−Rnaseh2bFL/FLCre-negative) littermates (n = 8) for Ki67. Representative image of Trp53−/−Rnaseh2EKO skin. Scale bar, 50 μm. Graph shows numbers of Ki67+ epidermal cells per randomly placed high power (40×) field. Mean and SD. ***, P < 0.0001. E, Quantification of γH2AX repair foci in immunostained back skin of 10–13 week-old p53−/−Rnaseh2EKO (EKO) mice (n = 4) and control (Trp53−/−Rnaseh2bFL/FLCre-negative) littermates (n = 4). Number of foci was counted in at least 220 nuclei (interfollicular epidermis and follicles) per animal. See Supplementary Fig. S4B for images. Mean and SD. ***, P < 0.0001 (one-way ANOVA). F, Flow cytometric quantification of epithelial stem cell populations in back skin of 10–12 week-old Trp53−/−Rnaseh2EKO mice (n = 6) and control (Trp53−/−Rnaseh2bFL/FLCre-negative) littermates (n = 5). See Supplementary Fig. S1F for gating. Mean and SD. ***, P = 0.0001. G, Comparison of ISG transcript levels determined by qRT-PCR in RNA extracted from total epidermis (left) or FACS-purified keratinocytes (right) of 10–12 week-old Trp53+/+Rnaseh2EKO and Trp53−/−Rnaseh2EKO mice with the mean obtained for the respective control group, which was set to 1. Left, ***, P = 0.0005 (Oasl1), ***, P < 0.0001 (Ifi44), **, P = 0.0066 (Viperin); right, **, P = 0.0042 (Oasl1; p53 deficiency alone did not cause a spontaneous IFN response; see Supplementary Fig. S4D.)

Figure 4.

Additional loss of p53 in Rnaseh2EKO mice enhances the epidermal IFN response and accelerates hyperproliferation and cancerogenesis. A, Incidence of invasive SCC (left) and of cancer of all stages (KIN and/or SCC; right) in Rnaseh2EKO mice that were Trp53WT/− or Trp53−/−. Cancer incidence of Trp53WT/WT Rnaseh2EKO mice from Fig. 3E is shown in light gray for comparison. ***, P < 0.0001, both (log-rank test). B, Macroscopic phenotype of a 10-week-old Trp53−/−Rnaseh2EKO mouse and an 11-week-old control (Trp53−/−Rnaseh2bFL/FLCre-negative). Note normal pigmentation (tail and paws) of the mutant, in contrast to Trp53wtwtRnaseh2EKO mice (Fig. 1A). Except for an inflamed, hyperkeratotic lesion in the neck with circumscribed hair loss, the fur is largely normal. C, Representative hematoxylin and eosin–stained neck skin section of a Trp53−/−Rnaseh2EKO mouse sacrificed for severe itching and erosions. Note typical features of KIN, including massive thickening of epidermis and severe dysplasia of the epithelium (irregular nuclear shape; numerous mitotic figures also in upper strata). See Supplementary Materials and Methods for detailed histopathologic findings. Scale bar, left, 50 μm; inset, 20 μm. D, Immunostaining of neck skin sections of 10–12 week-old Trp53−/−Rnaseh2EKO mice (n = 4) and control (Trp53−/−Rnaseh2bFL/FLCre-negative) littermates (n = 8) for Ki67. Representative image of Trp53−/−Rnaseh2EKO skin. Scale bar, 50 μm. Graph shows numbers of Ki67+ epidermal cells per randomly placed high power (40×) field. Mean and SD. ***, P < 0.0001. E, Quantification of γH2AX repair foci in immunostained back skin of 10–13 week-old p53−/−Rnaseh2EKO (EKO) mice (n = 4) and control (Trp53−/−Rnaseh2bFL/FLCre-negative) littermates (n = 4). Number of foci was counted in at least 220 nuclei (interfollicular epidermis and follicles) per animal. See Supplementary Fig. S4B for images. Mean and SD. ***, P < 0.0001 (one-way ANOVA). F, Flow cytometric quantification of epithelial stem cell populations in back skin of 10–12 week-old Trp53−/−Rnaseh2EKO mice (n = 6) and control (Trp53−/−Rnaseh2bFL/FLCre-negative) littermates (n = 5). See Supplementary Fig. S1F for gating. Mean and SD. ***, P = 0.0001. G, Comparison of ISG transcript levels determined by qRT-PCR in RNA extracted from total epidermis (left) or FACS-purified keratinocytes (right) of 10–12 week-old Trp53+/+Rnaseh2EKO and Trp53−/−Rnaseh2EKO mice with the mean obtained for the respective control group, which was set to 1. Left, ***, P = 0.0005 (Oasl1), ***, P < 0.0001 (Ifi44), **, P = 0.0066 (Viperin); right, **, P = 0.0042 (Oasl1; p53 deficiency alone did not cause a spontaneous IFN response; see Supplementary Fig. S4D.)

Close modal

Biallelic loss of Trp53 completely reverted the hyperpigmentation of Rnaseh2EKO (Fig. 4B; Supplementary Fig. S4A). Although no generalized alopecia was observed, all of these animals developed inflamed lesions associated with circumscribed loss of hair and extensive scratching in the neck area, starting at about 12 weeks of age (Fig. 4B). Massive pruritus and erosions of the neck skin invariably required euthanasia few weeks later. In all cases, histology identified KIN of different stages as the cause of the lesions (Fig. 4C; Supplementary Tables S5 and S6). Ki67 immunostaining showed the massive proliferation of cells in all strata of the epithelium typical of KIN (Fig. 4D). Inactivation of p53 signaling resulted in survival of more cells with a higher damage load as demonstrated by immunostaining for repair foci (Fig. 4E; Supplementary Fig. S4B) as compared with p53-competent Rnaseh2EKO skin (Fig. 3A). Flow cytometric analysis of Rnaseh2EKOTrp53−/− skin revealed that the loss of hair follicle stem cells caused by epidermal RNase H2 deficiency was prevented by lack of p53. Interestingly, hair follicle stem cell numbers of Rnaseh2EKOTrp53−/− epidermis ranged about 3-fold higher compared with control numbers (Fig. 4F).

Compared with p53-competent Rnaseh2EKO mice, leukocyte infiltration was more pronounced in Rnaseh2EKOTrp53−/− skin (Supplementary Fig. S4C). A potential reason for enhanced skin inflammation of these animals might be potentiated cytokine responses as we found strongly increased ISG mRNA levels in total epidermis and in keratinocytes of Trp53−/− as compared with p53-proficient Rnaseh2EKO mice (Fig. 4G), demonstrating that p53-dependent DNA damage responses control the intensity of the STING-induced IFN response mounted by RNase H2-deficient keratinocytes. Although we have shown that type I IFN was not causing the inflammation of Rnaseh2EKO skin (Fig. 2F), STING-mediated, NFκB-dependent proinflammatory cytokine expression, triggered in parallel to the IFN response, could be responsible for the enhanced inflammation of Ifnar1-deficient Rnaseh2EKO skin.

In summary, cancer development in Rnaseh2EKO mice is accelerated by additional absence of p53. This finding demonstrates that p53-dependent elimination of epithelial cells with high damage load potently antagonizes oncogenic transformation of RNase H2-deficient epidermis.

DNA polymerases do not replicate genomic DNA without mistakes. They incorporate nucleotides carrying the wrong base or wrong sugar at substantial rates. To ensure genome integrity, two DNA repair pathways operate on the newly synthesized strand to correct these replication errors. MMR, which removes nucleotides that do not base-pair correctly with the template strand, and RER, which removes rNTPs, that are incorporated into the genomic DNA at a rate of 1 every 7,600 nucleotides (12) and threaten genome stability. Humans and mice with genetic defects causing loss of MMR activity are viable but exhibit strongly increased cancer risk (6, 7, 43). Genetic defects of RER leading to inability to remove rNTPs from genomic DNA result in embryonic lethality in the mouse (12, 23) and most likely also in humans, as partial loss-of-function mutations of the genes encoding RNase H2, the enzyme essential for initiation of RER, can result in severe disease and biallelic complete null alleles have not been found (44).

We bypassed embryonic lethality of global RER deficiency by conditional inactivation of the Rnaseh2b gene only in the epidermis of the skin and observed spontaneous DNA damage and epithelial hyperproliferation, resulting in spontaneous cancer development in 100% of the animals within the first year of life. Thus, we showed that RER is essential to prevent malignant transformation.

Although the epithelium proliferated rapidly, it featured a high rate of keratinocyte apoptosis induced by p53 activation. Cancer development driven by epidermal RER deficiency was accelerated upon additional loss of one allele of the Trp53 gene. In Rnaseh2EKO mice harboring two Trp53 null alleles, the skin of neck and back transformed into one large confluent carcinoma in situ between 10 and 20 weeks of age. It seems likely that this process would have affected the entire skin and would have rapidly progressed to invasive SCC with increasing age, if not for the fact that the animals had to be euthanized for severe itch and erosions.

Carcinogenesis triggered by defective RNase H2 seems to require reduction of RNase H2 activity to low levels. Mice in which RNase H2 activity was reduced to 30% of control levels because of a homozygous Rnaseh2b partial loss-of-function mutation showed spontaneous activation of type I IFN responses but no sign of cancer by the age of 1 year (32). Likewise, mice that carried a biallelic knock in of a Rnaseh2a partial loss-of-function mutation, reducing RNase H2 activity to few percent of normal levels, and that were rescued from perinatal lethality by additional inactivation of STING were not reported to develop neoplastic disease (31). However, investigation of cancer incidence in a larger cohort of older animals was hampered by the fact that only a small fraction of these mice survived into adulthood (31).

We observed accelerated proliferation of Rnaseh2EKO epidermis resulting in epidermal hyperplasia and thickening of the cornified layer (hyperkeratosis) despite greatly increased rates of keratinocyte apoptosis. An important inducer of apoptosis in Rnaseh2EKO skin was p53, because additional loss of p53 resulted in significant further acceleration of epithelial hyperproliferation. These rapid cell divisions in response to keratinocyte genome damage ensuing from RER deficiency likely reflect a response program that the epidermis executes upon exposure to multiple forms of DNA damage, aiming at improved light protection by increasing thickness of epidermis and cornified layer, associated with enhanced melanin production, as we also observed in our animals. This response is activated, for example, upon overexposure to UV irradiation (45). Enhanced epidermal proliferation upon genome damage also results in elimination of damaged cells through rapid epithelial turnover, which represents an alternative keratinocyte disposal pathway in addition to apoptosis (46). Damaged basal cells arrested in late S-phase, detach from the basement membrane and get eliminated into the cornified layer (46).

Accelerated epithelial turnover in Rnaseh2EKO skin was associated with disturbed hair cycle regulation with progressive loss of hair follicle stem cells and hair follicle function. This might be a result of exhaustion of this stem cell population due to increased demand for differentiated keratinocytes, but may also reflect loss of hair follicle stem cells through apoptosis or terminal differentiation and transepidermal elimination. The latter was described to be the major cause for physiologic hair follicle aging and for premature hair follicle deterioration enforced by ionizing irradiation or genetic defects of nucleotide excision repair (47). Other epidermal stem cell populations were not reduced in numbers, likely reflecting differential responses to genome damage caused by defective RER. Specific responses of different stem cell populations to genome damage are documented in various tissues (48). Hair follicle stem cells were found to be more resistant to cell death induction by ionizing irradiation compared with other epidermal stem cell populations because of their capability to rapidly repair double-strand breaks by nonhomologous end joining and high level expression of the antiapoptotic protein Bcl2 (49). Epidermis-specific inactivation of homologous recombination repair by conditional knock out of BRCA1, however, resulted in selective depletion of hair follicle stem cells (50), similar to the loss of hair follicle stem cells in our Rnaseh2EKO animals. In both mouse lines, loss of hair and hair follicle stem cells are largely rescued by additional inactivation of the Trp53 gene (50), indicating p53-dependent cell death induction by RER deficiency in these cells. Additional absence of p53 in Rnaseh2EKO mice not only restored stem cell numbers to control levels, but resulted in three-fold higher hair follicle stem cell numbers compared with RER competent controls, likely reflecting preneoplastic expansion of this population.

Rapid progression to cancer in 100% of the animals with epidermis-specific loss of RER activity, even despite intact p53 signaling, raises the question of whether mutations impairing RER are relevant for human cancer. However, the studies of cancer incidence in patients with AGS carrying defects of RNASEH2 genes are hampered by their reduced life span (51).

Germline mutations compromising RER as well as somatic defects of RER acquired by growing tumor clones could affect tumor biology in two opposite ways. On one hand, DNA damage ensuing from unrepaired rNTPs causes genomic instability, thus promoting carcinogenesis and tumor progression. On the other hand, high rates of tumor cell death or senescence induction by persistent DNA lesions will antagonize malignant growth and can confer a competitive disadvantage to tumor cells. The net effect, tumor promotion or suppression, likely depends on the degree to which RER activity is reduced. A window of low RER activity advantageous for tumor growth may exist. The net effect of RER deficiency on cancer biology likely also depends on tumor type, as different cell types and cells of the same type but different differentiation stage show specific responses to DNA damage ranging from apoptosis and senescence to enforced differentiation (48, 52). In addition to detection of damage by nuclear sensors likely occurring in RNase H2-deficient cells, micronucleus formation associated with RER deficiency and subsequent micronuclear envelope rupture were shown to lead to frequent exposure of chromosomal DNA to the cytosolic DNA sensor cGAS, resulting in the activation of STING (24). STING is also activated upon DNA damage by herniation of chromatin into the cytosol (26). The potent type I IFN response we observed in Rnaseh2EKO skin most likely reflects activation of the cGAS-STING axis. Potentiation of this IFN response in the absence of functional p53 strongly suggests that genome damage is responsible for STING activation in our animals, as lack of p53 allows cells with high damage load producing IFN to survive longer. STING can trigger differential responses, including senescence or apoptosis as well as antiviral and inflammatory cytokine production stimulating immune activation, depending on signal strength and cell type. RER-deficient epidermis invariably develops cancer despite robust STING activation and production of type IFN, which can strongly boost antitumor immunity (53). Intriguingly, STING-mediated innate immune responses can also promote neoplastic growth, as STING activation by micronuclear DNA in chromosomally unstable cancer cells was recently shown to drive metastasis by activation of noncanonical NF-kB signaling (54).

Assessing contributions of RER defects to human cancer requires studies correlating frequencies of partial loss-of-function mutations in RNASEH2 genes, either in the germline or the tumor tissue of the patients, in particular cancer entities. Increased frequencies of RNASEH2B partial loss-of-function mutations were found in the germline of patients with sporadic prostate cancer compared with control cohorts (55). RNASEH2B partial loss-of-function mutations were also detected at higher frequencies in nontumor DNA of patients with glioma with familial cancer predisposition and cosegregated with manifestation of cancer in some of the families (55). These findings provide first hints that compromised RER activity can be associated with human cancer.

No potential conflicts of interest were disclosed.

Conception and design: W. Müller, T.A. Kunkel, R. Behrendt, A. Roers

Development of methodology: B. Hiller, A. Hoppe, C. Haase, C. Hiller

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): B. Hiller, A. Hoppe, C. Haase, C. Hiller, M.A.M. Reijns, A.P. Jackson, J. Wenzel, R. Behrendt

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): B. Hiller, A. Hoppe, C. Haase, M.A.M. Reijns, T.A. Kunkel, J. Wenzel, R. Behrendt

Writing, review, and/or revision of the manuscript: N. Schubert, W. Müller, M.A.M. Reijns, T.A. Kunkel, J. Wenzel, R. Behrendt, A. Roers

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases):

Study supervision: R. Behrendt, A. Roers

Others: Assisted with enzymatic assays (A.P. Jackson)

Discussion of the experiments at planning stage and discussions of the results (W. Müller)

We thank Tobias Häring for expert technical assistance.The study was supported by the Else Kröner-Fresenius Foundation (grant no. 2015_A100 to A. Roers and R. Behrendt), German Research Foundation (DFG; grant no. Ro2133/6-1 in the setting of Clinical Research Unit KFO249 to A. Roers), and CRTD seed grant (CRTD-FZ 111 to A. Roers). R. Behrendt is supported by the Aicardi-Goutières Syndrome Americas Association. Work in the laboratory of A.P. Jackson was supported by the Medical Research Council (MRC, U127580972).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Hoeijmakers
JH
. 
DNA damage, aging, and cancer
.
N Engl J Med
2009
;
361
:
1475
85
.
2.
Lindahl
T
,
Wood
RD
. 
Quality control by DNA repair
.
Science
1999
;
286
:
1897
905
.
3.
Ganai
RA
,
Johansson
E
. 
DNA Replication-A Matter of Fidelity
.
Mol Cell
2016
;
62
:
745
55
.
4.
Kunkel
TA
. 
Evolving views of DNA replication (in)fidelity
.
Cold Spring Harb Symp Quant Biol
2009
;
74
:
91
101
.
5.
Kunkel
TA
,
Erie
DA
. 
Eukaryotic mismatch repair in relation to DNA replication
.
Annu Rev Genet
2015
;
49
:
291
313
.
6.
Lynch
HT
,
Snyder
CL
,
Shaw
TG
,
Heinen
CD
,
Hitchins
MP
. 
Milestones of Lynch syndrome: 1895–2015
.
Nat Rev Cancer
2015
;
15
:
181
94
.
7.
Wimmer
K
,
Etzler
J
. 
Constitutional mismatch repair-deficiency syndrome: have we so far seen only the tip of an iceberg? Hum Genet
2008
;
124
:
105
22
.
8.
Nick McElhinny
SA
,
Watts
BE
,
Kumar
D
,
Watt
DL
,
Lundström
EB
,
Burgers
PM
, et al
Abundant ribonucleotide incorporation into DNA by yeast replicative polymerases
.
Proc Natl Acad Sci U S A
2010
;
107
:
4949
54
.
9.
Joyce
CM
. 
Choosing the right sugar: how polymerases select a nucleotide substrate
.
Proc Natl Acad Sci U S A
1997
;
94
:
1619
22
.
10.
Nick McElhinny
SA
,
Kumar
D
,
Clark
AB
,
Watt
DL
,
Watts
BE
,
Lundström
EB
, et al
Genome instability due to ribonucleotide incorporation into DNA
.
Nat Chem Biol
2010
;
6
:
774
81
.
11.
Traut
TW
. 
Physiological concentrations of purines and pyrimidines
.
Mol Cell Biochem
1994
;
140
:
1
22
.
12.
Reijns
MA
,
Rabe
B
,
Rigby
RE
,
Mill
P
,
Astell
KR
,
Lettice
LA
, et al
Enzymatic removal of ribonucleotides from DNA is essential for mammalian genome integrity and development
.
Cell
2012
;
149
:
1008
22
.
13.
Williams
JS
,
Lujan
SA
,
Kunkel
TA
. 
Processing ribonucleotides incorporated during eukaryotic DNA replication
.
Nat Rev Mol Cell Biol
2016
;
17
:
350
63
.
14.
Sparks
JL
,
Chon
H
,
Cerritelli
SM
,
Kunkel
TA
,
Johansson
E
,
Crouch
RJ
, et al
RNase H2-initiated ribonucleotide excision repair
.
Mol Cell
2012
;
47
:
980
6
.
15.
Cerritelli
SM
,
Crouch
RJ
. 
The balancing act of ribonucleotides in DNA
.
Trends Biochem Sci
2016
;
41
:
434
45
.
16.
Lujan
SA
,
Williams
JS
,
Clausen
AR
,
Clark
AB
,
Kunkel
TA
. 
Ribonucleotides are signals for mismatch repair of leading-strand replication errors
.
Mol Cell
2013
;
50
:
437
43
.
17.
Ghodgaonkar
MM
,
Lazzaro
F
,
Olivera-Pimentel
M
,
Artola-Borán
M
,
Cejka
P
,
Reijns
MA
, et al
Ribonucleotides misincorporated into DNA act as strand-discrimination signals in eukaryotic mismatch repair
.
Mol Cell
2013
;
50
:
323
32
.
18.
Williams
JS
,
Clausen
AR
,
Lujan
SA
,
Marjavaara
L
,
Clark
AB
,
Burgers
PM
, et al
Evidence that processing of ribonucleotides in DNA by topoisomerase 1 is leading-strand specific
.
Nat Struct Mol Biol
2015
;
22
:
291
7
.
19.
Caldecott
KW
. 
Molecular biology. Ribose–an internal threat to DNA
.
Science
2014
;
343
:
260
1
.
20.
Kim
N
,
Huang
SN
,
Williams
JS
,
Li
YC
,
Clark
AB
,
Cho
JE
, et al
Mutagenic processing of ribonucleotides in DNA by yeast topoisomerase I
.
Science
2011
;
332
:
1561
4
.
21.
Huang
SN
,
Williams
JS
,
Arana
ME
,
Kunkel
TA
,
Pommier
Y
. 
Topoisomerase I-mediated cleavage at unrepaired ribonucleotides generates DNA double-strand breaks
.
EMBO J
2017
;
36
:
361
73
.
22.
Chon
H
,
Sparks
JL
,
Rychlik
M
,
Nowotny
M
,
Burgers
PM
,
Crouch
RJ
, et al
RNase H2 roles in genome integrity revealed by unlinking its activities
.
Nucleic Acids Res
2013
;
41
:
3130
43
.
23.
Hiller
B
,
Achleitner
M
,
Glage
S
,
Naumann
R
,
Behrendt
R
,
Roers
A
. 
Mammalian RNase H2 removes ribonucleotides from DNA to maintain genome integrity
.
J Exp Med
2012
;
209
:
1419
26
.
24.
Mackenzie
KJ
,
Carroll
P
,
Martin
CA
,
Murina
O
,
Fluteau
A
,
Simpson
DJ
, et al
cGAS surveillance of micronuclei links genome instability to innate immunity
.
Nature
2017
;
548
:
461
5
.
25.
Bartsch
K
,
Knittler
K
,
Borowski
C
,
Rudnik
S
,
Damme
M
,
Aden
K
, et al
Absence of RNase H2 triggers generation of immunogenic micronuclei removed by autophagy
.
Hum Mol Genet
2017
;
26
:
3960
72
.
26.
Gluck
S
,
Guey
B
,
Gulen
MF
,
Wolter
K
,
Kang
TW
,
Schmacke
NA
, et al
Innate immune sensing of cytosolic chromatin fragments through cGAS promotes senescence
.
Nat Cell Biol
2017
;
19
:
1061
70
.
27.
Crow
YJ
,
Leitch
A
,
Hayward
BE
,
Garner
A
,
Parmar
R
,
Griffith
E
, et al
Mutations in genes encoding ribonuclease H2 subunits cause Aicardi-Goutieres syndrome and mimic congenital viral brain infection
.
Nat Genet
2006
;
38
:
910
6
.
28.
Crow
YJ
,
Manel
N
. 
Aicardi-Goutieres syndrome and the type I interferonopathies
.
Nat Rev Immunol
2015
;
15
:
429
40
.
29.
Reijns
MA
,
Bubeck
D
,
Gibson
LC
,
Graham
SC
,
Baillie
GS
,
Jones
EY
, et al
The structure of the human RNase H2 complex defines key interaction interfaces relevant to enzyme function and human disease
.
J Biol Chem
2011
;
286
:
10530
9
.
30.
Chon
H
,
Vassilev
A
,
DePamphilis
ML
,
Zhao
Y
,
Zhang
J
,
Burgers
PM
, et al
Contributions of the two accessory subunits, RNASEH2B and RNASEH2C, to the activity and properties of the human RNase H2 complex
.
Nucleic Acids Res
2009
;
37
:
96
110
.
31.
Pokatayev
V
,
Hasin
N
,
Chon
H
,
Cerritelli
SM
,
Sakhuja
K
,
Ward
JM
, et al
RNase H2 catalytic core Aicardi-Goutieres syndrome-related mutant invokes cGAS-STING innate immune-sensing pathway in mice
.
J Exp Med
2016
;
213
:
329
36
.
32.
Mackenzie
KJ
,
Carroll
P
,
Lettice
L
,
Tarnauskaitė
Ž
,
Reddy
K
,
Dix
F
, et al
Ribonuclease H2 mutations induce a cGAS/STING-dependent innate immune response
.
EMBO J
2016
;
35
:
831
44
.
33.
Hafner
M
,
Wenk
J
,
Nenci
A
,
Pasparakis
M
,
Scharffetter-Kochanek
K
,
Smyth
N
, et al
Keratin 14 Cre transgenic mice authenticate keratin 14 as an oocyte-expressed protein
.
Genesis
2004
;
38
:
176
81
.
34.
Jacks
T
,
Remington
L
,
Williams
BO
,
Schmitt
EM
,
Halachmi
S
,
Bronson
RT
, et al
Tumor spectrum analysis in p53-mutant mice
.
Curr Biol
1994
;
4
:
1
7
.
35.
Jensen
KB
,
Driskell
RR
,
Watt
FM
. 
Assaying proliferation and differentiation capacity of stem cells using disaggregated adult mouse epidermis
.
Nat Protoc
2010
;
5
:
898
911
.
36.
Love
MI
,
Huber
W
,
Anders
S
. 
Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
.
Genome Biol
2014
;
15
:
550
.
37.
Chang
CH
,
Kuo
CJ
,
Ito
T
,
Su
YY
,
Jiang
ST
,
Chiu
MH
, et al
CK1alpha ablation in keratinocytes induces p53-dependent, sunburn-protective skin hyperpigmentation
.
Proc Natl Acad Sci U S A
2017
;
114
:
E8035
44
.
38.
Cui
R
,
Widlund
HR
,
Feige
E
,
Lin
JY
,
Wilensky
DL
,
Igras
VE
, et al
Central role of p53 in the suntan response and pathologic hyperpigmentation
.
Cell
2007
;
128
:
853
64
.
39.
Wenzel
J
,
Tuting
T
. 
An IFN-associated cytotoxic cellular immune response against viral, self-, or tumor antigens is a common pathogenetic feature in "interface dermatitis"
.
J Invest Dermatol
2008
;
128
:
2392
402
.
40.
Essers
MA
,
Offner
S
,
Blanco-Bose
WE
,
Waibler
Z
,
Kalinke
U
,
Duchosal
MA
, et al
IFNalpha activates dormant haematopoietic stem cells in vivo
.
Nature
2009
;
458
:
904
8
.
41.
Kamphuis
E
,
Junt
T
,
Waibler
Z
,
Forster
R
,
Kalinke
U
. 
Type I interferons directly regulate lymphocyte recirculation and cause transient blood lymphopenia
.
Blood
2006
;
108
:
3253
61
.
42.
Fernandez Figueras
MT
. 
From actinic keratosis to squamous cell carcinoma: pathophysiology revisited
.
J Eur Acad Dermatol Venereol
2017
;
31
Suppl 2
:
5
7
.
43.
Lee
K
,
Tosti
E
,
Edelmann
W
. 
Mouse models of DNA mismatch repair in cancer research
.
DNA Repair (Amst)
2016
;
38
:
140
6
.
44.
Rice
G
,
Patrick
T
,
Parmar
R
,
Taylor
CF
,
Aeby
A
,
Aicardi
J
, et al
Clinical and molecular phenotype of Aicardi-Goutieres syndrome
.
Am J Hum Genet
2007
;
81
:
713
25
.
45.
Sterenborg
HJ
,
de Gruijl
FR
,
van der Leun
JC
. 
UV-induced epidermal hyperplasia in hairless mice
.
Photodermatol
1986
;
3
:
206
14
.
46.
Stout
GJ
,
Westdijk
D
,
Calkhoven
DM
,
Pijper
O
,
Backendorf
CM
,
Willemze
R
, et al
Epidermal transit of replication-arrested, undifferentiated keratinocytes in UV-exposed XPC mice: an alternative to in situ apoptosis
.
Proc Natl Acad Sci U S A
2005
;
102
:
18980
5
.
47.
Matsumura
H
,
Mohri
Y
,
Binh
NT
,
Morinaga
H
,
Fukuda
M
,
Ito
M
, et al
Hair follicle aging is driven by transepidermal elimination of stem cells via COL17A1 proteolysis
.
Science
2016
;
351
:
aad4395
.
48.
Blanpain
C
,
Mohrin
M
,
Sotiropoulou
PA
,
Passegué
E
. 
DNA-damage response in tissue-specific and cancer stem cells
.
Cell Stem Cell
2011
;
8
:
16
29
.
49.
Sotiropoulou
PA
,
Candi
A
,
Mascré
G
,
De Clercq
S
,
Youssef
KK
,
Lapouge
G
, et al
Bcl-2 and accelerated DNA repair mediates resistance of hair follicle bulge stem cells to DNA-damage-induced cell death
.
Nat Cell Biol
2010
;
12
:
572
82
.
50.
Sotiropoulou
PA
,
Karambelas
AE
,
Debaugnies
M
,
Candi
A
,
Bouwman
P
,
Moers
V
, et al
BRCA1 deficiency in skin epidermis leads to selective loss of hair follicle stem cells and their progeny
.
Genes Dev
2013
;
27
:
39
51
.
51.
Crow
YJ
,
Chase
DS
,
Lowenstein Schmidt
J
,
Szynkiewicz
M
,
Forte
GM
,
Gornall
HL
, et al
Characterization of human disease phenotypes associated with mutations in TREX1, RNASEH2A, RNASEH2B, RNASEH2C, SAMHD1, ADAR, and IFIH1
.
Am J Med Genet A
2015
;
167A
:
296
312
.
52.
Lans
H
,
Vermeulen
W
. 
Tissue specific response to DNA damage: C. elegans as role model
.
DNA Repair (Amst)
2015
;
32
:
141
8
.
53.
Zitvogel
L
,
Galluzzi
L
,
Kepp
O
,
Smyth
MJ
,
Kroemer
G
. 
Type I interferons in anticancer immunity
.
Nat Rev Immunol
2015
;
15
:
405
14
.
54.
Bakhoum
SF
,
Ngo
B
,
Laughney
AM
,
Cavallo
JA
,
Murphy
CJ
,
Ly
P
, et al
Chromosomal instability drives metastasis through a cytosolic DNA response
.
Nature
2018
;
553
:
467
72
.
55.
Beyer
U
,
Brand
F
,
Martens
H
,
Weder
J
,
Christians
A
,
Elyan
N
, et al
Rare ADAR and RNASEH2B variants and a type I interferon signature in glioma and prostate carcinoma risk and tumorigenesis
.
Acta Neuropathol
2017
;
134
:
905
22
.