The transcriptome of the adenovirus infected cell

Differential innate immune responses to fowl adenovirus serotype 4 infection in Leghorn male hepatocellular and chicken embryo fibroblast cells

Fowl adenovirus serotype 4 (FAdV-4) infections result in substantial economic losses in the poultry industry. Recent findings have revealed that FAdV-4 significantly suppresses the host immune response upon infection; however, the specific viral and host factors contributing to this immunomodulatory activity remain poorly characterized. Moreover, diverse cell types exhibit differential immune responses to FAdV-4 infection. To elucidate cell-specific host responses, we performed transcriptomic analysis of FAdV-4 infected leghorn male hepatocellular (LMH) and chicken embryo fibroblast (CEF) cells. Although FAdV-4 replicated more efficiently in LMH cells, it provoked limited interferon-stimulated gene induction. In contrast, FAdV-4 infection triggered robust antiviral responses in CEF cells, including upregulation of cytosolic DNA sensing and interferon-stimulated genes. Knockdown of key cytosolic DNA sensing molecules enhanced FAdV-4 replication in LMH cells while reducing interferon-stimulated gene expression. Our findings reveal cell-specific virus-host interactions that provide insight into FAdV-4 pathogenesis while identifying factors that mediate antiviral immunity against FAdV-4.

Identification of functional pathways and potential genes associated with interferon signaling during human adenovirus type 7 infection by weighted gene coexpression network analysis

Human adenovirus type 7 (HAdV-7) can cause severe pneumonia and complications in children. However, the mechanism of pathogenesis and the genes involved remain largely unknown. We collected HAdV-7-infected and mock-infected A549 cells at 24, 48, and 72 hours postinfection (hpi) for RNA sequencing (RNA-Seq) and identified potential genes and functional pathways associated with HAdV-7 infection using weighted gene coexpression network analysis (WGCNA). Based on bioinformatics analysis, 12 coexpression modules were constructed by WGCNA, with the blue, tan, and brown modules significantly positively correlated with adenovirus infection at 24, 48, and 72 hpi, respectively. Functional enrichment analysis indicated that the blue module was mainly enriched in DNA replication and viral processes, the tan module was largely enriched in metabolic pathways and regulation of superoxide radical removal, and the brown module was predominantly enriched in regulation of cell death. qPCR was used to determine transcript abundance of some identified hub genes, and the results were consistent with those from RNA-Seq. Comprehensively analyzing hub genes and differentially expressed genes in the GSE68004 dataset, we identified SOCS3, OASL, ISG15, and IFIT1 as potential candidate genes for use as biomarkers or drug targets in HAdV-7 infection. We propose a multi-target inhibition of the interferon signaling mechanism to explain the association of HAdV-7 infection with the severity of clinical consequences. This study has allowed us to construct a framework of coexpression gene modules in A549 cells infected with HAdV-7, thus providing a basis for identifying potential genes and pathways involved in adenovirus infection and for investigating the pathogenesis of adenovirus-associated diseases.

Novel viral splicing events and open reading frames revealed by long-read direct RNA sequencing of adenovirus transcripts

Adenovirus is a common human pathogen that relies on host cell processes for transcription and processing of viral RNA and protein production. Although adenoviral promoters, splice junctions, and polyadenylation sites have been characterized using low-throughput biochemical techniques or short read cDNA-based sequencing, these technologies do not fully capture the complexity of the adenoviral transcriptome. By combining Illumina short-read and nanopore long-read direct RNA sequencing approaches, we mapped transcription start sites and RNA cleavage and polyadenylation sites across the adenovirus genome. In addition to confirming the known canonical viral early and late RNA cassettes, our analysis of splice junctions within long RNA reads revealed an additional 35 novel viral transcripts that meet stringent criteria for expression. These RNAs include fourteen new splice junctions which lead to expression of canonical open reading frames (ORFs), six novel ORF-containing transcripts, and 15 transcripts encoding for messages that could alter protein functions through truncation or fusion of canonical ORFs. In addition, we detect RNAs that bypass canonical cleavage sites and generate potential chimeric proteins by linking distinct gene transcription units. Among these chimeric proteins we detected an evolutionarily conserved protein containing the N-terminus of E4orf6 fused to the downstream DBP/E2A ORF. Loss of this novel protein, E4orf6/DBP, was associated with aberrant viral replication center morphology and poor viral spread. Our work highlights how long-read sequencing technologies combined with mass spectrometry can reveal further complexity within viral transcriptomes and resulting proteomes.

Transcriptome Analysis Reveals the Potential Role of Long Noncoding RNAs in Regulating Fowl Adenovirus Serotype 4-Induced Apoptosis in Leghorn Male Hepatocellular Cells

Hepatitis-hydropericardium syndrome (HHS) is caused by fowl adenovirus serotype 4 (FAdV-4) and has resulted in considerable economic losses to the poultry industry globally. FAdV-4 elicits apoptosis in host cells. Long noncoding RNAs (lncRNAs) have emerged as important regulatory RNAs with profound effects on various biological processes, including apoptosis. However, it remains unknown whether lncRNAs participate in FAdV-4-induced apoptosis. In this study, RNA sequencing was applied to determine the transcription of cellular lncRNA in leghorn male hepatocellular (LMH) cells infected with FAdV-4. Cellular RNA transcription analysis demonstrated that FAdV-4 infection elicited 1798 significantly differentially expressed (DE) lncRNAs in infected LMH cells at 24 h post-infection (hpi) compared to mock control infection. In addition, 2873 DE mRNAs were also found. Target prediction and analyses revealed that 775 DE lncRNAs whose 671 target mRNAs were among the DE mRNAs were involved in several signaling pathways, including the AMPK signaling pathway, p53 signaling pathway and insulin signaling pathway. From these 775 DE lncRNAs, we identified 71 DE lncRNAs related to apoptosis based on their target gene functions. Subsequently, lncRNA 54128 was selected from the 71 identified DE lncRNAs, and its role in FAdV-4-induced apoptosis was verified. LncRNA 54128 interference significantly suppressed the rate of apoptosis, which was accompanied by reduced BMP4 transcription levels. To the best of our knowledge, this is the first study to analyze host lncRNA transcription during FAdV-4 infection. Our findings provide a better understanding of host responses to FAdV-4 infection and provide new directions for understanding the potential association between lncRNAs and FAdV-4 pathogenesis.

