Search

Kaposi's sarcoma-associated herpesvirus is the etiologic agent of Kaposi's sarcoma and two B-cell malignancies. Recent advancements in sequencing technologies have led to high resolution transcriptomes for several human herpesviruses that densely encode genes on both strands. However, for KSHV progress remained limited due to the overall low percentage of KSHV transcripts, even during lytic replication. To address this challenge, we have developed a target enrichment method to increase the KSHV-specific reads for both short- and long-read sequencing platforms. Furthermore, we combined this approach with the Transcriptome Resolution through Integration of Multi-platform Data (TRIMD) pipeline developed previously to annotate transcript structures. TRIMD first builds a scaffold based on long-read sequencing and validates each transcript feature with supporting evidence from Illumina RNA-Seq and deepCAGE sequencing data. Our stringent innovative approach identified 994 unique KSHV transcripts, thus providing the first high-density KSHV lytic transcriptome. We describe a plethora of novel coding and non-coding KSHV transcript isoforms with alternative untranslated regions, splice junctions and open-reading frames, thus providing deeper insights on gene expression regulation of KSHV. Interestingly, as described for Epstein-Barr virus, we identified transcription start sites that augment long-range transcription and may increase the number of latency-associated genes potentially expressed in KS tumors., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

B-cell acute lymphoblastic leukemia (B-ALL) remains a hard-to-treat disease with a poor prognosis in adults. Mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1) is a para-caspase required for B-cell receptor (BCR)-mediated NF-κB activation. Inhibition of MALT1 in preclinical models has proven efficacious in many B-cell malignancies including chronic lymphocytic leukemia, mantle cell lymphoma and diffuse large B-cell lymphoma. We sought to examine the role of MALT1 in B-ALL and determine the biological consequences of its inhibition. Targeting MALT1 with both Z-VRPR-fmk and MI-2 efficiently kills B-ALL cells independent of the cell-of-origin (pro, pre, mature) or the presence of the Philadelphia chromosome, and spares normal B cells. The mechanism of cell death was through apoptotic induction, mostly in cycling cells. The proteolytic activity of MALT1 can be studied by measuring its ability to cleave its substrates. Surprisingly, with the exception of mature B-ALL, we did not detect cleavage of MALT1 substrates at baseline, nor after proteasomal inhibition or following activation of pre-BCR. To explore the possibility of a distinct role for MALT1 in B-ALL, independent of signaling through BCR, we studied the changes in gene expression profiling following a 24-hour treatment with MI-2 in 12 B-ALL cell lines. Our transcriptome analysis revealed a strong inhibitory effect on MYC-regulated gene signatures, further confirmed by Myc protein downregulation, concomitant with an increase in the Myc degrader FBXW7. In conclusion, our evidence suggests a novel role for MALT1 in B-ALL through Myc regulation and provides support for clinical testing of MALT1 inhibitors in B-ALL.

The calcium-sensing receptor (CaSR), a G protein-coupled receptor, regulates Ca 2+ concentration in plasma by regulating parathyroid hormone secretion. In other tissues, it is reported to play roles in cellular differentiation and migration and in secretion and absorption. We reported previously that CaSR can be conditionally deleted in the mouse esophagus. This conditional knockout (KO) ( Eso CaSR -/- ) model showed a significant reduction in the levels of adherens and tight junction proteins and had a marked buildup of bacteria on the luminal esophageal surface. To further examine the role of CaSR, we used RNA sequencing to determine gene expression profiles in esophageal epithelia of control and Eso CaSR -/- mice RNA Seq data indicated upregulation of gene sets involved in DNA replication and cell cycle in Eso CaSR -/- . This is accompanied by the downregulation of gene sets involved in the innate immune response and protein homeostasis including peptide elongation and protein trafficking. Ingenuity pathway analysis (IPA) demonstrated that these genes are mapped to important biological networks including calcium and Ras homologus A (RhoA) signaling pathways. To further explore the bacterial buildup in Eso CaSR -/- esophageal tissue, 16S sequencing of the mucosal-associated bacterial microbiome was performed. Three bacterial species, g_Rodentibacter , s _Rodentibacter_unclassified , and s_Lactobacillus_hilgardi were significantly increased in Eso CaSR -/- . Furthermore, metagenomic analysis of 16S sequences indicated that pathways related to oxidative phosphorylation and metabolism were downregulated in Eso CaSR -/- tissues. These data demonstrate that CaSR impacts major pathways of cell proliferation, differentiation, cell cycle, and innate immune response in esophageal epithelium. The disruption of these pathways causes inflammation and significant modifications of the microbiome. NEW & NOTEWORTHY Calcium-sensing receptor (CaSR) plays a significant role in maintaining the barrier function of esophageal epithelium. Using RNA sequencing, we show that conditional deletion of CaSR from mouse esophagus causes upregulation of genes involved in DNA replication and cell cycle and downregulation of genes involved in the innate immune response, protein translation, and cellular protein synthesis. Pathway analysis shows disruption of signaling pathways of calcium and actin cytoskeleton. These changes caused inflammation and esophageal dysbiosis.

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., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

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., (© 2024. The Author(s).)

