RNaseH2A is involved in human gliomagenesis through the regulation of cell proliferation and apoptosis

Loss of the DNA Repair Gene RNase H2 Identifies a Unique Subset of DDR-Deficient Leiomyosarcomas

Targeting the DNA damage response (DDR) pathway is an emerging therapeutic approach for leiomyosarcoma (LMS), and loss of RNase H2, a DDR pathway member, is a potentially actionable alteration for DDR-targeted treatments. Therefore, we designed a protein- and genomic-based RNase H2 screening assay to determine its prevalence and prognostic significance. Using a selective RNase H2 antibody on a pan-tumor microarray (TMA), RNase H2 loss was more common in LMS (11.5%, 9/78) than across all tumors (3.8%, 32/843). In a separate LMS cohort, RNase H2 deficiency was confirmed in uterine LMS (U-LMS, 21%, 23/108) and soft-tissue LMS (ST-LMS; 30%, 39/102). In the TCGA database, RNASEH2B homozygous deletions (HomDels) were found in 6% (5/80) of LMS cases, with a higher proportion in U-LMS (15%; 4/27) compared with ST-LMS (2%; 1/53). Using the SNiPDx targeted-NGS sequencing assay to detect biallelic loss of function in select DDR-related genes, we found RNASEH2B HomDels in 54% (19/35) of U-LMS cases with RNase H2 loss by IHC, and 7% (3/43) HomDels in RNase H2 intact cases. No RNASEH2B HomDels were detected in ST-LMS. In U-LMS patient cohort (n = 109), no significant overall survival difference was seen in patients with RNase H2 loss versus intact, or RNASEH2B HomDel (n = 12) versus Non-HomDel (n = 37). The overall diagnostic accuracy, sensitivity, and specificity of RNase H2 IHC for detecting RNA-SEH2B HomDels in U-LMS was 76%, 93%, and 71%, respectively, and it is being developed for future predictive biomarker driven clinical trials targeting DDR in U-LMS.

Distinct features of ribonucleotides within genomic DNA in Aicardi-Goutières syndrome ortholog mutants of Saccharomyces cerevisiae

Ribonucleoside monophosphates (rNMPs) are abundantly found within genomic DNA of cells. The embedded rNMPs alter DNA properties and impact genome stability. Mutations in ribonuclease (RNase) H2, a key enzyme for rNMP removal, are associated with the Aicardi-Goutières syndrome (AGS), a severe neurological disorder. Here, we engineered orthologs of the human RNASEH2A-G37S and RNASEH2C-R69W AGS mutations in yeast Saccharomyces cerevisiae: rnh201-G42S and rnh203-K46W. Using the ribose-seq technique and the Ribose-Map bioinformatics toolkit, we unveiled rNMP abundance, composition, hotspots, and sequence context in these AGS-ortholog mutants. We found a high rNMP presence in the nuclear genome of rnh201-G42S-mutant cells, and an elevated rCMP content in both mutants, reflecting preferential cleavage of RNase H2 at rGMP. We discovered unique rNMP patterns in each mutant, showing differential activity of the AGS mutants on the leading or lagging replication strands. This study guides future research on rNMP characteristics in human genomes with AGS mutations.

Detection of ribonucleotides embedded in DNA by Nanopore sequencing

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes.

Distinct features of ribonucleotides within genomic DNA in Aicardi-Goutières syndrome (AGS)-ortholog mutants of Saccharomyces cerevisiae

Ribonucleoside monophosphates (rNMPs) are abundantly found within genomic DNA of cells. The embedded rNMPs alter DNA properties and impact genome stability. Mutations in ribonuclease (RNase) H2, a key enzyme for rNMP removal, are associated with the Aicardi-Goutières syndrome (AGS), a severe neurological disorder. Here, we engineered two AGS-ortholog mutations in Saccharomyces cerevisiae: rnh201-G42S and rnh203-K46W. Using the ribose-seq technique and the Ribose-Map bioinformatics toolkit, we unveiled rNMP abundance, composition, hotspots, and sequence context in these yeast AGS-ortholog mutants. We found higher rNMP incorporation in the nuclear genome of rnh201-G42S than in wild-type and rnh203-K46W-mutant cells, and an elevated rCMP content in both mutants. Moreover, we uncovered unique rNMP patterns in each mutant, highlighting a differential activity of the AGS mutants towards rNMPs embedded on the leading or on the lagging strand of DNA replication. This study guides future research on rNMP characteristics in human genomic samples carrying AGS mutations.

