Search

Most cancers are driven by somatic mutations in proto-oncogenes and tumor suppressor genes. Genetic changes in a tumor may accumulate in the tissue self-renewal phase prior to neoplasm. The risk of sporadic mutations increases with age. In this regard, a positive association between patient age and the accumulated mutation burden in tumors exists for many cancer types. However, the reported lines of evidence for such a connection are still limited. TP53 is the most frequently mutated cancer gene. The encoded p53 protein plays crucial roles in DNA repair. Hereby, we speculate that mutant TP53 can disrupt the age-related accumulation patterns of somatic mutations in tumors. We performed linear model analysis on the clinically-annotated genomic data published by TCGA. We found that there was a significant interaction between TP53 genotype (mutant versus wild-type) and patient age at the initial clinical date on somatic mutation burden for five cancers. That is, the regression coefficients of mutation burden on patient age were significant (p < 0.05) for TP53 wild-type tumors but not the mutant counterparts. This disparity was further verified by comparing the group-specific regression coefficients. This finding confirmed our hypothesis and provided unique insights into p53-related tumorigenesis such as the potential temporal order of driver mutations. [ABSTRACT FROM AUTHOR]

Next generation sequencing (NGS) can globally interrogate the genetic composition of biological samples in an unbiased yet sensitive manner. The objective of this study was to utilize the capabilities of NGS to investigate the reported association between glioblastoma multiforme (GBM) and human cytomegalovirus (HCMV). A large-scale comprehensive virome assessment was performed on publicly available sequencing datasets from the Cancer Genome Atlas (TCGA), including RNA-seq datasets from primary GBM (n = 157), recurrent GBM (n = 13), low-grade gliomas (n = 514), recurrent low-grade gliomas (n = 17), and normal brain (n = 5), and whole genome sequencing (WGS) datasets from primary GBM (n = 51), recurrent GBM (n = 10), and normal matched blood samples (n = 20). In addition, RNA-seq datasets from MRI-guided biopsies (n = 92) and glioma stem-like cell cultures (n = 9) were analyzed. Sixty-four DNA-seq datasets from 11 meningiomas and their corresponding blood control samples were also analyzed. Finally, three primary GBM tissue samples were obtained, sequenced using RNA-seq, and analyzed. After in-depth analysis, the most robust virus findings were the detection of papillomavirus (HPV) and hepatitis B reads in the occasional LGG sample (4 samples and 1 sample, respectively). In addition, low numbers of virus reads were detected in several datasets but detailed investigation of these reads suggest that these findings likely represent artifacts or non-pathological infections. For example, all of the sporadic low level HCMV reads were found to map to the immediate early promoter intimating that they likely originated from laboratory expression vector contamination. Despite the detection of low numbers of Epstein-Barr virus reads in some samples, these likely originated from infiltrating B-cells. Finally, human herpesvirus 6 and 7 aligned viral reads were identified in all DNA-seq and a few RNA-seq datasets but detailed analysis demonstrated that these were likely derived from the homologous human telomeric-like repeats. Other low abundance viral reads were detected in some samples but for most viruses, the reads likely represent artifacts or incidental infections. This analysis argues against associations between most known viruses and GBM or mengiomas. Nevertheless, there may be a low percentage association between HPV and/or hepatitis B and LGGs. [ABSTRACT FROM AUTHOR]

The high level of accuracy and sensitivity of next generation sequencing for quantifying genetic material across organismal boundaries gives it tremendous potential for pathogen discovery and diagnosis in human disease. Despite this promise, substantial bacterial contamination is routinely found in existing human-derived RNA-seq datasets that likely arises from environmental sources. This raises the need for stringent sequencing and analysis protocols for studies investigating sequence-based microbial signatures in clinical samples. [ABSTRACT FROM AUTHOR]

Background: The AKT/mammalian target of rapamycin (mTOR) signaling pathway is regulated by 17α-estradiol (E2) signaling and mediates E2-induced proliferation and progesterone receptor (PgR) expression in breast cancer. Methods and results: Here we use deep sequencing analysis of previously published data from The Cancer Genome Atlas to demonstrate that expression of a key component of mTOR signaling, rapamycin-insensitive companion of mTOR (Rictor), positively correlated with an estrogen receptor-α positive (ERα+) breast tumor signature. Through increased microRNA-155 (miR-155) expression in the ERα+ breast cancer cells we demonstrate repression of Rictor enhanced activation of mTOR complex 1 (mTORC1) signaling with both qPCR and western blot. miR-155-mediated mTOR signaling resulted in deregulated ERα signaling both in cultured cells in vitro and in xenografts in vivo in addition to repressed PgR expression and activity. Furthermore we observed that miR-155 enhanced mTORC1 signaling (observed through western blot for increased phosphorylation on mTOR S2448) and induced inhibition of mTORC2 signaling (evident through repressed Rictor and tuberous sclerosis 1 (TSC1) gene expression). mTORC1 induced deregulation of E2 signaling was confirmed using qPCR and the mTORC1-specific inhibitor RAD001. Co-treatment of MCF7 breast cancer cells stably overexpressing miR-155 with RAD001 and E2 restored E2-induced PgR gene expression. RAD001 treatment of SCID/CB17 mice inhibited E2-induced tumorigenesis of the MCF7 miR-155 overexpressing cell line. Finally we demonstrated a strong positive correlation between Rictor and PgR expression and a negative correlation with Raptor expression in Luminal B breast cancer samples, a breast cancer histological subtype known for having an altered ERα-signaling pathway. Conclusions: miRNA mediated alterations in mTOR and ERα signaling establishes a new mechanism for altered estrogen responses independent of growth factor stimulation. [ABSTRACT FROM AUTHOR]