As a fundamental aspect of normal cell signaling and disease states, there is great interest in determining alternative splicing (AS) changes in physiologic, pathologic, and pharmacologic settings. High throughput RNA sequencing and specialized software to detect AS has greatly enhanced our ability to determine transcriptome-wide splicing changes. Despite the richness of this data, deriving meaning from sometimes thousands of AS events is a substantial bottleneck for most investigators. We present SpliceTools, a suite of data processing modules that arms investigators with the ability to quickly produce summary statistics, mechanistic insights, and functional significance of AS changes through command line or through an online user interface. Utilizing RNA-seq datasets for 186 RNA binding protein knockdowns, nonsense mediated RNA decay inhibition, and pharmacologic splicing inhibition, we illustrate the utility of SpliceTools to distinguish splicing disruption from regulated transcript isoform changes, we show the broad transcriptome footprint of the pharmacologic splicing inhibitor, indisulam, we illustrate the utility in uncovering mechanistic underpinnings of splicing inhibition, we identify predicted neo-epitopes in pharmacologic splicing inhibition, and we show the impact of splicing alterations induced by indisulam on cell cycle progression. Together, SpliceTools puts rapid and easy downstream analysis at the fingertips of any investigator studying AS., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 O’Grady, Baddoo, Flemington, Ishaq, Ungerleider and Flemington.)

Cell lines represent an essential tool used in preclinical research. Most hematologic malignancies have a wide array of cell lines representing their respective molecular and pathologic spectra. In mantle cell lymphoma (MCL), cell lines become specifically valuable in view of the heterogeneity of this disease. Unfortunately, the number of MCL cell lines that are available for the research community remains small, with only nine cell lines available for purchase through the American Type Culture Collection (ATCC). We have established a novel blastoid MCL cell line, isolated from the malignant pleural effusion of a 69-year-old male with refractory MCL. Arbo was fully characterized with cytogenetics, immunophenotyping, whole exome sequencing and drug sensitivity assays. One of the most notable mutations identified in Arbo (but not in normal tissue) was the missense mutation NOTCH2 R2400*, which has been proposed as a clinically significant mutation in MCL seen in 5% of cases. NOTCH2 R2400* results in a truncated Notch2 protein, leading to a more stable and active protein. Using pharmacologic inhibition of Notch2, we showed a dependence of Arbo on NOTCH2 signaling, as well as a link between CD23 expression on Arbo and NOTCH2 activity. Arbo represents a NOTCH2 mutated model that is useful in MCL as well as other lymphomas with such mutation. We plan to deposit Arbo at the ATCC to be available for the research community., Competing Interests: The authors declare they have no conflicts of interest., (© 2022 The Authors. eJHaem published by British Society for Haematology and John Wiley & Sons Ltd.)

The oncogenic gammaherpesviruses, including human Epstein-Barr virus (EBV), human Kaposi's sarcoma-associated herpesvirus (KSHV), and murine gammaherpesvirus 68 (MHV68, γHV68, MuHV-4), are associated with numerous malignancies, including B cell lymphomas and nasopharyngeal carcinoma. These viruses employ numerous molecular strategies to colonize the host, including the expression of noncoding RNAs (ncRNAs). As the first viral ncRNAs identified, EBV-encoded RNA 1 and 2 (EBER1 and EBER2, respectively) have been investigated extensively for decades; however, their specific in vivo functions remain largely unknown. In work here, we used chimeric MHV68 viruses in an in vivo complementation system to test whether EBV EBER2 contributes to acute and/or chronic phases of infection. Expression of EBER2 derived from EBV strain B95-8 resulted in a significant expansion of latently infected B cells in vivo , which was accompanied by a decrease in virus-infected plasma cells. EBV strains typically carry one of two variants of EBER2, which differ primarily by a 5-nucleotide core polymorphism identified initially in the EBV strain M81. Strikingly, mutation of the 5 nucleotides that define this core polymorphism resulted in the loss of the infected B cell expansion and restored plasma cell infection. This work reveals that the B95-8 variant of EBER2 promotes the expansion of the latently infected B cell pool in vivo and may do so in part through inhibition of terminal differentiation. These findings provide new insight into mechanisms by which viral ncRNAs promote in vivo colonization and further and provide further evidence of the inherent tumorigenic risks associated with gammaherpesvirus manipulation of B cell differentiation. IMPORTANCE The oncogenic gammaherpesviruses, including human Epstein-Barr virus (EBV), human Kaposi's sarcoma-associated herpesvirus (KSHV), and murine gammaherpesvirus 68, employ numerous strategies to colonize the host, including expression of noncoding RNAs (ncRNAs). As the first viral ncRNAs ever identified, EBV-encoded RNA 1 and 2 (EBER1 and EBER2) have been investigated extensively for decades; however, their specific in vivo functions remain largely unknown. Work here reveals that an EBV EBER2 variant highly associated with B cell lymphoma promoted a significantly increased expansion of the infected B cell pool in vivo , which coincided with altered B cell differentiation. Mutation of the 5 nucleotides that define this EBER2 variant resulted in the loss of B cell expansion and normal B cell differentiation. These findings provide new insight into the mechanisms by which EBV manipulates B cells in vivo to retain infected cells in the high-risk B cell differentiation pathway where they are poised for tumorigenesis.

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., Competing Interests: Declaration of interests Dr. Garry is co-founder of Zalgen Labs, a biotechnology company that develops countermeasures for emerging viruses., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