Myosin 7b is a regulatory long noncoding RNA (lncMYH7b) in the human heart

Myosin heavy chain 7b (MYH7b) is an ancient member of the myosin heavy chain motor protein family that is expressed in striated muscles. In mammalian cardiac muscle, MYH7b RNA is expressed along with two other myosin heavy chains, β-myosin heavy chain (β-MyHC) and α-myosin heavy chain (α-MyHC). However, unlike β-MyHC and α-MyHC, which are maintained in a careful balance at the protein level, the MYH7b locus does not produce a full-length protein in the heart due to a posttranscriptional exon-skipping mechanism that occurs in a tissue-specific manner. Whether this locus has a role in the heart beyond producing its intronic microRNA, miR-499, was unclear. Using cardiomyocytes derived from human induced pluripotent stem cells as a model system, we found that the noncoding exon-skipped RNA (lncMYH7b) affects the transcriptional landscape of human cardiomyocytes, independent of miR-499. Specifically, lncMYH7b regulates the ratio of β-MyHC to α-MyHC, which is crucial for cardiac contractility. We also found that lncMYH7b regulates beat rate and sarcomere formation in cardiomyocytes. This regulation is likely achieved through control of a member of the TEA domain transcription factor family (TEAD3, which is known to regulate β-MyHC). Therefore, we conclude that this ancient gene has been repurposed by alternative splicing to produce a regulatory long-noncoding RNA in the human heart that affects cardiac myosin composition.

Cell-to-cell and genome-to-genome variability of adenovirus transcription tuned by the cell cycle

In clonal cultures, not all cells are equally susceptible to virus infection, and the mechanisms underlying this are poorly understood. Here, we developed image-based single-cell measurements to scrutinize the heterogeneity of adenovirus (AdV) infection. AdV delivers, transcribes and replicates a linear double-stranded DNA genome in the nucleus. We measured the abundance of viral transcripts using single-molecule RNA fluorescence in situ hybridization (FISH) and the incoming 5-ethynyl-2'-deoxycytidine (EdC)-tagged viral genomes using a copper(I)-catalyzed azide-alkyne cycloaddition (click) reaction. Surprisingly, expression of the immediate early gene E1A only moderately correlated with the number of viral genomes in the cell nucleus. Intranuclear genome-to-genome heterogeneity was found at the level of viral transcription and, in accordance, individual genomes exhibited heterogeneous replication activity. By analyzing the cell cycle state, we found that G1 cells exhibited the highest E1A gene expression and displayed increased correlation between E1A gene expression and viral genome copy numbers. The combined image-based single-molecule procedures described here are ideally suited to explore the cell-to-cell variability in viral gene expression in a range of different settings, including the innate immune response.

Cell Bank Origin of MDCK Parental Cells Shapes Adaptation to Serum-Free Suspension Culture and Canine Adenoviral Vector Production

Phenotypic variation in cultured mammalian cell lines is known to be induced by passaging and culture conditions. Yet, the effect these variations have on the production of viral vectors has been overlooked. In this work we evaluated the impact of using Madin-Darby canine kidney (MDCK) parental cells from American Type Culture Collection (ATCC) or European Collection of Authenticated Cell Cultures (ECACC) cell bank repositories in both adherent and suspension cultures for the production of canine adenoviral vectors type 2 (CAV-2). To further explore the differences between cells, we conducted whole-genome transcriptome analysis. ECACC's MDCK showed to be a less heterogeneous population, more difficult to adapt to suspension and serum-free culture conditions, but more permissive to CAV-2 replication progression, enabling higher yields. Transcriptome data indicated that this increased permissiveness is due to a general down-regulation of biological networks of innate immunity in ECACC cells, including apoptosis and death receptor signaling, Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling, toll-like receptors signaling and the canonical pathway of nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling. These results show the impact of MDCK source on the outcome of viral-based production processes further elucidating transcriptome signatures underlying enhanced adenoviral replication. Following functional validation, the genes and networks identified herein can be targeted in future engineering approaches aiming at improving the production of CAV-2 gene therapy vectors.

Efficient gene editing of human long-term hematopoietic stem cells validated by clonal tracking

Targeted gene editing in hematopoietic stem cells (HSCs) is a promising treatment for several diseases. However, the limited efficiency of homology-directed repair (HDR) in HSCs and the unknown impact of the procedure on clonal composition and dynamics of transplantation have hampered clinical translation. Here, we apply a barcoding strategy to clonal tracking of edited cells (BAR-Seq) and show that editing activates p53, which substantially shrinks the HSC clonal repertoire in hematochimeric mice, although engrafted edited clones preserve multilineage and self-renewing capacity. Transient p53 inhibition restored polyclonal graft composition. We increased HDR efficiency by forcing cell-cycle progression and upregulating components of the HDR machinery through transient expression of the adenovirus 5 E4orf6/7 protein, which recruits the cell-cycle controller E2F on its target genes. Combined E4orf6/7 expression and p53 inhibition resulted in HDR editing efficiencies of up to 50% in the long-term human graft, without perturbing repopulation and self-renewal of edited HSCs. This enhanced protocol should broaden applicability of HSC gene editing and pave its way to clinical translation.

Design, validation and evaluation of a SYBR green-based quantitative PCR array for comprehensive analysis of adenovirus type 5 transcriptional patterns

The adenoviral genome encodes coordinately expressed early and late gene transcriptional units that specify a complex collection of extensively spliced overlapping mRNAs. These complexities confound the generation of compatible, validated and optimized qPCR assays that permit comprehensive evaluation of adenoviral transcription. We have developed and evaluated a compilation of qPCR assays that represent the majority of the human adenovirus 5 (hAdV5) genome and allow for absolute and relative quantification of transcriptional activity. A panel of specific adenovirus gene primer pairs was designed through computational modeling to be compatible under a single reaction condition, precisely amplify spliced transcript products within each gene class, and not result in cellular or viral RNA/DNA background amplification. Primer pairs and reaction conditions were optimized to generate a single amplification product that was specific for its target amplicon with minimal intra-assay variability. The specificity of target amplicons was confirmed by dissociation curve analysis, gel electrophoresis and sequencing. In all, thirty-two primer sets representing specific gene products, as well as, pan early and late gene regions were validated under identical amplification conditions, thereby enabling a comprehensive assessment of adenoviral transcription within a single plate array. In order to generate positive control templates and to facilitate absolute quantification of gene expression, all target amplicons were cloned to create gene target-specific standards. These plasmid amplicon controls demonstrated that the SYBR qPCR assays exhibited optimal amplification efficiencies with a high sensitivity of detection to less than 10 copies and a linear amplification across at least eight orders of magnitude. The effectiveness and utility of the comprehensive adenoviral transcriptional array was assessed by investigating the changes in Ad5Wt gene expression at 72 versus 24 h post infection. Predictably, overall gene expression was globally increased at 72 h post infection; however, levels of E2 and Late transcripts exhibited the greatest increased expression, reflecting their necessity at this time point for genomic replication and virion assembly. Taken together, these data demonstrate that the adenoviral qPCR transcriptional array is a modular, scalable, and cost-effective method to comprehensively and accurately assess hAdV5 gene transcription. This array is broadly applicable to facilitate: adenoviral vector development; assessment of cell complementation of knockout viruses; antiviral mechanism of action evaluation; next-generation sequencing data validation.

