An efficient adenovirus infection results in high-level accumulation of viral DNA and mRNAs in the infected cell population. However, the average viral DNA and mRNA content in a heterogeneous cell population does not necessarily reflect the same abundance in individual cells. Here, we describe a novel padlock probe-based rolling-circle amplification technique that enables simultaneous detection and analysis of human adenovirus type 5 (HAdV-5) genomic DNA and virus-encoded mRNAs in individual infected cells. We demonstrate that the method is applicable for detection and quantification of HAdV-5 DNA and mRNAs in short-term infections in human epithelial cells and in long-term infections in human B lymphocytes. Single-cell evaluation of these infections revealed high heterogeneity and unique cell subpopulations defined by differential viral DNA content and mRNA expression. Further, our single-cell analysis shows that the specific expression pattern of viral and mRNA splice variants is linked to HAdV-5 DNA content in the individual cells. Furthermore, we show that expression of a mature form of the HAdV-5 histone-like protein VII affects virus genome detection in HAdV-5-infected cells. Collectively, padlock probes combined with rolling-circle amplification should be a welcome addition to the method repertoire for the characterization of the molecular details of the HAdV life cycle in individual infected cells. Human adenoviruses (HAdVs) have been extensively used as model systems to study various aspects of eukaryotic gene expression and genome organization. The vast majority of the HAdV studies are based on standard experimental procedures carried out using heterogeneous cell populations, where data averaging often masks biological differences. As every cell is unique, characteristics and efficiency of an HAdV infection can vary from cell to cell. Therefore, the analysis of HAdV gene expression and genome organization would benefit from a method that permits analysis of individual infected cells in the heterogeneous cell population. Here, we show that the padlock probe-based rolling-circle amplification method can be used to study concurrent viral DNA accumulation and mRNA expression patterns in individual HAdV-5-infected cells. Hence, this versatile method can be applied to detect the extent of infection and virus gene expression changes in different HAdV-5 infections. read more
We have used high-throughput small RNA sequencing to characterize viral small RNA expression in purified tonsillar B and T lymphocytes isolated from patients tested positive for Epstein-Barr virus (EBV) or human adenovirus (HAdV) infections, respectively. In the small set of patients analyzed, the expression profile of EBV and HAdV miRNAs could not distinguish between patients diagnosed with tonsillar hypertrophy or chronic/recurrent tonsillitis. The EBV miR-BART expression profile among the patients diagnosed with tonsillar diseases resembles most closely the pattern seen in EBV+ tumors (Latency II/I). The miR-BARTs that appear to be absent in normal EBV infected cells are essentially all detectable in the diseased tonsillar B lymphocytes. In the EBV+ B cells we detected 44 EBV miR-BARTs derived from the proposed BART precursor hairpins whereof five are not annotated in miRBase v21. One previously undetected miRNA, BART16b-5p, originates from the miR-BART16 precursor hairpin as an alternative 5´ miR-BART16 located precisely upstream of the annotated miR-BART16-5p. Further, our analysis revealed an extensive sequence variation among the EBV miRNAs with isomiRs having a constant 5´ end but alternative 3´ ends. A range of small RNAs was also detected from the terminal stem of the EBER RNAs and the 3´ part of v-snoRNA1. During a lytic HAdV infection in established cell lines the terminal stem of the viral non-coding VA RNAs are processed to highly abundant viral miRNAs (mivaRNAs). In contrast, mivaRNA expression in HAdV positive tonsillar T lymphocytes was very low. The small RNA profile further showed that the 5´ mivaRNA from VA RNAI and the 3´ mivaRNA from VA RNAII were as predicted, whereas the 3´ mivaRNA from VA RNAI showed an aberrant processing upstream of the expected Dicer cleavage site. read more