While dysregulation of RNA splicing has been recognized as an emerging target for cancer therapy, the functional significance of RNA splicing and individual splicing factors in brain tumors is poorly understood. Here, we identify SON as a master regulator that activates PTBP1-mediated oncogenic splicing while suppressing RBFOX2-mediated non-oncogenic neuronal splicing in glioblastoma multiforme (GBM). SON is overexpressed in GBM patients and SON knockdown causes failure in intron removal from the PTBP1 transcript, resulting in PTBP1 downregulation and inhibition of its downstream oncogenic splicing. Furthermore, SON forms a complex with hnRNP A2B1 and antagonizes RBFOX2, which leads to skipping of RBFOX2-targeted cassette exons, including the PTBP2 neuronal exon. SON knockdown inhibits proliferation and clonogenicity of GBM cells in vitro and significantly suppresses tumor growth in orthotopic xenografts in vivo. Collectively, our study reveals that SON-mediated RNA splicing is a GBM vulnerability, implicating SON as a potential therapeutic target in brain tumors., (© 2021. The Author(s).)

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., (© 2020. The Author(s), under exclusive licence to Springer Nature America, Inc. part of Springer Nature.)

Sex is an important determinant of brain microvessels (MVs) function and susceptibility to cerebrovascular and neurological diseases, but underlying mechanisms are unclear. Using high throughput RNA sequencing analysis, we examined differentially expressed (DE) genes in brain MVs from young, male, and female rats. Bioinformatics analysis of the 23,786 identified genes indicates that 298 (1.2%) genes were DE using False Discovery Rate criteria (FDR; p < 0.05), of which 119 (40%) and 179 (60%) genes were abundantly expressed in male and female MVs, respectively. Nucleic acid binding, enzyme modulator, and transcription factor were the top three DE genes, which were more highly expressed in male than female MVs. Synthesis of glycosylphosphatidylinositol (GPI), biosynthesis of GPI-anchored proteins, steroid and cholesterol synthesis, were the top three significantly enriched canonical pathways in male MVs. In contrast, respiratory chain, ribosome, and 3 ́-UTR-mediated translational regulation were the top three enriched canonical pathways in female MVs. Different gene functions of MVs were validated by proteomic analysis and western blotting. Our novel findings reveal major sex disparities in gene expression and canonical pathways of MVs and these differences provide a foundation to study the underlying mechanisms and consequences of sex-dependent differences in cerebrovascular and other neurological diseases.

The Epstein Barr virus (EBV) contributes to the tumor phenotype through a limited set of primarily non-coding viral RNAs, including 31 mature miRNAs. Here we investigated the impact of EBV miRNAs on remodeling the tumor cell transcriptome. Strikingly, EBV miRNAs displayed exceptionally abundant expression in primary EBV-associated Burkitt's Lymphomas (BLs) and Gastric Carcinomas (GCs). To investigate viral miRNA targeting, we used the high-resolution approach, CLASH in GC and BL cell models. Affinity constant calculations of targeting efficacies for CLASH hits showed that viral miRNAs bind their targets more effectively than their host counterparts, as did Kaposi's sarcoma-associated herpesvirus (KSHV) and murine gammaherpesvirus 68 (MHV68) miRNAs. Using public BL and GC RNA-seq datasets, we found that high EBV miRNA targeting efficacies translates to enhanced reduction of target expression. Pathway analysis of high efficacy EBV miRNA targets showed enrichment for innate and adaptive immune responses. Inhibition of the immune response by EBV miRNAs was functionally validated in vivo through the finding of inverse correlations between EBV miRNAs and immune cell infiltration and T-cell diversity in BL and GC datasets. Together, this study demonstrates that EBV miRNAs are potent effectors of the tumor transcriptome that play a role in suppressing host immune response., Competing Interests: The authors have declared that no competing interests exist.

Upregulation of androgen receptor splice variants (AR-Vs), especially AR-V7, is associated with castration resistance of prostate cancer. At the RNA level, AR-V7 upregulation is generally coupled with increased full-length AR (AR-FL); consequently, AR-V7 and AR-Vs collectively constitute a minority of the AR population. However, Western blotting showed that the relative abundance of AR-V proteins is much higher in many castration-resistant prostate cancers (CRPCs). To address the mechanism underlying this discrepancy, we analyzed RNA-seq data from ~350 CRPC samples and found a positive correlation between all canonical and alternative AR splicing. This indicates that increased alternative splicing is not at the expense of canonical splicing. Instead, androgen deprivation releases AR-FL from repressing the transcription of the AR gene to induce coordinated increase of AR-FL and AR-V mRNAs. At the protein level, however, androgen deprivation induces AR-FL, but not AR-V, degradation. Moreover, AR-V7 is translated much faster than AR-FL. Thus, androgen-deprivation-induced AR-gene transcription and AR-FL protein decay, together with efficient AR-V7 translation, explain the discrepancy between the relative AR-V mRNA and protein abundances in many CRPCs, highlighting the inevitability of AR-V induction after endocrine therapy., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Astrocytes are an early and important target of Zika virus (ZIKV) infection in the developing brain, but the impacts of infection on astrocyte function remain controversial. Given that nonhuman primate (NHP) models of ZIKV infection replicate aspects of neurologic disease seen in human infections, we cultured primary astrocytes from the brain tissue of infant rhesus macaques and then infected the cells with Asian or African lineage ZIKV to identify transcriptional patterns associated with infection in these cells. The African lineage virus appeared to have greater infectivity and promote stronger antiviral signaling, but infection by either strain ultimately produced typical virus response patterns. Both viruses induced hypoxic stress, but the Asian lineage strain additionally had an effect on metabolic and lipid biosynthesis pathways. Together, these findings describe an NHP astrocyte model that may be used to assess transcriptional signatures following ZIKV infection.