Phosphorylation Time-Course Study of the Response during Adenovirus Type 2 Infection

PTMs such as phosphorylations are usually involved in signal transduction pathways. To investigate the temporal dynamics of phosphoproteome changes upon viral infection, a model system of IMR-90 cells infected with human adenovirus type 2 (Ad2) is used in a time-course quantitative analysis combining titanium dioxide (TiO₂ ) particle enrichment and SILAC-MS. Quantitative data from 1552 phosphorylated sites clustered the highly altered phosphorylated sites to the signaling by rho family GTPases, the actin cytoskeleton signaling, and the cAMP-dependent protein kinase A signaling pathways. Their activation is especially pronounced at early time post-infection. Changes of several phosphorylated sites involved in the glycolysis pathway, related to the activation of the Warburg effect, point at virus-induced energy production. For Ad2 proteins, 32 novel phosphorylation sites are identified and as many as 52 phosphorylated sites on 17 different Ad2 proteins are quantified, most of them at late time post-infection. Kinase predictions highlighted activation of PKA, CDK1/2, MAPK, and CKII. Overlaps of kinase motif sequences for viral and human proteins are observed, stressing the importance of phosphorylation during Ad2 infection.

Adenovirus in the omics era - a multipronged strategy

Human adenoviruses (HAdVs) are common pathogens associated with a wide variety of respiratory, ocular, and gastrointestinal diseases. To achieve its effective lytic mode of replication, HAdVs have to reprogram host-cell gene expression and fine-tune viral gene expression in a temporal manner. In two decades, omics revolution has advanced our knowledge about the HAdV and host-cell interplay at the RNA and protein levels. This review summarizes the current knowledge from large-scale datasets on how HAdV infections adjust coding and noncoding RNA expression, as well as how they reprogram host-cell proteome during the lytic course of infection.

The oncolytic adenovirus VCN-01 promotes anti-tumor effect in primitive neuroectodermal tumor models

Last advances in the treatment of pediatric tumors has led to an increase of survival rates of children affected by primitive neuroectodermal tumors, however, still a significant amount of the patients do not overcome the disease. In addition, the survivors might suffer from severe side effects caused by the current standard treatments. Oncolytic virotherapy has emerged in the last years as a promising alternative for the treatment of solid tumors. In this work, we study the anti-tumor effect mediated by the oncolytic adenovirus VCN-01 in CNS-PNET models. VCN-01 is able to infect and replicate in PNET cell cultures, leading to a cytotoxicity and immunogenic cell death. In vivo, VCN-01 increased significantly the median survival of mice and led to long-term survivors in two orthotopic models of PNETs. In summary, these results underscore the therapeutic effect of VCN-01 for rare pediatric cancers such as PNETs, and warrants further exploration on the use of this virus to treat them.

Disparate Entry of Adenoviruses Dictates Differential Innate Immune Responses on the Ocular Surface

Human adenovirus infection of the ocular surface is associated with severe keratoconjunctivitis and the formation of subepithelial corneal infiltrates, which may persist and impair vision for months to years following infection. Long term pathology persists well beyond the resolution of viral replication, indicating that the prolonged immune response is not virus-mediated. However, it is not clear how these responses are sustained or even initiated following infection. This review discusses recent work from our laboratory and others which demonstrates different entry pathways specific to both adenovirus and cell type. These findings suggest that adenoviruses may stimulate specific pattern recognition receptors in an entry/trafficking-dependent manner, leading to distinct immune responses dependent on the virus/cell type combination. Additional work is needed to understand the specific connections between adenoviral entry and the stimulation of innate immune responses by the various cell types present on the ocular surface.

Encounters with adenovirus

In this paper I describe aspects of work on the human adenoviruses in which my laboratory has participated. It consists of two sections-one historic dealing with work performed in the previous century, and one dealing with the application of 'omics' technologies to understand how adenovirus-infected cells become reprogrammed to benefit virus multiplication.

RNA-Seq analysis of duck embryo fibroblast cell gene expression during the early stage of egg drop syndrome virus infection

Egg drop syndrome virus (EDSV), a member of the family Adenoviridae and an economically important pathogen with a broad host range, leads to markedly decreased egg production. However, the molecular mechanism underlying the host-EDSV interaction remains unclear. Here, we performed high-throughput RNA sequencing (RNA-Seq) to study the dynamic changes in host gene expression at 6, 12, and 24 hours post-infection in duck embryo fibroblasts (DEFs) infected with EDSV. Atotal of 441 differentially expressed genes (DEGs) were identified after EDSV infection. Gene Ontology category and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that these DEGs were associated with multiple biological functions, including signal transduction, host immunity, virus infection, cell apoptosis, cell proliferation, and pathogenicity-related and metabolic process signaling pathways. We screened and identified 12 DEGs for further examination by using qRT-PCR. The qRT-PCR and RNA-Seq results were highly consistent. This study analyzed viral infection and host immunity induced by EDSV infection from a novel perspective, and the results provide valuable information regarding the mechanisms underlying host-EDSV interactions, which will prove useful for the future development of antiviral drugs or vaccines for poultry, thus benefiting the entire poultry industry.

Temporal dynamics of adenovirus 5 gene expression in normal human cells

