Functional prediction of differentially expressed lncRNAs in HSV-1 infected human foreskin fibroblasts

Integrated transcriptome profiling in THP-1 macrophages infected with bunyavirus SFTSV

Severe fever with thrombocytopenia syndrome virus (SFTSV) is a tick-borne bunyavirus that causes an emerging hemorrhagic fever termed SFTS with high mortality. However, knowledge of SFTSV-host interactions is largely limited. Here, we performed a global transcriptome analysis of mRNAs and lncRNAs in THP-1 macrophages infected with SFTSV for 24 and 48 h. A total of 2,334 differentially expressed mRNAs and 154 differentially expressed lncRNAs were identified with 577 mRNAs and 31 lncRNAs commonly changed at both time points. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed mRNAs were mainly associated with innate immune, cytokine signaling, systemic lupus erythematosus, and alcoholism. Differentially expressed lncRNAs were enriched in systemic lupus erythematosus, alcoholism, and ribosome. Bioinformatic analysis also revealed hub regulatory mRNAs including IL6, TNF, UBA52, SRC, IL10, CXCL10, and CDK1 and core regulatory lncRNAs including XLOC_083027 and XLOC_113317. Transcription factor analysis of the differentially expressed mRNAs revealed that IRF1, SPI1, SPIB, ELF5, and FEV were enriched during SFTSV infection. Taken together, our studies illustrate the complex interaction between THP-1 macrophages and SFTSV.

Mechanism of interaction between virus and host is inferred from the changes of gene expression in macrophages infected with African swine fever virus CN/GS/2018 strain

BACKGROUND

African swine fever virus (ASFV) is a highly lethal virus that can infect porcine alveolar macrophages (PAMs). Since ASFV, China has dealt with a heavy blow to the pig industry. However, the effect of infection of ASFV strains isolated from China on PAM transcription level is not yet clarified.

METHODS

In this study, RNA sequencing (RNA-seq) was used to detect the differential expression of genes in PAMs at different time points after ASFV-CN/GS/2018 infection. The fluorescent quantitative polymerase chain reaction (qPCR) method was used to confirm the altered expression of related genes in PAMs infected with ASFV.

RESULTS

A total of 1154 differentially expressed genes were identified after ASFV-CN/GS/2018 infection, of which 816 were upregulated, and 338 were downregulated. GO and KEGG analysis showed that these genes were dynamically enriched in various biological processes, including innate immune response, inflammatory response, chemokines, and apoptosis. Furthermore, qPCR verified that the DEAD box polypeptide 58 (DDX58), Interferon-induced helicase C domain-containing protein 1 (IFIH1), Toll-like receptor 3 (TLR3), and TLR7 of PAMs were upregulated after ASFV infection, while TLR4 and TLR6 had a significant downward trend during ASFV infection. The expression of some factors related to antiviral and inflammation was altered significantly after ASFV infection, among which interferon-induced protein with tetratricopeptide repeats 1 (IFIT1), IFIT2, Interleukin-6 (IL-6) were upregulated, and Ewing's tumor-associated antigen 1 homolog (ETAA1) and Prosaposin receptor GPR37 (GPR37) were downregulated. In addition, we discovered that ASFV infection is involved in the regulation of chemokine expression in PAMs, and the chemokines, such as C-X-C motif chemokine 8 (CXCL8) and CXCL10, were upregulated after infection. However, the expression of chemokine receptor C-X-C chemokine receptor type 2 (CXCR2) is downregulated. Also, that the transcriptional levels of pro-apoptotic and anti-apoptotic factors changed after infection.

CONCLUSIONS

After ASFV-CN/GS/2018 infection, the expression of some antiviral and inflammatory factors in PAMs changed significantly. The ASFV infection may activates the RLR and TLR signaling pathways. In addition, ASFV infection is involved in regulating of chemokine expression in PAMs and host cell apoptosis.

Meta-analytic approach for transcriptome profiling of herpes simplex virus type 1

