Search

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]

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

Obesity, an immense epidemic affecting approximately half a billion adults, has doubled in prevalence in the last several decades. Epidemiological data support that obesity, due to intake of a high-fat, western diet, increases the risk of colon cancer; however, the mechanisms underlying this risk remain unclear. Here, utilizing next generation RNA sequencing, we aimed to determine the high-fat diet (HFD) mediated expression profile in mouse colon and the azoxymethane/ dextran sulfate sodium model of colon cancer. Mice on HFD had significantly higher colonic inflammation, tumor burden, and a number of differentially expressed transcripts compared to mice on regular diet (RD). We identified 721 transcripts differentially expressed in mouse HFD colon that were in a shared pattern with colonic tumors (RD and HFD). Importantly, in mouse colon, HFD stimulated an expression signature strikingly similar to human colon cancer, especially those with inflammatory microsatellite instability. Furthermore, pathway analysis of these transcripts demonstrated their association with active inflammation and colon cancer signaling, with leptin and Wnt as the top two transcripts elevated in mouse HFD colon shared with tumors. Moreover, in mouse colon, HFD-stimulated tumorigenic Wnt pathway activation was further validated by upregulation of β-catenin transcriptional targets. Finally, in human colon cancer, upregulation of leptin pathway members was shown with a large network of dysregulated transcripts being linked with worse overall survival. [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]

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]

Background Monocytes are an important cell type in chronic periodontitis (CP) by interacting with oral bacteria and mediating host immune response. The aim of this study was to reveal new functional genes and pathways for CP at monocyte transcriptomic level. Methods We performed an RNA-sequencing (RNA-seq) study of peripheral blood monocytes (PBMs) in 5 non-smoking moderate to severe CP (case) individuals vs. 5 controls. We took advantage of a microarray study of periodontitis to support our findings. We also performed pathway-based analysis on the identified differentially expressed (DEx) transcripts/isoforms using DAVID (Database for Annotation, Visualization and Integrated Discovery). Results Through differential expression analyses at both whole gene (or whole non-coding RNA) and isoform levels, we identified 380 DEx transcripts and 5955 DEx isoforms with a PPEE (posterior probability of equal expression) of < 0.05. Pervasive up-regulation of transcripts at isoform level in CP vs. control individuals was observed, suggesting a more functionally active monocyte transcriptome for CP. By comparing with the microarray dataset, we identified several CP-associated novel genes (e.g., FACR and CUX1) that have functions to interact with invading microorganisms or enhance TNF production on lipopolysaccharide stimulation. DAVID analysis of both the RNA-seq and the microarray datasets leads to converging evidence supporting “endocytosis”, “cytokine production” and “apoptosis” as significant biological processes in CP. Conclusions As the first RNA-seq study of PBMs for CP, this study provided novel findings at both gene (e.g., FCAR and CUX1) and biological process level. The findings will contribute to better understanding of CP disease mechanisms. [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]