J Gen Virol Email Content Delivery
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 87 (2006), 1961-1975; DOI 10.1099/vir.0.81756-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Versteeg, G. A.
Right arrow Articles by Spaan, W. J. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Versteeg, G. A.
Right arrow Articles by Spaan, W. J. M.
Agricola
Right arrow Articles by Versteeg, G. A.
Right arrow Articles by Spaan, W. J. M.
© 2006 Society for General Microbiology

Transcriptional profiling of acute cytopathic murine hepatitis virus infection in fibroblast-like cells

Gijs A. Versteeg{dagger}, Olga Slobodskaya{dagger} and Willy J. M. Spaan

Molecular Virology Laboratory, Department of Medical Microbiology, Center of Infectious Diseases, Leiden University Medical Center, E4P, PO Box 9600, 2300 RC Leiden, The Netherlands

Correspondence
Willy J. M. Spaan
w.j.m.spaan{at}lumc.nl


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Understanding the orchestrated genome-wide cellular responses is critical for comprehending the early events of coronavirus infection. Microarray analysis was applied to assess changes in cellular expression profiles during different stages of two independent, highly controlled murine hepatitis virus (MHV) infections in vitro. Fibroblast-like L cells were infected at high multiplicity in order to study the direct effects of a synchronized lytic coronavirus infection. Total RNA was harvested from MHV- or mock-infected L cells at 3, 5 and 6 h post-infection and hybridized to Affymetrix microarrays representing approximately 12 500 murine genes and expressed sequences. The expression data were compared to their respective mock-infected controls. Quantitative RT-PCR of selected transcripts was used to validate the differential expression of transcripts and inter-experiment reproducibility of microarray analysis. It was concluded that MHV-A59 infection in fibroblast-like cells triggers very few transcriptional cellular responses in the first 3 h of infection. Later, after having established a productive infection, a chemokine response is induced together with other cellular changes associated with RNA and protein metabolism, cell cycle and apoptosis. Interferon responses are not triggered during infection, although the L cells can be readily stimulated to produce interferon by dsRNA, a known potent inducer of interferon. Possibly, the interferon response is actively counteracted by a virus-encoded antagonist as has been described previously for other RNA viruses.

{dagger}These authors contributed equally to this work. Back


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Murine hepatitis virus (MHV) is one of the best-studied representatives of the family Coronaviridae and may serve as a valuable model for coronaviruses such as SARS-CoV and NL63 that cause human diseases (Fouchier et al., 2003Down; van der Hoek et al., 2004Down). Its genome consists of single-stranded RNA of positive polarity and encompasses more than 31 000 nt. Having such large genomes – and nearly all of it devoted to coding capacity – gives coronaviruses adaptive advantages compared with other viruses with smaller RNA genomes. MHV genome expression starts with the translation of two large open reading frames (ORF), which encode two polyproteins that are subsequently processed by viral proteinases into at least 16 non-structural proteins. These include an RNA-dependent RNA polymerase, three proteinases, a helicase and a number of RNA-processing enzymes (Snijder et al., 2003Down). The other ORFs, which are located at the 3'-third of the genome, are expressed from subgenomic RNAs. Most of these ORFs encode proteins that constitute the viral particles. The function of the majority of the viral proteins and their interaction with the host cell have not yet been fully characterized.

Coronavirus infection of established cell lines – including fibroblast cell lines – has been extensively used as a model for studying viral genome expression and virus–host interaction. During lytic MHV infection, levels of virus-encoded intracellular mRNAs rapidly increase and they actively compete for ribosomes with cellular mRNAs (Hilton et al., 1986Down). In this way MHV downregulates host cell translation, without having any pronounced effect on cellular transcription (Siddell et al., 1981Down). Like other positive-stranded RNA viruses, MHV uses cellular membranes as sites of viral RNA replication (van der Meer et al., 1999Down; Shi et al., 1999Down). Viral infection triggers the formation of double-membrane vesicles, which originate from intracellular membranes (Gosert et al., 2002Down; Prentice et al., 2004Down; Shi et al., 1999Down; van der Meer et al., 1999Down). Cell surface expression of the viral spike glycoprotein results in cell-to-cell fusion (Vennema et al., 1990Down), which is accompanied by fragmentation and rearrangement of the Golgi apparatus (Lavi et al., 1996Down). Accumulation of intracellular complexes of spike and the MHV receptor Ceacam1 has been implemented in the development of a cytopathic effect (CPE) (Rao & Gallagher, 1998Down). Profound effects of MHV infection on cellular functions and structure eventually result in cell death.

In response to viral infection, cells produce cytokines and chemokines (Li et al., 2004Down; Rempel et al., 2005Down), which play a key role in modulating innate as well as adaptive immune responses. In addition to an inflammatory response, other changes of host cell expression during MHV infection have been described previously (Cai et al., 2003Down; Kyuwa et al., 1994Down; Ning et al., 2003aDown, bDown); yet the understanding of MHV–host interactions on the cellular level has been limited.

Here, we applied microarray analysis to assess global changes in cellular expression profiles during acute cytopathic coronavirus infection. We demonstrate that MHV infection in fibroblast-like cells triggers very few changes in cellular mRNA levels in the first hours of infection. After having established exponential virus production, innate immune and inflammatory responses were elicited together with other cellular changes associated with RNA and protein metabolism, oxidative stress, cell cycle and apoptosis. The analysis of transcriptional changes in cells undergoing lytic MHV infection provides new clues into coronavirus–host interaction.


   METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Cells, virus, plaque assay.
Mouse Sac(–) and L cell lines (Jacobs et al., 1981Down) were grown as monolayers in Dulbecco's modified Eagle's medium supplemented with 100 U penicillin ml–1, 100 µg streptomycin ml–1 and 10 % fetal calf serum (FCS) (growth medium). After infection or during poly (I : C) treatment, cells were incubated in the same medium with antibiotics but with 3 % FCS (supporting medium). Strain A59 of MHV was obtained from ATCC. The working stock of the virus was prepared using Sac(–) cells as described previously (Spaan et al., 1981Down). Viral titres in p.f.u. ml–1 were determined by plaque assay on L cells (Spaan et al., 1981Down).