Adenovirus executes a finely tuned transcriptional program upon infection of a cell. To better understand the temporal dynamics of the viral transcriptional program we performed highly sensitive digital PCR on samples extracted from arrested human lung fibroblasts infected with human adenovirus 5 strain dl309. We show that the first transcript made from viral genomes is the virus associated non-coding RNA, in particular we detected abundant levels of virus associated RNA II four hours after infection. Activation of E1 and E4 occurred nearly simultaneously later in infection, followed by other early genes as well as late genes. Our study determined that genomes begin to replicate between 29 and 30 hours after infection. This study provides a comprehensive view of viral mRNA steady-state kinetics in arrested human cells using digital PCR.

Transcriptomic and proteomic analyses reveal new insights into the regulation of immune pathways during adenovirus type 2 infection

Background

Human adenovirus (Ad) infection leads to the changes of host cell gene expression and biosynthetic processes. Transcriptomics in adenovirus type 2 (Ad2)-infected lung fibroblasts (IMR-90) cells has previously been studied using RNA sequencing. However, this study included only two time points (12 and 24 hpi) using constrained 76 bp long sequencing reads. Therefore, a more detailed study of transcription at different phases of infection using an up-graded sequencing technique is recalled. Furthermore, the correlation between transcription and protein expression needs to be addressed.

Results

In total, 3556 unique cellular genes were identified as differentially expressed at the transcriptional level with more than 2-fold changes in Ad2-infected cells as compared to non-infected cells by using paired-end sequencing. Based on the kinetics of the gene expression changes at different times after infection, these RNAs fell into 20 clusters. Among them, cellular genes involved in immune response were highly up-regulated in the early phase before becoming down-regulated in the late phase. Comparison of differentially expressed genes at transcriptional and posttranscriptional levels revealed low correlation. Particularly genes involved in cellular immune pathways showed a negative correlation. Here, we highlight the genes which expose inconsistent expression profiles with an emphasis on key factors in cellular immune pathways including NFκB, JAK/STAT, caspases and MAVS. Different from their transcriptional profiles with up- and down-regulation in the early and late phase, respectively, these proteins were up-regulated in the early phase and were sustained in the late phase. A surprising finding was that the target genes of the sustained activators failed to show response.

Conclusion

There were features common to genes which play important roles in cellular immune pathways. Their expression was stimulated at both RNA and protein levels during the early phase. In the late phase however, their transcription was suppressed while protein levels remained stable. These results indicate that Ad2 and the host cell use different strategies to regulate cellular immune pathways. A control mechanism at the post-translational level must thus exist which is under the control of Ad2.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1375-5) contains supplementary material, which is available to authorized users.

Time-resolved proteomics of adenovirus infected cells

Viral infections cause large problems in the world and deeper understanding of the disease mechanisms is needed. Here we present an analytical strategy to investigate the host cell protein changes during human adenovirus type 2 (HAdV-C2 or Ad2) infection of lung fibroblasts by stable isotope labelling of amino acids in cell culture (SILAC) and nanoLC-MS/MS. This work focuses on early phase of infection (6 and 12 h post-infection (hpi)) but the data is combined with previously published late phase (24 and 36 hpi) proteomics data to produce a time series covering the complete infection. As many as 2169 proteins were quantitatively monitored from 6 to 36 hpi, while some proteins were time-specific. After applying different filter criteria, 2027 and 2150 proteins were quantified at 6 and 12 hpi and among them, 431 and 544 were significantly altered at the two time points. Pathway analysis showed that the De novo purine and pyrimidine biosynthesis, Glycolysis and Cytoskeletal regulation by Rho GTPase pathways were activated early during infection while inactivation of the Integrin signalling pathway started between 6 and 12 hpi. Moreover, upstream regulator analysis predicted MYC to be activated with time of infection and protein and RNA data for genes controlled by this transcription factor showed good correlation, which validated the use of protein data for this prediction. Among the identified phosphorylation sites, a group related to glycolysis and cytoskeletal reorganization were up-regulated during infection. The results show specific aspects on how the host cell proteins, the final products in the genetic information flow, are influenced by Ad2 infection, which would be overlooked if only knowledge derived from mRNA data is considered.

MicroRNA transcriptome analysis of porcine vital organ responses to immunosuppressive porcine cytomegalovirus infection

Background

Porcine cytomegalovirus (PCMV) is an immunosuppressive virus that mainly inhibits T-lymphocyte and macrophage immune functions; it has significantly damaged the farming industry. Although recent studies have shown that miRNAs play important roles in immune responses, the regulatory mechanisms of miRNAs during immunosuppressive virus infection remain unclear.

Methods

In this study, porcine small-RNA transcriptomes of PCMV-infected and uninfected vital organs were first characterised by high-throughput sequencing. miRDeep2 software was used to predict novel pig-encoded miRNAs. To verify the accuracy of the high-throughput sequencing results, stem-loop qRT-PCR was performed on 12 significantly DE miRNAs.

The physical and functional interactions between the immune-related target genes of the DE miRNAs in PCMV-infected organs were analysed using the STRING database.

Results

In total, 306 annotated and 295 novel miRNAs were identified from PCMV-infected and uninfected porcine organs, respectively, through alignment with known Sus scrofa pre-miRNAs. Overall, 92, 107, 95, 77 and 111 miRNAs were significantly differentially expressed in lung, liver, spleen, kidney and thymus after PCMV infection, respectively.

According to Gene Ontology enrichment analysis, target genes of the differentially expressed miRNAs associated with immune system processes, regulation of biological processes and metabolic processes were enriched in every sample. Integrated expression analysis of the differentially expressed miRNAs and their target mRNAs in PCMV-infected thymus showed that the significant differential expression of specific miRNAs under the pressure of PCMV infection in central immune organs interfered with the expression of genes involved in important immune-related signalling pathways, thus promoting the viral infection.

Conclusions

This is the first comprehensive analysis of the responses of host small-RNA transcriptomes to PCMV infection in vital porcine organs. It provides new insights into the regulatory mechanisms of miRNAs during infection by immunosuppressive viruses.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-0922-x) contains supplementary material, which is available to authorized users.

Temporal characterization of the non-structural Adenovirus type 2 proteome and phosphoproteome using high-resolving mass spectrometry

