Quantitative RNAseq analysis of Ugandan KS tumors reveals KSHV gene expression dominated by transcription from the LTd downstream latency promoter
Ontology highlight
ABSTRACT:
PROVIDER: PRJNA477513 | ENA |
REPOSITORIES: ENA
Dataset's files
| Action | DRS | |||
|---|---|---|---|---|
| SRR7412070.fastq.gz | Fastqsanger.gz | |||
| SRR7412080.fastq.gz | Fastqsanger.gz | |||
| SRR7412090.fastq.gz | Fastqsanger.gz | |||
| SRR7412100.fastq.gz | Fastqsanger.gz | |||
| SRR7412071.fastq.gz | Fastqsanger.gz |
Items per page: 5 1 - 5 of 41 |
Similar Datasets
Project description:Purpose: We used RNA deep sequencing to characterize KSHV expression in a large collection of KS biopsies (n=41) from HIV-infected Ugandans. Using a novel approach to quantitate expression in complex genomes like KSHV, we found that RNA from a single KSHV promoter within the latency region constituted the majority of KSHV transcripts in the KS tumors. Alternate RNA processing produced different spliced and un-spliced transcripts with different coding potentials. Differential expression of other KSHV genes was detected which segregated the tumors into three different types depending on their expression of lytic or latency genes. Quantitative analysis of KSHV expression in KS tumors provides an important basis for future studies on the role of KSHV in the development of KS Methods: Total nucleic acids were extracted using RLT buffer (Qiagen) and then RNA was isolated using the RNeasy mini kit with DNAse treatment step. Total RNA integrity was checked using an Agilent 2200 TapeStation (Agilent Technologies, Inc., Santa Clara, CA) and quantified using a Trinean DropSense96 spectrophotometer (Caliper Life Sciences, Hopkinton, MA). Unstranded RNA-seq libraries were prepared from 300 ng of total RNA using the TruSeq RNA Sample Prep Kit v2 (Illumina, Inc., San Diego, CA, USA). Four KS tumors libraries were prepared using the TruSeq Stranded mRNA Library Kit (Illumina) from 100 ng of total RNA. Library size distributions were validated using an Agilent 2200 TapeStation (Agilent Technologies, Santa Clara, CA, USA). Additional library QC, pooling of indexed libraries, and cluster optimization was performed using Life Technologies’ Invitrogen Qubit® 2.0 Fluorometer (Life Technologies-Invitrogen, Carlsbad, CA, USA). The unstranded RNA-seq libraries were pooled (5-plex) and the stranded libraries were pooled (4-plex) and each pool was clustered onto a flow cell lane. Sequencing was performed using an Illumina HiSeq 2500 in “High Output” mode with a paired-end, 50 base reads (PE50) sequencing strategy for the unstranded libraries and non-paired end for the stranded libraries. Image analysis and base calling was performed using Illumina's Real Time Analysis v1.18 software, followed by 'demultiplexing' of indexed reads and generation of FASTQ files, using Illumina's bcl2fastq Conversion Software v1.8.2. Analysis: Low-quality reads were removed prior to alignment to the KSHV (aka HHV-8) reference sequence (accession number NC_009333 for the KSHV GK18 strain) using TopHat2 (version 2.0.14) in a local instance of Galaxy. For quantitation purposes, the reads from paired-end libraries (non-stranded) were analyzed as unpaired (single-end data) to allow each read of the pair to map unambiguously to a single gene feature. The reads from both strands of the stranded libraries (non-paired end) were either concatenated and analyzed together (librarytype=unstranded) for comparison to the non-stranded paired end libraries or were analyzed separately (librarytype=FR) to show strand specificity. TopHat2 was used to detect splicing events ab initio. The default presets were used except that the maximum intron length was decreased to 10,000 and the maximum number of alignments allowed was decreased from 20 to 1, to avoid overcounting reads to repetitive regions. HTSeq (version 0.6) was used to quantitate the reads mapping to the unique set of UCDS gene features within the novel revised gene feature file "KSHV NC_009333 UCDS ver 121715.gff". The “intersection (non-empty)” setting in HTSeq was used to count all reads mapping completely or partially to a UCDS feature to maximize read count (featuretype=UCDS; IDattribute=gene). No reads were eliminated by ambiguity since the UCDS features were 50 bp apart, the length of a read. The read count was expressed as transcripts per million (TPM) by first normalizing the read count to reads per kilobase (RPK) by dividing the read counts by the length of the UCDS gene feature, in kilobases. The “per million” scaling factor was determined by summing all of the RPK values in a sample and dividing by 1,000,000. Each RPK value was then divided by the “per million” scaling factor to give TPM of mapped KSHV reads. Hierarchical clustering of TPM normalized expression levels was performed using the algorithm implemented in CIMminer. Hierarchical clustering of the gene correlation matrix was performed by calculating the Pearson correlation between the normalized transcript levels (TPM) associated with each pair of UCDS gene features, using a script in R to create and output the correlation matrix. Shiny web applications were developed for R-based principal component analysis (available at https://efg-ds.shinyapps.io/pcaApp/) and boxplot analysis of gene expression levels (available at https://efg-ds.shinyapps.io/boxplotApp/). Conclusions: We have used RNAseq to analyze and quantitate KSHV gene expression in a large collection of 41 KS tumor biopsies from HIV-infected individuals in Uganda that were naïve to ART. The RNAseq libraries were sequenced to an average depth of 100 million reads yielding up to 159,000 KSHV-mapped reads, and RNA reads mapping to non-overlapping UCDS features were quantitated using the new UCDS gene feature file. Phylogenetic analysis of the complete ORF75 sequence revealed the presence of at least six different KSHV strains in these tumor samples. We identified a set of transcripts from the latency region at the right end of the genome that was highly and consistently expressed in all the KS tumors. Variable expression of genes involved with transcription regulation and immune modulation was observed, with minimal expression of genes involved with viral replication, virion structure, or assembly. No correlation was observed between the KSHV transcript patterns detected in macular, nodular or fungating KS lesions. In contrast, the transcript patterns of different tumor lesions from the same individual showed the most similarity regardless of the tumor morphotype.
2018-10-30 | GSE116160 | GEO
Project description:KSHV is endemic in Uganda and the HIV epidemic has dramatically increased the incidence of Kaposi sarcoma (KS). To investigate the role of KSHV in the development of KS, we obtained KS biopsies from ART-naïve, HIV-positive individuals in Uganda and analyzed the tumors using RNAseq to globally characterize the KSHV transcriptome. Phylogenetic analysis of ORF75 sequences from 23 tumors revealed 6 distinct genetic clusters with KSHV strains exhibiting M, N or P alleles. RNA reads mapping to specific unique coding sequence (UCDS) features were quantitated using a gene feature file previously developed to globally analyze and quantitate KSHV transcription in infected endothelial cells. A pattern of high level expression was detected in the KSHV latency region that was common to all KS tumors. The clear majority of transcription was derived from the downstream latency transcript promoter P3(LTd) flanking ORF72, with little evidence of transcription from the P1(LTc) latency promoter, which is constitutive in KSHV-infected lymphomas and tissue-culture cells. RNAseq data provided evidence of alternate P3(LTd) transcript editing, splicing and termination resulting in multiple gene products, with 90% of the P3(LTd) transcripts spliced to release the intronic source of the microRNAs K1-9 and 11. The spliced transcripts encode a regulatory uORF upstream of Kaposin A with alterations in intervening repeat sequences yielding novel or deleted Kaposin B/C-like sequences. Hierarchical clustering and PCA analysis of KSHV transcripts revealed three clusters of tumors with different latent and lytic gene expression profiles. Paradoxically, tumors with a latent phenotype had high levels of total KSHV transcription, while tumors with a lytic phenotype had low levels of total KSHV transcription. Morphologically distinct KS tumors from the same individual showed similar KSHV gene expression profiles suggesting that the tumor microenvironment and host response play important roles in the activation level of KSHV within the infected tumor cells.
| S-EPMC6312348 | biostudies-literature
Project description:In the current study, we examined the cellular and viral miRNA expression profile in KSHV-infected patients without clinical signs of KSHV-associated disease. In addition, we characterized the viral miRNAs encoded by KSHV and EBV during natural infections in these patients. Finally, we show that malaria infection results in major alterations of the molecular content of plasma exosomes.
2018-02-16 | PXD008107 | Pride
Project description:Streptococcus sp. KS 6 strain:KS 6 Genome sequencing
| PRJNA556118 | ENA
Project description:Streptomyces sp. KS 21 strain:KS 21 Genome sequencing and assembly
| PRJNA501172 | ENA
Project description:Raoultella planticola strain:KS-EMBGN | isolate:KS-EMBGN Genome sequencing and assembly
| PRJNA534097 | ENA
Project description:Chromatin-organizing factors, like CTCF and cohesins, have been implicated in the control of complex viral regulatory programs. We investigated the role of CTCF and cohesin in the control of the latent to lytic switch for Kaposi's Sarcoma-Associated Herpesvirus (KSHV). We found that cohesin subunits, but not CTCF, were required for the repression of KSHV immediate early gene transcription. Depletion of cohesin subunits Rad21, SMC1, or SMC3 resulted in lytic cycle gene transcription and viral DNA replication. In contrast, depletion of CTCF failed to induce lytic transcription or DNA replication. ChiP-Seq analysis revealed that cohesins and CTCF bound to several sites within the immediate early control regions for ORF50 and more distal 5' sites that also regulate the divergently transcribed ORF45-46-47 gene cluster. Rad21 depletion led to a robust increase in ORF45 and ORF47 transcripts, with similar kinetics to that observed with chemical induction by sodium butyrate. During latency, the chromatin between the ORF45 and ORF50 transcription start sites was enriched in histone H3K4me3 with elevated H3K9ac at the ORF45 promoter and elevated H3K27me3 at the ORF50 promoter. A paused form of RNA pol II was loosely associated with the ORF45 promoter region during latency, but was converted to an active elongating form upon reactivation induced by Rad21 depletion. Butyrate-induced transcription of ORF45 and ORF47 was resistant to cyclohexamide, suggesting that these genes have immediate early features similar to ORF50. Butyrate-treatment caused the rapid dissociation of cohesins and loss of CTCF binding at the immediate early gene locus, suggesting that cohesins may be a direct target of butyrate-mediated lytic induction. Our findings implicate cohesins as a major repressor of KSHV lytic gene activation, and function coordinately with CTCF to regulate the switch between latent and lytic gene activity. Study of chromatin-organizing factors, like CTCF and cohesins.
2012-06-02 | E-GEOD-38411 | biostudies-arrayexpress
Project description:In sub-Saharan Africa, endemic Kaposi’s sarcoma (EnKS) is still prevalent despite high incidence of epidemic Kaposi’s sarcoma (EpKS) resulting from the on-going HIV-1 epidemic. While KSHV is clearly the etiologic agent of KS, the mechanisms underlying KS development are not fully understood. For example, HIV-1 co-infection and concomitant immune dysfunction have been associated with EpKS development. However, the direct or indirect role(s) of HIV-1, and therefore of immune suppression, in EpKS remains unclear. How, or whether, EpKS is mechanistically distinct from EnKS is unknown. Thus, the absence of HIV-1 co-infection in EnKS provides a unique control for investigating and deciphering whether HIV-1 plays a direct or indirect role in the EpKS tumor microenvironment. We hypothesized that HIV-1 co-infection would induce transcriptome changes that differentiate EpKS from EnKS, thereby defining the direct intra-tumor role of HIV-1 in KS. Comparison of ART-treated and -naïve patients would further define the impact of ART on the KS transcriptome. We utilized RNA-seq followed by multiparameter bioinformatics analysis to compare transcriptomes from KS lesions to uninvolved control skin. We provide the first transcriptomic comparison of EpKS versus EnKS, ART-treated vs. –naïve EpKS and male vs. female EpKS to define the roles of HIV-1 co-infection, the impact of ART, and gender on KS gene expression profiles. Our findings suggest that ART-use and gender have minimal impact on transcriptome profiles of KS lesions. Gene expression profiles strongly correlated between EpKS and EnKS patients (Spearman r=0.83, p<10-10). A subset of genes involved in tumorigenesis and inflammation/immune responses showed higher magnitude, but not unique dysregulation in EnKS compared to EpKS. While gender and ART had no detectable contribution, the trend toward higher magnitude of gene dysregulation in EnKS coupled with the absence of HIV-1 transcripts in EpKS may suggest an indirect or systemic effect of HIV-1 to promote KS tumorigenesis.
2020-06-26 | GSE147704 | GEO