In this meta-analysis, we re-analysed and compared herpes simplex virus type 1 transcriptomic data generated by eight studies using various short- and long-read sequencing techniques and different library preparation methods. We identified a large number of novel mRNAs, non-coding RNAs and transcript isoforms, and validated many previously published transcripts. Here, we present the most complete HSV-1 transcriptome to date. Furthermore, we also demonstrate that various sequencing techniques, including both cDNA and direct RNA sequencing approaches, are error-prone, which can be circumvented by using integrated approaches. This work draws attention to the need for using multiple sequencing approaches and meta-analyses in transcriptome profiling studies to obtain reliable results.

Single-cell RNA-sequencing analysis identifies host long noncoding RNA MAMDC2-AS1 as a co-factor for HSV-1 nuclear transport

Herpes simplex virus (HSV) type 1 (HSV-1) infection exhibited high heterogeneity at individual cells level, including the different gene expression patterns and varying amounts of progeny virus. However, the underlying mechanism of such variability remains obscure. The importance of host long noncoding RNAs (lncRNAs) in virus infection had been recognized, while the contribution of lncRNAs to the heterogeneous infection remains unknown. Herein, a prior single-cell RNA sequencing data using HSV-1 reporter strain expressing ICP4-YFP was re-analyzed to obtain the differentially expressed lncRNA between the successfully initiated viral gene expression (ICP4-YFP⁺) cells and the aborted infection cells (ICP4-YFP^-). The ICP4-YFP⁺ population show a higher abundance of MAMDC2 antisense 1 (MAMDC2-AS1) lncRNA than ICP4-YFP^- population. MAMDC2-AS1 silencing reduces the expression of HSV-1 immediate early (IE) genes and limit HSV-1 infection in human host cells. Consistently, ectopic expression of MAMDC2-AS1 enhances HSV-1 IE genes transcription and facilitates the formation of HSV-1-induced plaques. Mechanically, both RNA-pull down and RNA immunoprecipitation assays show that MAMDC2-AS1 interacts with the RNA binding protein heat shock protein 90α (Hsp90α), a molecular chaperone involving in the nuclear import of HSV-1. The MAMDC2-AS1-Hsp90α interaction facilitates the nuclear transport of viral tegument protein VP16, the core factor initiating the expression of HSV-1 IE genes. The transcription factor YY1 mediates the induction of MAMDC2-AS1 upon HSV-1 infection. Our study elucidates the contribution of lncRNA to HSV-1 infection susceptibility in human cells and the role of Hsp90α RNA binding activity in HSV-1 infection.

Multiple Long-Read Sequencing Survey of Herpes Simplex Virus Dynamic Transcriptome

Long-read sequencing (LRS) has become increasingly important in RNA research due to its strength in resolving complex transcriptomic architectures. In this regard, currently two LRS platforms have demonstrated adequate performance: the Single Molecule Real-Time Sequencing by Pacific Biosciences (PacBio) and the nanopore sequencing by Oxford Nanopore Technologies (ONT). Even though these techniques produce lower coverage and are more error prone than short-read sequencing, they continue to be more successful in identifying polycistronic RNAs, transcript isoforms including splice and transcript end variants, as well as transcript overlaps. Recent reports have successfully applied LRS for the investigation of the transcriptome of viruses belonging to various families. These studies have substantially increased the number of previously known viral RNA molecules. In this work, we used the Sequel and MinION technique from PacBio and ONT, respectively, to characterize the lytic transcriptome of the herpes simplex virus type 1 (HSV-1). In most samples, we analyzed the poly(A) fraction of the transcriptome, but we also performed random oligonucleotide-based sequencing. Besides cDNA sequencing, we also carried out native RNA sequencing. Our investigations identified more than 2,300 previously undetected transcripts, including coding, and non-coding RNAs, multi-splice transcripts, as well as polycistronic and complex transcripts. Furthermore, we found previously unsubstantiated transcriptional start sites, polyadenylation sites, and splice sites. A large number of novel transcriptional overlaps were also detected. Random-primed sequencing revealed that each convergent gene pair produces non-polyadenylated read-through RNAs overlapping the partner genes. Furthermore, we identified novel replication-associated transcripts overlapping the HSV-1 replication origins, and novel LAT variants with very long 5' regions, which are co-terminal with the LAT-0.7kb transcript. Overall, our results demonstrated that the HSV-1 transcripts form an extremely complex pattern of overlaps, and that entire viral genome is transcriptionally active. In most viral genes, if not in all, both DNA strands are expressed.

