Viral reprogramming of host transcription initiation.

Viruses are master remodelers of the host cell environment in support of infection and virus production. For example, viruses typically regulate cell gene expression through modulating canonical cell promoter activity. Here, we show that Epstein Barr virus (EBV) replication causes 'de novo' transcription initiation at 29674 new transcription start sites throughout the cell genome. De novo transcription initiation is facilitated in part by the unique properties of the viral pre-initiation complex (vPIC) that binds a TATT[T/A]AA, TATA box-like sequence and activates transcription with minimal support by additional transcription factors. Other de novo promoters are driven by the viral transcription factors, Zta and Rta and are influenced by directional proximity to existing canonical cell promoters, a configuration that fosters transcription through existing promoters and transcriptional interference. These studies reveal a new way that viruses interact with the host transcriptome to inhibit host gene expression and they shed light on primal features driving eukaryotic promoter function.

Transcriptomic analysis identifies B-lymphocyte kinase as a therapeutic target for desmoplastic small round cell tumor cancer stem cell-like cells.

Desmoplastic small round cell tumor (DSRCT) is an aggressive pediatric cancer caused by the EWSR1-WT1 fusion oncoprotein. The tumor is refractory to treatment with a 5-year survival rate of only 15-25%, necessitating the development of novel therapeutics, especially those able to target chemoresistant subpopulations. Novel in vitro cancer stem cell-like (CSC-like) culture conditions increase the expression of stemness markers (SOX2, NANOG) and reduce DSRCT cell line susceptibility to chemotherapy while maintaining the ability of DSRCT cells to form xenografts. To gain insights into this chemoresistant model, RNA-seq was performed to elucidate transcriptional alterations between DSRCT cells grown in CSC-like spheres and normal 2-dimensional adherent state. Commonly upregulated and downregulated genes were identified and utilized in pathway analysis revealing upregulation of pathways related to chromatin assembly and disassembly and downregulation of pathways including cell junction assembly and extracellular matrix organization. Alterations in chromatin assembly suggest a role for epigenetics in the DSRCT CSC-like state, which was further investigated with ATAC-seq, identifying over 10,000 differentially accessible peaks, including 4444 sphere accessible peaks and 6,120 adherent accessible peaks. Accessible regions were associated with higher gene expression, including increased accessibility of the CSC marker SOX2 in CSC-like culture conditions. These analyses were further utilized to identify potential CSC therapeutic targets, leading to the identification of B-lymphocyte kinase (BLK) as a CSC-enriched, EWSR1-WT1-regulated, druggable target. BLK inhibition and knockdown reduced CSC-like properties, including abrogation of tumorsphere formation and stemness marker expression. Importantly, BLK knockdown reduced DSRCT CSC-like cell chemoresistance, making its inhibition a promising target for future combination therapy.

Loss of feedback regulation between FAM3B and androgen receptor driving prostate cancer progression.

BACKGROUND: Although the fusion of the transmembrane serine protease 2 gene (TMPRSS2) with the erythroblast transformation-specific-related gene (ERG), or TMPRSS2-ERG, occurs frequently in prostate cancer, its impact on clinical outcomes remains controversial. Roughly half of TMPRSS2-ERG fusions occur through intrachromosomal deletion of interstitial genes and the remainder via insertional chromosomal rearrangements. Because prostate cancers with deletion-derived TMPRSS2-ERG fusions are more aggressive than those with insertional fusions, we investigated the impact of interstitial gene loss on prostate cancer progression. METHODS: We conducted an unbiased analysis of transcriptome data from large collections of prostate cancer samples and employed diverse in vitro and in vivo models combined with genetic approaches to characterize the interstitial gene loss that imposes the most important impact on clinical outcome. RESULTS: This analysis identified FAM3B as the top-ranked interstitial gene whose loss is associated with a poor prognosis. The association between FAM3B loss and poor clinical outcome extended to fusion-negative prostate cancers where FAM3B downregulation occurred through epigenetic imprinting. Importantly, FAM3B loss drives disease progression in prostate cancer. FAM3B acts as an intermediator of a self-governing androgen receptor feedback loop. Specifically, androgen receptor upregulates FAM3B expression by binding to an intronic enhancer to induce an enhancer RNA and facilitate enhancer-promoter looping. FAM3B, in turn, attenuates androgen receptor signaling. CONCLUSION: Loss of FAM3B in prostate cancer, whether through the TMPRSS2-ERG translocation or epigenetic imprinting, causes an exit from this autoregulatory loop to unleash androgen receptor activity and prostate cancer progression. These findings establish FAM3B loss as a new driver of prostate cancer progression and support the utility of FAM3B loss as a biomarker to better define aggressive prostate cancer.

