Abstract
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.
Introduction
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.
Materials and Methods
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 CD45−CD49f+ 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.
Results
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).
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).
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.
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).
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.
Discussion
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.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
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)
Acknowledgments
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.