High-throughput RNA sequencing (RNA-seq) has become an instrumental assay for the analysis of multiple aspects of an organism's transcriptome. Further, the analysis of a biological specimen's associated microbiome can also be performed using RNA-seq data and this application is gaining interest in the scientific community. There are many existing bioinformatics tools designed for analysis and visualization of transcriptome data. Despite the availability of an array of next generation sequencing (NGS) analysis tools, the analysis of RNA-seq data sets poses a challenge for many biomedical researchers who are not familiar with command-line tools. Here we present RNA CoMPASS, a comprehensive RNA-seq analysis pipeline for the simultaneous analysis of transcriptomes and metatranscriptomes from diverse biological specimens. RNA CoMPASS leverages existing tools and parallel computing technology to facilitate the analysis of even very large datasets. RNA CoMPASS has a web-based graphical user interface with intrinsic queuing to control a distributed computational pipeline. RNA CoMPASS was evaluated by analyzing RNA-seq data sets from 45 B-cell samples. Twenty-two of these samples were derived from lymphoblastoid cell lines (LCLs) generated by the infection of naïve B-cells with the Epstein Barr virus (EBV), while another 23 samples were derived from Burkitt's lymphomas (BL), some of which arose in part through infection with EBV. Appropriately, RNA CoMPASS identified EBV in all LCLs and in a fraction of the BLs. Cluster analysis of the human transcriptome component of the RNA CoMPASS output clearly separated the BLs (which have a germinal center-like phenotype) from the LCLs (which have a blast-like phenotype) with evidence of activated MYC signaling and lower interferon and NF-kB signaling in the BLs. Together, this analysis illustrates the utility of RNA CoMPASS in the simultaneous analysis of transcriptome and metatranscriptome data. RNA CoMPASS is freely available at http://rnacompass.sourceforge.net/. [ABSTRACT FROM AUTHOR]

Single-nucleotide polymorphisms (SNPs) contribute to the between-individual expression variation of many genes. A regulatory (trait-associated) SNP is usually located near or within a (host) gene, possibly influencing the gene’s transcription or/and post-transcriptional modification. But its targets may also include genes that are physically farther away from it. A heuristic explanation of such multiple-target interferences is that the host gene transfers the SNP genotypic effects to the distant gene(s) by a transcriptional or signaling cascade. These connections between the host genes (regulators) and the distant genes (targets) make the genetic analysis of gene expression traits a promising approach for identifying unknown regulatory relationships. In this study, through a mixed model analysis of multi-source digital expression profiling for 140 human lymphocyte cell lines (LCLs) and the genotypes distributed by the international HapMap project, we identified 45 thousands of potential SNP-induced regulatory relationships among genes (the significance level for the underlying associations between expression traits and SNP genotypes was set at FDR < 0.01). We grouped the identified relationships into four classes (paradigms) according to the two different mechanisms by which the regulatory SNPs affect their cis- and trans- regulated genes, modifying mRNA level or altering transcript splicing patterns. We further organized the relationships in each class into a set of network modules with the cis- regulated genes as hubs. We found that the target genes in a network module were often characterized by significant functional similarity, and the distributions of the target genes in three out of the four networks roughly resemble a power-law, a typical pattern of gene networks obtained from mutation experiments. By two case studies, we also demonstrated that significant biological insights can be inferred from the identified network modules. [ABSTRACT FROM AUTHOR]

Background: The tumor microenvironment is a crucial determinant in tumor progression. Interstitial extracellular matrix (ECM), such as type I collagen (Col-1), is aberrantly enriched in the tumor microenvironment and promotes tumor progression. Long intergenic non-coding RNAs (lincRNA) are a new family of regulatory RNAs that modulate fundamental cellular processes via diverse mechanisms. Findings: We investigated whether the expression of lincRNAs was regulated by the tumor promoting Col-1. In a three-dimensional organotypic culture model using the reconstituted basement membrane ECM Matrigel (rBM 3-D), supplementation of Col-1 disrupted acini, a differentiation feature of well-differentiated lung adenocarcinoma cells, and concurrently induced the expression of a tumor-promoting lincRNA, HOX transcript antisense RNA (HOTAIR). Induction of HOTAIR by Col-1 was diminished by a neutralizing antibody against the Col-1 receptor α2β1 integrin. Col-1 activates the expression of a reporter gene controlled by the human HOTAIR promoter. Moreover the expression of HOTAIR and Col-1 was concurrently up-regulated in human non-small cell lung cancer. Conclusions: Our findings indicate that tumor-promoting Col-1 up-regulates the expression of HOTAIR in NSCLC cells. These initial results warrant further investigation of HOTAIR and other lincRNA genes in lung tumorigenesis. [ABSTRACT FROM AUTHOR]

We have previously shown that disruption of promyelocytic leukemia nuclear bodies (PML NBs) is sufficient to activate the EBV lytic cycle thus making infected cells susceptible to ganciclovir (GCV) mediated killing in vitro. Here we show that co-administration of GCV and arsenic trioxide (ATO), a PML NB disruptor, reduces tumor volume in a xenograft model of nasopharyngeal carcinoma utilizing CNE1 cells. When administered at pharmacologic levels, both GCV and ATO reduced tumor growth while co-treatment with GCV + ATO resulted in a diminution of tumor volume. Treatment with GCV or ATO individually resulted in an increased number of apoptotic cells while co-treatment with GCV + ATO synergistically induced apoptosis. Treatment with ATO or co-treatment with GCV + ATO resulted in expression of EBV lytic proteins. These data suggest that co-treatment with GCV + ATO may provide an effective treatment for nasopharyngeal carcinoma patients. [ABSTRACT FROM AUTHOR]

Epstein-Barr virus (EBV) is associated with roughly 10% of gastric carcinomas worldwide (EBVaGC). Although previous investigations provide a strong link between EBV and gastric carcinomas, these studies were performed using selected EBV gene probes. Using a cohort of gastric carcinoma RNA-seq data sets from The Cancer Genome Atlas (TCGA), we performed a quantitative and global assessment of EBV gene expression in gastric carcinomas and assessed EBV associated cellular pathway alterations. EBV transcripts were detected in 17% of samples but these samples varied significantly in EBV coverage depth. In four samples with the highest EBV coverage (hiEBVaGC – high EBV associated gastric carcinoma), transcripts from the BamHI A region comprised the majority of EBV reads. Expression of LMP2, and to a lesser extent, LMP1 were also observed as was evidence of abortive lytic replication. Analysis of cellular gene expression indicated significant immune cell infiltration and a predominant IFNG response in samples expressing high levels of EBV transcripts relative to samples expressing low or no EBV transcripts. Despite the apparent immune cell infiltration, high levels of the cytotoxic T-cell (CTL) and natural killer (NK) cell inhibitor, IDO1, was observed in the hiEBVaGCs samples suggesting an active tolerance inducing pathway in this subgroup. These results were confirmed in a separate cohort of 21 Vietnamese gastric carcinoma samples using qRT-PCR and on tissue samples using in situ hybridization and immunohistochemistry. Lastly, a panel of tumor suppressors and candidate oncogenes were expressed at lower levels in hiEBVaGC versus EBV-low and EBV-negative gastric cancers suggesting the direct regulation of tumor pathways by EBV. [ABSTRACT FROM AUTHOR]