A novel role of TRIM28 B box domain in L1 retrotransposition and ORF2p-mediated cDNA synthesis.

The long interspersed element 1 (LINE-1 or L1) integration is affected by many cellular factors through various mechanisms. Some of these factors are required for L1 amplification, while others either suppress or enhance specific steps during L1 propagation. Previously, TRIM28 has been identified to suppress transposable elements, including L1 expression via its canonical role in chromatin remodeling. Here, we report that TRIM28 through its B box domain increases L1 retrotransposition and facilitates shorter cDNA and L1 insert generation in cultured cells. Consistent with the latter, we observe that tumor specific L1 inserts are shorter in endometrial, ovarian, and prostate tumors with higher TRIM28 mRNA expression than in those with lower TRIM28 expression. We determine that three amino acids in the B box domain that are involved in TRIM28 multimerization are critical for its effect on both L1 retrotransposition and cDNA synthesis. We provide evidence that B boxes from the other two members in the Class VI TRIM proteins, TRIM24 and TRIM33, also increase L1 retrotransposition. Our findings could lead to a better understanding of the host/L1 evolutionary arms race in the germline and their interplay during tumorigenesis.

Metabolic pathways enriched according to ERG status are associated with biochemical recurrence in Hispanic/Latino patients with prostate cancer.

BACKGROUND: The role of ERG-status molecular subtyping in prognosis of prostate cancer (PCa) is still under debate. In this study, we identified differentially expressed genes (DEGs) according to ERG-status to explore their enriched pathways and implications in prognosis in Hispanic/Latino PCa patients. METHODS: RNA from 78 Hispanic PCa tissues from radical prostatectomies (RP) were used for RNA-sequencing. ERGhigh /ERGlow tumor groups were determined based on the 1.5-fold change median expression in non-tumor samples. DEGs with a False Discovery Rate (FDR) < 0.01 and a fold change >2 were identified between ERGhigh and ERGlow tumors and submitted to enrichment analysis in MetaCore. Survival and association analyses were performed to evaluate biochemical recurrence (BCR)-free survival. RESULTS: The identification of 150 DEGs between ERGhigh and ERGlow tumors revealed clustering of most of the non-BCR cases (60%) into de ERGhigh group and most of the BCR cases (60.8%) in ERGlow group. Kaplan-Meier survival curves showed a worst BCR-free survival for ERGlow patients, and a significant reduced risk of BCR was observed for ERGhigh cases (OR = 0.29 (95%CI, 0.10-0.8)). Enrichment pathway analysis identified metabolic-related pathways, such as the renin-angiotensin system and angiotensin maturation system, the linoleic acid metabolism, and polyamines metabolism in these ERG groups. CONCLUSIONS: ERGlow tumor cases were associated with poor BCR-free survival in our Hispanic/Latino patients, with metabolism-related pathways altered in the BCR progression. IMPACT: Our findings suggest the need to dissect the role of diet, metabolism, and lifestyle as risk factors for more aggressive PCa subtypes.

Apremilast in Recalcitrant Cutaneous Dermatomyositis: A Nonrandomized Controlled Trial.