Transcription Factor ELK3 Promotes Stemness and Oxaliplatin Resistance of Glioma Cells by Regulating RNASEH2A

Oxaliplatin is a member of the platinum group that is often used to treat glioma, a common type of malignant brain tumor, though it does not come with desirable and notable effects. This study attempted to investigate how ELK3 impacts the oxaliplatin resistance of glioma cells and its molecular mechanism. Bioinformatics analysis was employed to screen mRNAs with differential expression in glioma cells and predict the possible regulator downstream. We used qRT-PCR to detect the expression of ELK3 and RNASEH2A. Dual-luciferase and ChIP assays were adopted to reassure the regulatory relationship between the two. We also evaluated cell viability and sphere formation efficiency through CCK-8 and sphere formation assay and calculated the IC50 value by using CCK-8 assay. The expression of stemness-related proteins (ALDH1 and Nanog) was assessed through western blot. Glioma cells and tissues presented a significantly high expression of ELK3, the knock-down of which would reduce the cell viability, stemness and oxaliplatin resistance dramatically. Bioinformatics analysis predicted RNASEH2A to be the downstream regulator of ELK3. RNASEH2A was remarkably upregulated in glioma tissue and cells. The results from dual luciferase assay and ChIP experiment verified the binding relationship between RNASEH2A promoter region and ELK3. Then through rescue experiments, we confirmed that overexpression of RNASEH2A could compensate for the inhibition of glioma cell progression resulting from the knock-down of ELK3. ELK3 could promote stemness and oxaliplatin resistance of glioma cells by upregulating RNASEH2A, indicating that targeting ELK3/RNASEH2A axis may be a possible solution to overcome oxaliplatin resistance of glioma cells.

Identification and verification of a 4-gene signature predicting the overall survival of cervical cancer

Cervical cancer (CC) is one of the most common gynecological malignancies, ranking fourth in both incidence and mortality in women worldwide. Early screening and treatment are of great significance in reducing the incidence and mortality of CC. Due to the complex molecular mechanisms of tumor progression, the predictive power of traditional clinical information is limited. In this study, an effective molecular model is established to assess prognosis of patients with CC and guide clinical treatment so as to improve their survival rate. Three high quality datasets (GSE138080, GSE52904, GSE67522) of expression profiling were obtained from gene expression omnibus (GEO) database. Another mRNA expression and clinicopathological data of CC were obtained from The Cancer Genome Atlas (TCGA) dataset. The bioinformatic analyses such as univariate analysis, multivariate Cox proportional-hazards model (Cox) analysis and lasso regression analysis were conducted to select survival-related differentially expressed genes (DEGs) and further establish a prognostic gene signature. Moreover, the performance of prognostic gene signature was evaluated based on Kaplan-Meier curve and receiver operating characteristic (ROC) curve. Gene set enrichment analysis (GSEA) and tumor immunity analysis were carried out to elucidate the molecular mechanisms and immune relevance. A 4-gene signature comprising procollagen-lysine, 2-oxoglutarate 5-dioxygenase 2 (PLOD2), spondin1 (SPON1), secreted phosphoprotein 1 (SPP1), ribonuclease H2 subunit A (RNASEH2A) was established to predict overall survival (OS) of CC. The ROC curve indicated good performance of the 4-gene signature in predicting OS of CC based on the TCGA dataset. The 4-gene signature classified the patients into high-risk and low-risk groups with distinct OS rates of CC. Univariate analysis and multivariate Cox regression analysis revealed that the 4-gene signature was an independent factor affecting the prognosis of patients with CC. Our study developed a 4-gene signature capable of predicting the OS of CC. The findings may be beneficial to individualized clinical treatment and timely follow-up for patients with CC.