These data include secondary analysis of publicly available RNA-seq data from castration-resistant prostate cancer (CRPC) patients as well as RT-qPCR and Western blotting analyses of patient-derived xenograft models and a CRPC cell line. We applied Spearman correlation analysis to assess the relationship between canonical androgen receptor (AR) splicing and alternative AR splicing. We also assessed the ratio of AR splice variants (AR-Vs) to the full-length AR (AR-FL) at the RNA and protein levels by absolute RT-qPCR and Western blotting, respectively. These data are critical for studying the mechanisms underlying upregulated expression of AR-Vs after AR-directed therapies and the importance of AR-Vs to castration-resistant progression of prostate cancer. Data presented here are related to the research article by Ma et al., "Increased transcription and high translation efficiency lead to accumulation of androgen receptor splice variant after androgen deprivation therapy", Cancer Lett. In Press [1]., Competing Interests: This work was supported by the National Institutes of Health [grant numbers R01CA188609, R01AI106676, R01CA243793]. The Richard M Lucas Foundation supported the development of the LuCaP models. The authors declare that they have no known competing financial interests or personal relationships that have or could be perceived to have influenced the work reported in this article., (© 2021 The Authors.)

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.

A major racial disparity in prostate cancer (PCa) is that African American (AA) patients have a higher mortality rate than European American (EA) patients. We filtered the SEER 2009-2011 records and divided them into four groups regarding patient races and cancer grades. On such a partition, we performed a series of statistical analyses to further clarify the aforementioned disparity. Molecular evidence for a primary result of the epidemiological analysis was obtained from gene expression data. The results include: (1) Based on the registry-specific measures, a significant linear regression of total mortality rate (as well as PCa specific mortality rate) on the percentage of (Gleason pattern-based) high-grade cancers (PHG) is demonstrated in EAs (p < 0.01) but not in AAs; (2) PHG and its racial disparity are differentiated across ages and the groups defined by patient outcomes; (3) For patients with cancers in the same grade category, i.e. the high or low grade, the survival stratification between races is not significant in most geographical areas; and (4) The genes differentially expressed between AAs' and EAs' tumors of the same grade category are relatively rare. The perception that prostate tumors are more lethal in AAs than in EAs is reasonable regarding AAs' higher PHG, while high grade alone could not imply aggressiveness. However, this perception is questionable when the comparison is focused on cases within the same grade category. Supporting observations for this conclusion hold a remarkable implication for erasing racial disparity in PCa. That is, "Equal grade, equal outcomes" is not only a verifiable hypothesis but also an achievable public health goal.

Background: Numerous publications suggest an association between herpes virus infection and idiopathic pulmonary fibrosis (IPF). These reports have employed immunohistochemistry, in situ hybridization and/or PCR, which are susceptible to specificity artifacts., Methods: We investigated the possible association between IPF and viral RNA expression using next-generation sequencing, which has the potential to provide a high degree of both sensitivity and specificity. We quantified viral RNA expression for 740 viruses in 28 IPF patient lung biopsy samples and 20 controls. Key RNA-seq results were confirmed using Real-time RT-PCR for select viruses (EBV, HCV, herpesvirus saimiri and HERV-K)., Results: We identified sporadic low-level evidence of viral infections in our lung tissue specimens, but did not find a statistical difference for expression of any virus, including EBV, herpesvirus saimiri and HERV-K, between IPF and control lungs., Conclusions: To the best of our knowledge, this is the first publication that employs RNA-seq to assess whether viral infections are linked to the pathogenesis of IPF. Our results do not address the role of viral infection in acute exacerbations of IPF, however, this analysis patently did not support an association between herpes virus detection and IPF.

We proposed a data cleaning pipeline for single cell (SC) RNA-seq data, where we first screen genes (gene-wise screening) followed by screening cell libraries (library-wise screening). Gene-wise screening is based on the expectation that for a gene with a low technical noise, a gene's count in a library will tend to increase with the increase of library size, which was tested using negative binomial regression of gene count (as dependent variable) against library size (as independent variable). Library-wise screening is based on the expectation that across-library correlations for housekeeping (HK) genes is expected to be higher than the correlations for non-housekeeping (NHK) genes in those libraries with low technical noise. We removed those libraries, whose mean pairwise correlation for HK genes is NOT significantly higher than that for NHK genes. We successfully applied the pipeline to two large SC RNA-seq datasets. The pipeline was also developed into an R package., (Copyright © 2019. Published by Elsevier Inc.)

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.

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.