Importance: Cutaneous disease in dermatomyositis has no standardized treatment approach and so presents a challenging task for patients and clinicians. Objective: To study the efficacy and safety of apremilast as an add-on therapy in patients with recalcitrant cutaneous dermatomyositis. Design, Setting, and Participants: This phase 2a, open-label, single-arm nonrandomized controlled trial was conducted at a single center from June 2018 to June 2021. Participants were 8 patients with recalcitrant cutaneous dermatomyositis, defined by a cutaneous disease activity severity index (CDASI) score greater than 5 despite treatment with steroids, steroid-sparing agents, or both. Data were analyzed from June 2018 to June 2021. Interventions: Apremilast 30 mg orally twice daily was added to ongoing treatment regimens. Main Outcomes and Measures: The primary outcome was the overall response rate (ORR) at 3 months. Key secondary outcomes were the safety and toxicity of apremilast and the durability of response at 6 months. The CDASI, muscle score, dermatology life quality index (DLQI), and depression assessments were performed at baseline and regularly until month 7. Skin biopsies were performed at baseline and 3 months after apremilast (defined as 3 months into active apremilast therapy) and tested for gene expression profiling and immunohistochemical stains. Adverse events were assessed using the Common Terminology Criteria for Adverse Events version 5.0. Results: Among 8 patients with recalcitrant cutaneous dermatomyositis (all women; mean [SD] age, 54 [15.9] years), a response was found at 3 months after apremilast among 7 patients (ORR, 87.5%). The mean (SD) decrease in CDASI was 12.9 (6.3) points at 3 months (P < .001). Apremilast was well tolerated, with no grade 3 or higher adverse events. Sequencing of RNA was performed on skin biopsies taken from 7 patients at baseline and at 3 months after therapy. Appropriate negative (ie, no primary antibody) and positive (ie, tonsil and spleen) controls were stained in parallel with each set of slides studied. Of 39 076 expressed genes, there were 195 whose expression changed 2-fold or more at P < .01 (123 downregulated and 72 upregulated genes). Gene set enrichment analysis identified 13 pathways in which apremilast was associated with downregulated expression, notably signal transducers and activators of transcription 1 (STAT1), STAT3, interleukin 4 (IL-4), IL-6, IL-12, IL-23, interferon γ (IFNγ), and tumor necrosis factor α (TNFα) pathways. In immunohistochemical staining, there was a mean (SD) decrease in phosphorylation levels STAT1 (22.3% [28.3%] positive cells) and STAT3 (13.4% [11.6%] positive cells) at the protein level, a downstream signaling pathway for the downregulated cytokines. Conclusions and Relevance: These findings suggest that apremilast was a safe and efficacious add-on treatment in recalcitrant dermatomyositis, with an overall response rate of 87.5% and associations with downregulation of multiple inflammatory pathways. Trial Registration: ClinicalTrials.gov Identifier: NCT03529955.

SCIFER: approach for analysis of LINE-1 mRNA expression in single cells at a single locus resolution.

BACKGROUND: Endogenous expression of L1 mRNA is the first step in an L1-initiated mutagenesis event. However, the contribution of individual cell types to patterns of organ-specific L1 mRNA expression remains poorly understood, especially at single-locus resolution. We introduce a method to quantify expression of mobile elements at the single-locus resolution in scRNA-Seq datasets called Single Cell Implementation to Find Expressed Retrotransposons (SCIFER). SCIFER aligns scRNA-Seq reads uniquely to the genome and extracts alignments from single cells by cell-specific barcodes. In contrast to the alignment performed using default parameters, this alignment strategy increases accuracy of L1 locus identification by retaining only reads that are uniquely mapped to individual L1 loci. L1 loci expressed in single cells are unambiguously identified using a list of L1 loci manually validated to be expressed in bulk RNA-Seq datasets generated from the same cell line or organ. RESULTS: Validation of SCIFER using MCF7 cells determined technical parameters needed for optimal detection of L1 expression in single cells. We show that unsupervised analysis of L1 expression in single cells exponentially inflates both the levels of L1 expression and the number of expressed L1 loci. Application of SCIFER to analysis of scRNA-Seq datasets generated from mouse and human testes identified that mouse Round Spermatids and human Spermatogonia, Spermatocytes, and Round Spermatids express the highest levels of L1 mRNA. Our analysis also determined that similar to mice, human testes from unrelated individuals share as much as 80% of expressed L1 loci. Additionally, SCIFER determined that individual mouse cells co-express different L1 sub-families and different families of transposable elements, experimentally validating their co-existence in the same cell. CONCLUSIONS: SCIFER detects mRNA expression of individual L1 loci in single cells. It is compatible with scRNA-Seq datasets prepared using traditional sequencing methods. Validated using a human cancer cell line, SCIFER analysis of mouse and human testes identified key cell types supporting L1 expression in these species. This will further our understanding of differences and similarities in endogenous L1 mRNA expression patterns in mice and humans.

Ebola virus delta peptide is an enterotoxin.