Exploring Potential Regulatory Anesthetic Drugs Based on RNA Binding Protein and Constructing CESC Prognosis Model: A Study Based on TCGA Database

Objective

To investigate the differential expression of RBPs in cervical squamous cell carcinoma (CESC), analyze the regulatory effect of narcotic drugs on RBPs, and establish the prognostic risk model of CESC patients.

Methods

RNA-SEQ data and clinical case data of cancer and normal samples from CESC patients were obtained from the Cancer Genome Atlas (TCGA) database and Genotype-Tissue Expression (GTEx) database. Differentially expressed RBPs were screened by R language and enriched. The CMAP database is used to predict the anesthetic drugs that regulate the differential expression of RBPs. The prognostic risk score model was constructed by COX regression analysis. Risk score of each CESC patient was calculated and divided into high-risk group and low-risk group according to the median risk score. The prediction efficiency of prognostic risk model was evaluated by Kaplan-Meier (KM) analysis and receiver operating characteristic (ROC) curve, and the correlation between prognostic risk model and clinical characteristics was analyzed. Immunohistochemistry was used to detect the expression of RNASEH2A and HENMT1 in tissues.

Results

There were 65 differentially expressed RBPs in CESC. Five anesthetics, including benzocaine, procaine, pentoxyverine, and tetracaine were obtained to regulate RBPs. Survival analysis showed that seven genes were related to the prognosis of patients, and the CESC risk score model was constructed by COX regression. The risk score can be used as an independent prognostic factor. RNASEH2A and HENMT1 are up-regulated in tumors, which can effectively distinguish normal tissues from tumor tissues.

Conclusion

It is found that different anesthetic drugs have different regulatory effects on the differential expression of RBPs. Based on the differentially expressed RBPs, the prognostic risk score model of CESC patients was constructed. To provide ideas for the formulation of individualized precise anesthesia scheme and cancer pain analgesia scheme, which is helpful to improve the perioperative survival rate of cancer patients.

Long noncoding RNA LINC01287 promotes proliferation and inhibits apoptosis of lung adenocarcinoma cells via the miR-3529-5p/RNASEH2A axis under the competitive endogenous RNA pattern

Lung adenocarcinoma (LUAD) is regarded as the most common type of lung cancer. The molecular targeted therapies for LUAD have being extensively studied. Ribonuclease H2 subunit A (RNASEH2A) is a nucleotide degrading enzyme gene that exerts great influence on cell proliferation, DNA replication and genomic stability. According to bioinformatics analysis, RNASEH2A expression in LUAD tissues is predicted to be upregulated and high expression of RNASEH2A might be related to lower survival rate in LUAD patients. In the present study, we investigated functions of RNASEH2A in LUAD. The mRNA RNASEH2A showed high expression in LUAD cells, and its knockdown inhibited proliferation and induced apoptosis in LUAD cells. RNASEH2A was found to be a target gene of microRNA miR-3529-5p after their expression levels and interaction being examined. Long noncoding RNA LINC01287 upregulated RNASEH2A expression in LUAD cells by combining with miR-3529-5p in a competitive way. Rescue assays revealed that the overexpression of RNASEH2A reversed the suppression of cell proliferation and the promotion of cell apoptosis induced by miR-3529-5p overexpression or LINC01287 knockdown. Finally, forkhead box A1 (FOXA1) interacted with RNASEH2A promoter and LINC01287 promoter to upregulate the expression levels of RNASEH2A and LINC01287 in LUAD cells. Overall, FOXA1-induced LINC01287 serves as a competing endogenous RNA to promote proliferation and inhibit apoptosis of LUAD cells via upregulation of RNASEH2A expression at the posttranscriptional level by competitively combining with miR-3529-5p.

Hypoxia-induced RNASEH2A limits activation of cGAS-STING signaling in HCC and predicts poor prognosis

BACKGROUND

Hypoxia is a hallmark of solid cancers, including hepatocellular carcinoma (HCC). There is scarce information about how hypoxia avoids immunologic stress and maintains a cancer-promoting microenvironment.