Variation in the Untranslated Genome and Susceptibility to Infections

The clinical outcomes of infections are highly variable among individuals and are determined by complex host-pathogen interactions. Genome-wide association studies (GWAS) are powerful tools to unravel common genetic variations that are associated with disease risk and clinical outcomes. However, GWAS has only rarely revealed information on the exact genetic elements and their effects underlying an association because the majority of the hits are within non-coding regions. Some of the variants or the linked polymorphisms are now being discovered to have functional significance, such as regulatory elements in the promoter and enhancer regions or the microRNA binding sites in the 3′untranslated region of the protein-coding genes, which influence transcription, RNA stability, and translation of the protein-coding genes. However, only 3% of the entire transcriptome is protein-coding, signifying that non-coding RNAs represent most of the transcripts. Thus, a large portion of previously identified intergenic GWAS single nucleotide polymorphisms (SNPs) is in the non-coding RNAs. The non-coding RNAs form a large-scale regulatory network across the transcriptome, greatly expanding the complexity of gene regulation. Accumulating evidence also suggests that the “non-coding” genome regions actively regulate the highly dynamic three dimensional (3D) chromatin structures, which are critical for genome function. Epigenetic modulation like DNA methylation and histone modifications further affect chromatin accessibility and gene expression adding another layer of complexity to the functional interpretation of genetic variation associated with disease outcomes. We provide an overview of the current information on the influence of variation in these “untranslated” regions of the human genome on infectious diseases. The focus of this review is infectious disease-associated polymorphisms and gene regulatory mechanisms of pathophysiological relevance.

Transcriptional Profile and Integrative Analyses of Long Noncoding RNAs in Primary Human Corneal Epithelial Cells in Response to HSV-1 Infection

PURPOSE

Long noncoding RNAs (lncRNAs) have been demonstrated to have important regulatory functions in diverse cellular processes; however, the role of lncRNAs in the pathogenesis of herpes simplex keratitis (HSK) remains poorly understood.

METHODS

Primary human corneal epithelial cells (HCECs) were infected with herpes simplex virus-1 (HSV-1) and the total RNAs extracted from both the infected group and the mock-infected group subjected to microarray analysis to identify the differential expression of lncRNAs and mRNAs. We also performed bioinformatic analysis including gene ontology (GO) analysis, pathway analysis and co-expression network analysis.

RESULTS

Compared with mock-infected group, the expression of thousands of lncRNAs and mRNAs were significantly changed, and the microarray results were validated by qRT-PCR. The most enriched GOs targeted by up-regulated transcripts were defense response, intrinsic component of plasma membrane and cytokine activity，and the most enriched GOs targeted by the down-regulated transcripts were cellular metabolic process, intracellular part and poly (A) RNA binding. Pathway analysis indicated that the most correlated pathways for up- and down-regulated transcripts were cytokine-cytokine receptor interaction and RNA transport, respectively.

CONCLUSIONS

Our study identified the genome-wide profile of lncRNAs and mRNAs expression in primary corneal epithelial cells with HSV-1 infection. These transcriptomic data together with subsequent bioinformatic analysis will provide us with novel clue to the insight into molecular mechanism and potential therapeutic targets of HSK. Further studies are expected to verify the potentially functional genes and pathways and explore the critical lncRNAs.

ABBREVIATIONS

Long noncoding RNAs: lncRNAs; herpes simplex virus-1: HSV-1; herpes simplex virus keratitis: HSK; human corneal epithelial cells: HCECs.

Expression of pseudorabies virus-encoded long noncoding RNAs in epithelial cells and neurons

Long noncoding RNAs (lncRNAs) play important roles in regulating eukaryotic genome replication and gene expression in diverse biological systems. Here, we identified lncRNAs transcribed from pseudorabies virus (PRV)-infected PK-15 cells. Based on high-throughput sequencing data, we obtained 87,263,926 and 93,947,628 clean reads from mock-infected and PRV-infected PK-15 cells, respectively. Through a normalized analytic protocol, we identified three novel viral lncRNAs. According to an analysis of differential expression between the mock-infected and PRV-infected cells, 4151 host lncRNAs were significantly upregulated and 2327 host lncRNAs were significantly downregulated in the latter group. Viral lncRNAs and several host lncRNAs were verified by northern blotting and real-time PCR. The findings showed that the viral lncRNA LDI might regulate the expression of IE180, a potent transcriptional activator of viral genes. Furthermore, we characterized the expression of viral lncRNAs in a culture of infected primary chicken dorsal root ganglia (DRG). Collectively, the obtained data suggest that PRV generates lncRNAs in both epithelial cells and chick DRG neurons.