During the 2013-2016 West African (WA) Ebola virus (EBOV) outbreak, severe gastrointestinal symptoms were common in patients and associated with poor outcome. Delta peptide is a conserved product of post-translational processing of the abundant EBOV soluble glycoprotein (sGP). The murine ligated ileal loop model was used to demonstrate that delta peptide is a potent enterotoxin. Dramatic intestinal fluid accumulation follows injection of biologically relevant amounts of delta peptide into ileal loops, along with gross alteration of villous architecture and loss of goblet cells. Transcriptomic analyses show that delta peptide triggers damage response and cell survival pathways and downregulates expression of transporters and exchangers. Induction of diarrhea by delta peptide occurs via cellular damage and regulation of genes that encode proteins involved in fluid secretion. While distinct differences exist between the ileal loop murine model and EBOV infection in humans, these results suggest that delta peptide may contribute to EBOV-induced gastrointestinal pathology.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells.

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of 'transitional' L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

Reversal of splicing infidelity is a pre-activation step in B cell differentiation.

Introduction: B cell activation and differentiation is central to the adaptive immune response. Changes in exon usage can have major impacts on cellular signaling and differentiation but have not been systematically explored in differentiating B cells. Methods: We analyzed exon usage and intron retention in RNA-Seq data from subsets of human B cells at various stages of differentiation, and in an in vitro laboratory model of B cell activation and differentiation (Epstein Barr virus infection). Results: Blood naïve B cells were found to have an unusual splicing profile, with unannotated splicing events in over 30% of expressed genes. Splicing changed substantially upon naïve B cell entry into secondary lymphoid tissue and before activation, involving significant increases in exon commitment and reductions in intron retention. These changes preferentially involved short introns with weak splice sites and were likely mediated by an overall increase in splicing efficiency induced by the lymphoid environment. The majority of transcripts affected by splicing changes showed restoration of encoded conserved protein domains and/or reduced targeting to the nonsense-mediated decay pathway. Affected genes were enriched in functionally important immune cell activation pathways such as antigen-mediated signaling, cell cycle control and mRNA processing and splicing. Discussion: Functional observations from donor B cell subsets in progressive states of differentiation and from timecourse experiments using the in vitro model suggest that these widespread changes in mRNA splicing play a role in preparing naïve B cells for the decisive step of antigen-mediated activation and differentiation.

Elevated ATGL in colon cancer cells and cancer stem cells promotes metabolic and tumorigenic reprogramming reinforced by obesity.

Obesity is a worldwide epidemic associated with increased risk and progression of colon cancer. Here, we aimed to determine the role of adipose triglyceride lipase (ATGL), responsible for intracellular lipid droplet (LD) utilization, in obesity-driven colonic tumorigenesis. In local colon cancer patients, significantly increased ATGL levels in tumor tissue, compared to controls, were augmented in obese individuals. Elevated ATGL levels in human colon cancer cells (CCC) relative to non-transformed were augmented by an obesity mediator, oleic acid (OA). In CCC and colonospheres, enriched in colon cancer stem cells (CCSC), inhibition of ATGL prevented LDs utilization and inhibited OA-stimulated growth through retinoblastoma-mediated cell cycle arrest. Further, transcriptomic analysis of CCC, with inhibited ATGL, revealed targeted pathways driving tumorigenesis, and high-fat-diet obesity facilitated tumorigenic pathways. Inhibition of ATGL in colonospheres revealed targeted pathways in human colonic tumor crypt base cells (enriched in CCSC) derived from colon cancer patients. In CCC and colonospheres, we validated selected transcripts targeted by ATGL inhibition, some with emerging roles in colonic tumorigeneses (ATG2B, PCK2, PGAM1, SPTLC2, IGFBP1, and ABCC3) and others with established roles (MYC and MUC2). These findings demonstrate obesity-promoted, ATGL-mediated colonic tumorigenesis and establish the therapeutic significance of ATGL in obesity-reinforced colon cancer progression.

SON inhibits megakaryocytic differentiation via repressing RUNX1 and the megakaryocytic gene expression program in acute megakaryoblastic leukemia.