The proteome and phosphoproteome of non-structural proteins of Adenovirus type 2 (Ad2) were time resolved using a developed mass spectrometry approach. These proteins are expressed by the viral genome and important for the infection process, but not part of the virus particle. We unambiguously confirm the existence of 95% of the viral proteins predicted to be encoded by the viral genome. Most non-structural proteins peaked in expression at late time post infection. We identified 27 non-redundant sites of phosphorylation on seven different non-structural proteins. The most heavily phosphorylated protein was the DNA binding protein (DBP) with 15 different sites. The phosphorylation occupancy rate could be calculated and monitored with time post infection for 15 phosphorylated sites on various proteins. In the DBP, phosphorylations with time-dependent relation were observed. The findings show the complexity of the Ad2 non-structural proteins and opens up a discussion for potential new drug targets.

MicroRNAs as Important Players in Host-Adenovirus Interactions

MicroRNAs (miRNAs) are powerful regulators of gene expression and fine-tuning genes in all tissues. Cellular miRNAs can control 100s of biologic processes (e.g., morphogenesis of embryonic structures, differentiation of tissue-specific cells, and metabolic control in specific cell types) and have been involved in the regulation of nearly all cellular pathways. Inherently to their involvement in different physiologic processes, miRNAs deregulation has been associated with several diseases. Moreover, several viruses have been described as either, avoid and block cellular miRNAs or synthesize their own miRNA to facilitate infection and pathogenesis. Adenoviruses genome encodes two non-coding RNAs, known as viral-associated (VA) RNA_I and VA RNA_II, which seem to play an important role either by blocking important proteins from miRNA pathway, such as Exportin-5 and Dicer, or by targeting relevant cellular factors. Drastic changes in cellular miRNA expression profile are also noticeable and several cellular functions are affected by these changes. This review focuses on the mechanisms underlying the biogenesis and molecular interactions of miRNAs providing basic concepts of their functions as well as in the interplay between miRNAs and human adenoviruses.

Implications of MicroRNAs in Oncolytic Virotherapy

MicroRNAs (miRNAs) are an abundant class of small non-coding RNA molecules (~22 nt) that can repress gene expression. Deregulation of certain miRNAs is widely recognized as a robust biomarker for many neoplasms, as well as an important player in tumorigenesis and the establishment of tumoral microenvironments. The downregulation of specific miRNAs in tumors has been exploited as a mechanism to provide selectivity to oncolytic viruses or gene-based therapies. miRNA response elements recognizing miRNAs expressed in specific tissues, but downregulated in tumors, have been inserted into the 3′UTR of viral genes to promote the degradation of these viral mRNAs in healthy tissue, but not in tumor cells. Consequently, oncolytic virotherapy-associated toxicities were diminished, while therapeutic activity in tumor cells was preserved. However, viral infections themselves can modulate the miRNome of the host cell, and such miRNA changes under infection impact the normal viral lifecycle. Thus, there is a miRNA-mediated interplay between virus and host cell, affecting both viral and cellular activities. Moreover, the outcome of such interactions may be cell type or condition specific, suggesting that the impact on normal and tumoral cells may differ. Here, we provide an insight into the latest developments in miRNA-based viral engineering for cancer therapy, following the most recent discoveries in miRNA biology. Furthermore, we report on the relevance of miRNAs in virus–host cell interaction, and how such knowledge can be exploited to improve the control of viral activity in tumor cells.

Normal human cell proteins that interact with the adenovirus type 5 E1B 55kDa protein

Several of the functions of the human adenovirus type 5 E1B 55kDa protein are fulfilled via the virus-specific E3 ubiquitin ligase it forms with the viral E4 Orf6 protein and several cellular proteins. Important substrates of this enzyme have not been identified, and other functions, including repression of transcription of interferon-sensitive genes, do not require the ligase. We therefore used immunoaffinity purification and liquid chromatography-mass spectrometry of lysates of normal human cells infected in parallel with HAdV-C5 and E1B 55kDa protein-null mutant viruses to identify specifically E1B 55kDa-associated proteins. The resulting set of >90 E1B-associated proteins contained the great majority identified previously, and was enriched for those associated with the ubiquitin-proteasome system, RNA metabolism and the cell cycle. We also report very severe inhibition of viral genome replication when cells were exposed to both specific or non-specific siRNAs and interferon prior to infection.

Posttranscriptional Regulation in Adenovirus Infected Cells

A deeper understanding of how viruses reprogram their hosts for production of progeny is needed to combat infections. Most knowledge on the regulation of cellular gene expression during adenovirus infection is derived from mRNA studies. Here, we investigated the changes in protein expression during the late phase of adenovirus type 2 (Ad2) infection of the IMR-90 cell line by stable isotope labeling in cell culture with subsequent liquid chromatography-high resolution tandem mass spectrometric analysis. Two biological replicates of samples collected at 24 and 36 h post-infection (hpi) were investigated using swapped labeling. In total, 2648 and 2394 proteins were quantified at 24 and 36 hpi, respectively. Among them, 659 and 645 were deregulated >1.6-fold at the two time points. The protein expression was compared with RNA expression using cDNA sequencing data. The correlation was surprisingly low (r = 0.3), and several examples of posttranscriptional regulation were observed; e.g., proteins related to carbohydrate metabolism were up-regulated at the protein level but unchanged at the RNA level, whereas histone proteins were down-regulated at the protein level but up-regulated at the RNA level. The deregulation of cellular gene expression by adenovirus is mediated at multiple levels and more complex than hitherto believed.

13 C-metabolic flux analysis of human adenovirus infection: Implications for viral vector production

Adenoviruses are human pathogens increasingly used as gene therapy and vaccination vectors. However, their impact on cell metabolism is poorly characterized. We performed carbon labeling experiments with [1,2-¹³ C]glucose or [U-¹³ C]glutamine to evaluate metabolic alterations in the amniocyte-derived, E1-transformed 1G3 cell line during production of a human adenovirus type 5 vector (AdV5). Nonstationary ¹³ C-metabolic flux analysis revealed increased fluxes of glycolysis (17%) and markedly PPP (over fourfold) and cytosolic AcCoA formation (nearly twofold) following infection of growing cells. Interestingly, infection of growth-arrested cells increased overall carbon flow even more, including glutamine anaplerosis and TCA cycle activity (both over 1.5-fold), but was unable to stimulate the PPP and was associated with a steep drop in AdV5 replication (almost 80%). Our results underscore the importance of nucleic and fatty acid biosynthesis for adenovirus replication. Overall, we portray a metabolic blueprint of human adenovirus infection, highlighting similarities with other viruses and cancer, and suggest strategies to improve AdV5 production. Biotechnol. Bioeng. 2017;114: 195-207. © 2016 Wiley Periodicals, Inc.