Viral infection and RNA isolation.
L cells were seeded at 7.5x105 into plates of 5 cm diameter. The next day subconfluent monolayers were washed once with PBS containing 40 µg DEAE-dextran ml–1 and incubated for 1 h at 37 °C with PBS containing 40 µg DEAE-dextran ml–1, 2 % FCS and MHV-A59 at an m.o.i. of 10. Mock-infected cells were treated in the same way with the omission of virus. After adsorption, the inoculum was removed, cells were washed twice, supporting medium was added and cells were incubated at 37 °C. In order to reduce the effect of biological variation in cellular gene expression, we performed two independent infection experiments for microarray analysis with a 3 week interval. For the two experiments, we used the same batches of all solutions, the same cell passage and the same viral stock and strove to reproduce identical conditions. Medium for virus titration and cells for RNA isolation were harvested at the indicated time points. Total RNA was isolated from the cells using TRIzol reagent (Gibco-BRL) as recommended by the manufacturer. RNA concentration was determined spectrophotometrically. For quality and quantity control, 4 µg total RNA was subjected to electrophoresis in 1 % native agarose gel and stained with ethidium bromide.

Sample preparation and chip hybridization.
Samples for hybridization to microarrays were prepared from total cellular RNA of MHV- [at 3, 5 and 6 h post-infection (p.i.)] and mock-infected (at 3 and 5 h p.i.) cells according to the Affymetrix GeneChip manual. Briefly, RNA was purified on RNeasy columns (Qiagen). A total of 12 µg RNA eluted from the columns was used for cDNA synthesis with SuperScript II reverse transcriptase (Invitrogen) and a primer containing a T7 RNA polymerase promoter site with a 5'-stretch of (dT)24. Second strand cDNA was synthesized in the presence of Escherichia coli DNA polymerase I, DNA ligase and RNase H (all enzymes from Invitrogen). Half of the dsDNA served as a template for cRNA synthesis that was performed with T7 RNA polymerase (Ambion), 100 mM of each NTP and 10 mM of each biotin-11-CTP and biotin-16-UTP (PerkinElmer). In vitro transcription yielded 15- to 20-fold amplification of the starting amount of RNA. RNA was purified on RNeasy columns and fragmented. Hybridization of RNA to Affymetrix murine GeneChip U74Av2, staining and scanning was performed in the Leiden Genome Technology Center following the manufacturer's recommendations.

Chip data analysis.
Affymetrix murine GeneChip U74Av2 has probe sets for 12 422 characterized genes and expressed sequences. Each sequence is represented by 16 probe oligonucleotide pairs, one oligonucleotide of the pair being a perfect match for the sequence and another having a single mismatch. The microarray scanning data were analysed with the use of Rosetta Resolver software, version 5.0 (Rosetta Inpharmatics LLC) according to its manual. Rosetta Resolver uses intensity signals of perfect match and mismatch probes and of background to calculate a P-value that reflects the level of confidence that a sequence is expressed. Ten transcriptional profiles of virus- and mock-infected cells were subjected to two-dimensional clustering. Normalization for clustering was achieved by transformation of intensity values into z-scores to emphasize the relative change of intensity of a gene across transcriptional profiles. An agglomerative clustering algorithm, error-weighted Cosine correlation as a similarity measure for profiles and genes, a complete linkage method for profiles and mean linkage method for genes were applied for clustering. A total of 938 gene clusters with the coefficient of variation more than 0.3 were included. In order to identify sequences with virus-induced changes in the expression level, we compared profiles of infected cells at 3 and 5 h p.i. with corresponding profiles of mock-infected cells. Profiles of MHV-infected cells at 6 h p.i. were compared with those of mock-infected at 5 h p.i. Combined replicate profiles from two independent experiments were used in the above comparisons. A positive fold-induction is the ratio of sequence intensity values in infection and mock-infection profiles. A negative fold-induction is the reverse ratio (mock- to MHV-infection) with the minus sign. The P-value of the ratio takes into account differences in intensities of infected and control probe sets as well as intensity variance within each probe set and reflects the confidence of differential expression. For all analyses we set the P-value to be less than 0.01.