A high incidence of acute megakaryoblastic leukemia (AMKL) in Down syndrome patients implies that chromosome 21 genes have a pivotal role in AMKL development, but the functional contribution of individual genes remains elusive. Here, we report that SON, a chromosome 21-encoded DNA- and RNA-binding protein, inhibits megakaryocytic differentiation by suppressing RUNX1 and the megakaryocytic gene expression program. As megakaryocytic progenitors differentiate, SON expression is drastically reduced, with mature megakaryocytes having the lowest levels. In contrast, AMKL cells express an aberrantly high level of SON, and knockdown of SON induced the onset of megakaryocytic differentiation in AMKL cell lines. Genome-wide transcriptome analyses revealed that SON knockdown turns on the expression of pro-megakaryocytic genes while reducing erythroid gene expression. Mechanistically, SON represses RUNX1 expression by directly binding to the proximal promoter and two enhancer regions, the known +23 kb enhancer and the novel +139 kb enhancer, at the RUNX1 locus to suppress H3K4 methylation. In addition, SON represses the expression of the AP-1 complex subunits JUN, JUNB, and FOSB which are required for late megakaryocytic gene expression. Our findings define SON as a negative regulator of RUNX1 and megakaryocytic differentiation, implicating SON overexpression in impaired differentiation during AMKL development.

Somatic mutations in the DNA repairome in prostate cancers in African Americans and Caucasians.

Most hereditary tumors show aberrations in DNA repair genes or their regulators. In contrast, only a minority of sporadic tumors show alterations in these genes. As a result, genomic instability is currently considered an enhancer of tumorigenesis rather than an obligatory event in this process. However, tumor heterogeneity presents a significant technical challenge for most cancer genomics studies performed at less than 100× mean resolution depth. To address the importance of genomic instability in prostate carcinogenesis and tumor progression, we performed ultrahigh depth exome sequencing of 124 DNA damage repair/response (repairome) genes in 63 tumors and matched normal tissue samples in African Americans and Caucasians. The average sequence depth was 712-fold for DNA isolated from normal tissue and 368-fold for FFPE tumors. We identified 671 somatic mutations in tumors from African Americans and 762 somatic mutations in tumors in Caucasians. The most frequently mutated DNA repairome genes were EXO1, ATR, POLQ, NEIL3, ERCC6, BRCA2, BRCA1, XPC, JAG1, RPA1, POLE, ATM, and LIG1 in African American men, and POLQ, NEIL3, POLB, BRCA2, EXO1, ERCC6, ATR, RBBP8, BRCA1, ATM, JAG1, XPC, and POLE in Caucasians. We found that 89% of tumors had at least one mutation in nucleotide excision repair pathway genes in African Americans, whereas >40% of tumors had mutations in base excision repair pathway genes in Caucasians. We further identified a marginal increase in mutation rate in tumors in African Americans with increasing age. Tumors in Caucasians did not show a correlation with age, but a progressive increase in the mutation rate was observed at higher Gleason scores. Our data reveal significant differences in the molecular signatures in the DNA repairome in prostate cancer between African Americans and Caucasians. These data also have substantial implications regarding the well-known health disparities in prostate cancer, such as the higher mortality in African Americans than Caucasians.

Bacterial TLR4 and NOD2 signaling linked to reduced mitochondrial energy function in active inflammatory bowel disease.

Inflammatory bowel disease (IBD) has been linked to active signaling with bacterial components and reduced mitochondrial ATP production; however, synergism between both of these disease characteristics remains unclear. We aimed to determine in human IBD transcriptomes the link between a transcriptional signature unique to intestinal cells (ICs) with reduced mitochondrial ATP production (Mito-0) and bacteria triggered signaling using a bioinformatics approach. We generated an IC Mito-0 panel comprised of 199 differentially expressed (DE) transcripts mediated by reduced mitochondrial ATP function (DEGseq, log2 fold-change > |2|, p < .001). Transcripts from this panel were involved in diverse biological functions including regulation of mitochondrial energy (lower ATP), extracellular matrix, cell-cell contact, cytoskeleton, growth, metabolism, and inflammation. Next, unsupervised hierarchical clustering showed that the Mito-0 panel distinctly separated inflamed IBD from non-inflamed transcriptomes, which was also supported by principal component analysis (PCA) revealing distinct variation between sample types based on presence of the Mito-0 signature (PCA, p = 8.77e-09). Utilizing three independent IBD cohorts, we validated that 60 novel transcripts from the Mito-0 panel were significantly increased in inflamed tissue. Subsequently, KEGG generated bacterial TLR4 and NOD2 transcriptional signatures strongly associated with inflamed IBD transcriptomes and with the Mito-0 signature as determined by Spearman's analysis (coefficient of correlation, r = 0.92, p < .05). Herein, using a comprehensive analysis we demonstrated existence of an axis between bacteria triggered signaling and reduced mitochondrial energy function. Furthermore, we identified and validated novel transcripts within this axis as potential drivers and therapeutic targets for human IBD.