The count of the nucleotides in a cloned, short genomic sequence has become an important criterion to annotate such a sequence as a miRNA molecule. While the majority of human mature miRNA sequences consist of 22 nucleotides, there exists discrepancy in the characteristic lengths of the miRNA sequences. There is also a lack of systematic studies on such length distribution and on the biological factors that are related to or may affect this length. In this paper, we intend to fill this gap by investigating the sequence structure of human miRNA molecules using statistics tools. We demonstrate that the traditional discrete probability distributions do not model the length distribution of the human mature miRNAs well, and we obtain the statistical distribution model with a decent fit. We observe that the four nucleotide bases in a miRNA sequence are not randomly distributed, implying that possible structural patterns such as dinucleotide (trinucleotide or higher order) may exist. Furthermore, we study the relationships of this length distribution to multiple important factors such as evolutionary conservation, tumorigenesis, the length of precursor loop structures, and the number of predicted targets. The association between the miRNA sequence length and the distributions of target site counts in corresponding predicted genes is also presented. This study results in several novel findings worthy of further investigation that include: (1) rapid evolution introduces variation to the miRNA sequence length distribution; (2) miRNAs with extreme sequence lengths are unlikely to be cancer-related; and (3) the miRNA sequence length is positively correlated to the precursor length and the number of predicted target genes. [ABSTRACT FROM AUTHOR]

Recent studies have shown the contribution of miRNAs to cancer pathogenesis. Prostate cancer is the most commonly diagnosed cancer in men. Unlike other major types of cancer, no single gene has been identified as being mutated in the majority of prostate tumors. This implies that the expression profiling of genes, including the non-coding miRNAs, may substantially vary across individual cases of this cancer. The within-class variability makes it possible to reconstruct or infer disease-specific miRNA-mRNA correlation and regulatory modular networks using high-dimensional microarray data of prostate tumor samples. Furthermore, since miRNAs and tumor suppressor genes are usually tissue specific, miRNA-mRNA modules could potentially differ between primary prostate cancer (PPC) and metastatic prostate cancer (MPC). We herein performed an in silico analysis to explore the miRNA-mRNA correlation network modules in the two tumor subtypes. Our analysis identified 5 miRNA-mRNA module pairs (MPs) for PPC and MPC, respectively. Each MP includes one positiveconnection (correlation) module and one negative-connection (correlation) module. The number of miRNAs or mRNAs (genes) in each module varies from 2 to 8 or from 6 to 622. The modules discovered for PPC are more informative than those for MPC in terms of the implicated biological insights. In particular, one negative-connection module in PPC fits well with the popularly recognized miRNA-mediated post-transcriptional regulation theory. That is, the 3'UTR sequences of the involved mRNAs (∼620) are enriched with the target site motifs of the 7 modular miRNAs, has-miR-106b, -191, -19b, -92a, - 92b, -93, and -141. About 330 GO terms and KEGG pathways, including TGF-beta signaling pathway that maintains tissue homeostasis and plays a crucial role in the suppression of the proliferation of cancer cells, are over-represented (adj.p<0.05) in the modular gene list. These computationally identified modules provide remarkable biological evidence for the interference of miRNAs in the development of prostate cancers and warrant additional follow-up in independent laboratory studies. [ABSTRACT FROM AUTHOR]

In metazoans, miRNAs regulate gene expression primarily through binding to target sites in the 39 UTRs (untranslated regions) of messenger RNAs (mRNAs). Cis-acting variants within, or close to, a gene are crucial in explaining the variability of gene expression measures. Single nucleotide polymorphisms (SNPs) in the 39 UTRs of genes can affect the base-pairing between miRNAs and mRNAs, and hence disrupt existing target sites (in the reference sequence) or create novel target sites, suggesting a possible mechanism for cis regulation of gene expression. Moreover, because the alleles of different SNPs within a DNA sequence of limited length tend to be in strong linkage disequilibrium (LD), we hypothesize the variants of miRNA target sites caused by SNPs potentially function as bridges linking the documented cis-SNP markers to the expression of the associated genes. A large-scale analysis was herein performed to test this hypothesis. By systematically integrating multiple latest information sources, we found 21 significant gene-level SNP-involved miRNA-mediated posttranscriptional regulation modules (SNP-MPRMs) in the form of SNP-miRNA-mRNA triplets in lymphocyte cell lines for the CEU and YRI populations. Among the cognate genes, six including ALG8, DGKE, GNA12, KLF11, LRPAP1, and MMAB are related to multiple genetic diseases such as depressive disorder and Type-II diabetes. Furthermore, we found that ∼35% of the documented transcript intensity-related cis-SNPs (∼950) in a recent publication are identical to, or in significant linkage disequilibrium (LD) (p<0.01) with, one or multiple SNPs located in miRNA target sites. Based on these associations (or identities), 69 significant exon-level SNP-MPRMs and 12 disease genes were further determined for two populations. These results provide concrete in silico evidence for the proposed hypothesis. The discovered modules warrant additional followup in independent laboratory studies. [ABSTRACT FROM AUTHOR]

Abstract: Tumor viruses are a class of pathogens with well established roles in the development of malignant diseases. Numerous bodies of work have highlighted miRNAs (microRNAs) as critical regulators of tumor pathways and it is clear that the dysregulation of cellular miRNA expression can promote tumor formation. Tumor viruses encode their own miRNAs and/or manipulate the expression of cellular miRNAs to modulate their host cell environment, thereby facilitating their respective infection cycles. The modulation of these miRNA responsive pathways, however, often influences certain signal transduction cascades in ways that favor tumorigenesis. In this review, we discuss the roles of virally-encoded and virally-regulated cellular miRNAs in the respective viral life cycles and in virus associated pathogenesis. [Copyright &y& Elsevier]

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. [Display omitted] • Delta peptide is a potent enterotoxin • Injection of delta peptide into murine ileal loops induces fluid accumulation • Delta peptide downregulates the expression of ion transporters and other exchangers • Delta peptide may present a potential therapeutic target Melnik et al. use the murine ligated ileal loop model to demonstrate that delta peptide is a potent enterotoxin. Intestinal fluid accumulation, alteration of villous architecture, and loss of goblet cells follow injection of delta peptide. Delta peptide induces diarrhea via direct cellular damage and regulation of fluid secretion. [ABSTRACT FROM AUTHOR]