METHODS

The Cancer Genome Atlas, RNA-seq data, and Oncomine database were used to discover the correlation of RNASEH2A with tumor progression; then expression of RNASEH2A mRNA and protein were detected in HCC tissues and cells subjected to hypoxia or with the treatment of CoCl₂ via real-time quantitative polymerase chain reaction and immunochemistry assays. Finally, the effect of RNASEH2A on cell proliferation and the involved signaling pathway was explored further.

RESULTS

RNASEH2A was positively correlated with tumor grade, size, vascular invasion, and poor prognosis. The expression of RNASEH2A mRNA and protein were increased and dependent on hypoxia-inducible factor 2α in HCC tissues and cell lines. Knockout of RNASEH2A in HCC cells greatly reduced cell proliferation and induced the transcription of multiple cGAS-STING (cyclic GMP-AMP synthase-stimulator of interferon genes) targeted type 1 interferon-related genes, including IFIT1, USP18, and CXCL10, which suggests knockout of RNASEH2A may produce immunologic stress and tumor suppressive effects.

CONCLUSIONS

RNASEH2A plays a critical role and potentially predicts patient outcomes in HCC, which uncovers a new mechanism that RNASEH2A contributes to limit immunologic stress of cancer cells in the context of hypoxia.

Aicardi-Goutières syndrome-associated gene SAMHD1 preserves genome integrity by preventing R-loop formation at transcription-replication conflict regions

The comorbid association of autoimmune diseases with cancers has been a major obstacle to successful anti-cancer treatment. Cancer survival rate decreases significantly in patients with preexisting autoimmunity. However, to date, the molecular and cellular profiles of such comorbidities are poorly understood. We used Aicardi-Goutières syndrome (AGS) as a model autoimmune disease and explored the underlying mechanisms of genome instability in AGS-associated-gene-deficient patient cells. We found that R-loops are highly enriched at transcription-replication conflict regions of the genome in fibroblast of patients bearing SAMHD1 mutation, which is the AGS-associated-gene mutation most frequently reported with tumor and malignancies. In SAMHD1-depleted cells, R-loops accumulated with the concomitant activation of DNA damage responses. Removal of R-loops in SAMHD1 deficiency reduced cellular responses to genome instability. Furthermore, downregulation of SAMHD1 expression is associated with various types of cancer and poor survival rate. Our findings suggest that SAMHD1 functions as a tumor suppressor by resolving R-loops, and thus, SAMHD1 and R-loop may be novel diagnostic markers and targets for patient stratification in anti-cancer therapy.

Mutations in SAMHD1 cause Aicardi-Goutières syndrome (AGS), a monogenic lupus-like autoimmune disease. Among AGS-associated genes, SAMHD1 is most frequently mutates in various types of tumors and malignancies, suggesting that it is biologically relevant to cancer development. Here, we show that SAMHD1 resolves R-loops induced by transcription-replication conflicts, thereby contributing to the maintenance of genome stability. Our findings provide insight into the molecular and mechanical understanding of the autoimmunity and cancer comorbidity, and suggest that SAMHD1 and R-loops are potential and reliable biomarkers in anti-cancer therapeutics.

Gene Co-Expression Analysis of Human RNASEH2A Reveals Functional Networks Associated with DNA Replication, DNA Damage Response, and Cell Cycle Regulation

Simple Summary

RNASEH2A is the catalytic subunit of the ribonuclease (RNase) H2 ternary complex that plays an important role in maintaining DNA stability in cells. Recent studies have shown that the RNASEH2A subunit alone is highly expressed in certain cancer cell types. Via a series of bioinformatics approaches, we found that RNASEH2A is highly expressed in human proliferative tissues and many cancers. Our analyses reveal a possible involvement of RNASEH2A in cell cycle regulation in addition to its well established role in DNA replication and DNA repair. Our findings underscore that RNASEH2A could serve as a biomarker for cancer diagnosis and a therapeutic target.

Abstract