Circular RNAs add diversity to androgen receptor isoform repertoire in castration-resistant prostate cancer.

Deregulated expression of circular RNAs (circRNAs) is associated with various human diseases, including many types of cancer. Despite their growing links to cancer, there has been limited characterization of circRNAs in metastatic castration-resistant prostate cancer, the major cause of prostate cancer mortality. Here, through the analysis of an exome-capture RNA-seq dataset from 47 metastatic castration-resistant prostate cancer samples and ribodepletion and RNase R RNA-sequencing of patient-derived xenografts (PDXs) and cell models, we identified 13 circRNAs generated from the key prostate cancer driver gene-androgen receptor (AR). We validated and characterized the top four most abundant, clinically relevant AR circRNAs. Expression of these AR circRNAs was upregulated during castration-resistant progression of PDXs. The upregulation was not due to global increase of circRNA formation in these tumors. Instead, the levels of AR circRNAs correlated strongly with that of the linear AR transcripts (both AR and AR variants) in clinical samples and PDXs, indicating a transcriptional mechanism of regulation. In cultured cells, androgen suppressed the expression of these AR circRNAs and the linear AR transcripts, and the suppression was attenuated by an antiandrogen. Using nuclear/cytoplasmic fractionation and RNA in-situ hybridization assays, we demonstrated predominant cytoplasmic localization of these AR circRNAs, indicating likely cytoplasmic functions. Overall, this is the first comprehensive characterization of circRNAs arising from the AR gene. With greater resistance to exoribonuclease compared to the linear AR transcripts and detectability of AR circRNAs in patient plasma, these AR circRNAs may serve as surrogate circulating markers for AR/AR-variant expression and castration-resistant prostate cancer progression.

RNA Next-Generation Sequencing and a Bioinformatics Pipeline to Identify Expressed LINE-1s at the Locus-Specific Level.

Long INterspersed Elements-1 (LINEs/L1s) are repetitive elements that can copy and randomly insert in the genome resulting in genomic instability and mutagenesis. Understanding the expression patterns of L1 loci at the individual level will lend to the understanding of the biology of this mutagenic element. This autonomous element makes up a significant portion of the human genome with over 500,000 copies, though 99% are truncated and defective. However, their abundance and dominant number of defective copies make it challenging to identify authentically expressed L1s from L1-related sequences expressed as part of other genes. It is also challenging to identify which specific L1 locus is expressed due to the repetitive nature of the elements. Overcoming these challenges, we present an RNA-Seq bioinformatic approach to identify L1 expression at the locus specific level. In summary, we collect cytoplasmic RNA, select for polyadenylated transcripts, and utilize strand-specific RNA-Seq analyses to uniquely map reads to L1 loci in the human reference genome. We visually curate each L1 locus with uniquely mapped reads to confirm transcription from its own promoter and adjust mapped transcript reads to account for mappability of each individual L1 locus. This approach was applied to a prostate tumor cell line, DU145, to demonstrate the ability of this protocol to detect expression from a small number of the full-length L1 elements.

A positive role of c-Myc in regulating androgen receptor and its splice variants in prostate cancer.

Increased expression of the full-length androgen receptor (AR-FL) and AR splice variants (AR-Vs) drives the progression of castration-resistant prostate cancer (CRPC). The levels of AR-FL and AR-V transcripts are often tightly correlated in individual CRPC samples, yet our understanding of how their expression is co-regulated is limited. Here, we report a role of c-Myc in accounting for coordinated AR-FL and AR-V expression. Analysis of gene-expression data from 159 metastatic CRPC samples and 2142 primary prostate tumors showed that the level of c-Myc is positively correlated with that of individual AR isoforms. A striking positive correlation also exists between the activity of the c-Myc pathway and the level of individual AR isoforms, between the level of c-Myc and the activity of the AR pathway, and between the activities of the two pathways. Moreover, the c-Myc signature is highly enriched in tumors expressing high levels of AR, as is the AR signature in c-Myc-high-expressing tumors. Using shRNA knockdown, we confirmed c-Myc regulation of expression and activity of AR-FL and AR-Vs in cell models and a patient-derived xenograft model. Mechanistically, c-Myc promotes the transcription of the AR gene and enhances the stability of the AR-FL and AR-V proteins without altering AR RNA splicing. Importantly, inhibiting c-Myc sensitizes enzalutamide-resistant cells to growth inhibition by enzalutamide. Overall, this study highlights a critical role of c-Myc in regulating the coordinated expression of AR-FL and AR-Vs that is commonly observed in CRPC and suggests the utility of targeting c-Myc as an adjuvant to AR-directed therapy.