Abstract: The Epstein Barr virus (EBV) plays a role in maintenance of the tumor phenotype in a number of human cancers. The EBV latency replication factor, EBNA1, is required for persistence of the EBV episome, is anti-apoptotic, and is universally expressed in all EBV-associated tumors. Here, we show that EBNA1-specific siRNAs can inhibit EBNA1 expression and function. siRNAs were generated against three target sites in the EBNA1 messenger RNA, and two of these were found to inhibit EBNA1 expression from an ectopic EBNA1 expression cassette. EBNA1 siRNAs also inhibit endogenously expressed EBNA1 in EBV-positive epithelial and B-cell lines. Using a mini-EBV replication model, siRNA-mediated inhibition of EBNA1 expression suppressed the episomal maintenance function of EBNA1. Lastly, introduction of an EBNA1 siRNA into an EBV-positive tumor cell line inhibited tumor cell growth/survival. These data suggest that siRNAs against EBNA1 may have therapeutic value in EBV-associated diseases. [Copyright &y& Elsevier]

The gammaherpesviruses, including Epstein-Barr virus (EBV), Kaposi's sarcoma-associated herpesvirus (KSHV), and murine gammaherpesvirus 68 (MHV68, MuHV-4, γHV68), are etiologic agents of a wide range of lymphomas and non-hematological malignancies. These viruses possess large and highly dense dsDNA genomes that feature >80 bidirectionally positioned open reading frames (ORFs). The abundance of overlapping transcripts and extensive splicing throughout these genomes have until now prohibited high throughput-based resolution of transcript structures. Here, we integrate the capabilities of long-read sequencing with the accuracy of short-read platforms to globally resolve MHV68 transcript structures using the transcript resolution through integration of multi-platform data (TRIMD) pipeline. This approach reveals highly complex features, including: (1) pervasive overlapping transcript structures; (2) transcripts containing intra-gene or trans-gene splices that yield chimeric ORFs; (3) antisense and intergenic transcripts containing ORFs; and (4) noncoding transcripts. This work sheds light on the underappreciated complexity of gammaherpesvirus transcription and provides an extensively revised annotation of the MHV68 transcriptome. • Global resolution of transcript isoforms during lytic gammaherpesvirus infection • Genome-wide alternate isoform usage and readthrough transcription • Resolution of 258 transcript structures, including overlapping isoforms • Identification of long noncoding RNAs and unexpected ORFs The highly dense dsDNA genomes of herpesviruses feature an abundance of overlapping transcripts and extensive splicing, greatly hindering the global identification of individual transcripts. O'Grady et al. use integrated multi-platform genomics to globally resolve transcript structures from murine gammaherpesvirus 68, providing an extensively revised genome annotation. [ABSTRACT FROM AUTHOR]

Previous investigations proposed a link between the Epstein-Barr virus (EBV) and lung cancer (LC), but the results are highly controversial largely due to the insufficient sample size and the inherent limitation of the traditional viral screening methods such as PCR. Unlike PCR, current next-generation sequencing (NGS) utilizes an unbiased method for the global assessment of all exogenous agents within a cancer sample with high sensitivity and specificity. In our current study, we aim to resolve this long-standing controversy by utilizing our unbiased NGS-based informatics approaches in conjunction with traditional molecular methods to investigate the role of EBV in a total of 1127 LC. In situ hybridization analysis of 110 LC and 10 normal lung samples detected EBV transcripts in 3 LC samples. Comprehensive virome analyses of RNA sequencing (RNA-seq) data sets from 1017 LC and 110 paired adjacent normal lung specimens revealed EBV transcripts in three lung squamous cell carcinoma and one lung adenocarcinoma samples. In the sample with the highest EBV coverage, transcripts from the BamHI A region accounted for the majority of EBV reads. Expression of EBNA-1, LMP-1 and LMP-2 was observed. A number of viral circular RNA candidates were also detected. Thus, we for the first time revealed a type II latency-like viral transcriptome in the setting of LC in vivo. The high-level expression of viral BamHI A transcripts in LC suggests a functional role of these transcripts, likely as long non-coding RNA. Analyses of cellular gene expression and stained tissue sections indicated an increased immune cell infiltration in the sample expressing high levels of EBV transcripts compared to samples expressing low EBV transcripts. Increased level of immune checkpoint blockade factors was also detected in the sample with higher levels of EBV transcripts, indicating an induced immune tolerance. Lastly, inhibition of immune pathways and activation of oncogenic pathways were detected in the sample with high EBV transcripts compared to the EBV-low LC indicating the direct regulation of cancer pathways by EBV. Taken together, our data support the notion that EBV likely plays a pathological role in a subset of LC. [ABSTRACT FROM AUTHOR]

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. • Single cell RNA sequencing measures gene expression at level of individual cells. • Single cell RNA sequencing data normally contain a lot of technical noise. • A pipeline is proposed to clean the technical noise based on housekeeping genes. • The pipeline is based on rigorous statistical analyses that avoid arbitrary cut off. • A software program is available to implement the pipeline. [ABSTRACT FROM AUTHOR]

Transcription of the immediate early gene BZLF1 is mediated initially through the activation of cellular transcription factors. Reporter-based studies have provided evidence that following this initial activation, the BZLF1 gene product Zta may be involved in an autoactivation loop through binding to its promoter Zp. In contrast, other reports have shown that transfection of a Zta expression vector in latently infected cells does not activate endogenous Zp. Using chromatin immunoprecipitation (ChIP) assays, we show here that Zta binds to endogenous Zp following induction of the lytic cycle by anti-Ig and TGF-beta1 and that binding occurs early enough to play a role in the activation of Zp. We have also generated a dominant-negative Zta and shown that it inhibits activation of endogenous Zp. These data support a two-step model for Zp activation during reactivation involving initial activation by cellular factors followed by an autoactivation step. [Copyright &y& Elsevier]

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. [ABSTRACT FROM AUTHOR]

