The Human Adenovirus Type 2 Transcriptome: An Amazing Complexity of Alternatively Spliced mRNAs.

We have used the Nanopore long-read sequencing platform to demonstrate how amazingly complex the human adenovirus type 2 (Ad2) transcriptome is with a flexible splicing machinery producing a range of novel mRNAs both from the early and late transcription units. In total we report more than 900 alternatively spliced mRNAs produced from the Ad2 transcriptome whereof more than 850 are novel mRNAs. A surprising finding was that more than 50% of all E1A transcripts extended upstream of the previously defined transcriptional start site. The novel start sites mapped close to the inverted terminal repeat (ITR) and within the E1A enhancer region. We speculate that novel promoters or enhancer driven transcription, so-called eRNA transcription, is responsible for producing these novel mRNAs. Their existence was verified by a peptide in the Ad2 proteome that was unique for the E1A ITR mRNA. Although we show a high complexity of alternative splicing from most early and late regions, the E3 region was by far the most complex when expressed at late times of infection. More than 400 alternatively spliced mRNAs were observed in this region alone. These mRNAs included extended L4 mRNAs containing E3 and L5 sequences and readthrough mRNAs combining E3 and L5 sequences. Our findings demonstrate that the virus has a remarkable capacity to produce novel exon combinations, which will offer the virus an evolutionary advantage to change the gene expression repertoire and protein production in an evolving environment.IMPORTANCE Work in the adenovirus system led to the groundbreaking discovery of RNA splicing and alternative RNA splicing in 1977. These mechanisms are essential in mammalian evolution by increasing the coding capacity of a genome. Here, we have used a long-read sequencing technology to characterize the complexity of human adenovirus pre-mRNA splicing in detail. It is mindboggling that the viral genome, which only houses around 36,000 bp, not being much larger than a single cellular gene, generates more than 900 alternatively spliced mRNAs. Recently, adenoviruses have been used as the backbone in several promising SARS-CoV-2 vaccines. Further improvement of adenovirus-based vaccines demands that the virus can be tamed into an innocent carrier of foreign genes. This requires a full understanding of the components that govern adenovirus replication and gene expression.

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 (TiO2 ) 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.

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.

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.

Tyrosine phosphorylation profiling via in situ proximity ligation assay.

BACKGROUND: Tyrosine phosphorylation (pTyr) is an important cancer relevant posttranslational modification since it regulates protein activity and cellular localization. By controlling cell growth and differentiation it plays an important role in tumor development. This paper describes a novel approach for detection and visualization of a panel of pTyr proteins in tumors using in situ proximity ligation assay. METHODS: K562 leukemia cells were treated with tyrosine kinase and/or phosphatase inhibitors to induce differences in pTyr levels and mimic cells with different malignant properties. Cells were then probed with one antibody against the pTyr modification and another probe against the detected protein, resulting in a detectable fluorescent signal once the probes were in proximity. RESULTS: Total and protein specific pTyr levels on ABL, SHC, ERK2 and PI3K proteins were detected and samples of control and treated cells were distinguished at the pTyr level using this novel approach. Promising results were also detected for formalin fixed and paraffin embedded cells in the micro array format. CONCLUSIONS: This application of in situ proximity ligation assay is valuable in order to study the pTyr modification of a panel of proteins in large data sets to validate mass spectrometric data and to be combined with tissue microarrays. The approach offers new opportunities to reveal the pTyr signatures in cells of different malignant properties that can be used as biomarker of disease in the future.

Tumour expression of bladder cancer-associated urinary proteins.

UNLABELLED: WHAT'S KNOWN ON THE SUBJECT? AND WHAT DOES THE STUDY ADD?: The current basis for diagnosis and prognosis in urinary bladder cancer is based on the pathologists' assessment of a biopsy of the tumour. Urinary biomarkers are preferable as they can be non-invasively sampled. Urinary cytology is the only test with widespread use but is hampered by poor reproducibility and low sensitivity. By studying the protein expression in bladder tumour tissue samples of proteins previously found in elevated levels in the urine of patients with bladder cancer, we have been able to show that these proteins originate from the tumour. The immunoreactivity of three of the investigated proteins increased with higher stage. Also a serine peptidase inhibitor was found to be predictive of progression from non-muscle-invasive to muscle-invasive tumours. OBJECTIVES: To analyse the expression of five bladder cancer-associated urinary proteins and investigate if expression is related to the malignant phenotype of the tumour. To explore the possible prognostic value of these proteins. PATIENTS AND METHODS: Urine samples, 16 from patients with bladder cancer and 26 from controls, were used in Western Blotting experiments. Tissue microarrays with bladder tissue from 344 patients diagnosed with bladder cancer between 1984 and 2005 was used in immunohistochemistry experiments. The proteins apolipoprotein E (APOE), fibrinogen ß chain precursor (FGB), leucine-rich α2-glycoprotein (LRG1), polymerase (RNA) I polypeptide E (POLR1E), α1-antitrypsin (SERPINA1) and topoisomerase 2A (TOP2A) were probed with antibodies validated by the Human Protein Atlas. RESULTS: Increased expressions of APOE, FGB and POLR1E were correlated with increased tumour stage (P < 0.001). Expression of SERPINA1 in Ta and T1 tumours was found to increase the risk of tumour progression (hazard ratio 2.57, 95% confidence interval 1.13-5.87; P = 0.025) CONCLUSIONS: All proteins previously detected in urine from patients with bladder cancer were also expressed in bladder cancer tissue. The expression of APOE, FGB and POLR1E increased with stage and they are potential diagnostic markers. SERPINA1 was identified as a prognostic marker candidate.

