Visualization of uracils created by APOBEC3A using UdgX shows colocalization with RPA at stalled replication forks

Differentiation signals induce APOBEC3A expression via GRHL3 in squamous epithelia and squamous cell carcinoma

Two APOBEC (apolipoprotein-B mRNA editing enzyme catalytic polypeptide-like) DNA cytosine deaminase enzymes (APOBEC3A and APOBEC3B) generate somatic mutations in cancer, driving tumour development and drug resistance. Here we used single cell RNA sequencing to study APOBEC3A and APOBEC3B expression in healthy and malignant mucosal epithelia, validating key observations with immunohistochemistry, spatial transcriptomics and functional experiments. Whereas APOBEC3B is expressed in keratinocytes entering mitosis, we show that APOBEC3A expression is confined largely to terminally differentiating cells and requires Grainyhead-like transcription factor 3 (GRHL3). Thus, in normal tissue, neither deaminase appears to be expressed at high levels during DNA replication, the cell cycle stage associated with APOBEC-mediated mutagenesis. In contrast, we show that in squamous cell carcinoma tissues, there is expansion of GRHL3 expression and activity to a subset of cells undergoing DNA replication and concomitant extension of APOBEC3A expression to proliferating cells. These findings indicate a mechanism for acquisition of APOBEC3A mutagenic activity in tumours.

Structural and functional coupling in cross-linking uracil-DNA glycosylase UDGX

Enzymes in uracil-DNA glycosylase (UDG) superfamily are involved in removal of deaminated nucleobases such as uracil, methylcytosine derivatives such as formylcytosine and carboxylcytosine, and other base damage in DNA repair. UDGX is the latest addition of a new class to the UDG superfamily with a sporadic distribution in bacteria. UDGX type enzymes have a distinct biochemical property of cross-linking itself to the resulting AP site after uracil removal. Built on previous biochemical and structural analyses, this work comprehensively investigated the kinetic and enzymatic properties of Mycobacterium smegmatis UDGX. Kinetics and mutational analyses, coupled with structural information, defined the roles of E52, D56, D59, F65 of motif 1, H178 of motif 2 and N91, K94, R107 and H109 of motif 3 play in uracil excision and cross-linking. More importantly, a series of quantitative analyses underscored the structural coupling through inter-motif and intra-motif interactions and subsequent functional coupling of the uracil excision and cross-linking reactions. A catalytic model is proposed, which underlies this catalytic feature unique to UDGX type enzymes. This study offers new insight on the catalytic mechanism of UDGX and provides a unique example of enzyme evolution.

Enzymatic Cleavage-Mediated Extension Stalling Enables Accurate Recognition and Quantification of Locus-Specific Uracil Modification in DNA

Chemical modifications in DNA have profound influences on the structures and functions of DNA. Uracil, a naturally occurring DNA modification, can originate from the deamination of cytosine or arise from misincorporation of dUTP into DNA during DNA replication. Uracil in DNA will imperil genomic stability due to their potential in producing detrimental mutations. An in-depth understanding of the functions of uracil modification requires the accurate determination of its site as well as content in genomes. Herein, we characterized that a new member of the uracil-DNA glycosylase (UDG) family enzyme (UdgX-H109S) could selectively cleave both uracil-containing single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA). Based on this unique property of UdgX-H109S, we developed an enzymatic cleavage-mediated extension stalling (ECES) method for the locus-specific detection and quantification of uracil in genomic DNA. In the ECES method, UdgX-H109S specifically recognizes and cleaves the N-glycosidic bond of uracil from dsDNA and generates an apurinic/apyrimidinic (AP) site, which could be broken by APE1 to form a one-nucleotide gap. The specific cleavage by UdgX-H109S is then evaluated and quantified by qPCR. With the developed ECES approach, we demonstrated that the level of uracil at position Chr4:50566961 in genomic DNA of breast cancer tissues was significantly decreased. Collectively, the ECES method has been proved to be accurate and reproducible in the locus-specific quantification of uracil in genomic DNA from biological and clinical samples.

Global DNA and protein interactomes of FLT1P1 (Fms-related tyrosine kinase 1 pseudogene 1) revealed its molecular regulatory functions associated with preeclampsia

BACKGROUND

Preeclampsia (PE) is one of the most serious pregnancy complications with unknown pathogenesis. Emerging evidence has demonstrated that Fms-related tyrosine kinase 1 (FLT1) is highly involved in PE development. As a pseudogene of FLT1, FLT1P1 increased in PE samples. However, its functions remain largely unknown.

METHODS AND RESULTS

In this study, co-expression analysis was performed to identify the potential target genes of FTL1P1. Then chromatin isolation using RNA purification (ChIRP) method was employed to explore the interactomes of FLT1P1, including interacting with DNA fragments and proteins. We found that in PE samples, both FLT1P1 and FLT1 were highly expressed and closely correlated. ChIRP-protein data revealed that FLT1P1 interacts with translation- and transcription-related proteins, including 4 transcription factors (TFs). ChIRP-DNA analysis revealed that FLT1P1 preferentially interacted with DNA fragments downstream of transcription start sites (TSSs). Functional analysis of its interacting genes revealed that they were enriched in transcriptional regulation and apoptosis-related pathways. Twenty-six TFs, including CREB1 and SRF, were extracted from the potential FLT1P1-interacting gene sets and were potential targets of FLT1P1. CREB1 could bind to FLT1 promoter, and was negatively correlated with FLT1 at the expression level, making it a potential regulator of FLT1.

CONCLUSIONS

Our study extensively investigated the interactome profiles of FLT1P1, especially the prompter region of TF gene CREB1, and revealed the potential molecular regulatory mechanisms of FLT1 expression in PE samples. Our results provide a novel view of PE pathogenesis, and suggest that FLT1P1 could serve as a potential therapeutic target in PE diagnosis and treatment.

Competition for DNA binding between the genome protector replication protein A and the genome modifying APOBEC3 single-stranded DNA deaminases

The human APOBEC family of eleven cytosine deaminases use RNA and single-stranded DNA (ssDNA) as substrates to deaminate cytosine to uracil. This deamination event has roles in lipid metabolism by altering mRNA coding, adaptive immunity by causing evolution of antibody genes, and innate immunity through inactivation of viral genomes. These benefits come at a cost where some family members, primarily from the APOBEC3 subfamily (APOBEC3A-H, excluding E), can cause off-target deaminations of cytosine to form uracil on transiently single-stranded genomic DNA, which induces mutations that are associated with cancer evolution. Since uracil is only promutagenic, the mutations observed in cancer genomes originate only when uracil is not removed by uracil DNA glycosylase (UNG) or when the UNG-induced abasic site is erroneously repaired. However, when ssDNA is present, replication protein A (RPA) binds and protects the DNA from nucleases or recruits DNA repair proteins, such as UNG. Thus, APOBEC enzymes must compete with RPA to access their substrate. Certain APOBEC enzymes can displace RPA, bind and scan ssDNA efficiently to search for cytosines, and can become highly overexpressed in tumor cells. Depending on the DNA replication conditions and DNA structure, RPA can either be in excess or deficient. Here we discuss the interplay between these factors and how despite RPA, multiple cancer genomes have a mutation bias at cytosines indicative of APOBEC activity.

A redox-sensitive iron-sulfur cluster in murine FAM72A controls its ability to degrade the nuclear form of uracil-DNA glycosylase

Murine FAM72A, mFAM72A, binds the nuclear form of uracil-DNA glycosylase, mUNG2, inhibits its activity and causes its degradation. In immunoprecipitation assays the human paralog, hFAM72A, binds hUNG2 and is a potential anti-cancer drug target because of its high expression in many cancers. Using purified mFAM72A, and mUNG2 proteins we show that mFAM72A binds mUNG2, and the N-terminal 25 amino acids of mUNG2 bind mFAM72A at a nanomolar dissociation constant. We also show that mFAM72A is present throughout the cells, and mUNG2 helps localize it to nuclei. Based on in silico models of mFAM72A-mUNG2 interactions, we constructed several mutants of mFAM72A and found that while they have reduced ability to deplete mUNG2, the mutations also destabilized the former protein. We confirmed that Withaferin A, a predicted lead molecule for the design of FAM72A inhibitors, binds mFAM72A with micromolar affinity but has little affinity to mUNG2. We identified two potential metal-binding sites in mFAM72A and show that one of the sites contains an Fe-S cluster. This redox-sensitive cluster is involved in the mFAM72A-mUNG2 interaction and modulates mFAM72A activity. Hydrogen peroxide treatment of cells increases mUNG2 depletion in a FAM72A-dependent fashion suggesting that mFAM72A activity is redox-sensitive.

Genome-Wide Analysis of Human Long Noncoding RNAs: A Provocative Review

Do long noncoding RNAs (lncRNAs) contribute little or substantively to human biology? To address how lncRNA loci and their transcripts, structures, interactions, and functions contribute to human traits and disease, we adopt a genome-wide perspective. We intend to provoke alternative interpretation of questionable evidence and thorough inquiry into unsubstantiated claims. We discuss pitfalls of lncRNA experimental and computational methods as well as opposing interpretations of their results. The majority of evidence, we argue, indicates that most lncRNA transcript models reflect transcriptional noise or provide minor regulatory roles, leaving relatively few human lncRNAs that contribute centrally to human development, physiology, or behavior. These important few tend to be spliced and better conserved but lack a simple syntax relating sequence to structure and mechanism, and so resist simple categorization. This genome-wide view should help investigators prioritize individual lncRNAs based on their likely contribution to human biology.

Comprehensive characterization of clonality of driver genes revealing their clinical relevance in colorectal cancer

Background

Genomic studies of colorectal cancer have revealed the complex genomic heterogeneity of the tumor. The acquisition and selection of genomic alterations may be critical to understanding the initiation and progression of this disease.

Methods

In this study, we have systematically characterized the clonal architecture of 97 driver genes in 536 colorectal cancer patients from TCGA.

Results

A high proportion of clonal mutations in 93 driver genes were observed. 40 genes showed significant associations between their clonality and multiple clinicopathologic factors. Kaplan–Meier analysis suggested that the mutation clonality of ANK1, CASP8, SMAD2, and ARID1A had a significant impact on the CRC patients' outcomes. Multivariable analysis revealed that subclonal ANK1 mutations, clonal CASP8 mutations, and clonal SMAD2 mutations independently predicted for shorter overall survival after adjusting for clinicopathological factors. The poor outcome of the subclonal ANK1 mutation may be caused by upregulation of IL4I1, IDO1, IFNG and MAPK12 which showed potential roles in tumor immune evasion through accumulation of immunosuppressive cells such as regulatory T cells and myeloid derived suppressor cells.

Conclusion

These results suggested that the clonality of driver genes could act as prognostic markers and potential therapeutic targets in human colorectal cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03529-x.

Cross Brain-Gut Analysis Highlighted Hub Genes and LncRNA Networks Differentially Modified During Leucine Consumption and Endurance Exercise in Mice with Depression-Like Behaviors

Depression is a frequent mood disorder that might impair the brain-gut axis. In this study, we divided 30 mice into five groups: untreated mice, mice with depression-like behaviors, mice with depression-like behaviors treated with consumed leucine, mice with depression-like behaviors treated with exercise training, mice with depression-like behaviors treated with exercise training along with consumed leucine. According to artificial intelligence biological analysis, we found some mediators such as lncRNAs profile and Kdr/Vegfα/Pten/Bdnf interactions network in the hippocampus region and ileum tissue which could be decisive molecules in the brain-gut axis. Moreover, KDR as a principal cutpoint protein in the network was identified as the pharmaceutical approach for major depressive ameliorating based on pharmacophore modeling and molecular docking outcomes. Furthermore, we indicated that the mRNA and protein level of the Pten enhanced and Vegfα/Kdr/Bdnf mRNAs, as well as the protein level of KDR, decreased in mice with depression-like behaviors. Moreover, exercise and leucine ameliorated the brain-gut axis in mice with depression-like behaviors. Exercise and leucine regulated the lncRNAs network in the hippocampus and ileum of mice with depression-like behaviors. We suggest that the lncRNAs profiles could be considered as diagnosis and prognosis biomarkers, and exercise + leucine might be a practical approach to improve depression.

UdgX-Mediated Uracil Sequencing at Single-Nucleotide Resolution

As an aberrant base in DNA, uracil is generated by either deoxyuridine (dU) misincorporation or cytosine deamination, and involved in multiple physiological and pathological processes. Genome-wide profiles of uracil are important for study of these processes. Current methods for whole-genome mapping of uracil all rely on uracil-DNA N-glycosylase (UNG) and are limited in resolution, specificity, and/or sensitivity. Here, we developed a UdgX cross-linking and polymerase stalling sequencing ("Ucaps-seq") method to detect dU at single-nucleotide resolution. First, the specificity of Ucaps-seq was confirmed on synthetic DNA. Then the effectiveness of the approach was verified on two genomes from different sources. Ucaps-seq not only identified the enrichment of dU at dT sites in pemetrexed-treated cancer cells with globally elevated uracil but also detected dU at dC sites within the "WRC" motif in activated B cells which have increased dU in specific regions. Finally, Ucaps-seq was utilized to detect dU introduced by the cytosine base editor (nCas9-APOBEC) and identified a novel off-target site in cellular context. In conclusion, Ucaps-seq is a powerful tool with many potential applications, especially in evaluation of base editing fidelity.

OUP accepted manuscript

Activation-induced deaminase (AID) is a DNA-cytosine deaminase that mediates maturation of antibodies through somatic hypermutation and class-switch recombination. While it causes mutations in immunoglobulin heavy and light chain genes and strand breaks in the switch regions of the immunoglobulin heavy chain gene, it largely avoids causing such damage in the rest of the genome. To help understand targeting by human AID, we expressed it in repair-deficient Escherichia coli and mapped the created uracils in the genomic DNA using uracil pull-down and sequencing, UPD-seq. We found that both AID and the human APOBEC3A preferentially target tRNA genes and transcription start sites, but do not show preference for highly transcribed genes. Unlike A3A, AID did not show a strong replicative strand bias or a preference for hairpin loops. Overlapping uracilation peaks between these enzymes contained binding sites for a protein, FIS, that helps create topological domains in the E. coli genome. To confirm whether these findings were relevant to B cells, we examined mutations from lymphoma and leukemia genomes within AID-preferred sequences. These mutations also lacked replicative strand bias or a hairpin loop preference. We propose here a model for how AID avoids causing mutations in the single-stranded DNA found within replication forks.

Single-stranded DNA binding proteins influence APOBEC3A substrate preference

The cytidine deaminase, APOBEC3A (A3A), is a prominent source of mutations in multiple cancer types. These APOBEC-signature mutations are non-uniformly distributed across cancer genomes, associating with single-stranded (ss) DNA formed during DNA replication and hairpin-forming sequences. The biochemical and cellular factors that influence these specificities are unclear. We measured A3A's cytidine deaminase activity in vitro on substrates that model potential sources of ssDNA in the cell and found that A3A is more active on hairpins containing 4 nt ssDNA loops compared to hairpins with larger loops, bubble structures, replication fork mimics, ssDNA gaps, or linear DNA. Despite pre-bent ssDNAs being expected to fit better in the A3A active site, we determined A3A favors a 4 nt hairpin substrate only 2- to fivefold over linear ssDNA substrates. Addition of whole cell lysates or purified RPA to cytidine deaminase assays more severely reduced A3A activity on linear ssDNA (45 nt) compared to hairpin substrates. These results indicate that the large enrichment of A3A-driven mutations in hairpin-forming sequences in tumor genomes is likely driven in part by other proteins that preferentially bind longer ssDNA regions, which limit A3A's access. Furthermore, A3A activity is reduced at ssDNA associated with a stalled T7 RNA polymerase, suggesting that potential protein occlusion by RNA polymerase also limits A3A activity. These results help explain the small transcriptional strand bias for APOBEC mutation signatures in cancer genomes and the general targeting of hairpin-forming sequences in the lagging strand template during DNA replication.

Non-coding RNAs as Novel Biomarkers in Cancer Drug Resistance

Chemotherapy is often the primary and most effective anticancer treatment; however, drug resistance remains a major obstacle to it being curative. Recent studies have demonstrated that non-coding RNAs (ncRNAs), especially microRNAs and long non-coding RNAs, are involved in drug resistance of tumor cells in many ways, such as modulation of apoptosis, drug efflux and metabolism, epithelial-to-mesenchymal transition, DNA repair, and cell cycle progression. Exploring the relationships between ncRNAs and drug resistance will not only contribute to our understanding of the mechanisms of drug resistance and provide ncRNA biomarkers of chemoresistance, but will also help realize personalized anticancer treatment regimens. Due to the high cost and low efficiency of biological experimentation, many researchers have opted to use computational methods to identify ncRNA biomarkers associated with drug resistance. In this review, we summarize recent discoveries related to ncRNA-mediated drug resistance and highlight the computational methods and resources available for ncRNA biomarkers involved in chemoresistance.

Base-Resolution Analysis of Deoxyuridine at the Genome Scale Based on the Artificial Incorporation Modified Nucleobase

Deamination of cytosine and dUMP misincorporation have been found to be capable of producing uracil in the genome. This study presents the AI-seq (artificial incorporation modified nucleobase for sequencing), a "base substitution", which not only is capable of profiling uracil at single-nucleotide resolution and showing its centromeric enrichment but could also reveal that the identified uracil sites are derived from cytosine deamination. All the results indicate the potential biological significance of uracil as the epigenetic modification.

Annotation of snoRNA abundance across human tissues reveals complex snoRNA-host gene relationships

BACKGROUND

Small nucleolar RNAs (snoRNAs) are mid-size non-coding RNAs required for ribosomal RNA modification, implying a ubiquitous tissue distribution linked to ribosome synthesis. However, increasing numbers of studies identify extra-ribosomal roles of snoRNAs in modulating gene expression, suggesting more complex snoRNA abundance patterns. Therefore, there is a great need for mapping the snoRNome in different human tissues as the blueprint for snoRNA functions.

RESULTS

We used a low structure bias RNA-Seq approach to accurately quantify snoRNAs and compare them to the entire transcriptome in seven healthy human tissues (breast, ovary, prostate, testis, skeletal muscle, liver, and brain). We identify 475 expressed snoRNAs categorized in two abundance classes that differ significantly in their function, conservation level, and correlation with their host gene: 390 snoRNAs are uniformly expressed and 85 are enriched in the brain or reproductive tissues. Most tissue-enriched snoRNAs are embedded in lncRNAs and display strong correlation of abundance with them, whereas uniformly expressed snoRNAs are mostly embedded in protein-coding host genes and are mainly non- or anticorrelated with them. Fifty-nine percent of the non-correlated or anticorrelated protein-coding host gene/snoRNA pairs feature dual-initiation promoters, compared to only 16% of the correlated non-coding host gene/snoRNA pairs.

CONCLUSIONS

Our results demonstrate that snoRNAs are not a single homogeneous group of housekeeping genes but include highly regulated tissue-enriched RNAs. Indeed, our work indicates that the architecture of snoRNA host genes varies to uncouple the host and snoRNA expressions in order to meet the different snoRNA abundance levels and functional needs of human tissues.

Preferential CEBP binding to T:G mismatches and increased C-to-T human somatic mutations

DNA cytosine methylation in mammals modulates gene expression and chromatin accessibility. It also impacts mutation rates, via spontaneous oxidative deamination of 5-methylcytosine (5mC) to thymine. In most cases the resulting T:G mismatches are repaired, following T excision by one of the thymine DNA glycosylases, TDG or MBD4. We found that C-to-T mutations are enriched in the binding sites of CCAAT/enhancer binding proteins (CEBP). Within a CEBP site, the presence of a T:G mismatch increased CEBPβ binding affinity by a factor of >60 relative to the normal C:G base pair. This enhanced binding to a mismatch inhibits its repair by both TDG and MBD4 in vitro. Furthermore, repair of the deamination product of unmethylated cytosine, which yields a U:G DNA mismatch that is normally repaired via uracil DNA glycosylase, is also inhibited by CEBPβ binding. Passage of a replication fork over either a T:G or U:G mismatch, before repair can occur, results in a C-to-T mutation in one of the daughter duplexes. Our study thus provides a plausible mechanism for accumulation of C-to-T human somatic mutations.

Comprehensive analyses of correlation and survival reveal informative lncRNA prognostic signatures in colon cancer

BACKGROUND

Colon cancer is a commonly worldwide cancer with high morbidity and mortality. Long non-coding RNAs (lncRNAs) are involved in many biological processes and are closely related to the occurrence of colon cancer. Identification of the prognostic signatures of lncRNAs in colon cancer has great significance for its treatment.

METHODS

We first identified the colon cancer-related mRNAs and lncRNAs according to the differential analysis methods using the expression data in TCGA. Then, we performed correlation analysis between the identified mRNAs and lncRNAs by integrating their expression values and secondary structure information to estimate the co-regulatory relationships between the cancer-related mRNAs and lncRNAs. Besides, the competing endogenous RNA regulation network based on co-regulatory relationships was constructed to reveal cancer-related regulatory patterns. Meanwhile, we used traditional regression analysis (univariate Cox analysis, random survival forest analysis, and lasso regression analysis) to screen the cancer-related lncRNAs. Finally, by combining the identified colon cancer-related lncRNAs according to the above analyses, we constructed a risk prognosis model for colon cancer through multivariate Cox analysis and also validated the model in the colon cancer dataset in TCGA cohorts.

RESULTS

Six lncRNAs were found highly correlated with the overall survival of colon cancer patients, and a risk prognosis model based on them was constructed to predict the overall survival of colon cancer patients. In particular, EVX1-AS, ZNF667-AS1, CTC-428G20.6, and CTC-297N7.9 were first reported to be related to colon cancer by using our model, among which EVX1-AS and ZNF667-AS1 have been predicted to be related to colon cancer in LncRNADisease database.

CONCLUSIONS

This study identified the potential regulatory relationships between lncRNAs and mRNAs by integrating their expression values and secondary structure information and presented a significant 6-lncRNA risk prognosis model to predict the overall survival of colon cancer patients.

Detection of Genomic Uracil Patterns

The appearance of uracil in the deoxyuridine moiety of DNA is among the most frequently occurring genomic modifications. Three different routes can result in genomic uracil, two of which do not require specific enzymes: spontaneous cytosine deamination due to the inherent chemical reactivity of living cells, and thymine-replacing incorporation upon nucleotide pool imbalances. There is also an enzymatic pathway of cytosine deamination with multiple DNA (cytosine) deaminases involved in this process. In order to describe potential roles of genomic uracil, it is of key importance to utilize efficient uracil-DNA detection methods. In this review, we provide a comprehensive and critical assessment of currently available uracil detection methods with special focus on genome-wide mapping solutions. Recent developments in PCR-based and in situ detection as well as the quantitation of genomic uracil are also discussed.

Integrated RNA Sequencing and Single-Cell Mass Cytometry Reveal a Novel Role of LncRNA HOXA-AS2 in Tumorigenesis and Stemness of Hepatocellular Carcinoma

Purpose

Accumulating evidence indicates that long non-coding RNAs (lncRNAs) play critical roles in the development of many cancer types. However, the changes of lncRNAs expression profiles in hepatocarcinogenesis remain largely unknown. Therefore, the purpose of this study was to identify the clinical significance, oncogenic functions, and potential mechanism of cancer-related lncRNAs in hepatocellular carcinoma (HCC).

Materials and Methods

An in vitro hepatocellular carcinoma model was established via oncogene-mediated transformation with a combination of three genetic alterations, including hTERT overexpression, inactivation of P53, and KRAS activation. Changes of biological function and transcriptome profile in these cell lines were determined by colony formation assay, MTT assay, wound-healing scratch assay, xenograft nude mice model, mass cytometry and RNA sequencing (RNA-Seq). Furthermore, 116 HCC tissues and its corresponding normal tumor-adjacent tissues were explored to validate the results of cell lines. Finally, RNA sequencing, single-cell mass cytometry and fluorescence-activated cell sorter were applied to evaluate the potential association between the expression of lncRNA and the stemness of HCC.

Results

LncRNA HOXA-AS2 was aberrantly upregulated and it may be involved in the regulation of cancer stem cells during oncogenic transformation. Consistently, lncRNA HOXA-AS2 expression was significantly upregulated in HCC and its higher expression positively correlated with poor prognosis and stem cell-related functions. Moreover, a specific cancer stem cell subpopulation with EPCAM⁺, C-MYC⁺ and CK19⁺ may exist in higher HOXA-AS2 expression HCC patients.

Conclusion

LncRNA HOXA-AS2 plays pivotal roles in the occurrence and progression of HCC, which may act as a therapeutic target for prognostic biomarker in hepatocellular carcinoma.