Recent studies have identified circular RNAs (circRNAs) expressed from the Epstein-Barr virus (EBV) and Kaposi's sarcoma herpesvirus (KSHV) human DNA tumor viruses. To gain initial insights into the potential relevance of EBV circRNAs in virus biology and disease, we assessed the circRNAome of the interspecies homologue rhesus macaque lymphocryptovirus (rLCV) in a naturally occurring lymphoma from a simian immunodeficiency virus (SIV)-infected rhesus macaque. This analysis revealed rLCV orthologues of the latency-associated EBV circular RNAs circRPMS1_E4_E3a and circEBNA_U. Also identified in two samples displaying unusually high lytic gene expression was a novel rLCV circRNA that contains both conserved and rLCV-specific RPMS1 exons and whose backsplice junctions flank an rLCV lytic origin of replication (OriLyt). Analysis of a lytic infection model for the murid herpesvirus 68 (MHV68) rhadinovirus identified a cluster of circRNAs near an MHV68 lytic origin of replication, with the most abundant of these, circM11_ORF69, spanning the OriLyt. Lastly, analysis of KSHV latency and reactivation models revealed the latency associated circRNA originating from the vIRF4 gene as the predominant viral circRNA. Together, the results of this study broaden our appreciation for circRNA repertoires in the Lymphocryptovirus and Rhadinovirus genera of gammaherpesviruses and provide evolutionary support for viral circRNA functions in latency and viral replication. [ABSTRACT FROM AUTHOR]

Nef-specific CD8+ T lymphocytes (CD8TL) are linked to extraordinary control of primate lentiviral replication, but the mechanisms underlying their efficacy remain largely unknown. The immunodominant, Mamu-B*017:01+-restricted Nef195-203MW9 epitope in SIVmac239 partially overlaps a sorting motif important for interactions with host AP-2 proteins and, hence, downmodulation of several host proteins, including Tetherin (CD317/BST-2), CD28, CD4, SERINC3, and SERINC5. We reasoned that CD8TL-driven evolution in this epitope might compromise Nef's ability to modulate these important molecules. Here, we used deep sequencing of SIV from nine B*017:01+ macaques throughout infection with SIVmac239 to characterize the patterns of viral escape in this epitope and then assayed the impacts of these variants on Nef-mediated modulation of multiple host molecules. Acute variation in multiple Nef195-203MW9 residues significantly compromised Nef's ability to downregulate surface Tetherin, CD4, and CD28 and reduced its ability to prevent SERINC5-mediated reduction in viral infectivity but did not impact downregulation of CD3 or major histocompatibility complex class I, suggesting the selective disruption of immunomodulatory pathways involving Nef AP-2 interactions. Together, our data illuminate a pattern of viral escape dictated by a selective balance to maintain AP-2-mediated downregulation while evading epitope-specific CD8TL responses. These data could shed light on mechanisms of both CD8TL-driven viral control generally and on Mamu-B*017:01-mediated viral control specifically. IMPORTANCE A rare subset of humans infected with HIV-1 and macaques infected with SIV can control the virus without aid of antiviral medications. A common feature of these individuals is the ability to mount unusually effective CD8 T lymphocyte responses against the virus. One of the most formidable aspects of HIV is its ability to evolve to evade immune responses, particularly CD8 T lymphocytes. We show that macaques that target a specific peptide in the SIV Nef protein are capable of better control of the virus and that, as the virus evolves to escape this response, it does so at a cost to specific functions performed by the Nef protein. Our results help show how the virus can be controlled by an immune response, which could help in designing effective vaccines. [ABSTRACT FROM AUTHOR]

The relationship between Epstein Barr Virus (EBV) and miR-155 is well established. EBV infection induces miR-155 expression, which is expressed at higher levels in EBV latency type III cells compared to EBV latency type I cells. However, the mechanism by which EBV latency genes activate miR-155 expression is still unclear. Here we present data showing that DNA methylation regulates miR-155 expression. We also provide evidence that the AP1 signaling pathway is involved in EBV-mediated miR-155 activation, and that Bay11 influences signaling of the miR-155 promoter AP1 element. Lastly, we show that LMP2A, LMP1 and EBNAs cannot activate miR-155 expression alone, indicating that the regulation of miR-155 by EBV is dependent on more than one EBV gene or cell signaling pathway. We conclude that the regulation of miR-155 in EBV-positive cells occurs through multiple cell signaling processes involving EBV-mediated chromatin remodeling, cell signaling regulation and transcription factor activation. [ABSTRACT FROM AUTHOR]

In the oral epithelium, peripheral stores of Epstein-Barr virus (EBV) are transmitted from infiltrating B cells to epithelial cells. Once the virus is transmitted to epithelial cells, the highly permissive nature of this cell type for lytic replication allows virus amplification and exchange to other hosts. Since the initial transfer of EBV from B cells to epithelial cells requires transitioning of the B-cell to a state that induces virus reactivation, we hypothesized that there might be epithelium-specific signals that allow the infiltrating B cells to sense the appropriate environment to initiate reactivation and begin this exchange process. We previously found that the epithelium-specific miR-200 family of microRNAs promotes EBV lytic replication. Here we show that there are high levels of miR-200 family members in oral and tonsillar epithelia and in saliva. Analysis of cultured oral epithelial cells (OKF6) showed that they actively secrete membrane vesicles (exosomes) that are enriched with miR-200 family members. Coculturing of EBV-positive B cells with OKF6 cells induced viral reactivation. Further, treatment of EBV-positive B cells with OKF6 cell-derived membrane vesicles promoted reactivation. Using a cell system that does not naturally express miR-200 family members, we found that enforced expression of a miR-200 family member produced membrane vesicles that were able to induce the lytic cascade in EBV-positive B cells. We propose that membrane vesicles secreted by oral and tonsillar epithelial cells may serve as a tissue-specific environmental cue that initiates reactivation in B cells, promoting the transfer of virus from peripheral B-cell stores to the oral epithelium to facilitate virus amplification and exchange to other hosts. [ABSTRACT FROM AUTHOR]