Ribonuclease (RNase) H2 is a key enzyme for the removal of RNA found in DNA-RNA hybrids, playing a fundamental role in biological processes such as DNA replication, telomere maintenance, and DNA damage repair. RNase H2 is a trimer composed of three subunits, RNASEH2A being the catalytic subunit. RNASEH2A expression levels have been shown to be upregulated in transformed and cancer cells. In this study, we used a bioinformatics approach to identify RNASEH2A co-expressed genes in different human tissues to underscore biological processes associated with RNASEH2A expression. Our analysis shows functional networks for RNASEH2A involvement such as DNA replication and DNA damage response and a novel putative functional network of cell cycle regulation. Further bioinformatics investigation showed increased gene expression in different types of actively cycling cells and tissues, particularly in several cancers, supporting a biological role for RNASEH2A but not for the other two subunits of RNase H2 in cell proliferation. Mass spectrometry analysis of RNASEH2A-bound proteins identified players functioning in cell cycle regulation. Additional bioinformatic analysis showed that RNASEH2A correlates with cancer progression and cell cycle related genes in Cancer Cell Line Encyclopedia (CCLE) and The Cancer Genome Atlas (TCGA) Pan Cancer datasets and supported our mass spectrometry findings.

An integrative multi-omics network-based approach identifies key regulators for breast cancer

Although genome-wide association studies (GWASs) have successfully identified thousands of risk variants for human complex diseases, understanding the biological function and molecular mechanisms of the associated SNPs involved in complex diseases is challenging. Here we developed a framework named integrative multi-omics network-based approach (IMNA), aiming to identify potential key genes in regulatory networks by integrating molecular interactions across multiple biological scales, including GWAS signals, gene expression-based signatures, chromatin interactions and protein interactions from the network topology. We applied this approach to breast cancer, and prioritized key genes involved in regulatory networks. We also developed an abnormal gene expression score (AGES) signature based on the gene expression deviation of the top 20 rank-ordered genes in breast cancer. The AGES values are associated with genetic variants, tumor properties and patient survival outcomes. Among the top 20 genes, RNASEH2A was identified as a new candidate gene for breast cancer. Thus, our integrative network-based approach provides a genetic-driven framework to unveil tissue-specific interactions from multiple biological scales and reveal potential key regulatory genes for breast cancer. This approach can also be applied in other complex diseases such as ovarian cancer to unravel underlying mechanisms and help for developing therapeutic targets.

Simple inheritance, complex regulation: Supergene-mediated fire ant queen polymorphism

The fire ant Solenopsis invicta exists in two alternate social forms: monogyne nests contain a single reproductive queen and polygyne nests contain multiple reproductive queens. This colony-level social polymorphism corresponds with individual differences in queen physiology, queen dispersal patterns and worker discrimination behaviours, all evidently regulated by an inversion-based supergene that spans more than 13 Mb of a "social chromosome," contains over 400 protein-coding genes and rarely undergoes recombination. The specific mechanisms by which this supergene influences expression of the many distinctive features that characterize the alternate forms remain almost wholly unknown. To advance our understanding of these mechanisms, we explore the effects of social chromosome genotype and natal colony social form on gene expression in queens sampled as they embarked on nuptial flights, using RNA-sequencing of brains and ovaries. We observe a large effect of natal social form, that is, of the social/developmental environment, on gene expression profiles, with similarly substantial effects of genotype, including: (a) supergene-associated gene upregulation, (b) allele-specific expression and (c) pronounced extra-supergene trans-regulatory effects. These findings, along with observed spatial variation in differential and allele-specific expression within the supergene region, highlight the complex gene regulatory landscape that emerged following divergence of the inversion-mediated Sb haplotype from its homologue, which presumably largely retained the ancestral gene order. The distinctive supergene-associated gene expression trajectories we document at the onset of a queen's reproductive life expand the known record of relevant molecular correlates of a complex social polymorphism and point to putative genetic factors underpinning the alternate social syndromes.

A multidimensional systems biology analysis of cellular senescence in aging and disease

BACKGROUND