After an adaptive immune response is mounted, gammaherpesviruses achieve persistence through the utilization of viral noncoding RNAs to craft a suitable host cell environment in an immunologically transparent manner. While gammaherpesvirus long noncoding RNAs (lncRNAs) and microRNAs have been recognized for some time and have been actively investigated, a recent spate of reports have now identified repertoires of the circular RNA (circRNA) class of noncoding RNAs in both the lymphocryptovirus and rhadinovirus genera of gammaherpesviruses. Despite the recent nature of these findings, the detection of circRNAs across viruses and viral gene expression programs, the conservation of some viral circRNAs, and their detection in the clinical setting already raises the spectrum of functional importance in gammaherpesvirus biology and associated malignancies. Here, we provide an overview of currently known gammaherpesvirus circular RNAs and discuss reported physical and contextual properties that may be germane to future functional studies. With the Epstein-Barr virus (EBV) circRNAome being the most extensively studied to date, our discussions will be weighted toward EBV circRNAs while also addressing circRNAs discovered in the rhesus macaque lymphocryptovirus (rLCV), the Kaposi's sarcoma herpesvirus (KSHV), and the murid gammaherpesvirus 68 (MHV68). We hope that this will help set the stage for future investigations into the functions and relevance of this new class of viral noncoding RNAs in infection and disease., (Copyright © 2019 Ungerleider et al.)

EBV and KSHV are etiologic agents of multiple types of lymphomas and carcinomas. The frequency of EBV + or KSHV + malignancies arising in immunocompromised individuals reflects the intricate evolutionary balance established between these viruses and their immunocompetent hosts. However, the specific mechanisms by which these pathogens drive tumorigenesis remain poorly understood. In recent years an enormous array of cellular and viral noncoding RNAs (ncRNAs) have been discovered, and host ncRNAs have been revealed as contributory factors to every single cancer hallmark cellular process. As new evidence emerges that gammaherpesvirus ncRNAs also alter host processes and viral factors dysregulate host ncRNA expression, and as novel viral ncRNAs continue to be discovered, we examine the contribution of small, non-miRNA ncRNAs and long ncRNAs to gammaherpesvirus tumorigenesis., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Kaposi's sarcoma-associated herpesvirus (KSHV) is the causative agent of several human cancers such as Kaposi's sarcoma (KS), which represents the most common AIDS-associated malignancy that lacks effective treatment options. Despite its clear role in AIDS malignancies, the fact that only a small set of KSHV-infected patients will eventually develop these tumors implies that additional co-factors are required for the development of KSHV-related cancers. In the current study, we demonstrate for the first time that KSHV de novo infection or viral latent proteins are able to transactivate human endogenous retrovirus K (HERV-K) through a variety of cellular signaling pathways and transcriptional factors. Moreover, we found that HERV-K transactivation, particularly activation of its encoded oncogenic NP9 protein, plays an important role in KSHV pathogenesis and tumorigenesis in vitro and in vivo. Our data provide innovative insights into the mechanisms of HERV-K transactivation contributing to viral oncogenesis, which may represent a promising target for KS treatment.

Background: Bacterial invasive infection and host immune response is fundamental to the understanding of pathogen pathogenesis and the discovery of effective therapeutic drugs. However, there are very few experimental studies on the signaling cross-talks between bacteria and human host to date., Methods: In this work, taking M. tuberculosis H37Rv (MTB) that is co-evolving with its human host as an example, we propose a general computational framework that exploits the known bacterial pathogen protein interaction networks in STRING database to predict pathogen-host protein interactions and their signaling cross-talks. In this framework, significant interlogs are derived from the known pathogen protein interaction networks to train a predictive l 2 -regularized logistic regression model., Results: The computational results show that the proposed method achieves excellent performance of cross validation as well as low predicted positive rates on the less significant interlogs and non-interlogs, indicating a low risk of false discovery. We further conduct gene ontology (GO) and pathway enrichment analyses of the predicted pathogen-host protein interaction networks, which potentially provides insights into the machinery that M. tuberculosis H37Rv targets human genes and signaling pathways. In addition, we analyse the pathogen-host protein interactions related to drug resistance, inhibition of which potentially provides an alternative solution to M. tuberculosis H37Rv drug resistance., Conclusions: The proposed machine learning framework has been verified effective for predicting bacteria-host protein interactions via known bacterial protein interaction networks. For a vast majority of bacterial pathogens that lacks experimental studies of bacteria-host protein interactions, this framework is supposed to achieve a general-purpose applicability. The predicted protein interaction networks between M. tuberculosis H37Rv and Homo sapiens, provided in the Additional files, promise to gain applications in the two fields: (1) providing an alternative solution to drug resistance; (2) revealing the patterns that M. tuberculosis H37Rv genes target human immune signaling pathways.