A strategy for identification of protein tyrosine phosphorylation.

To develop methods for studying phosphorylation of protein tyrosine residues is an important task since this protein modification regulates many cellular functions and often is involved in oncogenesis. An optimal protocol includes enrichment of tyrosine phosphorylated (pTyr) peptides or proteins, followed by a high resolving analytical method for identification of the enriched components. In this Methods paper, we describe a working strategy on how immunoaffinity enrichments, using anti-pTyr antibodies, combined with mass spectrometric (MS) analysis can be used to study the pTyr proteome. We describe in detail how our procedure was used to characterize the pTyr proteome of K562 leukemia cells. Important questions concerning the use of different anti-pTyr antibodies, enrichments performed at the peptide and/or the protein level, pooling of enrichments and requirements for the MS characterization are discussed.

Proteomic analysis of urinary biomarker candidates for nonmuscle invasive bladder cancer.

Nonmuscle invasive tumors of the bladder often recur and thereby bladder cancer patients need regular re-examinations which are invasive, unpleasant, and expensive. A noninvasive and less expensive method, e.g. a urine dipstick test, for monitoring recurrence would thus be advantageous. In this study, the complementary techniques mass spectrometry (MS) and Western blotting (WB)/dot blot (DB) were used to screen the urine samples from bladder cancer patients. High resolving MS was used to analyze and quantify the urinary proteome and 29 proteins had a significantly higher abundance (p<0.05) in bladder cancer samples compared with control urine samples. The increased abundance found in urine from bladder cancer patients compared with controls was confirmed with Western blot for four selected proteins; fibrinogen ß chain precursor, apolipoprotein E, α-1-antitrypsin, and leucine-rich α-2-glycoprotein 1. Dot blot analysis of an independent urine sample set pointed out fibrinogen ß chain and α-1-antitrypsin as most interesting biomarkers having sensitivity and specificity values in the range of 66-85%. Exploring the Human Protein Atlas (HPA) also revealed that bladder cancer tumors are the likely source of these proteins. They have the potential of being useful in diagnosis, monitoring of recurrence and thus may improve the treatment of bladder tumors, especially nonmuscle invasive tumors.

Optimization of immunoaffinity enrichment and detection: toward a comprehensive characterization of the phosphotyrosine proteome of K562 cells by liquid chromatography-mass spectrometry.

Phosphorylation of protein tyrosine residues regulates many cell functions and has also been proved to be involved in oncogenesis. Thus, the identification of the phosphotyrosine (pTyr) proteome of cells is a very important task. Since tyrosine phosphorylation represents only around 1% of the total human phosphoproteome, the study of pTyr proteins is rather challenging. Here we report the optimization study of the phosphotyrosine proteome using K562 cells as a model system. A substantial segment of the phosphotyrosine proteome of K562 cells was characterized by immunoaffinity enrichment with 4G10 and PYKD1 antibodies followed by LC-MS/MS analysis. 480 non-redundant pTyr peptides corresponding to 342 pTyr proteins were found. 141 pTyr peptides were not described elsewhere. The mass spectrometry approach involving high-resolving FTMS analysis of precursor ions and subsequent detection of CID fragments in a linear ion trap was considered as optimal. For detection of low abundant pTyr peptides pooling of individual immunoaffinity enrichments for one LC-MS/MS analysis was crucial. The enrichment properties of the monoclonal PYKD1 antibody were presented for the first time, also in comparison to the 4G10 antibody. PYKD1 was found to be more effective for protein enrichment (1.2 and 5% efficiency at peptide and protein level correspondingly), while 4G10 showed better results when peptide enrichment was performed (15% efficiency versus 3.6% at protein level). Substantially different subsets of the phosphoproteome were enriched by these antibodies. This finding together with previous studies demonstrates that comprehensive pTyr proteome characterization by immunoprecipitation requires multiple antibodies to be used for the affinity enrichment.

Detection of tyrosine phosphorylated proteins by combination of immunoaffinity enrichment, two-dimensional difference gel electrophoresis and fluorescent Western blotting.

