Overexpression of human DNA polymerase theta is a biomarker of aggressive and DNA repair-deficient papillary thyroid cancers.

BACKGROUND: DNA polymerase theta (POLQ) is an enzyme that repairs double-strand DNA breaks. POLQ is overexpressed in several cancer types, and increased expression is associated with a poor prognosis. Ablating POLQ function in vitro increases drug sensitivity to agents that cause double-strand DNA breaks, including chemotherapies and ionizing radiation. POLQ's role in thyroid cancer remains poorly understood. METHODS: Expression of POLQ and other genes of interest were analyzed in 513 papillary thyroid cancers (505 primary tumors and 8 metastatic lesions) and 59 normal thyroid samples available in the Cancer Genome Atlas. The Cancer Genome Atlas RNA and DNA sequencing data were queried with the Xena platform. The Recombination Proficiency Score was calculated to assess DNA repair efficiency. Other signaling events associated with thyroid tumorigenesis and clinical outcomes were analyzed. Univariate and multivariate analyses were performed. Treatment with the POLQ inhibitors ART558 and Novobiocin tested the effect of POLQ inhibition on in vitro thyroid cancer growth. RESULTS: POLQ expression was increased in papillary thyroid cancers compared to normal thyroid tissue (P < .05). POLQ expression levels were inversely correlated with Recombination Proficiency Score levels (P < .05). POLQ expression was highest in tall cell papillary thyroid cancers and in metastases. Higher POLQ expression was also associated with dedifferentiation, BRAF signaling, and shorter progression-free intervals (P < .05). Treatment with POLQ inhibitors decreased in vitro thyroid cancer growth (P < .05). CONCLUSION: These findings suggest that increased POLQ expression could serve as a valuable clinical marker and a potential therapeutic target in the treatment of thyroid cancer.

CrypticProteinDB: an integrated database of proteome and immunopeptidome derived non-canonical cancer proteins.

Translated non-canonical proteins derived from noncoding regions or alternative open reading frames (ORFs) can contribute to critical and diverse cellular processes. In the context of cancer, they also represent an under-appreciated source of targets for cancer immunotherapy through their tumor-enriched expression or by harboring somatic mutations that produce neoantigens. Here, we introduce the largest integration and proteogenomic analysis of novel peptides to assess the prevalence of non-canonical ORFs (ncORFs) in more than 900 patient proteomes and 26 immunopeptidome datasets across 14 cancer types. The integrative proteogenomic analysis of whole-cell proteomes and immunopeptidomes revealed peptide support for a nonredundant set of 9760 upstream, downstream, and out-of-frame ncORFs in protein coding genes and 12811 in noncoding RNAs. Notably, 6486 ncORFs were derived from differentially expressed genes and 340 were ubiquitously translated across eight or more cancers. The analysis also led to the discovery of thirty-four epitopes and eight neoantigens from non-canonical proteins in two cohorts as novel cancer immunotargets. Collectively, our analysis integrated both bottom-up proteogenomic and targeted peptide validation to illustrate the prevalence of translated non-canonical proteins in cancer and to provide a resource for the prioritization of novel proteins supported by proteomic, immunopeptidomic, genomic and transcriptomic data, available at https://www.maherlab.com/crypticproteindb.

Characterization of cell-type specific circular RNAs associated with colorectal cancer metastasis.

Colorectal cancer (CRC) is the most common gastrointestinal malignancy and a leading cause of cancer deaths in the United States. More than half of CRC patients develop metastatic disease (mCRC) with an average 5-year survival rate of 13%. Circular RNAs (circRNAs) have recently emerged as important tumorigenesis regulators; however, their role in mCRC progression remains poorly characterized. Further, little is known about their cell-type specificity to elucidate their functions in the tumor microenvironment (TME). To address this, we performed total RNA sequencing (RNA-seq) on 30 matched normal, primary and metastatic samples from 14 mCRC patients. Additionally, five CRC cell lines were sequenced to construct a circRNA catalog in CRC. We detected 47 869 circRNAs, with 51% previously unannotated in CRC and 14% novel candidates when compared to existing circRNA databases. We identified 362 circRNAs differentially expressed in primary and/or metastatic tissues, termed circular RNAs associated with metastasis (CRAMS). We performed cell-type deconvolution using published single-cell RNA-seq datasets and applied a non-negative least squares statistical model to estimate cell-type specific circRNA expression. This predicted 667 circRNAs as exclusively expressed in a single cell type. Collectively, this serves as a valuable resource, TMECircDB (accessible at https://www.maherlab.com/tmecircdb-overview), for functional characterization of circRNAs in mCRC, specifically in the TME.

Pan-Cancer Analysis Reveals Recurrent BCAR4 Gene Fusions across Solid Tumors.

Chromosomal rearrangements often result in active regulatory regions juxtaposed upstream of an oncogene to generate an expressed gene fusion. Repeated activation of a common downstream partner-with differing upstream regions across a patient cohort-suggests a conserved oncogenic role. Analysis of 9,638 patients across 32 solid tumor types revealed an annotated long noncoding RNA (lncRNA), Breast Cancer Anti-Estrogen Resistance 4 (BCAR4), was the most prevalent, uncharacterized, downstream gene fusion partner occurring in 11 cancers. Its oncogenic role was confirmed using multiple cell lines with endogenous BCAR4 gene fusions. Furthermore, overexpressing clinically prevalent BCAR4 gene fusions in untransformed cell lines was sufficient to induce an oncogenic phenotype. We show that the minimum common region to all gene fusions harbors an open reading frame that is necessary to drive proliferation. IMPLICATIONS: BCAR4 gene fusions represent an underappreciated class of gene fusions that may have biological and clinical implications across solid tumors.

Pan-cancer proteogenomic analysis reveals long and circular noncoding RNAs encoding peptides.

Recent studies show that annotated long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) encode for stable, functional peptides that contribute to human development and disease. To systematically discover lncRNAs and circRNAs encoding peptides, we performed a comprehensive integrative analysis of mass spectrometry-based proteomic and transcriptomic sequencing data from >900 patients across nine cancer types. This enabled us to identify 19,871 novel peptides derived from 8,903 lncRNAs. Further, we exploited open reading frames overlapping the backspliced region of circRNAs to identify 3,238 peptides that are uniquely derived from 2,834 circRNAs and not their corresponding linear RNAs. Collectively, our pan-cancer proteogenomic analysis will serve as a resource for evaluating the coding potential of lncRNAs and circRNAs that could aid future mechanistic studies exploring their function in cancer.

Mining biosynthetic gene clusters in Virgibacillus genomes.

BACKGROUND: Biosynthetic gene clusters produce a wide range of metabolites with activities that are of interest to the pharmaceutical industry. Specific interest is shown towards those metabolites that exhibit antimicrobial activities against multidrug-resistant bacteria that have become a global health threat. Genera of the phylum Firmicutes are frequently identified as sources of such metabolites, but the biosynthetic potential of its Virgibacillus genus is not known. Here, we used comparative genomic analysis to determine whether Virgibacillus strains isolated from the Red Sea mangrove mud in Rabigh Harbor Lagoon, Saudi Arabia, may be an attractive source of such novel antimicrobial agents. RESULTS: A comparative genomics analysis based on Virgibacillus dokdonensis Bac330, Virgibacillus sp. Bac332 and Virgibacillus halodenitrificans Bac324 (isolated from the Red Sea) and six other previously reported Virgibacillus strains was performed. Orthology analysis was used to determine the core genomes as well as the accessory genome of the nine Virgibacillus strains. The analysis shows that the Red Sea strain Virgibacillus sp. Bac332 has the highest number of unique genes and genomic islands compared to other genomes included in this study. Focusing on biosynthetic gene clusters, we show how marine isolates, including those from the Red Sea, are more enriched with nonribosomal peptides compared to the other Virgibacillus species. We also found that most nonribosomal peptide synthases identified in the Virgibacillus strains are part of genomic regions that are potentially horizontally transferred. CONCLUSIONS: The Red Sea Virgibacillus strains have a large number of biosynthetic genes in clusters that are not assigned to known products, indicating significant potential for the discovery of novel bioactive compounds. Also, having more modular synthetase units suggests that these strains are good candidates for experimental characterization of previously identified bioactive compounds as well. Future efforts will be directed towards establishing the properties of the potentially novel compounds encoded by the Red Sea specific trans-AT PKS/NRPS cluster and the type III PKS/NRPS cluster.

In silico exploration of Red Sea Bacillus genomes for natural product biosynthetic gene clusters.

BACKGROUND: The increasing spectrum of multidrug-resistant bacteria is a major global public health concern, necessitating discovery of novel antimicrobial agents. Here, members of the genus Bacillus are investigated as a potentially attractive source of novel antibiotics due to their broad spectrum of antimicrobial activities. We specifically focus on a computational analysis of the distinctive biosynthetic potential of Bacillus paralicheniformis strains isolated from the Red Sea, an ecosystem exposed to adverse, highly saline and hot conditions. RESULTS: We report the complete circular and annotated genomes of two Red Sea strains, B. paralicheniformis Bac48 isolated from mangrove mud and B. paralicheniformis Bac84 isolated from microbial mat collected from Rabigh Harbor Lagoon in Saudi Arabia. Comparing the genomes of B. paralicheniformis Bac48 and B. paralicheniformis Bac84 with nine publicly available complete genomes of B. licheniformis and three genomes of B. paralicheniformis, revealed that all of the B. paralicheniformis strains in this study are more enriched in nonribosomal peptides (NRPs). We further report the first computationally identified trans-acyltransferase (trans-AT) nonribosomal peptide synthetase/polyketide synthase (PKS/ NRPS) cluster in strains of this species. CONCLUSIONS: B. paralicheniformis species have more genes associated with biosynthesis of antimicrobial bioactive compounds than other previously characterized species of B. licheniformis, which suggests that these species are better potential sources for novel antibiotics. Moreover, the genome of the Red Sea strain B. paralicheniformis Bac48 is more enriched in modular PKS genes compared to B. licheniformis strains and other B. paralicheniformis strains. This may be linked to adaptations that strains surviving in the Red Sea underwent to survive in the relatively hot and saline ecosystems.

Complete genome sequence of Mycobacterium fortuitum subsp. fortuitum type strain DSM46621.

Mycobacterium fortuitum is a member of the rapidly growing nontuberculous mycobacteria (NTM). It is ubiquitous in water and soil habitats, including hospital environments. M. fortuitum is increasingly recognized as an opportunistic nosocomial pathogen causing disseminated infection. Here we report the genome sequence of M. fortuitum subsp. fortuitum type strain DSM46621.