Distinct temporal changes in host cell lncRNA expression during the course of an adenovirus infection

The deregulation of cellular long non-coding RNA (lncRNA) expression during a human adenovirus infection was studied by deep sequencing. Expression of lncRNAs increased substantially following the progression of the infection. Among 645 significantly expressed lncRNAs, the expression of 398 was changed more than 2-fold. More than 80% of them were up-regulated and 80% of them were detected during the late phase. Based on the genomic locations of the deregulated lncRNAs in relation to known mRNAs and miRNAs, they were predicted to be involved in growth, structure, apoptosis and wound healing in the early phase, cell proliferation in the intermediate phase and protein synthesis, modification and transport in the late phase. The most significant functions of cellular RNA-binding proteins, previously shown to interact with the deregulated lncRNAs identified here, are involved in RNA splicing, nuclear export and translation events. We hypothesize that adenoviruses exploit the lncRNA network to optimize their reproduction.

•

The expression of 398 lncRNAs showed a distinct temporal pattern during Ad2 infection.

•

80% of the deregulated lncRNAs were up-regulated during the late phase of infection.

•

The deregulated lncRNAs potentiallyinteract with 33 cellular RNA binding proteins.

•

These RBPs are involved in RNA splicing, nuclear export and translation.

•

Adenovirus exploits the cellular lncRNA network to optimize its replication.

Expression

Tight control of gene expression is highly characteristic of virus replication, although the mechanisms used to achieve this vary depending on the host cell. There are considerable differences in patterns of gene expression in prokaryotic and eukaryotic cells. Control of gene expression is also achieved at a posttranscriptional level by control over translation of proteins using a number of different mechanisms. Finally, epigenetics plays a significant role in controlling virus gene expression, just as it does for cellular genes. The net results of all these interlocking processes are a complex web of interactions between the virus genome and the machinery of the host cell.

Comparison of the Life Cycles of Genetically Distant Species C and Species D Human Adenoviruses Ad6 and Ad26 in Human Cells

Our understanding of adenovirus (Ad) biology is largely extrapolated from human species C Ad5. Most humans are immune to Ad5, so lower-seroprevalence viruses like human Ad6 and Ad26 are being tested as therapeutic vectors. Ad6 and Ad26 differ at the DNA level by 34%. To better understand how this might impact their biology, we examined the life cycle of the two viruses in human lung cells in vitro. Both viruses infected A549 cells with similar efficiencies, executed DNA replication with identical kinetics within 12 h, and began killing cells within 72 h. While Ad6-infected cells remained adherent until death, Ad26-infected cells detached within 12 h of infection but remained viable. Next-generation sequencing (NGS) of mRNA from infected cells demonstrated that viral transcripts constituted 1% of cellular mRNAs within 6 h and 8 to 16% within 12 h. Quantitative PCR and NGS revealed the activation of key early genes at 6 h and transition to late gene activation by 12 h by both viruses. There were marked differences in the balance of E1A and E1B activation by the two viruses and in the expression of E3 immune evasion mRNAs. Ad6 was markedly more effective at suppressing major histocompatibility complex class I (MHC I) display on the cell surface and in evading TRAIL-mediated apoptosis than was Ad26. These data demonstrate shared as well as divergent life cycles in these genetically distant human adenoviruses. An understanding of these differences expands the knowledge of alternative Ad species and may inform the selection of related Ads for therapeutic development.

IMPORTANCE

A burgeoning number of adenoviruses (Ads) are being harnessed as therapeutics, yet the biology of these viruses is generally extrapolated from Ad2 and Ad5. Here, we are the first to compare the transcriptional programs of two genetically distant Ads by mRNA next-generation sequencing (NGS). Species C Ad6 and Ad26 are being pursued as lower-seroprevalence Ad vectors but differ at the DNA level by 34%. Head-to-head comparison in human lung cells by NGS revealed that the two viruses generally conform to our general understanding of the Ad transcriptional program. However, fine mapping revealed subtle and strong differences in how these two viruses execute these programs, including differences in the balance of E1A and E1B mRNAs and in E3 immune evasion genes. This suggests that not all adenoviruses behave like Ad2 and Ad5 and that they may have unique strategies to infect cells and evade the immune system.

Transcriptome sequencing and development of an expression microarray platform for liver infection in adenovirus type 5-infected Syrian golden hamsters

The Syrian golden hamster is an attractive animal for research on infectious diseases and other diseases. We report here the sequencing, assembly, and annotation of the Syrian hamster transcriptome. We include transcripts from ten pooled tissues from a naïve hamster and one stimulated with lipopolysaccharide. Our data set identified 42,707 non-redundant transcripts, representing 34,191 unique genes. Based on the transcriptome data, we generated a custom microarray and used this new platform to investigate the transcriptional response in the Syrian hamster liver following intravenous adenovirus type 5 (Ad5) infection. We found that Ad5 infection caused a massive change in regulation of liver transcripts, with robust up-regulation of genes involved in the antiviral response, indicating that the innate immune response functions in the host defense against Ad5 infection of the liver. The data and novel platforms developed in this study will facilitate further development of this important animal model.

•

Syrian hamster transcriptome; 42,707 transcripts representing 34,191 unique genes

•

Syrian hamster custom expression microarray platform

•

Ad5 intravenous infection of the Syrian hamster liver

•

Ad5 upregulation of hamster liver genes involved in innate antiviral response.

Fluctuating expression of microRNAs in adenovirus infected cells

The changes in cellular microRNA (miRNA) expression during the course of an adenovirus type 2 infection in human lung fibroblast were studied by deep RNA sequencing. Expressions of 175 miRNAs with over 100 transcripts per million nucleotides were changed more than 1.5-fold. The expression patterns of these miRNAs changed dramatically during the course of the infection, from upregulation of the miRNAs known as tumor suppressors (such as miR-22, miR-320, let-7, miR-181b, and miR-155) and down-regulation of oncogenic miRNAs (such as miR-21 and miR-31) early to downregulation of tumor suppressor miRNAs (such as let-7 family, mir-30 family, 23/27 cluster) and upregulation of oncogenic miRNAs (include miR-125, miR-27, miR-191) late after infection. The switch in miRNA expression pattern occurred when adenovirus DNA replication started. Furthermore, deregulation of cellular miRNA expression was a step-wise and special sets of miRNAs were deregulated in different phases of infection.

An activating transcription factor of Litopenaeus vannamei involved in WSSV genes Wsv059 and Wsv166 regulation