Summary Dysregulation of MLL complex-mediated histone methylation plays a pivotal role in gene expression associated with diseases, but little is known about cellular factors modulating MLL complex activity. Here, we report that SON, previously known as an RNA splicing factor, controls MLL complex-mediated transcriptional initiation. SON binds to DNA near transcription start sites, interacts with menin, and inhibits MLL complex assembly, resulting in decreased H3K4me3 and transcriptional repression. Importantly, alternatively spliced short isoforms of SON are markedly upregulated in acute myeloid leukemia. The short isoforms compete with full-length SON for chromatin occupancy but lack the menin-binding ability, thereby antagonizing full-length SON function in transcriptional repression while not impairing full-length SON-mediated RNA splicing. Furthermore, overexpression of a short isoform of SON enhances replating potential of hematopoietic progenitors. Our findings define SON as a fine-tuner of the MLL-menin interaction and reveal short SON overexpression as a marker indicating aberrant transcriptional initiation in leukemia. [ABSTRACT FROM AUTHOR]

The Deepwater Horizon oil spill (BP oil spill) in the Gulf of Mexico was a unique disaster event, where a huge amount of oil spilled from the sea bed and a large volume of dispersants were applied to clean the spill. The operation lasted for almost 3 months and involved > 50,000 workers. The potential health hazards to these workers may be significant as previous research suggested an association of persistent respiratory symptoms with exposure to oil and oil dispersants. To reveal the potential effects of oil and oil dispersants on the respiratory system at the molecular level, we evaluated the transcriptomic profile of human airway epithelial cells grown under treatment of crude oil, the dispersants Corexit 9500 and Corexit 9527, and oil-dispersant mixtures. We identified a very strong effect of Corexit 9500 treatment, with 84 genes (response genes) differentially expressed in treatment vs. control samples. We discovered an interactive effect of oil-dispersant mixtures; while no response gene was found for Corexit 9527 treatment alone, cells treated with Corexit 9527 + oil mixture showed an increased number of response genes (46 response genes), suggesting a synergic effect of 9527 with oil on airway epithelial cells. Through GO (gene ontology) functional term and pathway-based analysis, we identified upregulation of gene sets involved in angiogenesis and immune responses and downregulation of gene sets involved in cell junctions and steroid synthesis as the prevailing transcriptomic signatures in the cells treated with Corexit 9500, oil, or Corexit 9500 + oil mixture. Interestingly, these key molecular signatures coincide with important pathological features observed in common lung diseases, such as asthma, cystic fibrosis and chronic obstructive pulmonary disease. Our study provides mechanistic insights into the detrimental effects of oil and oil dispersants to the respiratory system and suggests significant health impacts of the recent BP oil spill to those people involved in the cleaning operation. [ABSTRACT FROM AUTHOR]

The Epstein-Barr virus (EBV) BNLF2a gene product provides immune evasion properties to infected cells through inhibition of transporter associated with antigen processing (TAP)-mediated transport of antigen peptides. Although BNLF2a is considered to be a lytic gene, we demonstrate that it is expressed in nearly half of the EBV-associated gastric carcinomas analyzed. Further, we show that BNLF2a expression is dissociated from lytic gene expression. BNLF2a is therefore expressed in this latency setting, potentially helping protect the infected tumor cells from immunosurveillance. [ABSTRACT FROM AUTHOR]

Constitutively active androgen receptor splice variants (AR-V) lacking the ligand-binding domain have been implicated in the pathogenesis of castration-resistant prostate cancer and in mediating resistance to newer drugs that target the androgen axis. AR-V regulates expression of both canonical AR targets and a unique set of cancer-specific targets that are enriched for cell-cycle functions. However, little is known about how AR-V controls gene expression. Here, we report that two major AR-Vs, termed AR-V7 and ARv567es, not only homodimerize and heterodimerize with each other but also heterodimerize with full-length androgen receptor (AR-FL) in an androgen-independent manner. We found that heterodimerization of AR-V and AR-FL was mediated by N- and C-terminal interactions and by the DNA-binding domain of each molecule, whereas AR-V homodimerization was mediated only by DNA-binding domain interactions. Notably, AR-V dimerization was required to transactivate target genes and to confer castrationresistant cell growth. Our results clarify the mechanism by which AR-Vs mediate gene regulation and provide a pivotal pathway for rational drug design to disrupt AR-V signaling as a rational strategy for the effective treatment of advanced prostate cancer. [ABSTRACT FROM AUTHOR]

We have previously shown that the Epstein-Barr virus (EBV) likely encodes hundreds of viral long noncoding RNAs (vlncRNAs) that are expressed during reactivation. Here we show that the EBV latency origin of replication (oriP) is transcribed bi-directionally during reactivation and that both leftward (oriPtLs) and rightward (oriPtRs) transcripts are largely localized in the nucleus. While the oriPtLs are most likely noncoding, at least some of the oriPtRs contain the BCRF1/vIL10 open reading frame. Nonetheless, oriPtR transcripts with long 5= untranslated regions may partially serve noncoding functions. Both oriPtL and oriPtR transcripts are expressed with late kinetics, and their expression is inhibited by phosphonoacetic acid.RNAsequencing (RNA-seq) analysis showed that oriPtLs and oriPtRs exhibited extensive "hyperediting" at their Family of Repeat (FR) regions.RNAsecondary structure prediction revealed that the FR region of both oriPtLs and oriPtRs may form large evolutionarily conserved and thermodynamically stable hairpins. The doublestranded RNA-binding protein and RNA-editing enzymeADARwas found to bind to oriPtLs, likely facilitating editing of the FR hairpin. Further, the multifunctional paraspeckle protein, NONO, was found to bind to oriPt transcripts, suggesting that oriPts interact with the paraspeckle-based innate antiviral immune pathway. Knockdown and ectopic expression of oriPtLs showed that it contributes to global viral lytic gene expression and viralDNAreplication. Together, these results show that these new vlncRNAs interact with cellular innate immune pathways and that they help facilitate progression of the viral lytic cascade. [ABSTRACT FROM AUTHOR]