Cellular senescence, a permanent state of replicative arrest in otherwise proliferating cells, is a hallmark of aging and has been linked to aging-related diseases. Many genes play a role in cellular senescence, yet a comprehensive understanding of its pathways is still lacking.

RESULTS

We develop CellAge (http://genomics.senescence.info/cells), a manually curated database of 279 human genes driving cellular senescence, and perform various integrative analyses. Genes inducing cellular senescence tend to be overexpressed with age in human tissues and are significantly overrepresented in anti-longevity and tumor-suppressor genes, while genes inhibiting cellular senescence overlap with pro-longevity and oncogenes. Furthermore, cellular senescence genes are strongly conserved in mammals but not in invertebrates. We also build cellular senescence protein-protein interaction and co-expression networks. Clusters in the networks are enriched for cell cycle and immunological processes. Network topological parameters also reveal novel potential cellular senescence regulators. Using siRNAs, we observe that all 26 candidates tested induce at least one marker of senescence with 13 genes (C9orf40, CDC25A, CDCA4, CKAP2, GTF3C4, HAUS4, IMMT, MCM7, MTHFD2, MYBL2, NEK2, NIPA2, and TCEB3) decreasing cell number, activating p16/p21, and undergoing morphological changes that resemble cellular senescence.

CONCLUSIONS

Overall, our work provides a benchmark resource for researchers to study cellular senescence, and our systems biology analyses reveal new insights and gene regulators of cellular senescence.

One, No One, and One Hundred Thousand: The Many Forms of Ribonucleotides in DNA

In the last decade, it has become evident that RNA is frequently found in DNA. It is now well established that single embedded ribonucleoside monophosphates (rNMPs) are primarily introduced by DNA polymerases and that longer stretches of RNA can anneal to DNA, generating RNA:DNA hybrids. Among them, the most studied are R-loops, peculiar three-stranded nucleic acid structures formed upon the re-hybridization of a transcript to its template DNA. In addition, polyribonucleotide chains are synthesized to allow DNA replication priming, double-strand breaks repair, and may as well result from the direct incorporation of consecutive rNMPs by DNA polymerases. The bright side of RNA into DNA is that it contributes to regulating different physiological functions. The dark side, however, is that persistent RNA compromises genome integrity and genome stability. For these reasons, the characterization of all these structures has been under growing investigation. In this review, we discussed the origin of single and multiple ribonucleotides in the genome and in the DNA of organelles, focusing on situations where the aberrant processing of RNA:DNA hybrids may result in multiple rNMPs embedded in DNA. We concluded by providing an overview of the currently available strategies to study the presence of single and multiple ribonucleotides in DNA in vivo.

Single-Gene Deletions Contributing to Loss of Heterozygosity in Saccharomyces cerevisiae: Genome-Wide Screens and Reproducibility

Loss of heterozygosity (LOH) is a phenomenon commonly observed in cancers; the loss of chromosomal regions can be both causal and indicative of underlying genome instability. Yeast has long been used as a model organism to study genetic mechanisms difficult to study in mammalian cells. Studying gene deletions leading to increased LOH in yeast aids our understanding of the processes involved, and guides exploration into the etiology of LOH in cancers. Yet, before in-depth mechanistic studies can occur, candidate genes of interest must be identified. Utilizing the heterozygous Saccharomyces cerevisiae deletion collection (≈ 6500 strains), 217 genes whose disruption leads to increased LOH events at the endogenously heterozygous mating type locus were identified. Our investigation to refine this list of genes to candidates with the most definite impact on LOH includes: secondary testing for LOH impact at an additional locus, gene ontology analysis to determine common gene characteristics, and positional gene enrichment studies to identify chromosomal regions important in LOH events. Further, we conducted extensive comparisons of our data to screens with similar, but distinct methodologies, to further distinguish genes that are more likely to be true contributors to instability due to their reproducibility, and not just identified due to the stochastic nature of LOH. Finally, we selected nine candidate genes and quantitatively measured their impact on LOH as a benchmark for the impact of genes identified in our study. Our data add to the existing body of work and strengthen the evidence of single-gene knockdowns contributing to genome instability.