Background: Mitochondrial reprogramming has emerged as a hallmark of cancer pathobiology. Although it is believed this reprogramming is essential for cancer cells to thrive, how it supports cancer pathobiology is unclear. We previously generated colonic ρ0 (rho0) cells with reduced mitochondrial energy function and acquired their transcriptional signature. Here, we utilized a bioinformatics approach to identify their changes linked to cancer pathobiology., Methods: Human colon cancer HCT116 cells, control and ρ0, were used for qPCR. Bioinformatics analysis: GeneCards, Kaplan-Meier Survival, GENT, cBioPortal., Results: The colonic ρ0 transcriptome was linked with proliferation, DNA replication, survival, tumor morphology, and cancer. Among differentially expressed transcripts, 281 were regulators or biomarkers of human colon cancer especially those with inflammatory microsatellite instability (MSI). We identified and validated novel transcripts in ρ0 cells with altered expression in human colon cancer. Among them DGK1, HTR7, FLRT3, and ZBTB18 co-occurred with established regulators of human colon cancer pathobiology. Also, increased levels of DGKI, FLRT3, ZBTB18, and YPEL1 as well as decreased levels of HTR7, and CALML6 were linked to substantially poorer patient survival., Conclusion: We identified established and novel regulators in colon cancer pathobiology that are dependent on mitochondrial energy reprogramming and linked to poorer patient survival., Competing Interests: CONFLICTS OF INTERESTS Suzana D. Savkovic wishes to disclose ownership in Pegasus Biosolution, LLC. No other conflicts of interest, financial or otherwise, are declared by the authors.

Elevated overexpression of the lncRNA HOTAIR mediates invasion and metastasis in breast cancer. In an apparent paradox, we observed low expression of HOTAIR in the invasive Claudin-low MDA-MB-231 and Hs578T cells in two-dimensional culture (2D). However, HOTAIR expression exhibited robust induction in laminin-rich extracellular matrix-based three-dimensional organotypic culture (lrECM 3D) over that in 2D culture. Induction of HOTAIR required intact ECM signaling, namely integrin α2 and SRC kinase activity. Moreover, invasive growth was suppressed by HOTAIR-specific siRNA. Induction of HOTAIR in lrECM 3D culture resulted from the activation of a novel isoform of HOTAIR (HOTAIR-N) whose transcription is started from the first intron of the HOXC11 gene. The HOTAIR-N promoter exhibited increased trimethylation of histone H3 lysine 4, a histone marker of active transcription, and binding of BRD4, a reader of transcriptionally active histone markers. Inhibition of BRD4 substantially reduced the expression of HOTAIR in lrECM 3D culture. In summary, our results indicate that HOTAIR expression is activated by BRD4 binding to a novel HOTAIR-N promoter in Claudin-low breast cancer cells that are attached to ECM. Induction of HOTAIR is required for invasive growth of Claudin-low breast cancer cells in lrECM 3D culture., (© 2017 The Authors. Published by FEBS Press and John Wiley & Sons Ltd.)

Lipids are essential for mammalian cells to maintain many physiological functions. Emerging evidence has shown that cancer cells can develop specific alterations in lipid biosynthesis and metabolism to facilitate their survival and various malignant behaviors. To date, the precise role of cellular lipids and lipid metabolism in viral oncogenesis is still largely unclear with only a handful of literature covering this topic to implicate lipid metabolism in oncogenic virus associated pathogenesis. In this review, we focus on the role of lipid biosynthesis and metabolism in the pathogenesis of the Kaposi's sarcoma-associated herpesvirus, a common causative factor for cancers arising in the immunocompromised settings.

High-throughput sequencing of RNA is used to analyze the transcriptomes of viruses and cells, providing information about transcript structure and abundance. A wide array of programs and pipelines has been created to manage and interpret the abundance of data generated from high-throughput RNA sequencing experiments. This protocol details the use of free and open-source programs to align RNA-Seq reads to a reference genome, visualize read coverage and splice junctions, estimate transcript abundance, and evaluate differential expression of transcripts in different conditions. Particular concerns related to EBV and viral transcriptomics are addressed and access to EBV reference files is provided.