Using high-throughput RNA sequencing data from 50 common lymphoma cell culture models from the Cancer Cell Line Encyclopedia project, we performed an unbiased global interrogation for the presence of a panel of 740 viruses and strains known to infect human and other mammalian cells. This led to the findings of previously identified infections by Epstein-Barr virus (EBV), Kaposi's sarcoma herpesvirus (KSHV), and human T-lymphotropic virus type 1 (HTLV-1). In addition, we also found a previously unreported infection of one cell line (DEL) with a murine leukemia virus. High expression of murine leukemia virus (MuLV) transcripts was observed in DEL cells, and we identified four transcriptionally active integration sites, one being in the TNFRSF6B gene. We also found low levels of MuLV reads in a number of other cell lines and provided evidence suggesting crosscontamination during sequencing. Analysis of HTLV-1 integrations in two cell lines, HuT 102 and MJ, identified 14 and 66 transcriptionally active integration sites with potentially activating integrations in immune regulatory genes, including interleukin- 15 (IL-15), IL-6ST, STAT5B, HIVEP1, and IL-9R. Although KSHV and EBV do not typically integrate into the genome, we investigated a previously identified integration of EBV into the BACH2 locus in Raji cells. This analysis identified a BACH2 disruption mechanism involving splice donor sequestration. Through viral gene expression analysis, we detected expression of stable intronic RNAs from the EBV BamHIWrepeats that may be part of long transcripts spanning the repeat region. We also observed transcripts at the EBV vIL-10 locus exclusively in the Hodgkin's lymphoma cell line, Hs 611.T, the expression of which were uncoupled from other lytic genes. Assessment of the KSHV viral transcriptome in BCP-1 cells showed expression of the viral immune regulators, K2/vIL-6, K4/vIL-8-like vCCL1, and K5/E2-ubiquitin ligase 1 that was significantly higher than expression of the latency-associated nuclear antigen. Together, this investigation sheds light into the virus composition across these lymphoma model systems and provides insights into common viral mechanistic principles. [ABSTRACT FROM AUTHOR]

Epstein-Barr virus (EBV) reactivation involves the ordered induction of approximately 90 viral genes that participate in the generation of infectious virions. Using strand-specific RNA-seq to assess the EBV transcriptome during reactivation, we found extensive bidirectional transcription extending across nearly the entire genome. In contrast, only 4% of the EBV genome is currently bidirectionally annotated. Most of the newly identified transcribed regions show little evidence of coding potential, supporting noncoding roles for most of these RNAs. Based on previous cellular long noncoding RNA size calculations, we estimate that there are likely hundreds more EBV genes expressed during reactivation than was previously known. Limited 5= and 3= rapid amplification of cDNA ends (RACE) experiments and findings of novel splicing events by RNA-seq suggest that the complexity of the viral genome during reactivation may be even greater. Further analysis of antisense transcripts at some of the EBV latency gene loci showed that they are "late" genes, they are nuclear, and they tend to localize in areas of the nucleus where others find newly synthesized viral genomes. This raises the possibility that these transcripts perform functions such as new genome processing, stabilization, organization, etc. The finding of a significantly more complex EBV transcriptome during reactivation changes our view of the viral production process from one that is facilitated and regulated almost entirely by previously identified viral proteins to a process that also involves the contribution of a wide array of virus encoded noncoding RNAs. [ABSTRACT FROM AUTHOR]

Abstract: Three-dimensional organotypic culture using reconstituted basement membrane matrix (rBM 3-D) is an invaluable tool to characterize morphogenesis of epithelial cells and to elucidate the tumor-modulating actions of extracellular matrix. microRNAs (miRNA) are a novel class of tumor modulating genes. A substantial amount of investigation of miRNAs in cancer is carried out using monolayer 2-D culture on plastic substratum, which lacks a consideration of the matrix-mediated regulation of miRNAs. In the current study we compared the expression of miRNAs in rBM 3-D and 2-D cultures of two lung adenocarcinoma cell lines. Our findings revealed a profound difference in miRNA profiles between 2-D and rBM 3-D cultures of lung adenocarcinoma cells. The rBM 3-D culture-specific miRNA profile was highlighted with higher expression of the tumor suppressive miRNAs (i.e., miR-200 family) and lower expression of the oncogenic miRNAs (i.e., miR-17–92 cluster and miR-21) than that of 2-D culture. Moreover, the expression pattern of miR-17, miR-21, and miR-200a in rBM 3-D culture correlated with the expression of their targets and acinar morphogenesis, a differentiation behavior of lung epithelial cells in rBM 3-D culture. Over-expression of miR-21 suppressed its target PTEN and disrupted acinar morphogenesis. In summary, we provide the first miRNA profile of lung adenocarcinoma cells in rBM 3-D culture with respect to acinar morphogenesis. These results indicate that rBM 3-D culture is essential to a comprehensive understanding of the miRNA biology in lung epithelial cells pertinent to lung adenocarcinoma. [Copyright &y& Elsevier]

Abstract: Promyelocytic leukemia protein nuclear bodies (PML NBs) have been implicated in host immune response to viral infection. PML NBs are targeted for degradation during reactivation of herpes viruses, suggesting that disruption of PML NB function supports this aspect of the viral life cycle. The Epstein–Barr virus (EBV) Latent Membrane Protein 1 (LMP1) has been shown to suppress EBV reactivation. Our finding that LMP1 induces PML NB immunofluorescence intensity led to the hypothesis that LMP1 may modulate PML NBs as a means of maintaining EBV latency. Increased PML protein and morphometric changes in PML NBs were observed in EBV infected alveolar epithelial cells and nasopharyngeal carcinoma cells. Treatment with low dose arsenic trioxide disrupted PML NBs, induced expression of EBV lytic proteins, and conferred ganciclovir susceptibility. This study introduces an effective modality to induce susceptibility to ganciclovir in epithelial cells with implications for the treatment of EBV associated pathologies. [Copyright &y& Elsevier]

Abstract: The Epstein–Barr virus (EBV) is associated with lymphoid and epithelial cancers. Initial EBV infection alters lymphocyte gene expression, inducing cellular proliferation and differentiation as the virus transitions through consecutive latency transcription programs. Cellular microRNAs (miRNAs) are important regulators of signaling pathways and are implicated in carcinogenesis. The extent to which EBV exploits cellular miRNAs is unknown. Using micro-array analysis and quantitative PCR, we demonstrate differential expression of cellular miRNAs in type III versus type I EBV latency including elevated expression of miR-21, miR-23a, miR-24, miR-27a, miR-34a, miR-146a and b, and miR-155. In contrast, miR-28 expression was found to be lower in type III latency. The EBV-mediated regulation of cellular miRNAs may contribute to EBV signaling and associated cancers. [Copyright &y& Elsevier]