Human foreskin fibroblasts: from waste bag to important biomedical applications

Circumcision is one of the most performed surgical procedures worldwide, and it is estimated that one in three men worldwide is circumcised, which makes the preputial skin removed after surgery an abundant material for possible applications. In particular, it is possible efficiently to isolate the cells of the foreskin, with fibroblasts being the most abundant cells of the dermis and the most used in biomedical research. This work aimed to review the knowledge and obtain a broad view of the main applications of human foreskin fibroblast cell culture. A literature search was conducted, including clinical trials, preclinical basic research studies, reviews and experimental studies. Several medical and laboratory applications of human foreskin fibroblast cell culture have been described, especially when it comes to the use of human foreskin fibroblasts as feeder cells for the cultivation of human embryonic stem cells, in addition to co-culture with other cell types. The culture of foreskin fibroblasts has also been used to: obtain induced pluripotent stem cells; the diagnosis of Clostridium difficile; to test the toxicity and effect of substances on normal cells, especially the toxicity of possible antineoplastic drugs; in viral culture, mainly of the human cytomegalovirus, study of the pathogenesis of other microorganisms; varied studies of cellular physiology and cellular interactions. Fibroblasts are important for cell models for varied application cultures, demonstrating how the preputial material can be reused, making possible new applications.

Level of evidence: Not applicable for this multicentre audit.

ZIKV infection effects changes in gene splicing, isoform composition and lncRNA expression in human neural progenitor cells

BACKGROUND

The Zika virus (ZIKV) is a mosquito-borne flavivirus that causes microcephaly and Guillain-Barré syndrome in infected individuals. To obtain insights into the mechanism of ZIKV infection and pathogenesis, we analyzed the transcriptome of ZIKV infected human neural progenitor cells (hNPCs) for changes in alternative splicing (AS), gene isoform (ISO) composition and long noncoding RNAs (lncRNAs) expression.

METHODS

We analyzed differentially expressed lncRNAs, AS, ISO from RNA-seq data in ZIKV infected hNPCs.

RESULTS

We obtained 149 differentially expressed lncRNAs, including potential viral targets to modulate cellular processes such as cell cycle, apoptosis and immune response. The infection induced 262 cases of AS occurring in 229 genes, which were enriched in cell death, RNA processing, transport, and neuron development. Among 691 differentially expressed ISOs, upregulated ISOs were enriched in signaling, regulation of transcription, and amino acid biosynthesis, while downregulated ISOs were mostly enriched in cell cycle. Importantly, these analyses revealed specific links between ZIKV induced changes in cellular pathways and the type of changes in the host transcriptome, suggesting important regulatory mechanisms.

CONCLUSIONS

Our analyses revealed candidate lncRNAs, AS events and ISOs which may function in ZIKV infection induced cell cycle disruption, apoptosis and attenuation of neurogenesis, and shed light on the roles of lncRNAs, AS and ISOs in virus-host interactions, and would facilitate future studies of ZIKV infection and pathogenesis.

Widespread activation of antisense transcription of the host genome during herpes simplex virus 1 infection

Background

Herpesviruses can infect a wide range of animal species. Herpes simplex virus 1 (HSV-1) is one of the eight herpesviruses that can infect humans and is prevalent worldwide. Herpesviruses have evolved multiple ways to adapt the infected cells to their needs, but knowledge about these transcriptional and post-transcriptional modifications is sparse.

Results