Ovarian carcinoma is the fifth-leading cause of cancer death among women in the United States. Major reasons for this persistent mortality include the poor understanding of the underlying biology and a lack of reliable biomarkers. Previous studies have shown that aberrantly expressed MicroRNAs (miRNAs) are involved in carcinogenesis and tumor progression by post-transcriptionally regulating gene expression. However, the interference of miRNAs in tumorigenesis is quite complicated and far from being fully understood. In this work, by an integrative analysis of mRNA expression, miRNA expression and clinical data published by The Cancer Genome Atlas (TCGA), we studied the modularity and dynamicity of miRNA-mRNA interactions and the prognostic implications in high-grade serous ovarian carcinomas. With the top transcriptional correlations (Bonferroni-adjusted p-value<0.01) as inputs, we identified five miRNA-mRNA module pairs (MPs), each of which included one positive-connection (correlation) module and one negative-connection (correlation) module. The number of miRNAs or mRNAs in each module varied from 3 to 7 or from 2 to 873. Among the four major negative-connection modules, three fit well with the widely accepted miRNA-mediated post-transcriptional regulation theory. These modules were enriched with the genes relevant to cell cycle and immune response. Moreover, we proposed two novel algorithms to reveal the group or sample specific dynamic regulations between these two RNA classes. The obtained miRNA-mRNA dynamic network contains 3350 interactions captured across different cancer progression stages or tumor grades. We found that those dynamic interactions tended to concentrate on a few miRNAs (e.g. miRNA-936), and were more likely present on the miRNA-mRNA pairs outside the discovered modules. In addition, we also pinpointed a robust prognostic signature consisting of 56 modular protein-coding genes, whose co-expression patterns were predictive for the survival time of ovarian cancer patients in multiple independent cohorts., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Epigenetic regulation of gene expression is critical to phenotypic maintenance and transition of human breast cancer cells. HOX antisense intergenic RNA (HOTAIR) is a long intergenic non-coding RNA that epigenetically represses gene expression via recruitment of enhancer of zeste homolog 2 (EZH2), a histone methyltransferase. Elevated expression of HOTAIR promotes progression of breast cancer. In the current study we examined the expression and function of HOTAIR in MCF-7-TNR cells, a derivative of the luminal-like breast cancer cell line MCF-7 that acquired resistance to TNF-α-induced cell death. The expression of HOTAIR, markers of the luminal-like and basal-like subtypes, and growth were compared between MCF-7 and MCF-7-TNR cells. These variables were further assessed upon inhibition of HOTAIR, EZH2, p38 MAPK, and SRC kinase in MCF-7-TNR cells. When compared with MCF-7 cells, MCF-7-TNR cells exhibited an increase in the expression of HOTAIR, which correlated with characteristics of a luminal-like to basal-like transition as evidenced by dysregulated gene expression and accelerated growth. MCF-7-TNR cells exhibited reduced suppressive histone H3 lysine27 trimethylation on the HOTAIR promoter. Inhibition of HOTAIR and EZH2 attenuated the luminal-like to basal-like transition in terms of gene expression and growth in MCF-7-TNR cells. Inhibition of p38 and SRC diminished HOTAIR expression and the basal-like phenotype in MCF-7-TNR cells. HOTAIR was robustly expressed in the native basal-like breast cancer cells and inhibition of HOTAIR reduced the basal-like gene expression and growth. Our findings suggest HOTAIR-mediated regulation of gene expression and growth associated with the basal-like phenotype of breast cancer cells., (© 2014 Wiley Periodicals, Inc.)

microRNAs (miRNA) are regulators of cellular pathways and alterations of normal miRNA expression levels have been shown to increase tumorigenesis. miR-24 has been demonstrated as having both tumor suppressive and oncogenic properties depending on cell context. Here, we demonstrate a possible role for pre-miR-24-2 as a tumor suppressor in the MCF-7 breast cancer cell line through the preferential processing of mature miR-24-2* over miR-24. Specifically, we show that the ectopic expression of miR-24-2* in MCF-7 breast cancer cells results in a suppression of cellular survival both in vivo and in vitro. Notably, the overexpression of miR-24-2* results in a dampening of cell survival through the targeted suppression of PKCα. In addition, a similar biological change is observed in vivo where MCF-7 cells overexpressing pre-miR-24-2 have decreased tumorigenicity and tumor incidence. Taken together our data demonstrate that when overexpressed biogenesis of the pre-miR-24-2 favors miR-24-2* in the MCF-7 breast cancer cell line and suggests a tumor suppressive role for miR-24-2* observed through the inhibition of PKCα-mediated cellular survival., (© 2013 Wiley Periodicals, Inc.)

Three-dimensional organotypic culture using reconstituted basement membrane matrix Matrigel (rBM 3-D) is an indispensable tool to characterize morphogenesis of mammary epithelial cells and to elucidate the tumor-modulating actions of extracellular matrix (ECM). microRNAs (miRNAs) are a novel class of oncogenes and tumor suppressors. The majority of our current knowledge of miRNA expression and function in cancer cells is derived from monolayer 2-D culture on plastic substratum, which lacks consideration of the influence of ECM-mediated morphogenesis on miRNAs. In the present study, we compared the expression of miRNAs in rBM 3-D and 2-D cultures of the non-invasive MCF-7 and the invasive MDA-MB231 cells. Our findings revealed a profound difference in miRNA profiles between 2-D and rBM 3-D cultures within each cell type. Moreover, rBM 3-D culture exhibited greater discrimination in miRNA profiles between MCF-7 and MDA-MB231 cells than 2-D culture. The disparate miRNA profiles correlated with distinct mass morphogenesis of MCF-7 and invasive stellate morphogenesis of MDA-MB231 cells in rBM 3-D culture. Supplementation of the tumor promoting type I collagen in rBM 3-D culture substantially altered the miRNA signature of mass morphologenesis of MCF-7 cells in rBM 3-D culture. Overexpression of the differentially expressed miR-200 family member miR429 in MDA-MB231 cells attenuated their invasive stellate morphogenesis in rBM 3-D culture. In summary, we provide the first miRNA signatures of morphogenesis of human breast cancer cells in rBM 3-D culture and warrant further utilization of rBM 3-D culture in investigation of miRNAs in breast cancer.