The cellular microRNA miR-155 has been shown to be involved in lymphocyte activation and is expressed in Epstein-Barr virus (EBV)-infected cells displaying type III latency gene expression but not type I latency gene expression. We show here that the elevated levels of miR-155 in type III latency cells is due to EBV gene expression and not epigenetic differences in cell lines tested, and we show that expression in EBV-infected cells requires a conserved AP-1 element in the miR-155 promoter. Gene expression analysis was carried out in a type I latency cell line transduced with an miR-155-expressing retrovirus. This analysis identified both miR-155-suppressed and -induced cellular mRNAs and suggested that in addition to direct targeting of 3' untranslated regions (UTRs), miR-155 alters gene expression in part through the alteration of signal transduction pathways. 3' UTR reporter analysis of predicted miR-155 target genes identified the transcriptional regulatory genes encoding BACH1, ZIC3, HIVEP2, CEBPB, ZNF652, ARID2, and SMAD5 as miR-155 targets. Western blot analysis of the most highly suppressed of these, BACH1, showed lower expression in cells transduced with a miR-155 retrovirus. Inspection of the promoters from genes regulated in EBV-infected cells and in cells infected with an miR-155 retrovirus identified potential binding sequences for BACH1 and ZIC3. Together, these experiments suggest that the induction of miR-155 by EBV contributes to EBV-mediated signaling in part through the modulation of transcriptional regulatory factors. [ABSTRACT FROM AUTHOR]

Many cell lines commonly used for biological studies have been found to harbor exogenous agents such as the human tumor viruses Epstein-Barr virus (EBV) and human papillomavirus. Nevertheless, broad-based, unbiased approaches to globally assess the presence of ectopic organisms within cell model systems have not previously been available. We reasoned that high-throughput sequencing should provide unparalleled insights into the microbiomes of tissue culture cell systems. Here we have used our RNA-seq analysis pipeline, PARSES (Pipeline for Analysis of RNA-Seq Exogenous Sequences), to investigate the presence of ectopic organisms within two EBV-positive B-cell lines commonly used by EBV researchers. Sequencing data sets from both the Akata and JY B-cell lines were found to contain reads for EBV, and the JY data set was found to also contain reads from the murine leukemia virus (MuLV). Further investigation revealed that MuLV transcription in JY cells is highly active. We also identified a number of MuLV alternative splicing events, and we uncovered evidence of APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G)-dependent DNA editing. Finally, reverse transcription-PCR analysis showed the presence of MuLV in three other human B-cell lines (DG75, Ramos, and P3HR1 Cl.13) commonly used by investigators in the Epstein-Barr virus field. We believe that a thorough examination of tissue culture microbiomes using RNA-seq/PARSES-like approaches is critical for the appropriate utilization of these systems in biological studies. [ABSTRACT FROM AUTHOR]

The miR-200 microRNA family is important for maintaining the epithelial phenotype, partially through suppressing ZEB1 and ZEB2. Since ZEB1 inhibits Epstein-Barr virus (EBV) reactivation, we hypothesized that expression of miR-200 family members in epithelial cells may partly account for higher levels of EBV reactivation in this tissue (relative to nonplasma B cells). Here we show that, whereas miR-200 family members are expressed in epithelial cells, their expression is low in latently infected B cells. Furthermore, the miR-200 family member miR-429 shows elevated expression in plasma cell lines and is induced by B-cell-receptor activation in Akata cells. Lastly, expression of miR-429 can break latency. [ABSTRACT FROM AUTHOR]

MicroRNA miR-155 is expressed at elevated levels in human cancers including cancers of the lung, breast, colon, and a subset of lymphoid malignancies. In B cells, miR-155 is induced by the oncogenic latency gene expression program of the human herpesvirus Epstein-Barr virus (EBV). Two other oncogenic herpesviruses, Kaposi's sarcoma-associated herpesvirus and Marek's disease virus, encode functional homologues of miR-155, suggesting a role for this microRNA in the biology and pathogenesis of these viruses. Bone morphogenetic protein (BMP) signaling is involved in an array of cellular processes, including differentiation, growth inhibition, and senescence, through context-dependent interactions with multiple signaling pathways. Alteration of this pathway contributes to a number of disease states including cancer. Here, we show that miR-155 targets the 3' untranslated region of multiple components of the BMP signaling cascade, including SMAD1, SMAD5, HIVEP2, CEBPB, RUNX2, and MYO10. Targeting of these mediators results in the inhibition of BMP2-, BMP6-, and BMP7-induced ID3 expression as well as BMP-mediated EBV reactivation in the EBV-positive B-cell line, Mutu I. Further, miR-155 inhibits SMAD1 and SMAD5 expression in the lung epithelial cell line A549, it inhibits BMP-mediated induction of the cyclin-dependent kinase inhibitor p21, and it reverses BMP-mediated cell growth inhibition. These results suggest a role for miR-155 in controlling BMP-mediated cellular processes, in regulating BMPinduced EBV reactivation, and in the inhibition of antitumor effects of BMP signaling in normal and virus-infected cells. [ABSTRACT FROM AUTHOR]

The Epstein-Barr virus (EBV)-encoded latent membrane protein 1 (LMP1) is a functional homologue of the tumor necrosis factor receptor family and contributes substantially to the oncogenic potential of EBV through activation of nuclear factor κB (NF-κB). MicroRNAs (miRNAs) are a class of small RNA molecules that are involved in the regulation of cellular processes such as growth, development, and apoptosis and have recently been linked to cancer phenotypes. Through miRNA microarray analysis, we demonstrate that LMP1 dysregulates the expression of several cellular miRNAs, including the most highly regulated of these, miR-146a. Quantitative reverse transcription-PCR analysis confirmed induced expression of miR-146a by LMP1. Analysis of miR-146a expression in EBV latency type III and type I cell lines revealed substantial expression of miR-146a in type III (which express LMP1) but not in type I cell lines. Reporter studies demonstrated that LMP1 induces miR-146a predominantly through two NF-κB binding sites in the miR-146a promoter and identified a role for an Oct-1 site in conferring basal and induced expression. Array analysis of cellular mRNAs expressed in Akata cells transduced with an miR-146a-expressing retrovirus identified genes that are directly or indirectly regulated by miR-146a, including a group of interferon-responsive genes that are inhibited by miR-146a. Since miR-146a is known to be induced by agents that activate the interferon response pathway (including LMP1), these results suggest that miR-146a functions in a negative feedback loop to modulate the intensity and/or duration of the interferon response. [ABSTRACT FROM AUTHOR]