Here, we show that HSV-1 induces the expression of about 1000 antisense transcripts from the human host cell genome. A subset of these is also activated by the closely related varicella zoster virus. Antisense transcripts originate either at gene promoters or within the gene body, and they show different susceptibility to the inhibition of early and immediate early viral gene expression. Overexpression of the major viral transcription factor ICP4 is sufficient to turn on a subset of antisense transcripts. Histone marks around transcription start sites of HSV-1-induced and constitutively transcribed antisense transcripts are highly similar, indicating that the genetic loci are already poised to transcribe these novel RNAs. Furthermore, an antisense transcript overlapping with the BBC3 gene (also known as PUMA) transcriptionally silences this potent inducer of apoptosis in cis.

Conclusions

We show for the first time that a virus induces widespread antisense transcription of the host cell genome. We provide evidence that HSV-1 uses this to downregulate a strong inducer of apoptosis. Our findings open new perspectives on global and specific alterations of host cell transcription by viruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1329-5) contains supplementary material, which is available to authorized users.

Identification and characterization of a long non-coding RNA up-regulated during HIV-1 infection

Long non-coding RNAs (lncRNAs) are rapidly emerging as important regulators of a diverse array of cellular functions. Here, we describe a meta-analysis of two independent RNA-seq studies to identify lncRNAs that are differentially expressed upon HIV-1 infection. Only three lncRNA genes exhibited altered expression of ≥ 2-fold in HIV-1-infected cells. Of these, the uncharacterized lncRNA LINC00173 was chosen for further study. Both transcript variants of LINC00173 (lnc173 TSV1 and 2) could be detected by qPCR, localized predominantly to the nucleus and were reproducibly up-regulated during infection. Knock-out of the LINC00173 locus did not have detectable effects on HIV-1 replication. Interestingly, however, stimulation of Jurkat T cells with PMA/ionomycin resulted in a decrease of lnc173 expression, and Jurkat cells deficient for lnc173 on average expressed higher levels of specific cytokines than control cells. These data suggest that lnc173 may have a role in the regulation of cytokines in T cells.

Transcriptomic profiling of long non-coding RNAs in hepatitis B virus-related hepatocellular carcinoma

Long non-coding RNAs (lncRNAs) have been reported to be involved in the development and progression of hepatocellular carcinoma (HCC). However, few studies have focus on the dyregulation and the role of lncRNAs in HBV-related HCC. We performed a comprehensive analysis of lncRNAs expression profile in HBV-related HCC tissues samples using deep sequencing. We revealed that a total of 1242 lncRNA transcripts (983 up-regulated and 259 down-regulated) and 1841 mRNA transcripts were significantly differentially expressed in HBV-related HCC patients. Pathway and gene ontology analysis showed that they are involved in the biological process related to HCC development by cis-regulation of co-expressed protein-coding genes. 10 candidate lncRNAs were selected and validated with quantitative real-time PCR analysis. Furthermore, we found that one of most down-regulated lncRNAs, n346077, could suppress HCC cells invasion and migration in vitro. Our findings provide an overview of aberrantly expressed lncRNAs in HBV-related HCC and will be useful for further functional studies of lncRNAs in HBV-related pathogenesis.

Long-Read Isoform Sequencing Reveals a Hidden Complexity of the Transcriptional Landscape of Herpes Simplex Virus Type 1

In this study, we used the amplified isoform sequencing technique from Pacific Biosciences to characterize the poly(A)⁺ fraction of the lytic transcriptome of the herpes simplex virus type 1 (HSV-1). Our analysis detected 34 formerly unidentified protein-coding genes, 10 non-coding RNAs, as well as 17 polycistronic and complex transcripts. This work also led us to identify many transcript isoforms, including 13 splice and 68 transcript end variants, as well as several transcript overlaps. Additionally, we determined previously unascertained transcriptional start and polyadenylation sites. We analyzed the transcriptional activity from the complementary DNA strand in five convergent HSV gene pairs with quantitative RT-PCR and detected antisense RNAs in each gene. This part of the study revealed an inverse correlation between the expressions of convergent partners. Our work adds new insights for understanding the complexity of the pervasive transcriptional overlaps by suggesting that there is a crosstalk between adjacent and distal genes through interaction between their transcription apparatuses. We also identified transcripts overlapping the HSV replication origins, which may indicate an interplay between the transcription and replication machineries. The relative abundance of HSV-1 transcripts has also been established by using a novel method based on the calculation of sequencing reads for the analysis.