Composition of extracellular matrix (ECM) is crucial to the establishment and maintenance of epithelial apical-basolateral polarity. Increased ECM rigidity caused by deposition of fibrillar collagen, for example, collagen type I (Col-1), promotes loss of epithelial polarity and tumor progression. microRNAs are small non-coding RNAs that regulate gene expression and fundamental cellular processes. The current study explored a link between microRNAs and Col-1 using organotypic three-dimensional culture in which epithelial cells are embedded within Matrigel, a mimic of basement membrane matrix (Matrigel 3-D). Matrigel 3-D culture of A549, MCF-7, and mK-ras-LE cells (lung and mammary epithelial cell lines) gave rise to acinus, an in vitro equivalent of apical-basolateral polarity that consists of a polarized monolayer of epithelial cells facing a central lumen. Supplementation of Col-1 disrupted acinus. Moreover, Col-1 up-regulated the expression of miR-21, a well-documented oncogenic microRNA, via a post-transcriptional mechanism. Similar post-transcriptional up-regulation of miR-21 correlated with deposition of Col-1 in a murine model of lung fibrogenesis. In summary, our findings link altered ECM composition/rigidity and the expression of oncogenic microRNAs. The current study also suggests a novel post-transcriptional mechanism for regulation of miR-21 expression at maturation from pre-miR-21 to mature miR-21., (Copyright © 2011 Wiley-Liss, Inc.)

The histopathology of idiopathic pulmonary fibrosis (IPF) includes the presence of myofibroblasts within so-called fibroblastic foci, and studies suggest that lung myofibroblasts may be derived from epithelial cells through epithelial--mesenchymal transition (EMT). Transforming growth factor (TGF)-β1 is expressed and/or activated in fibrogenesis, and induces EMT in lung epithelial cells in a dose-dependent manner. A higher occurrence of Epstein-Barr virus (EBV) has been reported in the lung tissue of patients with IPF. EBV expresses latent membrane protein (LMP) 1 during the latent phase of infection, and may play a role in the pathogenesis of pulmonary fibrosis inasmuch as LMP-1 may act as a constitutively active TNF-α receptor. Our data show a remarkable increase in mesenchymal cell markers, along with a concurrent reduction in the expression of epithelial cell markers in lung epithelial cells cotreated with LMP-1, and very low doses of TGF-β1. This effect was mirrored in lung epithelial cells infected with EBV expressing LMP1 and cotreated with TGF-β1. LMP1 pro-EMT signaling was identified, and occurs primarily through the nuclear factor-κB pathway and secondarily through the extracellular signal--regulated kinase (ERK) pathway. Activation of the ERK pathway was shown to be critical for aspects of TGF-β1-induced EMT. LMP1 accentuates the TGF-β1 activation of ERK. Together, these data demonstrate that the presence of EBV-LMP1 in lung epithelial cells synergizes with TGF-β1 to induce EMT. Our in vitro data may help to explain the observation that patients with IPF demonstrating positive staining for LMP1 in lung epithelial cells have a more rapid demise than patients in whom LMP1 is not detected.

Computational prediction of microRNA targets remains a challenging problem. The existing rule-based, data-driven and expression profiling approaches to target prediction are mostly approached from the gene-level. The increasing availability of RNA-seq data provides a new perspective for microRNA target prediction on the isoform-level. We hypothesize that the splicing isoform is the ultimate effector in microRNA targeting and that the proposed isoform-level approach is capable of predicting non-dominant isoform targets as well as their targeting regions that are otherwise invisible to many existing approaches. To test the hypothesis, we used an iterative expectation maximization (EM) algorithm to quantify transcriptomes at the isoform-level. The performance of the EM algorithm in transcriptome quantification was examined in simulation studies using FluxSimulator. We used joint evidence from isoform-level down-regulation and seed enrichment to predict microRNA-155 targets. We validated our computational approach using results from 149 in-house performed in vitro 3'-UTR assays. We also augmented the splicing database using exon-exon junction evidence, and applied the EM algorithm to predict and quantify 1572 cell line specific novel isoforms. Combined with seed enrichment analysis, we predicted 51 novel microRNA-155 isoform targets. Our work is among the first computational studies advocating the isoform-level microRNA target prediction., (© The Author(s) 2011. Published by Oxford University Press.)

Although C-MYC is overexpressed in a number of tumors, the mechanisms governing its expression in normal or tumor cells are not completely understood. Recruitment of the Retinoblastoma protein family members to gene promoters by E2F factors has a dominant negative effect on their activity during the G(0) and G(1) phases of the cell cycle. Despite the presence of an E2F-binding site on the C-MYC promoter, it escapes the repressive effect of E2F-Retinoblastoma complexes through unknown mechanisms during exit from quiescence. We hypothesized that occupancy of E2F elements by factors distinct from E2F might account for this escape. To test this hypothesis, we investigated whether the E2F element in the C-MYC promoter is regulated differently than E2F elements in promoters that are activated at the G(1)-S transition. Employing gel shift analysis, the E2F element from the C-MYC promoter was found to form a unique non-E2F complex, referred to as E2F C-MYC Specific (EMYCS), which is not observed with E2F elements from several other promoters. The DNA contact residues for EMYCS are distinct but overlapping with residues required for binding of E2F proteins. Finally, the approximate estimated molecular weight of the DNA-binding component of EMCYS is 105 kDa. Functional studies indicate that EMYCS has transcriptional transactivation capacity and suggest that it is required to activate the C-MYC promoter during exit from quiescence.