Sequences that showed differential expression in infected cells were grouped according to Gene Ontology terms of biological processes or cellular components, which were available on the website of the National Center for Biotechnology Information in March, 2005. Sequences not yet annotated by Gene Ontology were not analysed further. Genes that could be placed into more than one group according to their annotation, were put arbitrarily into a single group. Here, we present differentially expressed genes that belong to the selected 14 groups and that were significantly changed by 2.5-fold or more (P-value<0.01). We used gene symbols defined by the rules of the International Committee on Standardized Genetic Nomenclature for Mice (http://www.informatics.jax.org/mgihome/nomen/gene.shtml).

Quantitative RT-PCR.
cDNA was synthesized with 1 µg total cellular RNA as template, 0.5 µg random hexamers (Promega) and SuperScript II reverse transcriptase (Invitrogen). Random hexamers were chosen as primers in order to represent the relative amount of each RNA species in the cells. Aliquots of the same cDNA solution were used for quantitative PCR (qPCR) on different templates. Gene-specific primers (nucleotide sequences are available on request) for qPCR were selected with the help of Primer3 software (Rozen & Skaletsky, 2000Down). A primer set for the MHV-A59 3'-untranslated region, annealing to genome-length and subgenomic viral RNA, was used for the detection of total viral RNA. Primers for MHV-A59 ORF1b were used for the detection of genome-length viral RNA. qPCR was performed with HotStar Taq polymerase (Qiagen) and SYBR Green I (Molecular Probes) in an iCycler PCR machine (Bio-Rad) for 40 cycles with three steps in a cycle, each step for 20 s, a denaturing step at 94 °C, an annealing step at a temperature specific for a given pair of primers (58–63 °C) and an extension step at 72 °C. Data were analysed with iCycler software. Specificity of qPCR was confirmed by the melting curve of amplified products. Ten-fold dilutions of a cDNA sample with the highest concentration of a given template were used to build standard curves. In a similar way, standard curves for total and genome-length viral RNA were built using cDNA of a mixture of 0.9 µg total RNA from uninfected cells and 0.1 µg virion RNA. Virion RNA was isolated from virions purified by isokinetic sedimentation in a sucrose gradient as described before (Spaan et al., 1981Down). The amount of RNA in a sample was determined with respect to a standard sample and expressed in relative units. Results of quantification were normalized to the amount of Gapdh mRNA in the same sample in order to correct for efficiency of reverse transcription. Finally, results were arbitrarily normalized to fit a convenient scale. Each qPCR was performed in triplicate, means and standard deviation (SD) were calculated. The level of interferon {alpha}4 and beta (Ifna4 and Ifnb) mRNA in infected and mock-infected cells was below the detection limit. Ifn fold-induction after poly (I : C) treatment was calculated as a ratio of the IFN RNA level at a given time point and the detection limit.

Ifn induction by poly (I : C).
L cells were seeded as for an infection experiment. The next day cells were washed with PBS and incubated for 2 h at 37 °C in supporting medium containing 400 µg DEAE-dextran ml–1. After that, cells were washed with PBS and incubated for 32 h at 37 °C in supporting medium with 100 µg poly (I : C) ml–1 (Sigma). At the indicated time after poly (I : C) addition, cells were dissolved in TRIzol for total cellular RNA isolation.


   RESULTS AND DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Virus replication
We set out to characterize changes in host cell expression in fibroblast-like L cells during lytic infection with MHV. The course of a robust cytopathic MHV infection was described more than 30 years ago (Sturman & Takemoto, 1972Down). In order to be able to correlate changes in cellular expression with parameters of viral infection, we determined the kinetics of virus production, viral RNA content and CPE development in L cells that are maintained in our lab. We used an m.o.i. of 10 in order to infect all cells simultaneously. During the first 4 h p.i., no visible changes in the cell morphology were observed and no virus was secreted into the medium (Fig. 1aDown). The end of the lag period and the beginning of the exponential phase of virus production coincided with the appearance of CPE, i.e. the formation of multinuclear giant cells, at 5 h p.i. Viral production reached a maximum level at 9 h p.i. when nearly all cells were fused into syncytia. At this stage virus yield was 1000 p.f.u. cell–1 (Fig. 1aDown). The use of a sensitive RT-qPCR allowed us to measure the levels of viral RNA from the beginning of infection. The increase of viral RNA levels preceded that of viral production by 2 h (Fig. 1bDown). With a 1 h interval between time points, the accumulation of both genome-length and subgenomic RNA started simultaneously at 2 h p.i. The molar ratio of total viral RNA to genome-size viral RNA reached 10 : 1 at 4 h p.i. Later in infection this ratio remained constant, which was similar to that described for bovine coronavirus infection (Hofmann et al., 1990Down).


Figure 1
View larger version (12K):
[in this window]
[in a new window]
 
Fig. 1. Kinetics of viral production (a) and viral RNA content (b) in L cells infected with MHV A59. Virus was added at time point 0 h at an m.o.i. of 10 p.f.u. per cell. Reference experiment (dashed line) and microarray experiments (solid line) are shown. Microarray experiments values are the mean of two replicate values. Bars show variation between results from two infections. Arrows indicate time points used for microarray analysis. RNA content (b) was determined by RT-qPCR and expressed in relative units (RU). Values for total viral RNA ({blacktriangleup}) and genome-size viral RNA ({blacksquare}) are shown.

 
For analysis of cellular expression, we repeated the infection of L cells twice, harvested medium and isolated RNA from MHV- and mock-infected cells at hourly intervals. The kinetics of viral production and RNA synthesis were very similar in both experiments (Fig. 1a and bUp). We chose 3, 5 and 6 h p.i. as time points for expression profiling. At 3 h p.i., the concentration of viral products within infected cells was apparently sufficient to support exponential viral RNA accumulation, but there was no virus secretion and no obvious CPE yet. At 5 and 6 h p.i., viral production was in the exponential phase. The levels of viral RNA increased by 100-fold (genome-size) and 200-fold (total) at 5 h p.i. as compared with 3 h p.i. Multinuclear cells started to appear at 5 h p.i. and harboured approximately 25 % of the nuclei at 6 h p.i. Expression profiles in mock-infected cells at 3 and 5 h p.i. served as controls.

Overall cellular gene expression during MHV infection
Approximately 44–48 % of 12 422 sequences represented on the chips were expressed in mock-infected L cells and early in infection. Later in infection the number of expressed genes decreased to 41 and 38 % at 5 and 6 h p.i., respectively. Clustering analysis allowed us to visualize differences in cellular gene expression. Profiles of infected L cells at 3 h p.i. clustered together with those of mock-infected cells, which indicated a very low impact of viral infection on the steady-state levels of cellular RNAs at this time point (Fig. 2Down). A pronounced effect of infection became obvious at 5 and especially 6 h p.i. (Fig. 2Down). Clustering illustrated biological variance in two independent experiments. Despite our attempts to reproduce the experimental conditions, profiles of mock-infected cells and of MHV-infected cells at 3 h p.i. formed two separate clusters according to the experiment (Fig. 2Down). Combining two replicates for further analysis increased the statistical significance of the analysis and filtered out non-consistent changes that were present only in one experiment. Reliability of the analysis was confirmed further by the consistency of many changes throughout the course of infection (see Tables 1–4DownDownDownDown).


Figure 2
View larger version (43K):
[in this window]
[in a new window]
 
Fig. 2. Clustering of expression profiles of MHV- and mock-infected L cells. Gene expression is shown as a matrix with rows representing profiles and columns representing genes. Gene expression in an individual profile was compared to a mean value for this gene based on 10 profiles. Red colour indicates upregulation, green indicates downregulation and black is used for nochange in gene expression. The names of the profiles include experiment number (I or II), infection (MHV or mock) and time p.i. The dendrogram on the left illustrates the structure of clusters of the profiles.

 

View this table:
[in this window]
[in a new window]
 
Table 1. Innate immunity, inflammatory and defence responses and signalling

 

View this table:
[in this window]
[in a new window]
 
Table 2. RNA metabolism

 

View this table:
[in this window]
[in a new window]
 
Table 3. Protein metabolism

 

View this table:
[in this window]
[in a new window]
 
Table 4. Oxidative stress, cell cycle, DNA repair and apoptosis

 
Validation of microarray hybridization results
The level of expression of several genes was measured by RT-qPCR to validate the microarray data by an independent method. mRNA levels of Fos, macrophage inflammatory protein 2 (MIP-2 or Cxcl2), beta-2 microglobulin (B2m), non-POU-domain-containing octamer-binding protein (Nono) and Ifnb and Ifna4 were determined in MHV- and mock-infected cells not only at the time points of the microarray analysis but throughout a 6 h time-course. qPCR and microarray results for these genes were in agreement (Fig. 3Down, inset tables). Levels of Fos and Cxcl2 transcripts increased, B2m and Nono transcript content dropped (Fig. 3Down), and Ifnb and Ifna4 mRNA was undetectable during MHV infection (data not shown).


Figure 3
View larger version (10K):
[in this window]
[in a new window]
 
Fig. 3. Kinetics of mRNA content of Fos, Cxcl2, B2m and Nono inMHV- ({blacksquare}) and mock-infected ({square}) cells. mRNA content in the cells was determined with the use of RT-qPCR and expressed in relative units (RU). Inset tables show fold-induction as determined by hybridization to microarrays (MA) and by RT-qPCR (qPCR). Microarray data for Nono are from two probe sets. Values marked with an asterisk had a P-value more than 0.01.

 
Cellular genes affected by MHV infection
We grouped differentially expressed genes according to the Gene Ontology annotation of biological processes or cellular components. Here, we present genes (see Tables 1–4UpUpUpUp) that belonged to selected groups and had more than 2.5-fold changes in expression in comparison with control samples.

Innate immunity, inflammatory and defence responses, and signalling.
MHV infection in L cells resulted in an induction of three inflammatory chemokines: Cxcl1, Cxcl2 and Cx3cl1 (Table 1Up). An increase in Cxcl2 RNA concentration was detected at 3 h p.i. and reached a sevenfold increase at 6 h p.i. Cxcl1 and Cxcl2 are immediate-early chemokines that attract neutrophils to sites of injury in vivo. Cx3cl1 functions as a chemokine as well as an adhesion molecule. Together with Cxcl1 and Cxcl2, induction of interleukin-6 (Il6) and tumour necrosis factor alpha (Tnfa) mRNA were found to be prominent responses to MHV infection in neural cells (Rempel et al., 2005Down). However, MHV did not induce expression of Il6 and Tnfa in L cells. Differences in cytokine expression patterns between MHV-infected neural and fibroblast-like cells suggest different mechanisms of induction of these cytokines. An anti-inflammatory response in virus-infected L cells was illustrated by suppression of Jak3 and elevation of suppressor of cytokine signalling 3 (Socs3) expression (Table 1Up).

Despite the fact that L cells were able to sense the viral infection and respond to it with chemokine induction, the interferon pathway was not activated at any time point during MHV infection. The absence of transcripts of immediate-early interferons (Ifna4 and Ifnb) was confirmed by RT-qPCR (data not shown). In order to exclude the possibility that our line of L cells had a defect in the interferon induction pathway, we treated the cells with a known interferon inducer, poly (I : C). Immediate-early Ifns mRNA levels increased more than 10-fold after 3 h of treatment (Fig. 4Down). Maximum induction of approximately 1000-fold was achieved at 22 h. Induction was also obtained with Sendai virus infection in L cells (data not shown). It was shown previously that various strains of MHV, including A59, are poor inducers of interferon in cultured cells (Garlinghouse et al., 1984Down; Sturman & Takemoto, 1972Down). Whether MHV hides from recognition by the interferon induction system or actively suppresses it remains a question for future research.


Figure 4
View larger version (7K):
[in this window]
[in a new window]
 
Fig. 4. Induction of immediate-early interferons-{alpha}4 ({blacktriangleup}) and -beta ({triangleup}) mRNA in L cells by poly (I : C) treatment. L cells were pre-treated for 2 h with DEAE-dextran and at time point 0 h100 µg poly (I : C) ml–1 was added. Ifn RNA content was determined by RT-qPCR and fold induction was calculated as described in Methods.

 
MHV infection led to downregulation of genes encoding major histocompatibility complex (MHC)-I proteins (Table 1Up) that present intracellular antigens to specialized immune cells. Tripeptidyl peptidase II (Tpp2), which plays a crucial role in MHC-I presentation (Reits et al., 2004Down), was also downregulated (Table 3Up). Reduction in mRNA levels and surface expression of MHC-I was found previously in a macrophage cell line infected with MHV (Kyuwa et al., 1994Down). However, in non-lytically infected cells MHC-I has been shown to be upregulated, probably by paracrine action (Lavi et al., 1989Down; Rempel et al., 2005Down). Mouse cells ubiquitously express Crry, a cell surface-bound complement regulatory protein, capable of inhibiting C3 convertase activity and thereby protecting cells from homologous complement (Nomura et al., 2002Down). RNA levels for this gene were greatly reduced in MHV-infected cells (Table 1Up), potentially rendering infected cells vulnerable to destructive complement action.

MHV infection affected expression of a number of genes associated with signalling cascades that regulate homeostasis, growth, differentiation, response to stress and death. Mapk signalling mediates inflammatory responses and it might be involved in induction of the chemokines described above. Dual specificity phosphatases (Dusp) modulate intracellular signalling by dephosporylation of Mapks thus limiting an inflammatory reaction. During MHV infection, Dusp1, 2 and 8 were upregulated (Table 1Up). This result supports the observation that both inflammatory and anti-inflammatory signals were activated in MHV-infected cells.

RNA metabolism.
The largest group of genes differentially expressed in MHV-infected cells was associated with cellular transcription. In this group, transcription factor Fos showed the strongest induction. Fos, as well as Jun, Junb and Egr1 (Krox-24), which were also upregulated in infected cells (Table 2Up), are induced in response to diverse signals. These transcription cofactors initiate ordered expression of target genes thus regulating cell division, apoptosis and inflammatory responses. Fos, Jun and Junb are constituents of the activator protein-1 (AP-1) transcription factor complex that was shown to play a key role in the induction of IL-8 by the spike glycoprotein of SARS-CoV (Chang et al., 2004Down). In the same way, upregulation of AP-1 genes in MHV-infected mouse cells could lead to an induction of chemokine Cxcl2, a mouse functional counterpart of human IL-8. On the other hand, downregulation of some other transcription factors implied suppressive action of the virus. Activation of Nfkb (NF-{kappa}B) by its release from a complex with inhibitor Nfkbia (I{kappa}B{alpha}) or by transcriptional induction is a common feature of a response to pathogens (Jenner & Young, 2005Down). However, in MHV-A59-infected L cells, the level of Nfkb mRNA was decreasing (–2.3-fold at 6 h p.i.), while its inhibitor Nfkbia was induced (Table 2Up).

Differentially expressed genes participating in RNA processing and splicing and in export of RNA from the nucleus were prevalently downregulated (Table 2Up). The earliest change in expression in the RNA metabolism group was detected for the Nono gene that was downregulated from 3 h p.i. onwards (Table 2Up). Nono in complex with proline/glutamine-rich splicing factor Sfpq, the expression of which was also downregulated, participates in a variety of nuclear processes including transcription and pre-mRNA processing (Shav-Tal & Zipori, 2002Down). The finding that the human homologue of Nono constantly shuttles between the nucleus and cytoplasm (Zolotukhin et al., 2003Down) suggested that this protein may also have a cytoplasmic function.

Protein metabolism.
Infection by various viruses results in endoplasmic reticulum (ER) stress that can be caused by either synthesis of proteins foreign to the host cells (Dimcheff et al., 2004Down) or rapid accumulation of substantial amounts of viral, especially misfolded, proteins in the ER (Li et al., 2005Down; Netherton et al., 2004Down; Su et al., 2002Down). MHV replication in its exponential phase led to induction of three genes that are considered to be key markers of the unfolded protein response: homocysteine-responsive endoplasmic reticulum resident protein 1 (Herpud1 or Herp), DNA-damage inducible transcript 3 (Ddit3 or Chop) (Table 3Up) and heat shock 70 kDa protein 5 (Hspa5 or Bip; 1.5-fold at 6 h p.i.).

A number of genes encoding ubiquitin-specific proteases (Usp) were downregulated in infected cells (Table 3Up). Suppression of Usp genes could indicate enhanced ubiquitin-mediated protein degradation as part of the unfolded protein response and could also affect other possible functions of Usps, such as regulation of transcription, membrane protein trafficking and signal transduction (Amerik & Hochstrasser, 2004Down). It was recently shown that coronavirus papain-like proteinase PLpro can act as a deubiquitinating enzyme (Barretto et al., 2005Down; Lindner et al., 2005Down). Knowledge of changes in the cellular environment might help in identifying the function and specificity of the virus-coded enzyme.

Oxidative stress, cell cycle, DNA repair and apoptosis.
Additional evidence for ER-overload in MHV-infected cells was supplied by detection of upregulated expression of the myeloperoxidase (Mpo) and thioredoxin interacting protein (Txnip) genes, which are characteristic for oxidative stress (Table 4Up). Stress conditions can eventually lead to cell cycle arrest and apoptosis. Genes with a strongly changed expression later in infection were consistent with an anti-apoptotic programme (Table 4Up). Caspase 7 (Casp7) and death-associated kinase 3 (Dapk3) (Kawai et al., 1998Down) participate in the induction of apoptosis and their mRNA levels were significantly decreased in infected cells. Other genes involved in apoptosis had smaller changes in expression. Taken together, both pro- as well as anti-apoptotic genes were differentially expressed as a result of lytic infection. The overall effect on the induction of apoptosis in coronavirus-infected L cells remains yet inconclusive.

Concluding remarks
We applied microarray analysis to study changes in the steady-state levels of cellular mRNAs during a one-step growth curve of MHV. Viral infection can influence the transcriptional profile of infected cells by affecting rates of mRNA synthesis and processing, transport from the nucleus to cytoplasm, and decay in the cytoplasm. Coronavirus proteins including RNA modifying enzymes (Snijder et al., 2003Down) and proteinases can act directly on cellular mRNA or on cellular enzymes involved in RNA metabolism. In addition, transcript levels can be modulated as a cellular response to viral products. Ongoing improvement of nucleotide sequence annotations and microarray analysis tools might support new findings in the same datasets. Future studies are needed to reveal whether protein levels of the differentially expressed genes are also changed in MHV-infected cells and what role these proteins play in coronavirus replication.

Recently, Rempel et al. (2005)Down reported microarray analysis of primary culture of a mixed population of neural cells infected with MHV. Infection in the neural culture was slow and non-lytic with only 10 % of the cells infected by JHM strain of MHV at the time of transcription profiling on day 3 p.i. (Rempel et al., 2005Down). Interestingly, despite big differences in experimental setup, some cellular genes showed the same trends in expression in two different models of MHV–host interaction. Among the common responses were upregulated genes encoding chemokines, Junb and Nfkbia. Probably, mechanisms of induction of these genes are common for these two models and can therefore be conveniently studied in an established cell line in the future. The relatively small overlap of two expression datasets from neural and fibroblast-like cells can be explained by reasoning that the majority of differentially expressed genes in the primary neural culture represented paracrine effects – changes in uninfected cells that constituted 90 % of the cells in the culture.

In the course of MHV infection in L cells, chemokine induction became gradually more prominent. In contrast to chemokine induction, no increase in mRNA levels of interferons or any of their known downstream effectors was detected at any of the time points investigated. These data are in agreement with the long-standing observation that many coronavirus strains are sensitive to interferon pre-treatment, but poor interferon inducers in infected cells (Aurisicchio et al., 2000Down; Garlinghouse et al., 1984Down; Matsuyama et al., 2000Down; Taguchi & Siddell, 1985Down; Virelizier & Gresser, 1978Down). Recently, Spiegel et al. (2005)Down concluded that SARS-CoV interferes with the key transcription factor, interferon regulatory factor 3 (IRF-3), involved in interferon induction, which suggests that interferon response can be counteracted by a coronavirus-encoded antagonist as has previously been described for several other (+) RNA-containing viruses (Frolova et al., 2002Down; van Pesch et al., 2001Down; Zoll et al., 2002Down).

Our analysis provided data on the cellular environment of coronavirus replication. On the basis of our current knowledge of the strategy of coronavirus genome expression, some of the genes with changed expression (Nono, Sfpq, Usps) might play a direct role in the viral life cycle. Genes encoding chemokines, constituents of transcription factor complex AP-1, proteins involved in stress responses (to unfolded protein and oxidative stress) are likely to represent cellular responses to viral infection. Newly identified cellular genes from either group can serve as potential targets for the development of antivirals if these responses are common to coronavirus infections. Current research in our laboratory is aimed at identification of differences in expression profiles between acutely and persistently infected cells as well as unravelling of induction mechanisms of the chemokine responses described in this work.


   ACKNOWLEDGEMENTS
 
We are grateful to Peter Bredenbeek for helpful suggestions during the manuscript preparation and critical reading of the manuscript. Thanks are due to Alexander Gorbalenya for help in selecting PCR primers and stimulating discussions. We acknowledge Evelyne Mank and Peter Taschner from the Leiden Genome Technology Center for their assistance during the experiments involving the Affymetrix gene arrays and the use of Rosetta Resolver.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Amerik, A. Y. & Hochstrasser, M. (2004). Mechanism and function of deubiquitinating enzymes. Biochim Biophys Acta 1695, 189–207.[Medline]

Aurisicchio, L., Delmastro, P., Salucci, V., Paz, O. G., Rovere, P., Ciliberto, G., La Monica, N. & Palombo, F. (2000). Liver-specific alpha 2 interferon gene expression results in protection from induced hepatitis. J Virol 74, 4816–4823.[Abstract/Free Full Text]

Barretto, N., Jukneliene, D., Ratia, K., Chen, Z., Mesecar, A. D. & Baker, S. C. (2005). The papain-like protease of severe acute respiratory syndrome coronavirus has deubiquitinating activity. J Virol 79, 15189–15198.[Abstract/Free Full Text]

Cai, Y., Liu, Y., Yu, D. & Zhang, X. (2003). Down-regulation of transcription of the proapoptotic gene BNip3 in cultured astrocytes by murine coronavirus infection. Virology 316, 104–115.[CrossRef][Medline]

Chang, Y. J., Liu, C. Y., Chiang, B. L., Chao, Y. C. & Chen, C. C. (2004). Induction of IL-8 release in lung cells via activator protein-1 by recombinant baculovirus displaying severe acute respiratory syndrome-coronavirus spike proteins: identification of two functional regions. J Immunol 173, 7602–7614.[Abstract/Free Full Text]

Dimcheff, D. E., Faasse, M. A., McAtee, F. J. & Portis, J. L. (2004). Endoplasmic reticulum (ER) stress induced by a neurovirulent mouse retrovirus is associated with prolonged BiP binding and retention of a viral protein in the ER. J Biol Chem 279, 33782–33790.[Abstract/Free Full Text]

Fouchier, R. A., Kuiken, T., Schutten, M. & 7 other authors (2003). Aetiology: Koch's postulates fulfilled for SARS virus. Nature 423, 240.[CrossRef][Medline]

Frolova, E. I., Fayzulin, R. Z., Cook, S. H., Griffin, D. E., Rice, C. M. & Frolov, I. (2002). Roles of nonstructural protein nsP2 and Alpha/Beta interferons in determining the outcome of Sindbis virus infection. J Virol 76, 11254–11264.[Abstract/Free Full Text]

Garlinghouse, L. E., Jr, Smith, A. L. & Holford, T. (1984). The biological relationship of mouse hepatitis virus (MHV) strains and interferon: in vitro induction and sensitivities. Arch Virol 82, 19–29.[CrossRef][Medline]

Gosert, R., Kanjanahaluethai, A., Egger, D., Bienz, K. & Baker, S. C. (2002). RNA replication of mouse hepatitis virus takes place at double-membrane vesicles. J Virol 76, 3697–3708.[Abstract/Free Full Text]

Hilton, A., Mizzen, L., MacIntyre, G., Cheley, S. & Anderson, R. (1986). Translational control in murine hepatitis virus infection. J Gen Virol 67, 923–932.[Abstract/Free Full Text]

Hofmann, M. A., Sethna, P. B. & Brian, D. A. (1990). Bovine coronavirus mRNA replication continues throughout persistent infection in cell culture. J Virol 64, 4108–4114.[Abstract/Free Full Text]

Jacobs, L., Spaan, W. J., Horzinek, M. C. & van der Zeijst, B. A. (1981). Synthesis of subgenomic mRNAs of mouse hepatitis virus is initiated independently: evidence from UV transcription mapping. J Virol 39, 401–406.[Abstract/Free Full Text]

Jenner, R. G. & Young, R. A. (2005). Insights into host responses against pathogens from transcriptional profiling. Nat Rev Microbiol 3, 281–294.[CrossRef][Medline]

Kawai, T., Matsumoto, M., Takeda, K., Sanjo, H. & Akira, S. (1998). ZIP kinase, a novel serine/threonine kinase which mediates apoptosis. Mol Cell Biol 18, 1642–1651.[Abstract/Free Full Text]

Kyuwa, S., Cohen, M., Nelson, G., Tahara, S. M. & Stohlman, S. A. (1994). Modulation of cellular macromolecular synthesis by coronavirus: implication for pathogenesis. J Virol 68, 6815–6819.[Abstract/Free Full Text]

Lavi, E., Suzumura, A., Murray, E. M., Silberberg, D. H. & Weiss, S. R. (1989). Induction of MHC class I antigens on glial cells is dependent on persistent mouse hepatitis virus infection. J Neuroimmunol 22, 107–111.[CrossRef][Medline]

Lavi, E., Wang, Q., Weiss, S. R. & Gonatas, N. K. (1996). Syncytia formation induced by coronavirus infection is associated with fragmentation and rearrangement of the Golgi apparatus. Virology 221, 325–334.[CrossRef][Medline]

Li, Y., Fu, L., Gonzales, D. M. & Lavi, E. (2004). Coronavirus neurovirulence correlates with the ability of the virus to induce proinflammatory cytokine signals from astrocytes and microglia. J Virol 78, 3398–3406.[Abstract/Free Full Text]

Li, X. D., Lankinen, H., Putkuri, N., Vapalahti, O. & Vaheri, A. (2005). Tula hantavirus triggers pro-apoptotic signals of ER stress in Vero E6 cells. Virology 333, 180–189.[CrossRef][Medline]

Lindner, H. A., Fotouhi-Ardakani, N., Lytvyn, V., Lachance, P., Sulea, T. & Menard, R. (2005). The papain-like protease from the severe acute respiratory syndrome coronavirus is a deubiquitinating enzyme. J Virol 79, 15199–15208.[Abstract/Free Full Text]

Matsuyama, S., Henmi, S., Ichihara, N., Sone, S., Kikuchi, T., Ariga, T. & Taguchi, F. (2000). Protective effects of murine recombinant interferon-beta administered by intravenous, intramuscular or subcutaneous route on mouse hepatitis virus infection. Antiviral Res 47, 131–137.[CrossRef][Medline]

Netherton, C. L., Parsley, J. C. & Wileman, T. (2004). African swine fever virus inhibits induction of the stress-induced proapoptotic transcription factor CHOP/GADD153. J Virol 78, 10825–10828.[Abstract/Free Full Text]

Ning, Q., Berger, L., Luo, X., Yan, W., Gong, F., Dennis, J. & Levy, G. (2003a). STAT1 and STAT3 alpha/beta splice form activation predicts host responses in mouse hepatitis virus type 3 infection. J Med Virol 69, 306–312.[CrossRef][Medline]

Ning, Q., Lakatoo, S., Liu, M., Yang, W., Wang, Z., Phillips, M. J. & Levy, G. A. (2003b). Induction of prothrombinase fgl2 by the nucleocapsid protein of virulent mouse hepatitis virus is dependent on host hepatic nuclear factor-4 alpha. J Biol Chem 278, 15541–15549.[Abstract/Free Full Text]

Nomura, M., Kurita-Taniguchi, M., Kondo, K., Inoue, N., Matsumoto, M., Yamanishi, K., Okabe, M. & Seya, T. (2002). Mechanism of host cell protection from complement in murine cytomegalovirus (CMV) infection: identification of a CMV-responsive element in the CD46 promoter region. Eur J Immunol 32, 2954–2964.[CrossRef][Medline]

Prentice, E., Jerome, W. G., Yoshimori, T., Mizushima, N. & Denison, M. R. (2004). Coronavirus replication complex formation utilizes components of cellular autophagy. J Biol Chem 279, 10136–10141.[Abstract/Free Full Text]

Rao, P. V. & Gallagher, T. M. (1998). Intracellular complexes of viral spike and cellular receptor accumulate during cytopathic murine coronavirus infections. J Virol 72, 3278–3288.[Abstract/Free Full Text]

Reits, E., Neijssen, J., Herberts, C., Benckhuijsen, W., Janssen, L., Drijfhout, J. W. & Neefjes, J. (2004). A major role for TPPII in trimming proteasomal degradation products for MHC class I antigen presentation. Immunity 20, 495–506.[CrossRef][Medline]

Rempel, J. D., Quina, L. A., Blakely-Gonzales, P. K., Buchmeier, M. J. & Gruol, D. L. (2005). Viral induction of central nervous system innate immune responses. J Virol 79, 4369–4381.[Abstract/Free Full Text]

Rozen, S. & Skaletsky, H. (2000). Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol 132, 365–386.[Medline]

Shav-Tal, Y. & Zipori, D. (2002). PSF and p54(nrb)/NonO – multi-functional nuclear proteins. FEBS Lett 531, 109–114.[CrossRef][Medline]

Shi, S. T., Schiller, J. J., Kanjanahaluethai, A., Baker, S. C., Oh, J. W. & Lai, M. M. (1999). Colocalization and membrane association of murine hepatitis virus gene 1 products and de novo-synthesized viral RNA in infected cells. J Virol 73, 5957–5969.[Abstract/Free Full Text]

Siddell, S., Wege, H., Barthel, A. & ter Meulen, V. (1981). Coronavirus JHM: intracellular protein synthesis. J Gen Virol 53, 145–155.[Abstract/Free Full Text]

Snijder, E. J., Bredenbeek, P. J., Dobbe, J. C. & 7 other authors (2003). Unique and conserved features of genome and proteome of SARS-coronavirus, an early split-off from the coronavirus group 2 lineage. J Mol Biol 331, 991–1004.[CrossRef][Medline]

Spaan, W. J., Rottier, P. J., Horzinek, M. C. & van der Zeijst, B. A. (1981). Isolation and identification of virus-specific mRNAs in cells infected with mouse hepatitis virus (MHV-A59). Virology 108, 424–434.[CrossRef][Medline]

Spiegel, M., Pichlmair, A., Martinez-Sobrido, L., Cros, J., Garcia-Sastre, A., Haller, O. & Weber, F. (2005). Inhibition of Beta interferon induction by severe acute respiratory syndrome coronavirus suggests a two-step model for activation of interferon regulatory factor 3. J Virol 79, 2079–2086.[Abstract/Free Full Text]

Sturman, L. S. & Takemoto, K. K. (1972). Enhanced growth of a murine coronavirus in transformed mouse cells. Infect Immun 6, 501–507.[Abstract/Free Full Text]

Su, H. L., Liao, C. L. & Lin, Y. L. (2002). Japanese encephalitis virus infection initiates endoplasmic reticulum stress and an unfolded protein response. J Virol 76, 4162–4171.[Abstract/Free Full Text]

Taguchi, F. & Siddell, S. G. (1985). Difference in sensitivity to interferon among mouse hepatitis viruses with high and low virulence for mice. Virology 147, 41–48.[CrossRef][Medline]

van der Hoek, L., Pyrc, K., Jebbink, M. F. & 7 other authors (2004). Identification of a new human coronavirus. Nat Med 10, 368–373.[CrossRef][Medline]

van der Meer, Y., Snijder, E. J., Dobbe, J. C., Schleich, S., Denison, M. R., Spaan, W. J. & Locker, J. K. (1999). Localization of mouse hepatitis virus nonstructural proteins and RNA synthesis indicates a role for late endosomes in viral replication. J Virol 73, 7641–7657.[Abstract/Free Full Text]

van Pesch, V., van Eyll, O. & Michiels, T. (2001). The leader protein of Theiler's virus inhibits immediate-early alpha/beta interferon production. J Virol 75, 7811–7817.[Abstract/Free Full Text]

Vennema, H., Heijnen, L., Zijderveld, A., Horzinek, M. C. & Spaan, W. J. (1990). Intracellular transport of recombinant coronavirus spike proteins: implications for virus assembly. J Virol 64, 339–346.[Abstract/Free Full Text]

Virelizier, J. L. & Gresser, I. (1978). Role of interferon in the pathogenesis of viral diseases of mice as demonstrated by the use of anti-interferon serum. V. Protective role in mouse hepatitis virus type 3 infection of susceptible and resistant strains of mice. J Immunol 120, 1616–1619.[Abstract/Free Full Text]

Zoll, J., Melchers, W. J., Galama, J. M. & van Kuppeveld, F. J. (2002). The mengovirus leader protein suppresses alpha/beta interferon production by inhibition of the iron/ferritin-mediated activation of NF-{kappa}B. J Virol 76, 9664–9672.[Abstract/Free Full Text]

Zolotukhin, A. S., Michalowski, D., Bear, J., Smulevitch, S. V., Traish, A. M., Peng, R., Patton, J., Shatsky, I. N. & Felber, B. K. (2003). PSF acts through the human immunodeficiency virus type 1 mRNA instability elements to regulate virus expression. Mol Cell Biol 23, 6618–6630.[Abstract/Free Full Text]

Received 9 December 2005; accepted 2 March 2006.


This article has been cited by other articles:


Home page
J. Virol.Home page
J. K. Roth-Cross, S. J. Bender, and S. R. Weiss
Murine Coronavirus Mouse Hepatitis Virus Is Recognized by MDA5 and Induces Type I Interferon in Brain Macrophages/Microglia
J. Virol., October 15, 2008; 82(20): 9829 - 9838.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
J. Bechill, Z. Chen, J. W. Brewer, and S. C. Baker
Coronavirus Infection Modulates the Unfolded Protein Response and Mediates Sustained Translational Repression
J. Virol., May 1, 2008; 82(9): 4492 - 4501.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
D. D. C. Ireland, S. A. Stohlman, D. R. Hinton, R. Atkinson, and C. C. Bergmann
Type I Interferons Are Essential in Controlling Neurotropic Coronavirus Infection Irrespective of Functional CD8 T Cells
J. Virol., January 1, 2008; 82(1): 300 - 310.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
O. Slobodskaya, A. Laarman, and W. J. M. Spaan
Intracellular Restriction of a Productive Noncytopathic Coronavirus Infection
J. Virol., January 1, 2008; 82(1): 451 - 460.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
G. A. Versteeg, P. S. van de Nes, P. J. Bredenbeek, and W. J. M. Spaan
The Coronavirus Spike Protein Induces Endoplasmic Reticulum Stress and Upregulation of Intracellular Chemokine mRNA Concentrations
J. Virol., October 15, 2007; 81(20): 10981 - 10990.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
J. K. Roth-Cross, L. Martinez-Sobrido, E. P. Scott, A. Garcia-Sastre, and S. R. Weiss
Inhibition of the Alpha/Beta Interferon Response by Mouse Hepatitis Virus at Multiple Levels
J. Virol., July 1, 2007; 81(13): 7189 - 7199.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
Y. Ye, K. Hauns, J. O. Langland, B. L. Jacobs, and B. G. Hogue
Mouse Hepatitis Coronavirus A59 Nucleocapsid Protein Is a Type I Interferon Antagonist
J. Virol., March 15, 2007; 81(6): 2554 - 2563.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
H. Zhou and S. Perlman
Mouse Hepatitis Virus Does Not Induce Beta Interferon Synthesis and Does Not Inhibit Its Induction by Double-Stranded RNA
J. Virol., January 15, 2007; 81(2): 568 - 574.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Versteeg, G. A.
Right arrow Articles by Spaan, W. J. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Versteeg, G. A.
Right arrow Articles by Spaan, W. J. M.
Agricola
Right arrow Articles by Versteeg, G. A.
Right arrow Articles by Spaan, W. J. M.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS