TPMCalculator: one-step software to quantify mRNA abundance of genomic features.

The quantification of RNA sequencing (RNA-seq) abundance using a normalization method that calculates transcripts per million (TPM) is a key step to compare multiple samples from different experiments. TPMCalculator is a one-step software to process RNA-seq alignments in BAM format and reports TPM values, raw read counts and feature lengths for genes, transcripts, exons and introns. The program describes the genomic features through a model generated from the gene transfer format file used during alignments reporting of the TPM values and the raw read counts for each feature. In this paper, we show the correlation for 1256 samples from the TCGA-BRCA project between TPM and FPKM reported by TPMCalculator and RSeQC. We also show the correlation for raw read counts reported by TPMCalculator, HTSeq and featureCounts.TPMCalculator is freely available at https://github.com/ncbi/TPMCalculator. It is implemented in C++14 and supported on Mac OS X, Linux and MS Windows.Supplementary data are available at Bioinformatics online

DNA and RNA Cleavage Complexes and Repair Pathway for TOP3B RNA- and DNA-Protein Crosslinks

The present study demonstrates that topoisomerase 3B (TOP3B) forms both RNA and DNA cleavage complexes (TOP3Bccs) in vivo and reveals a pathway for repairing TOP3Bccs. For inducing and detecting cellular TOP3Bccs, we engineer a self-trapping mutant of TOP3B (R338W-TOP3B). Transfection with R338W-TOP3B induces R-loops, genomic damage, and growth defect, which highlights the importance of TOP3Bcc repair mechanisms. To determine how cells repair TOP3Bccs, we deplete tyrosyl-DNA phosphodiesterases (TDP1 and TDP2). TDP2-deficient cells show elevated TOP3Bccs both in DNA and RNA. Conversely, overexpression of TDP2 lowers cellular TOP3Bccs. Using recombinant human TDP2, we demonstrate that TDP2 can process both denatured and proteolyzed TOP3Bccs. We also show that cellular TOP3Bccs are ubiquitinated by the E3 ligase TRIM41 before undergoing proteasomal processing and excision by TDP2

Genome-Wide Analysis Unveils DNA Helicase RECQ1 as a Regulator of Estrogen Response Pathway in Breast Cancer Cells

Susceptibility to breast cancer is significantly increased in individuals with germ line mutations in RECQ1 (also known as RECQL or RECQL1), a gene encoding a DNA helicase essential for genome maintenance. We previously reported that RECQ1 expression predicts clinical outcomes for sporadic breast cancer patients stratified by estrogen receptor (ER) status. Here, we utilized an unbiased integrative genomics approach to delineate a cross talk between RECQ1 and ERα, a known master regulatory transcription factor in breast cancer. We found that expression of ESR1, the gene encoding ERα, is directly activated by RECQ1. More than 35% of RECQ1 binding sites were cobound by ERα genome-wide. Mechanistically, RECQ1 cooperates with FOXA1, the pioneer transcription factor for ERα, to enhance chromatin accessibility at the ESR1 regulatory regions in a helicase activity-dependent manner. In clinical ERα-positive breast cancers treated with endocrine therapy, high RECQ1 and high FOXA1 coexpressing tumors were associated with better survival. Collectively, these results identify RECQ1 as a novel cofactor for ERα and uncover a previously unknown mechanism by which RECQ1 regulates disease-driving gene expression in ER-positive breast cancer cells

Dynamics of genomic and immune responses during primary immunotherapy resistance in mismatch repair–deficient tumors