Members of activating transcription factor/cyclic adenosine 3', 5'-monophosphate response element binding protein (ATF/CREB) family are induced by various stress signals and function as effector molecules. Consequently, cellular changes occur in response to discrete sets of instructions. In this work, we found an ATF transcription factor in Litopenaeus vannamei designated as LvATFβ. The full-length cDNA of LvATFβ was 1388 bp long with an open reading frame of 939 bp that encoded a putative 313 amino acid protein. The protein contained a basic region-leucine zipper (bZip) domain that was a common feature among ATF/CREB transcription factors. LvATFβ was highly expressed in intestines, gills, and heart. LvATFβ expression was dramatically upregulated by white spot syndrome virus (WSSV) infection. Pull-down assay revealed that LvATFβ had strong affinity to promoters of WSSV genes, namely, wsv059 and wsv166. Dual-luciferase reporter assay showed that LvATFβ could upregulate the expression of wsv059 and wsv166. Knocked down LvATFβ resulted in decreased expression of wsv059 and wsv166 in WSSV-challenged L. vannamei. Knocked down expression of wsv059 and wsv166 by RNA interference inhibited the replication and reduce the mortality of L. vannamei during WSSV challenge inoculation. The copy numbers of WSSV in wsv059 and wsv166 knocked down group were significant lower than in the control. These results suggested that LvATFβ may be involved in WSSV replication by regulating the expression of wsv059 and wsv166.

Transcriptome analysis of porcine thymus following porcine cytomegalovirus infection

Porcine cytomegalovirus (PCMV) is a major immunosuppressive virus that mainly affects the immune function of T lymphocytes and macrophages. Despite being widely distributed around the world, no significantly different PCMV serotypes have been found. Moreover, the molecular immunosuppressive mechanisms of PCMV, along with the host antiviral mechanisms, are still not well characterized. To understand the potential impact of PCMV on the function of immune organs, we examined the transcriptome of PCMV-infected thymuses by microarray analysis. We identified 5,582 genes that were differentially expressed as a result of PCMV infection. Of these, 2,161 were upregulated and 3,421 were downregulated compared with the uninfected group. We confirmed the expression of 13 differentially expressed immune-related genes using quantitative real-time RT-PCR, and further confirmed the expression of six of those cytokines by western blot. Gene ontology, gene interaction networks, and KEGG pathway analysis of our results indicated that PCMV regulates multiple functional pathways, including the immune system, cellular and metabolic processes, networks of cytokine-cytokine receptor interactions, the TGF-β signaling pathway, the lymphocyte receptor signaling pathway, and the TNF-α signaling pathway. Our study is the first comprehensive attempt to explore the host transcriptional response to PCMV infection in the porcine immune system. It provides new insights into the immunosuppressive molecular mechanisms and pathogenesis of PCMV. This previously unrecognized endogenous antiviral mechanism has implications for the development of host-directed strategies for the prevention and treatment of immunosuppressive viral diseases.

Adenovirus small E1A employs the lysine acetylases p300/CBP and tumor suppressor Rb to repress select host genes and promote productive virus infection

Oncogenic transformation by adenovirus small e1a depends on simultaneous interactions with the host lysine acetylases p300/CBP and the tumor suppressor RB. How these interactions influence cellular gene expression remains unclear. We find that e1a displaces RBs from E2F transcription factors and promotes p300 acetylation of RB1 K873/K874 to lock it into a repressing conformation that interacts with repressive chromatin-modifying enzymes. These repressing p300-e1a-RB1 complexes specifically interact with host genes that have unusually high p300 association within the gene body. The TGF-β, TNF-, and interleukin-signaling pathway components are enriched among such p300-targeted genes. The p300-e1a-RB1 complex condenses chromatin in a manner dependent on HDAC activity, p300 lysine acetylase activity, the p300 bromodomain, and RB K873/K874 and e1a K239 acetylation to repress host genes that would otherwise inhibit productive virus infection. Thus, adenovirus employs e1a to repress host genes that interfere with viral replication.

Unabated adenovirus replication following activation of the cGAS/STING-dependent antiviral response in human cells

The cGAS/STING DNA sensing complex has recently been established as a predominant pathogen recognition receptor (PRR) for DNA-directed type I interferon (IFN) innate immune activation. Using replication-defective adenovirus vectors and replication-competent wild-type adenovirus, we have modeled the influence of the cGAS/STING cascade in permissive human cell lines (A549, HeLa, ARPE19, and THP1). Wild-type adenovirus induced efficient early activation of the cGAS/STING cascade in a cell-specific manner. In all responsive cell lines, cGAS/STING short hairpin RNA (shRNA) knockdown resulted in a loss of TBK1 and interferon response factor 3 (IRF3) activation, a lack of beta interferon transcript induction, loss of interferon-dependent STAT1 activation, and diminished induction of interferon-stimulated genes (ISGs). Adenoviruses that infect through the coxsackievirus-adenovirus receptor (CAR) (Ad2 and Ad5) and the CD46 (Ad35) and desmoglein-2 (Ad7) viral receptors all induce the cGAS/STING/TBK1/IRF3 cascade. The magnitude of the IRF3/IFN/ISG antiviral response was strongly influenced by serotype, with Ad35>Ad7>Ad2. For each serotype, no enhancement of viral DNA replication or virus production occurred in cGAS or STING shRNA-targeted cell line pools. We found no replication advantage in permissive cell lines that do not trigger the cGAS/STING cascade following infection. The cGAS/STING/TBK1/IRF3 cascade was not a direct target of viral antihost strategies, and we found no evidence that Ad stimulation of the cGAS/STING DNA response had an impact on viral replication efficiency.

IMPORTANCE

This study shows for the first time that the cGAS DNA sensor directs a dominant IRF3/IFN/ISG antiviral response to adenovirus in human cell lines. Activation of cGAS occurs with viruses that infect through different high-affinity receptors (CAR, CD46, and desmoglein-2), and the magnitude of the cGAS/STING DNA response cascade is influenced by serotype-specific functions. Furthermore, activation of the cGAS cascade occurred in a cell-specific manner. Activation of the cGAS/STING response did not impact viral replication, and viral immune evasion strategies did not target the cGAS/STING/TBK1/IRF3 cascade. These studies provide novel insight into the early innate recognition response to adenovirus.

Identification of adenovirus-encoded small RNAs by deep RNA sequencing