We describe fluorescence-based 2-D gel electrophoresis methods for visualization of low abundant, cancer relevant tyrosine phosphorylated (pTyr) proteins. The methods investigated were fluorescent Western blotting and two-dimensional difference gel electrophoresis (2-D DIGE) for detection of non-enriched and immunoaffinity enriched pTyr protein patterns. The same anti-phosphotyrosine specific antibody, 4G10, was used for both approaches. The results from fluorescent Western blotting of total proteins and from enriched CyDye DIGE pre-labeled pTyr proteins showed similar down regulation of phosphorylation upon treating of cells from a cancer model system (K562 chronic myeloid leukemia cells) with imatinib. This treatment introduced a known perturbation of phosphorylation that enabled testing of these new approaches to analyze variations in tyrosine phosphorylation levels. Enrichment of pTyr proteins was found highly advantageous for the outcome. Out of a simplified 2-D DIGE experiment of immunoaffinity enriched control and treated pTyr proteins, differential analysis as well as protein identification by mass spectrometry (MS) was possible.

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.

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.

Support for schizophrenia susceptibility locus on chromosome 2q detected in a Swedish isolate using a dense map of microsatellites and SNPs.

Extended pedigrees are not only very useful to identify disease genes for rare Mendelian conditions, but they may also help unravel the genetics of complex diseases such as schizophrenia. In this study we performed genome-wide multipoint non-parametric linkage (NPL) score calculations using 825 microsatellites and 5,366 single nucleotide polymorphisms (SNPs), respectively, and searched for haplotypes shared by affected individuals, in three multiplex families including 29 genotyped affected individuals which in total contains 49 relative pairs useful for linkage studies. The most consistent results for microsatellites and SNPs were observed on 2q12.3-q14.1 (NPL scores 2.0, empirical P-value 0.009). However, the overall highest NPL score was observed on chromosome 2q33.3 using SNPs (NPL score 2.2, empirical P-value 0.007). Other chromosomal regions were detected on 5q15-q22.1, with microsatellites (NPL scores 1.7, empirical P-value 0.021) and with SNPs (NPL scores 2.0, empirical P-value 0.010) and on 5q23.1 (NPL score 1.9, empirical P-value 0.012) and 8q24.1-q24.2 (NPL score 2.1, empirical P-value 0.009) when using SNPs. The analysis of extended pedigrees allowed the search for haplotypes inherited identical by decent (IBD) by affected individuals. In all regions with NPL score >1.9 we found haplotypes inherited IBD by multiple cases. However, no common haplotypes were found for affected individuals in all families. In conclusion our NPL results support earlier findings suggesting that 2q and possibly 5q and 8q contain susceptibility loci for schizophrenia. Haplotype sharing in families helped to delimit the detected regions that potentially are susceptibility loci for schizophrenia.

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.

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.

Identification of the adenovirus type 2 C-168 protein.

A hitherto predicted but undetected protein, C-168, in adenovirus type 2 (Ad2) has been identified using mass spectrometry (MS) based proteomics. The gene of this 17.7kDa protein is located on the forward strand in the major late transcription unit between base pairs 9294 and 9797. A tryptic peptide, derived from the C-terminal part of the protein, was identified with high amino acid sequence coverage. A candidate splice site for the corresponding mRNA is also presented. The protein sequence is unusual with repeats of serine, glycine and arginine. A bioinformatics prediction of protein function and localization is presented.

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.

Data on the expression of cellular lncRNAs in human adenovirus infected cells.

Expression of cellular long non-coding RNAs (lncRNAs) in human primary lung fibroblasts (IMR-90) during the course of adenovirus type 2 (Ad2) infection was studied by strand-specific whole transcriptome sequencing. In total, 645 cellular lncRNAs were expressed at a significant level and 398 of them were changed more than 2-fold. The changes in expression followed a distinct temporal pattern. Significantly, 80% of the changes occurred at the late phase and 80% of the de-regulated lncRNAs were up-regulated. The three largest groups of deregulated lncRNAs were 125 antisense RNAs, 111 pseudogenes and 85 long intergenic non-coding RNAs (lincRNAs). Lastly, more than 36% of lncRNAs have been shown to interact with RNA binding proteins.

Impact of the interaction between adenovirus E1A and CtBP on host cell gene expression.

In cell lines harbouring inducible adenovirus E1A genes, the cytotoxicity of wild type E1A was manifested by poor and subsiding expression of the E1A protein during prolonged induction. In contrast, cells expressing E1A deleted in the C-terminal binding protein (CtBP)-interaction domain (E1ADeltaCID) demonstrated high levels of expression for extended time. Microarray analyses of host cell gene expression demonstrated that approximately 70% of the regulated genes were increased upon E1A induction and that the majority of E1A-regulated genes were similarly regulated by wild type E1A and E1ADeltaCID. However, for 29 genes, regulation by wild type E1A and E1ADeltaCID were different. Consistent with the altered transforming capacity of E1A unable to bind CtBP, genes involved in tumour cell progression and growth suppression were found among the differently regulated genes. Moreover, promoter sequences of genes up regulated by wild type E1A and/or repressed by E1ADeltaCID demonstrated a higher prevalence of potential binding sites for the CtBP-targeted transcription factors Ets, Ikaros and/or partial differentialEF1/ZEB, suggesting that the failure to block CtBP-repression contributed to the "hyper-transforming" phenotype of E1ADeltaCID. Since E1ADeltaCID also specifically activated host cell gene expression, we find it likely that additional, possibly CtBP-independent, mechanisms contribute to the altered phenotype of E1ADeltaCID-expressing cells.