Mismatch repair–deficient (dMMR) cancers generate a substantial number of immunogenic neoantigens, rendering them sensitive to immunotherapy. Yet, there is considerable variability in responses, and roughly one-half of dMMR cancers are refractory to immunotherapy. Here we study a patient with dMMR lung cancer refractory to immunotherapy. The tumor exhibited typical dMMR molecular features, including exceptionally high frameshift insertions and deletions (indels). Despite the treatment inducing abundant intratumoral T-cell infiltrates, it failed to elicit tumor regression, pointing to the T cells lacking cytotoxic activity. A post-treatment tumor demonstrated compound heterozygous frameshift deletions located upstream of the kinase domain in the gene encoding JAK1 protein, down-regulation of JAK1 and mediators of its signal transduction, and total loss of JAK1 phosphorylation. Importantly, one of the JAK1 mutations, despite not being detected in the pretreatment tumor, was found at low variant allele frequency in the pretreatment circulating tumor DNA, suggesting clonal selection of the mutation. To our knowledge, this report provides the most detailed look yet at defective JAK1 signaling in the context of dMMR and immunotherapy resistance. Together with observations of JAK1 frameshift indels being enriched in dMMR compared with MMR-proficient tumors, our findings demonstrate the critical function of JAK1 in immunological surveillance of dMMR cancer

Sarcoma_CellminerCDB: A tool to interrogate the genomic and functional characteristics of a comprehensive collection of sarcoma cell lines

Summary: Sarcomas are a diverse group of rare malignancies composed of multiple different clinical and molecular subtypes. Due to their rarity and heterogeneity, basic, translational, and clinical research in sarcoma has trailed behind that of other cancers. Outcomes for patients remain generally poor due to an incomplete understanding of disease biology and a lack of novel therapies. To address some of the limitations impeding preclinical sarcoma research, we have developed Sarcoma_CellMinerCDB, a publicly available interactive tool that merges publicly available sarcoma cell line data and newly generated omics data to create a comprehensive database of genomic, transcriptomic, methylomic, proteomic, metabolic, and pharmacologic data on 133 annotated sarcoma cell lines. The reproducibility, functionality, biological relevance, and therapeutic applications of Sarcoma_CellMinerCDB described herein are powerful tools to address and generate biological questions and test hypotheses for translational research. Sarcoma_CellMinerCDB (https://discover.nci.nih.gov/SarcomaCellMinerCDB) aims to contribute to advancing the preclinical study of sarcoma

SIRT1 Prevents R-Loops during Chronological Aging by Modulating DNA Replication at rDNA Loci

In metazoans, the largest sirtuin, SIRT1, is a nuclear protein implicated in epigenetic modifications, circadian signaling, DNA recombination, replication, and repair. Our previous studies have demonstrated that SIRT1 binds replication origins and inhibits replication initiation from a group of potential initiation sites (dormant origins). We studied the effects of aging and SIRT1 activity on replication origin usage and the incidence of transcription–replication collisions (creating R-loop structures) in adult human cells obtained at different time points during chronological aging and in cancer cells. In primary, untransformed cells, SIRT1 activity declined and the prevalence of R-loops rose with chronological aging. Both the reduction in SIRT1 activity and the increased abundance of R-loops were also observed during the passage of primary cells in culture. All cells, regardless of donor age or transformation status, reacted to the short-term, acute chemical inhibition of SIRT1 with the activation of excessive replication initiation events coincident with an increased prevalence of R-loops. However, cancer cells activated dormant replication origins, genome-wide, during long-term proliferation with mutated or depleted SIRT1, whereas, in primary cells, the aging-associated SIRT1-mediated activation of dormant origins was restricted to rDNA loci. These observations suggest that chronological aging and the associated decline in SIRT1 activity relax the regulatory networks that protect cells against excess replication and that the mechanisms protecting from replication–transcription collisions at the rDNA loci manifest as differentially enhanced sensitivities to SIRT1 decline and chronological aging

Dynamics of replication origin over-activation

DNA replication processes are often dysregulated in cancer. Here the authors analyse DNA synthesis patterns in cancer cells undergoing partial genome re-replication to reveal that re-replication exhibits aberrant replication fork dynamics and a skewed distribution of replication initiation that over-duplicates early-replicating genomic regions