Using deep RNA sequencing, we have studied the expression of adenovirus-encoded small RNAs at different times after infection. Nineteen small RNAs which comprised more than 1% of the total pool of small RNAs at least one time point were identified. These small RNAs were between 25 and 35 nucleotides long and mapped in the region of the VA RNAI and RNAII genes. However, the overlap was incomplete and some contained a few extra nucleotides at the 3' end. This finding together with the observation that some of the small RNAs were detected before VA RNA expression had started might indicate that they are derived from other precursors than VA RNAI and II. Interestingly, the small RNAs displayed different expression profiles during the course of the infection suggesting that they have different functions. An effort was made to identify their mRNA targets by using computer prediction and deep cDNA sequencing. The most significant targets for the earliest small RNAs were genes involved in signaling pathways.

[Adenovirus infections in immunocompromised patients]

Human adenoviruses function as genetic models and vectors for gene therapy. Upper respiratory, gastrointestinal or ocular infections usually have mild course without any major complication in immunocompetent individuals. However, reactivation from latency in immunocompromised patients may lead to death. Depending on the underlying diseases, different adenovirus serotypes damage different organs. In children with severe combined immunodeficiency syndrome, serotypes of species A and C induce lung, liver or bladder inflammation. Paediatric hematopoietic stem cell transplantation is frequently followed by serotype 31-induced pneumonia, enteritis, cystitis. B serotypes can destroy transplanted organs. In AIDS patients, D and novel F serotypes cause enteritis. Recombinants of B serotypes induce urinary tract infections. Progression of lymphomas, tumours, and systemic lupus erythematosus might be facilitated by immunosuppressive effects of adenoviruses. As far as the diagnostic work-up of adenoviruses, detection of viral DNA and virus copy number is predictive, while serology testing is quite unreliable. For treatment, cidofovir derivates, ribavirin, ganciclovir, vidarabine and microRNA have been used.

De novo derivation of proteomes from transcriptomes for transcript and protein identification

Identification of proteins by tandem mass spectrometry requires a reference protein database, but these are only available for model species. Here we demonstrate that, for a non-model species, the sequencing of expressed mRNA can generate a protein database for mass spectrometry-based identification. This combination of high-throughput sequencing and protein identification technologies allows detection of genes and proteins. We use human cells infected with human adenovirus as a complex and dynamic model to demonstrate the robustness of this approach. Our proteomics informed by transcriptomics (PIT) technique identifies >99% of over 3,700 distinct proteins identified using traditional analysis that relies on comprehensive human and adenovirus protein lists. We show that this approach can also be used to highlight genes and proteins undergoing dynamic changes in post-transcriptional protein stability.

Arenavirus infection induces discrete cytosolic structures for RNA replication

Arenaviruses are responsible for acute hemorrhagic fevers with high mortality and pose significant threats to public health and biodefense. These enveloped negative-sense RNA viruses replicate in the cell cytoplasm and express four proteins. To better understand how these proteins insinuate themselves into cellular processes to orchestrate productive viral replication, we have identified and characterized novel cytosolic structures involved in arenavirus replication and transcription. In cells infected with the nonpathogenic Tacaribe virus or the attenuated Candid#1 strain of Junín virus, we find that newly synthesized viral RNAs localize to cytosolic puncta containing the nucleoprotein (N) of the virus. Density gradient centrifugation studies reveal that these replication-transcription complexes (RTCs) are associated with cellular membranes and contain full-length genomic- and antigenomic-sense RNAs. Viral mRNAs segregate at a higher buoyant density and are likewise scant in immunopurified RTCs, consistent with their translation on bulk cellular ribosomes. In addition, confocal microscopy analysis reveals that RTCs contain the lipid phosphatidylinositol-4-phosphate and proteins involved in cellular mRNA metabolism, including the large and small ribosomal subunit proteins L10a and S6, the stress granule protein G3BP1, and a subset of translation initiation factors. Elucidating the structure and function of RTCs will enhance our understanding of virus-cell interactions that promote arenavirus replication and mitigate against host cell immunity. This knowledge may lead to novel intervention strategies to limit viral virulence and pathogenesis.

Pervasive initiation and 3'-end formation of poxvirus postreplicative RNAs

Poxviruses are large DNA viruses that replicate within the cytoplasm and encode a complete transcription system, including a multisubunit RNA polymerase, stage-specific transcription factors, capping and methylating enzymes, and a poly(A) polymerase. Expression of the more than 200 open reading frames by vaccinia virus, the prototype poxvirus, is temporally regulated: early mRNAs are synthesized immediately after infection, whereas intermediate and late mRNAs are synthesized following genome replication. The postreplicative transcripts are heterogeneous in length and overlap the entire genome, which pose obstacles for high resolution mapping. We used tag-based methods in conjunction with high throughput cDNA sequencing to determine the precise 5'-capped and 3'-polyadenylated ends of postreplicative RNAs. Polymerase slippage during initiation of intermediate and late RNA synthesis results in a 5'-poly(A) leader that allowed the unambiguous identification of true transcription start sites. Ninety RNA start sites were located just upstream of intermediate and late open reading frames, but many more appeared anomalous, occurring within coding and non-coding regions, indicating pervasive transcription initiation. We confirmed the presence of functional promoter sequences upstream of representative anomalous start sites and demonstrated that alternative start sites within open reading frames could generate truncated isoforms of proteins. In an analogous manner, poly(A) sequences allowed accurate mapping of the numerous 3'-ends of postreplicative RNAs, which were preceded by a pyrimidine-rich sequence in the DNA coding strand. The distribution of postreplicative promoter sequences throughout the genome provides enormous transcriptional complexity, and the large number of previously unmapped RNAs may have novel functions.

Next generation deep sequencing and vaccine design: today and tomorrow

Next generation sequencing (NGS) technologies have redefined the modus operandi in both human and microbial genetics research, allowing the unprecedented generation of very large sequencing datasets on a short time scale and at affordable costs. Vaccine development research is rapidly taking full advantage of the advent of NGS. This review provides a concise summary of the current applications of NGS in relation to research seeking to develop vaccines for human infectious diseases, incorporating studies of both the pathogen and the host. We focus on rapidly mutating viral pathogens, which are major targets in current vaccine research. NGS is unraveling the complex dynamics of viral evolution and host responses against these viruses, thus contributing substantially to the likelihood of successful vaccine development.