Loss of Forkhead Box O3 Facilitates Inflammatory Colon Cancer: Transcriptome Profiling of the Immune Landscape and Novel Targets.

BACKGROUND & AIMS: Diminished forkhead box O3 (FOXO3) function drives inflammation and cancer growth; however, mechanisms fostering these pathobiologies are unclear. Here, we aimed to identify in colon loss of FOXO3-dependent cellular and molecular changes that facilitate inflammation-mediated tumor growth. METHODS: FOXO3 knockout (KO) and wild-type (WT) mice were used in the AOM/DSS model of inflammation-mediated colon cancer. Bioinformatics were used for profiling of mRNA sequencing data from human and mouse colon and tumors; specific targets were validated in human colon cancer cells (shFOXO3). RESULTS: In mice, FOXO3 deficiency led to significantly elevated colonic tumor burden (incidence and size) compared with WT (P < .05). In FOXO3 KO colon, activated molecular pathways overlapped with those associated with mouse and human colonic inflammation and cancer, especially human colonic tumors with inflammatory microsatellite instability (false discovery rate < 0.05). FOXO3 KO colon, similar to tumors, had increased neutrophils, macrophages, B cells, T cells, and decreased natural killer cells (false discovery rate < 0.05). Moreover, in KO colon differentially expressed transcripts were linked to activation of inflammatory nuclear factor kappa B, tumorigenic cMyc, and bacterial Toll-like receptor signaling. Among differentially expressed transcripts, we validated altered expression of integrin subunit alpha 2 (ITGA2), ADAM metallopeptidase with thrombospondin type 1 motif 12, and ST8 alpha-N-acetyl-neuraminide alpha-2,8-sialyltransferase 5 in mouse WT and FOXO3 KO colon and tumors (P < .05). Similarly, their altered expression was found in human inflammatory bowel disease and colon cancer tissues and linked to poor patient survival. Ultimately, in human colon cancer cells, FOXO3 knockdown (shFOXO3) led to significantly increased ITGA2, and silencing ITGA2 (siRNA) alone diminished cell growth. CONCLUSIONS: We identified the loss of FOXO3-mediated immune landscape, pathways, and transcripts that could serve as biomarkers and new targets for inflammatory colon cancer treatment.

The Epstein Barr virus circRNAome.

Our appreciation for the extent of Epstein Barr virus (EBV) transcriptome complexity continues to grow through findings of EBV encoded microRNAs, new long non-coding RNAs as well as the more recent discovery of over a hundred new polyadenylated lytic transcripts. Here we report an additional layer to the EBV transcriptome through the identification of a repertoire of latent and lytic viral circular RNAs. Utilizing RNase R-sequencing with cell models representing latency types I, II, and III, we identified EBV encoded circular RNAs expressed from the latency Cp promoter involving backsplicing from the W1 and W2 exons to the C1 exon, from the EBNA BamHI U fragment exon, and from the latency long non-coding RPMS1 locus. In addition, we identified circular RNAs expressed during reactivation including backsplicing from exon 8 to exon 2 of the LMP2 gene and a highly expressed circular RNA derived from intra-exonic backsplicing within the BHLF1 gene. While expression of most of these circular RNAs was found to depend on the EBV transcriptional program utilized and the transcription levels of the associated loci, expression of LMP2 exon 8 to exon 2 circular RNA was found to be cell model specific. Altogether we identified over 30 unique EBV circRNAs candidates and we validated and determined the structural features, expression profiles and nuclear/cytoplasmic distributions of several predominant and notable viral circRNAs. Further, we show that two of the EBV circular RNAs derived from the RPMS1 locus are detected in EBV positive clinical stomach cancer specimens. This study increases the known EBV latency and lytic transcriptome repertoires to include viral circular RNAs and it provides an essential foundation and resource for investigations into the functions and roles of this new class of EBV transcripts in EBV biology and diseases.