key: cord-0948804-71trm8pq
authors: Dar, Arshud; Munir, Shirin; Vishwanathan, Satya; Manuja, Anju; Griebel, Philip; Tikoo, Suresh; Townsend, Hugh; Potter, Andrew; Kapur, Vivek; Babiuk, Lorne A.
title: Transcriptional analysis of avian embryonic tissues following infection with avian infectious bronchitis virus
date: 2005-03-03
journal: Virus Res
DOI: 10.1016/j.virusres.2005.01.006
sha: b73388cfd1c1f5cafc537b3e41d53d64b2570ce5
doc_id: 948804
cord_uid: 71trm8pq

Avian infectious bronchitis virus (IBV) infection is one of the major viral respiratory diseases of chickens. Better understanding of the molecular basis of viral pathogenesis should contribute significantly towards the development of improved prophylactic, therapeutic and diagnostic reagents to control infections. In the present investigation, transcriptional profiles were analyzed by using RNA recovered from the lung tissue of IBV infected 18-day-old chicken embryos at 6, 24, 48 and 72 h post IBV infection. This microarray analysis was completed using avian cDNA arrays comprised of fragments of 1191 unique chicken and turkey gene transcripts. These arrays were generated from normalized cDNA subtraction libraries that were derived from avian pneumovirus (APV) infected chicken embryo fibroblast (CEF) cultures and tissues obtained from APV infected turkeys subtracted with their respective uninfected cultures and tissues. Of the 1191 unique genes represented on the array, the expression of a total of 327 genes (27% of total) were altered by two-fold or more from 6 through 72 h post-infection. A comparative analysis of IBV regulated genes with genes previously reported to change in expression following infection with other avian respiratory viruses revealed both conserved and unique changes. Real-time qRT-PCR was used to confirm the regulated expression of genes related to several functional classes including kinases, interferon induced genes, chemokines and adhesion molecules, vesicular trafficking and fusion protein genes, extracellular matrix protein genes, cell cycle, metabolism, cell physiology and development, translation, RNA binding, lysosomal, protein degradation and ubiquitination related genes. Microarray analysis served as an efficient tool in facilitating a comparative analysis of avian respiratory viral infections and provided insight into host transcriptional changes that were conserved as well as those which were unique to individual pathogens.

Avian viral respiratory infections caused by avian influenza virus, avian infectious bronchitis virus (IBV), avian pneumovirus (APV) and Newcastle disease virus (NDV) result in considerable economic losses to the poultry industry worldwide (Cavanagh, 2003) . IBV is one of the major respiratory viruses of chickens that is endemic in all infection in kidney and caecal tonsils. The endemic nature of infection, lack of effective vaccines due to the absence of inter-strain cross protection and circulation of various IBV strains and genotypes in poultry flocks makes the control of IBV infection a challenging task (Farsang et al., 2002) . It is well recognized that characterization of in vivo survival mechanisms of the invading organism and identification of the components of host response leading to the elimination of the invading pathogen may provide the basis of pathogenesis and immunity and may lead to improvements in control strategies (Qureshi et al., 1999) . The recent introduction of microarray technology has facilitated the large scale analysis of the molecular events occurring during avian respiratory viral infections (Munir and Kapur, 2003; Munir et al., 2005) . In order to understand the host cellular determinants involved in the pathogenesis and immunity to IBV, we initiated experiments involving the study of post-IBV infection gene transcriptional responses. Microarray analyses were performed with RNA isolated from lung tissue at various times following IBV infection of 18-day-old chick embryos. We selected 18-day-old embryos for this study because they proved an excellent model to investigate the molecular elements governing host responses to a single pathogen while excluding the potential for co-infection by other pathogens. Furthermore, chicken embryos are immunologically competent at this age and support IBV replication with pathological lesions in epithelial cells of various tissues including lungs, trachea, kidney and bursa of fabricius following IBV inoculation (Kapczynski et al., 2002; Lee et al., 2002) . Recently, even younger chicken embryos (8 and 14 days old) have successfully been used as models to study the tissue tropism and kinetics of other avian viral infections (Qureshi et al., 1999; Worthington et al., 2003) . Thus, chicken embryos provide a physiologically relevant model system for studying host responses to virus infection. In the current study, we identified several unique genes and pathways that may be involved in IBV infection. Further characterization of the role of these host factors in IBV pathogenesis will help reveal the commonalities and unique differences in the host response to avian viral pathogens and will offer better strategies for disease control.

The clinical IBV isolate (B8358) used in these studies was kindly provided by Dr. Keith West of the Prairie Diagnostic Services (PDS), University of Saskatchewan, Canada. The IBV isolate B8358 was isolated from 28-week-old layers during an IBV infection outbreak at a local poultry farm. The clinical signs associated with virus infection include subdued activity, hoarse sounding noises and squeak and mouth breathing of birds in barn. Sick birds submitted for examination and diagnosis 24 h after the appearance of clinical signs, showed mild catarrhal exudates and mildly hyperemic areas in trachea. Histopathological examination of tissues showed loss of cilia, sloughing of epithelial cells, mucosal congestion, edema with mild hemorrhage and infiltration of heterophils and lymphocytes and severe cuboidal metaplasia in trachea. No significant changes were seen in lung, heart, spleen, and kidney tissues. The virus was isolated from tracheal tissues, identified by electron microscopy and propagated in embryonated chicken eggs at the virology laboratory of PDS. The phylogenetic analysis of the S1 gene of IBV B8358 showed 80.3% sequence homology to the IBV serotype California (CA/633/85) (PDS unpublished data). Further, intra-tracheal inoculation of IBV (B8358) in 4 weeks old layers at VIDO showed similar disease symptoms and pathology as described above. Two hundred microliters of choroallontoic fluid containing IBV was inoculated into 18-day-old chick embryos via the chorioallontoic route [175 embryo lethal dose 50 (ELD 50 ) ]. This dose was selected as it was known to be sufficient to produce moderate clinical symptoms in 4-week-old chicks. IBV genes were not present on the array but to verify the presence of productive infection of embryonic lung tissue we quantified IBV N gene transcripts in RNA used for microarray analysis and gene expression. The virus present in IBV infected lung tissue collected at 6 h PI was not detectable by real-time PCR. However, 4.6 × 10 5 and 5.3 × 10 6 -fold increase in N gene transcript was detected in RNA samples derived from IBV infected embryonic lung tissue collected at 24 and 72 h PI, respectively. Moreover, pathological lesions including curling of toes, and congestion of embryos and embryonic lung tissue was observed at 24 h PI. Furthermore, at 72 h PI there was substantial difference in growth and development of virus infected and non-infected embryos. These data confirm a progressive increase in IBV viral gene expression within infected embryonic lung tissue.

In each set of experiment (n = 2), the lung tissues from virus infected (n = 3) and sham inoculated (n = 3) (sham inoculum was 200 l of sterile chorioallontoic fluid) embryos were pooled and collected in ice cold TRIZOL (Invitrogen, Life technologies, Carlsbad, CA) at 6, 24, 48 and 72 h post infection (PI). Within 30 min of collection, the tissues were homogenized (Polytron, PT MR 3100, Kinematica, AG, Switzerland) and stored at −80 • C until used.

The chicken-turkey array used for hybridization in these studies was an extended version of the avian array used by Munir and Kapur (2003) . In brief, the array was generated from normalized cDNA subtraction libraries derived from avian pneumovirus (APV, sub type C) infected chicken embryo fibroblast (CEF) cultures and bursa and thymus tissues obtained from APV infected turkeys subtracted with uninfected cultures and tissues, respectively. The chicken-turkey cDNA array is composed of 1247 chicken and turkey genes (with 3741 spots on array, i.e. each transcript was spotted in triplicate). After subtracting the redundancy of cDNAs derived from 56 genes, these arrays represent 1191 unique chicken and turkey genes. As per a recent public database blast search of the 1191 unique expressed sequence tags (ESTs) spotted onto the array there are 753 genes with no known orthologs in public DNA databases (Munir, unpublished data) .

Total cellular RNA was extracted from tissue homogenized in TRIZOL according to the manufacturer's instruction. The RNA was quantified and checked for quality, using an Agilent RNA 6000 nano kit 5065-4476 (Agilent Technologies, 7637 Waldbronn, Germany). Five to ten micrograms of good quality total RNA (with 28s/18s ratio > 1.7) was used to synthesize cDNA probes for hybridization onto a chicken array. The cDNA was prepared by using oligo dT primers and superscript II (Invitrogen, Life technologies, Carlsbad, CA) reverse transcriptase. The cDNA from infected and control samples were labeled with the monofunctional dyes Cy3 and Cy5 (Amersham, Piscataway, NJ). Probes were hybridized onto arrays at 65 • C for 8 h. The buffers and the protocols for hybridization and the post-hybridization processing were reported previously (http://www.agac.umn.edu/microarray/protocols/protocols.htm, Munir and Kapur, 2003) . All hybridizations were performed under optimized and uniform conditions. The complete experiment including embryo infection and RNA isolation was performed twice. Probe synthesis and microarray hybridization was repeated at least three to four times per time point within each replicate experiment.

The spot intensities of the hybridized images were obtained by using an Affymetrix scanner (428 Array scanner) and the Jaguar software version 2.0. Before importing the gene identifiers and the gene description, the scanned data files were analyzed for their quality by observing individual image spot quality and by subtracting the local background fluorescence from fluorescent intensity of each spot for both the Cy3 and Cy5 channels. Spots with background intensity higher than either dye were flagged as absent and discarded. The replicate spots that had infected/control intensity ratios of 2 or more standard deviations higher than the mean intensity ratio were also discarded. For quantitative analysis and further normalization, the data sets were analyzed by GeneSpring (Silicon Genetics, Redwood city, CA) software version 6.1. The entire data set was normalized by intensity dependent (Lowess) normalization.

Genes with known critical functions in host pathogen interaction and genes previously characterized in similar studies (Munir and Kapur, 2003) were selected for further validation. Differential expression of selected genes was validated for all four time points PI by real-time quantitative RT-PCR (qRT-PCR) using the SYBR Green based detection system. The specificity of the amplified products was determined by visualizing on 2% agarose gel.

Of the 1191 genes on the array, a total of 327 genes (27%) were altered by two-fold or more at one or more time points between 6 and 72 h PI. Of 200 transcripts that belonged to the category of unknown genes, 177 (88%) and 23 (12%) were up and down-regulated by two-fold or more, respectively. In contrast, of the 127 annotated genes (39% of differentially expressed genes), 111 (87%) and 16 (13%) were two-fold or more up and down-regulation, respectively. The analysis of altered genes revealed expression changes in many vital functional classes including kinase signaling, interferon stimulated genes (ISGs), chemokines and adhesion molecules, vesicular protein trafficking and fusion proteins genes, host genes involved in viral RNA synthesis, extracellular matrix protein genes, and genes belonging to cell cycle, metabolism, physiology and development, RNA binding, transcription and translation regulation, and ubiquitination. A partial list of differentially expressed genes and a graphical representation of kinase gene expression is shown in Table 1 and Fig. 1 , respectively. A comparison of the gene expression profile for IBV infection in ovo (Table 2 ) and APV infection of chicken embryo fibroblast (CEF) cells in vitro (Munir and Kapur, 2003 ) revealed a list of 14 genes (with >2-fold change in expression) that were common to both viral infections (Table 2 ). This comparison was completed for genes differentially expressed at 24 and 48 h PI. Gene expression changes at 6 and 72 h PI in IBV infection were not compared with APV since similar time points were not available for APV infection.

To confirm results obtained in array experiment by another independent method, 14 genes were selected for realtime RT-PCR assay at all time points PI. The results obtained are shown in Table 3 . The genes selected for real-time PCR belonged to several important functional groups including interferon induced genes, antigen presentation, chemokine and adhesion molecules and transcription factors protein genes. Moreover, the expression of these genes was either up-regulated, down-regulated or unchanged from 6 to 72 h PI. RT-PCR data show 85-90% agreement (except at 72 h PI) with the expression patterns seen with microarray analysis. Thus, confirmation rate is similar to previous reports (Munir and Kapur, 2003) . Though the trends of increase and decrease in gene expression at different time points for all genes are consistent for both the assay systems. We noticed a wide gap in fold changes expression of interferon induced protein 6-16 precursor and alpha interferon induciable protein (fragnment) p27-h genes in microarray and RT-PCR assay. We assume that variation in expression levels detected by microarray and RT-PCR may be attributed to differences in the sensitivity and the functional kinetics of both assay systems.

The set of kinases and phosphatase genes showing significant and stable changes in expression during IBV infection of embryonic lungs, include 14-3-3 protein zeta/delta, the protein kinase C inhibitor protein-1 (KCIP-1) or fac- tor activating exoenzyme S (Fas), protein kinase (similar to serine/threonine protein kinase (PKR), and protein-tyrosinephosphatase alpha (Table 1 and Fig. 1 ). Amongst the eight kinase genes spotted on the array, the gene coding for Fas showed the highest and most stable elevation in expression. The gene was 70-fold up regulated at 6 h PI and gradually declined to 22-fold at 72 h PI. 14-3-3 protein represent a novel class of kinases that recognize phosphoserine motifs on their ligand and play an important role in regulating various cellular processes including signal transduction, cell cycle, apoptosis, stress responses and viral and bacterial pathogenesis. These cellular activities are regulated through an interaction of over100 cellular proteins with Fas protein (van Hemert et al., 2001) . Proteins binding to 14-3-3 protein are very diverse and include kinases, phophatases, receptors, structural proteins and transcription factors. In general, most 14-3-3 protein ligands require a phosphorylation step before they can interact with their receptor and thus intracellular kinase/phosphatase signaling networks play a key role in Fas ligand interaction (Fu et al., 2000) . The role of 14-3-3 protein in IBV infection is not known. However, recently it has been reported that during hepatitis C virus (HCV) infection, viral C protein directly interacts with 14-3-3 and leads to activation of Raf-1 protein kinase signaling pathway (Aoki et al., 2000) . Activation of serine/threonine kinase and phosphatase genes along with the Fas gene in our data set may suggest their important role during IBV infection. The induction of dsRNA dependent PKR as a mediator of antiproliferative and antiviral activity of IFN is well known. However, recently the enhanced expression and increased activity of PKR in certain viral infections has been associated with the activation of specific viral proteins that interact with 14-3-3 protein (Schulze zur Wiesch et al., 2003) . Another kinase related gene that showed 61-fold upregulation at 6 h and remained up-regulated up to 48 h PI was lymphocyte specific protein 1 (LSP1) gene. Since the regulatory activity of KCIP-1 is mediated through PKC, it is interesting to note that LSP1 acts as a substrate for PKC and MAPKAP kinase 2 (Carballo et al., 1996) . LSP1 is an intracellular phosphoprotein that is expressed in B-lineage cells, functional T cells, thymocytes, macrophages and neutrophils derived from fetal liver and bone marrow. LSP1 is an important signaling molecule, regulating the cytoskeletal architecture, cell motility, receptor-induced apoptosis and a negative regulator of neutrophil chemotaxis (Jongstra-Bilen et al., 2000) . Alterations in expression of LSP-1, and KCIP-1 genes in our data may suggest their involvement in IBV infection and modulation of immune cell function. However, the precise interaction of these factors with each other or with viral proteins in IBV infection needs further investigation perhaps by array based studies with a larger number of gene sets specifically comprising all or the majority of genes involved in kinase and phosphatase signaling cascades. Further studies are anticipated to provide more information regarding the role of these signaling factors in IBV pathogenesis

RNA synthesis in positive strand RNA viruses is associated with viral replication complexes and is mediated through viral RNA-dependent RNA polymerase (RdRP). RdRP is the key enzyme for coronavirus replication and its sequences is conserved among coronavirus isolates (Lu et al., 2004) . However, the formation of the replication complex and recruitment of host cell components to complete viral replication are not fully understood (Brockway et al., 2003) . In general, the specificity of RdRP to select viral RNA as a template is dictated by cellular factors including tubulin, actin and heat shock proteins (Oglesbee et al., 1996) . The analysis of cellular factors involved in replication of RNA viruses such as Q␤ phage, brome mosaic virus (BMV) and tobacco mosaic virus (TMV) has shown that the RdRP holoenzyme can direct the synthesis of positive strand RNA from negative strand RNA. However, the recruitment of RdRP to the transcriptional promoters of positive strand viral RNA to synthesize negative RNA strand requires the interaction of cellular factors like TATA-binding protein (TBP), other general transcription factors and eukaryotic translation initiation factor 3 (eIF3) (Novak and Kirkegaard, 1994) . Moreover, the interaction of host elongation factor EF-1 and ribosomal S1 protein with RdRP complex to initiate transcription of polio virus RNA is also well documented (Lai, 1998) . Increased gene expression for some of these host factors following IBV infection suggests their involvement in IBV infection. For instance, chicken elongation factor 1 (EF-I) and TBP genes were two-and seven-fold up-regulated, respectively, at different time points following IBV infection (Table 1 ). In contrast, the eIF3 subunit 6 gene showed marked down regulation (>3-fold) at 6 h PI and then gradual up-regulation (two-fold) at 24 h PI. However, other subunits of eIF3 present on the array (e.g. homolog of mouse eIF3 subunit 7) showed no change of expression during IBV infection. In accordance with findings reported elsewhere (Guo et al., 2000) the down regulation of the eIF3 gene with simultaneous activation of the IFI 56 K gene in our data (Table 1 ) may possibly be a part of the interferon-induced antiviral pathway. Moreover, four-and five-fold upregulation of the 60s ribosomal protein 13 and ribosomal protein s6, respectively, at 6 h PI in IBV infection may also play a role in viral pathogenesis. There are examples where the antiviral effect of certain plant proteins is mediated through inhibition of ribosomal protein l3 (Rajamohan et al., 2001) .

Heterogeneous nuclear ribo-nucleoprotein 1 (hnRNP1) is another cellular protein that is reported to be involved in corona virus replication. Biologically, hnRNP1 is involved in pre-mRNA splicing in the nucleus and translational regulation in the cytoplasm (Shi et al., 2000) . The eight-fold up regulation of polypyrimidine tract-binding protein (PTB) or heterogeneous nuclear ribo-nucleoprotein 1 (hnRNP1) gene at 24 h PI suggests a possible role for hnRNP 1 in IBV transcription. Several lines of evidence are consistent with this hypothesis. For instance, the interaction between viral 3 -UTR and cellular hnRNP1 is a necessary requirement for the replication of the viral genome in bovine coronavirus (BCV) and mouse hepatitis virus (MHV) transcription (Lai, 1997; Shi et al., 2000; Spagnolo and Hogue, 2000) . Moreover, during MHV infection hnRNPA1 mediates the leader-body fusion event in synthesis of viral sub genomic RNA transcripts (Shen and Masters, 2001) . The binding of hnRNP-1 to the internal ribosome entry site on viral RNA of poliovirus, hepatitis A and hepatitis C virus is required for viral mRNA translation (Chung and Kaplan, 1999) . However, in a nonsegmented negative strand RNA virus infection (e.g. APV) the up-regulation of hnRNP1 gene was observed at later stages (96 h PI) of infection (Munir and Kapur, 2003) . We hypothesize that this difference in hnRNP1 gene expression profile for IBV and APV may be due to the different transcriptional strategy (synthesis of nested sets of viral mRNA) of coronavirus.

The induction of chicken interferon (IFN) by different IBV strains greatly varies in different culture environments (Holmes and Darbyshire, 1978) . It has been suggested that IBV's ability to induce IFN is linked to the virulence and adaptability of IBV strain for a particular host system (Holmes and Darbyshire, 1978; Otsuki et al., 1987) . In the present investigation we were unable to detect the induction of IFN genes due to the absence of IFN coding cDNA on the array. Interferons mediate their effect by binding to cell surface receptors which activate members of the JAK kinase family of proteins. Activated JAK kinases phosphorylate STAT factors. The STAT proteins homo-or heterodimerize and form complexes with other transcription factors to activate transcription of interferon stimulated genes (ISGs) (Haque and Williams, 1998) . Amongst the STAT group of ISGs, IBV infection of chicken embryos induced up regulation (two-to three-fold) of STAT 1 and another similar gene in lung tissue at 6 h PI. However, at 48 h PI only the STAT1 gene showed three-fold up-regulation. Similar to many other viral infections the activation of STAT1 seems to be an important component of host innate antiviral activity in coronavirus infection. For instance, STAT-1 null mice show enhanced viral susceptibility, greatly reduced responses to IFN-␣ and IFN-␥ and compromised innate immunity. Moreover, inactivation of the STAT-1 gene following Sendai virus infection (through an interaction to the viral C protein) is believed to be an important viral strategy to evade host innate immune responses (Garcin et al., 2002) .

In many viral infections IFI-56K, a member of the interferon induced proteins with tetratricopeptide repeats, is considered to be the most strongly induced gene amongst the known ISGs (Sen, 2001) . The marked up regulation of IFI-56K gene, as much as 22-fold in IBV infected chicken embryonic lung (Table 1) , is similar to what was observed in many other viral infections including APV infection of CEF cells, HCV and cytomegalovirus infection (Bigger et al., 2001; Munir and Kapur, 2003) . The inhibition of protein synthesis by IFI-56K protein is mediated through its interaction with the translation initiation factor eIF3 making it unavailable for protein translation (Guo et al., 2000) . Interestingly, the up-regulation of IFI-56K gene and a simultaneous reduction in eIF3 protein gene expression during IBV infection at 6 h PI may implicate strong activation of host innate immune response during the early phase of IBV infection (Table 1) . However, subsequent up-regulation of the eIF3 gene with its involvement in viral RNA synthesis, may reveal that effective viral strategies exist for escape from interferon's antiviral effects mediated through activation of IFI-56K gene expression. The overall ISG expression profile in this study and previous studies indicate that stable and marked activation of IFI-56K may be crucial for effective host protection in many viral infections including IBV, APV and HCV (Zhu et al., 2003) .

Other ISGs related to IFI-56K, including genes encoding for IFN-induced 58 kDa protein and IFN-induced TPR4 protein, did not show changes in expression following IBV infection. However, these ISGs showed marked up-regulation in other avian respiratory viral infections (Munir and Kapur, 2003) . Interferons regulate gene expression in various functional classes including signal transduction proteins, chemokines and adhesion molecules, antigen processing and presentation, signal transduction and protein degradation systems (de Veer et al., 2001) .

The G proteins are the class of signaling proteins that mediate many intracellular functions including cytoskeletal remodeling, vesicle transport and growth (Evers et al., 2000; Willard and Crouch, 2000) . Activation of G-protein genes and genes related to G-protein signaling in response to IFNs treatment of host cells have been reported in other microarray studies (de Veer et al., 2001) . However, the precise mechanism of gene activation is not clear. Moreover, the activation of G proteins in response to the binding of chemokines to transmembrane G-protein receptors and the requirement of JAK/STAT and G-protein interaction for chemotaxis and Ca flux responses in host cells is well characterized (Soriano et al., 2003) . In our studies, the up-regulation (17-fold at 6 h PI) of a gene similar to the G protein gene may occur as part of the IFN or chemokine induction pathways mentioned above (Table 1) . Further analysis of genes that were altered over two-fold and were common in IBV and other avian viral respiratory infection including APV showed marked differences in the magnitude of expression in the 6-16 precursor protein gene and STAT1 gene (Table 2 ). Microarray analysis of ISGs in non-infected cells have shown 31-and 21-fold increased expression of MxA protein in response to IFN-␤ and IFN-␣ treatment, respectively (Der et al., 1998) . Moreover, the Mx protein gene is considered to be a predominant ISG in the majority of negative strand RNA virus infections such as APV, influenza virus, vesicular stomatitis virus and measles virus infection (Atreya and Kulkarni, 1999; Munir and Kapur, 2003; Pavlovic et al., 1990; Schnorr et al., 1993 ) and a few positive strand RNA viruses such as Toga viruses (Landis et al., 1998) . However, there was no change in expression of the Mx protein gene in IBV infected embryos. The absence of change in Mx gene expression is supported by recent findings in other coronavirus studies indicating that the antiviral effect of IFNs in SARS corona virus infection in vitro is not mediated by MxA protein (Spiegel et al., 2004) . The induction of STAT-1 and IFI-56K gene expression supports the conclusion that IFN pathway was activated. However, no change in Mx gene expression is intriguing and may possibly be an outcome of an IBV anti-IFN virulence mechanisms.

IBV infection of embryonic chicken lungs resulted in a five-fold increase in 9E3 gene expression at 24 h PI ( Table 1 ). The 9E3 and K60 genes are two of four chicken chemokines known to belong to the CXC family. The 9E3 gene is a chicken homolog of human interleukin 8 (IL8). The role of IL8 gene induction in viral infections is well known. For example, cytomegalovirus (CMV) and Epstein Bar virus (EBV) infections activate the transcription factors NF-kB and AP1 which in turn induce interleukin 8 (IL8) expression and thus enhance the attraction of neutrophils,which serve as the target cell for CMV infection (Murayama et al., 1997; Penfold et al., 1999) . In addition to the attraction of heterophils (avian neutrophils), monocytes and lymphocytes, 9E3 is involved in growth regulation, wound healing and angiogenesis (Sick et al., 2000) . A five-fold increase in expression of the 9E3 gene in IBV infected embryonic lungs at 24 h PI (Table 1) may be linked to the common patho-physiological and microbiological states observed in IBV infection. For instance, in many respiratory viral infections including IBV, APV and human respiratory syncytial virus, peak accumulation of heterophils in lung lavage has been reported at 24 h PI (Kotani et al., 2000) . Similarly, peak IBV replication in the lungs of 18-day-old chicken embryos inoculated with IBV via the chorioallontoic route was observed at 24 h PI (Kapczynski et al., 2002) . Elevated 9E3 expression at 24 h post IBV infection in ovo may be correlated to the occurrence of peak IBV titers in lungs at 24 h PI and represent a conserved host response to viral infection. We did not detect significant alterations in other chemokine genes including K60, complement factor C3F and ornithokinin receptor genes which were present on the array. Unaltered expression of these chemokine genes may be an IBV strategy to escape possible adverse effects of chemokines. For example, it is well known for influenza virus infection, that an interaction between regulatory sequences in the viral nucleoprotein gene and host cell factors inhibits neutrophil chemotaxis and activation (Cooper et al., 2001) .

Differences in chemokine expression pattern between IBV (Table 1) and APV infection (Munir and Kapur, 2003 ) may contribute to differences in the patho-physiology of APV (a negative strand RNA virus) and IBV (a positive strand RNA virus) infection. The type of facial edema, infiltration of inflammatory heterophils, exudation and edema of sinuses is much more severe during APV infection of turkeys than during IBV infection. These differences may be the outcome of more extensive chemokine gene upregulation during APV infection. Moreover, unaltered expression of some ISGs, such as RI58, TPR-4, and Mx, and down-regulation of STAT I gene at 24 h post-IBV infection may be attributed to the antagonistic effect of IL-8 on IFN-␣ induction. Such an antagonistic activity for IL-8 has been reported for a majority of positive stranded RNA viral infections including polio, encephlomyocarditis and some DNA viruses such as herpes simplex virus 1 (HSV-I) infections. Similar inhibitory activity for IL-8 has not been observed in negative strand RNA viruses such as vesicular stomatitis virus and APV (Khabar et al., 1997a,b; Munir and Kapur, 2003) .

Another gene (LSP1) coding for a negative regulator of neutrophil chemotaxis also displayed marked and stable upregulation during IBV infection. This gene is believed to participate in the inhibition of chemokine and complement gene induction (Jongstra-Bilen et al., 2000) . Interestingly, chemokines other than 9E3 and complement factor genes (CF3) did not show any change in expression in our data. Further studies are needed to evaluate the impact of LSP1 on cytokine production and virus replication to determine if LSP1 plays a critical role in IBV pathogenesis

In contrast to many well studied enveloped RNA viruses, coronaviruses have a different budding mechanism. Coronaviruses acquire their membrane envelop by budding into the lumen of the Golgi and pre-Golgi compartments with virus release restricted to certain areas of the cell (Tooze et al., 1987) . The reason for this intracellular budding and selection of specific sites for budding of the virus is not known (Corse and Machamer, 2003) . Viral genomic factors, including Golgi targeting signal sequences on envelop proteins of IBV, involved in the process of virus egress have been well characterized. However, the cellular counterparts interacting in this process are not well known (Corse and Machamer, 2003) . In our experiments, alterations in the expression of seven genes with functions in the vesicular proteins (Vps) trafficking including clathrin light chain and transitional endoplasmic reticulum atpase genes, showed greater than twofold increased expression within 6 h PI. Signal recognition particle 54 kDa protein gene and Vps45 showed six-and twofold up regulation at 24 PI, respectively, but Vps45 gene showed two-fold down-regulation at 48 h post-IBV infection (Table 1 ). In addition, the TC10 protein gene and the COP1 gene related to the Vps pathway also showed two-and eight-fold enhanced expression at 48 and 72 h PI, respectively (Table 1 ). The Vps 45 and signal recognition particle 54 kDa proteins are thought to play important roles in vac-uolar protein sorting from golgi to the vacuole (Miller et al., 1995) .

The TC 10 protein is a member of the Rho family of GTP binding proteins that is closely related to Cdc42. Recent studies have shown that these proteins function as molecular switches linking cell surface receptors and extracellular signals to the regulation of actin dynamics (Hall, 1998) . The interaction between TC10 and coatmer protein COP1 causes massive Golgi membrane actin polymerization and thus play an important role in the control of perinuclear actin dynamics (Kanzaki et al., 2002) . The up-regulation of both of these genes (TC10 and COP1) in IBV (Table 1) infection suggests a possible role in the transport of coronavirus proteins through the golgi network (Munir and Kapur, 2003) . The involvement of different COP proteins in vacuolar transport of viral proteins of different viruses has been previously reported. For example, upregulation of coatamer complex gene COP9 is reported in APV infection (Munir and Kapur, 2003) . The M glycoprotein of coronavirus is one of the viral proteins that interacts with cellular factors and is accumulated in the golgi region and is believed to possess a signal for its retention at this location. The induction of cellular genes including clathrin light polypeptide, COP1 and GTP binding protein TC10 genes in IBV infection may suggest the interaction of these cellular proteins with viral proteins (like the M glycoprotein of the coronavirus) in the process of virus assembly and exocytic pathway (Gottesman and Maurizi, 1992) .

The genes related to ubiquitin and the proteasomal degradation pathway that showed over two fold alteration in expression included deubiquitinating enzyme 18 (hubp43) and a gene similar to the 26s proteasome non-ATPase regulatory subunit gene (Table 1 ). Both genes are over-expressed within 6 h PI. Protein ubiquitinylation is a classical pathway to target proteins for degradation by the 26S proteasome (Khor et al., 2003) . The up-regulation of the ubp18 gene and its role in viral replication and evasion of host immunity has also been reported in other viral infections including porcine reproductive and respiratory syndrome (PRRS) virus, CMV and APV. It is believed that the induction of ubiquitin carboxylterminal hydrolase 18 (Ubp18) gene in the early phase of infection may be a viral strategy to prevent degradation of its nascent proteins (Munir and Kapur, 2003) . However, unchanged expression (during IBV infection) of genes encoding the set of enzymes including E1 (Ub activating), E2 (Ub conjugating) and E3 (Ub substrate ligase), which are found to be regulated in certain other viral infections including APV (Munir and Kapur, 2003) , may be consistent with the substrate specificity of these enzymes and may suggest the activity of some other ubiquition pathway during IBV infection. Recently, a number of ubiquitin-like modifier (Ubls) proteins were discovered to have a variety of ubiquitin-like functions, including substrate conjugation (Haas and Siepmann, 1997 ). In addition, there also exists a class of ubiquitin related proteins in which an ubiquitin domain is built into a larger polypeptide and is not excised or attached to other proteins. These domains may impart properties anal-ogous to those conferred by a transferable Ubl called UbD (Hochstrasser, 2000; Munir and Kapur, 2003) . The absence of expression changes in other ubiquitin genes and activation of ubp18 gene and gene similar to 26s proteosome may suggest the involvement of ubiquitin related proteins like Ubl or UbD (that have not yet been identified in avian species) during IBV infection.

Over two-fold alterations in six genes belonging to the regulation of cellular transcription family were observed. These included transcriptional repressor protein Yin Yang-1 (YY1), TAF12 RNA pol II, human nuclear ribonuclear protein 1 (hnRNP1), similar to enhancer factor I (EF-1), E1a binding protein p300 and PCAF-associated factor 400 (Table 1) . Our data indicate a greater than seven-fold up regulation of the YY1 gene at 24 h PI. The YY1 protein is a multifunctional zinc finger protein that regulates the expression of various cellular as well as viral transcription factors including TBP, TFIIB, sequence specific DNA binding transcriptional activator (Sp1), c-myc, ATF/CREB, transcriptional co-regulators such as E1A, CBP, p300, PCAF, histone deacetylase 1 (HDAC1), HDAC2, and HDAC3, c-Fos, p53, ␣ actin, ␤ and ␥ interferon, P5 of adeno-associated virus, E6 and E7 genes of human papillomavirus, and a number of other viral long terminal repeats including Moloney murine leukemia virus and human immunodeficiency virus type I (Austen et al., 1997; Yao et al., 2001) . The role of YY1 protein in inhibition or propagation of virus infection is documented in various virus infections including adeno-associated virus replication where interaction of YY1 protein with E1A at promoter P5 relieves YY1 dependent suppression of the promoter. However, in the absence of E1A, the binding of YY1 to its cis elements represses transcription (Shrivastava and Calame, 1994) . The binding of YY1 to the upstream region of the MuLV promoter down-regulates MuLV promoter activity (Flanagan et al., 1992) . Moreover, it has been shown that an interaction between c-myc and YY1 protein suppressess the activity of the CXCR4 promoter where, CXCR4 serves as a major co-factor for T-cell tropism of HIV-I (Bleul et al., 1996) . However, mutations in the YY1 binding site removes suppression and facilitates HIV propagation (Moriuchi et al., 1999) . The role of the YY1 protein in coronavirus infection remains to be investigated. On the basis of alterations in expression of YY1 gene (>7-fold upregulation at 24 h PI) and other transcription factors including p300, PCAF, TBP and hnRNP we hypothesize that YY1 protein interacts with different cellular and viral elements during IBV infection. In addition to the role of hnRNP1 protein discussed earlier the alteration in hnRNP1 gene expression is reported to be involved in cellular immunobiology. It represses the expression of c-fos and zif268 genes and activates transcription of Ha-ras and SP6 kappa promoters that participate in host immune responses (Bemark et al., 1998; Mikheev et al., 2000) .

Transcription initiation by RNA polymerase II largely depends on TBP and associated factors (TAFIIs). Proteins like p300 and PCAF are transcriptional co-activators that are said to be involved in various cellular activities including cell growth, cellular transformation and differentiation, transcription initiation and developmental processes (Mitsiou and Stunnenberg, 2003) . On the basis of p300 interaction with cellular transcription factor TFIIB and RNA pol II, this protein is thought to play a direct role in transcription initiation thus serving as a molecular bridge between upstream transcriptional activators and basal machinery (Felzien et al., 1999) . The CH3 region of p300 has been reported to interact with many viral and cellular transcription factors like adenovirus E1A, SV40 large antigen and transcription co-activator PCAF (Arany et al., 1995; Bhattacharya et al., 1996) . The interaction between the CHI region of p300 and or CREB-binding proteins and the carboxy terminal region of STAT2 (an ISG) is supposed to be essential for the antiviral activity of STAT2. Hence, the targeting of p300 by EIA of adenovirus abolishes this antiviral effect of STAT2 and indicates a significant role for this protein in cellular antiviral responses. Moreover, it was recently reported that interaction between P300 and TBP plays an important role in transcription and embryogenesis in yeast (Walker et al., 2004) . In some recent studies p300 and PCAF have been found to have intrinsic histone acetyltransferase (HAT) activity and the ability to recruit additional HAT factors (Ji et al., 2003; Schiltz and Nakatani, 2000) . HATs have important role in chromatin remodeling and transcription (Baluchamy et al., 2003) . The activation of p300, TBP and PCAF genes at 6, 24 and 48 h PI during IBV infection (Table 1 ) may play such a role in the infection process.

Alteration of host translational machinery to optimize viral protein synthesis is an important viral strategy to establish a productive viral infection. For this purpose multiple biochemical mechanisms operating at the initiation phase of mRNA translation are altered in virus infected cells. For example, eukaryotic initiation factor 4a (eIF4a) is one among the DEAD box family of proteins that has shown marked induction during IBV infection (Table 1 ). The DEAD box proteins are highly conserved from bacteria to human and are thought to be involved in a wide variety of cellular processes including translation, ribosome biogenesis, RNA splicing, RNA maturation and degradation (Owsianka and Patel, 1999; Wassarman and Steitz, 1991) . The eIF4A is an ATP dependent RNA helicase that facilitates binding of the 40S ribosome to mRNA by melting 5 proximal secondary structure on cellular mRNA (Pause and Sonenberg, 1992) . We found seven-and three-fold up regulation of the eIF4a and DDX gene at 24 h PI, respectively during IBV infection. The AT-Pase and helicase activities of these proteins are essential for replication of some viruses and thus represent attractive targets for intervention (Du et al., 2002) . The role of eIF3 down regulation in conjunction to the upregulation of IFI-56K has already been discussed with reference to our data and have also been well studied in other viral infections (Guo et al., 2000) .

Substantial changes in expression profiles of a number of unique cellular pathways and genes including the activation of 14-3-3 protein kinase mediated signaling pathway, cellular factors like hNRNP1, TBP, eIF3, and PTB that are supposed to interact with RdRP enzyme gene, exclusive and stable alteration of IFI-56 k gene expression in interferon induced genes, enhanced expression of 9E3 and LSP1 genes as chemoattractant regulators, up-regulation of YY1, P300, PCAF associated factor 400, TBP, TAFII genes as transcriptionnal regulators, EF1a, ribosomal proteins l3, ribosomal protein s6, eIF4a and DDX 3 genes in translation regulation suggest the pathobiological significances of these genes and pathways during IBV infection process. However, the biological involvement of these cellular factors during IBV infection has been based on assumptions derived from previously known cellular functions for these genes. Further detailed characterization of these factors, using other molecular approaches will validate their role in IBV infection and will enhance our understanding of respiratory viral pathogenesis and immunity. In addition, a search of regulatory motifs (such as TPR and GTP binding motifs) in unknown but significantly altered genes identified in this study (data not shown here) may lead to the identification of more IFN related genes in the chicken. Moreover, precise determination of the role of viral cis acting elements (like UTR and 3 -poly A tail) and their interaction with cellular factors such as hnRNP1 and kinase C may be helpful in identifying potential targets to interrupt the viral replication cycle.

Hepatitis C virus core protein interacts with 14-3-3 protein and activates the kinase Raf-1

A family of transcriptional adaptor proteins targeted by the E1A oncoprotein

Respiratory syncytial virus strain A2 is resistant to the antiviral effects of type I interferons and human MxA

Characterization of the transcriptional regulator YY1. The bipartite transactivation domain is independent of interaction with the TATA box-binding protein, transcription factor IIB, TAFII55, or cAMP-responsive element-binding protein (CPB)-binding protein

Repression of c-Myc and inhibition of G1 exit in cells conditionally overexpressing p300 that is not dependent on its histone acetyltransferase activity

Purification and characterization of a protein binding to the SP6 kappa promoter A potential role for CArG-box binding factor-A in kappa transcription

Cooperation of Stat2 and p300/CBP in signalling induced by interferon-alpha

DNA microarray analysis of chimpanzee liver during acute resolving hepatitis C virus infection

The lymphocyte chemoattractant SDF-1 is a ligand for LESTR/fusin and blocks HIV-1 entry

Characterization of the expression, intracellular localization, and replication complex association of the putative mouse hepatitis virus RNA-dependent RNA polymerase

Characterization and purification of a protein kinase C substrate in human B cells Identification as lymphocyte-specific protein 1 (LSP1)

Severe acute respiratory syndrome vaccine development: experiences of vaccination against avian infectious bronchitis coronavirus

Heterogeneous nuclear ribonucleoprotein I (hnRNP-I/PTB) selectively binds the conserved 3' terminus of hepatitis C viral RNA

Attenuation of interleukin 8-induced nasal inflammation by an inhibitor peptide

The cytoplasmic tails of infectious bronchitis virus E and M proteins mediate their interaction

Functional classification of interferon-stimulated genes identified using microarrays

Identification of genes differentially regulated by interferon alpha, beta, or gamma using oligonucleotide arrays

Comparative characterization of two DEAD-box RNA helicases in superfamily II: human translation-initiation factor 4A and hepatitis C virus non-structural protein 3 (NS3) helicase

Rho family proteins in cell adhesion and cell migration

Molecular epizootiology of infectious bronchitis virus in Sweden indicating the involvement of a vaccine strain

Specificity of cyclin E-Cdk2, TFIIB, and E1A interactions with a common domain of the p300 coactivator

Cloning of a negative transcription factor that binds to the upstream conserved region of Moloney murine leukemia virus

14-3-3 proteins: structure, function, and regulation

All four Sendai Virus C proteins bind Stat1, but only the larger forms also induce its mono-ubiquitination and degradation

Regulation by proteolysis: energydependent proteases and their targets

A new pathway of translational regulation mediated by eukaryotic initiation factor 3

Pathways of ubiquitin conjugation

Rho GTPases and the actin cytoskeleton

Signal transduction in the interferon system

Evolution and function of ubiquitin-like proteinconjugation systems

Induction of chicken interferon by avian infectious bronchitis virus

EID-2, a novel member of the EID family of p300-binding proteins inhibits transactivation by MyoD

LSP1 modulates leukocyte populations in resting and inflamed peritoneum

Small GTP-binding protein TC10 differentially regulates two distinct populations of filamentous actin in 3T3L1 adipocytes

Detection of in ovo-inoculated infectious bronchitis virus by immunohistochemistry and in situ hybridization with a riboprobe in epithelial cells of the lung and cloacal bursa

The alpha chemokine, interleukin 8, inhibits the antiviral action of interferon alpha

Interleukin-8 selectively enhances cytopathic effect (CPE) induced by positive-strand RNA viruses in the human WISH cell line

The ubiquitin-vacuolar protein sorting system is selectively required during entry of influenza virus into host cells

Epithelial cell kinetics in the inflammatory process of chicken trachea infected with infectious bronchitis virus

RNA-protein interactions in the regulation of coronavirus RNA replication and transcription

Cellular factors in the transcription and replication of viral RNA genomes: a parallel to DNA-dependent RNA transcription

Human MxA protein confers resistance to Semliki Forest virus and inhibits the amplification of a Semliki Forest virus-based replicon in the absence of viral structural proteins

Tissue distribution of avian infectious bronchitis virus following in ovo inoculation of chicken embryos examined by in situ hybridization with antisense digoxigeninlabeled universal riboprobe

Attenuation of SARS coronavirus by a short hairpin RNA expression plasmid targeting RNA-dependent RNA polymerase

CArG binding factor A (CBF-A) is involved in transcriptional regulation of the rat Ha-ras promoter

The beta subunit of the signal recognition particle receptor is a transmembrane GTPase that anchors the alpha subunit, a peripheral membrane GTPase, to the endoplasmic reticulum membrane

p300 is involved in formation of the TBP-TFIIA-containing basal transcription complex

USF/c-Myc enhances, while Yin-Yang 1 suppresses, the promoter activity of CXCR4, a coreceptor for HIV-1 entry

Regulation of host cell transcriptional physiology by the avian pneumovirus provides key insights into hostpathogen interactions

Transcriptional response of avian cells to infection with Newcastle disease virus

Human cytomegalovirus induces interleukin-8 production by a human monocytic cell line, THP-1, through acting concurrently on AP-1-and NF-B-binding sites of the interleukin-8 gene

Coupling between genome translation and replication in an RNA virus

The highly inducible member of the 70 kDa family of heat shock proteins increases canine distemper virus polymerase activity

Serological relationship among ten strains of avian infectious bronchitis virus

Hepatitis C virus core protein interacts with a human DEAD box protein DDX3

Mutational analysis of a DEAD box RNA helicase: the mammalian translation initiation factor eIF-4A

Resistance to influenza virus and vesicular stomatitis virus conferred by expression of human MxA protein

Cytomegalovirus encodes a potent alpha chemokine

Comparative genomics and host resistance against infectious diseases

Active center cleft residues of pokeweed antiviral protein mediate its high-affinity binding to the ribosomal protein L3

The PCAF acetylase complex as a potential tumor suppressor

MxA-dependent inhibition of measles virus glycoprotein synthesis in a stably transfected human monocytic cell line

The proteins of the Hepatitis C virus: their features and interactions with intracellular protein phosphorylation

Viruses and interferons

Evaluation of the role of heterogeneous nuclear ribonucleoprotein A1 as a host factor in murine coronavirus discontinuous transcription and genome replication

Heterogeneous nuclear ribonucleoprotein A1 regulates RNA synthesis of a cytoplasmic virus

An analysis of genes regulated by the multi-functional transcriptional regulator Yin Yang-1

Novel chicken CXC and CC chemokines

Chemokines integrate JAK/STAT and G-protein pathways during chemotaxis and calcium flux responses

Host protein interactions with the 3' end of bovine coronavirus RNA and the requirement of the poly(A) tail for coronavirus defective genome replication

The antiviral effect of interferon-beta against SARS-coronavirus is not mediated by MxA protein

Sorting of progeny coronavirus from condensed secretory proteins at the exit from the trans-Golgi network of AtT20 cells

14-3-3 proteins: key regulators of cell division, signalling and apoptosis

An extensive requirement for transcription factor IID-specific TAF-1 in Caenorhabditis elegans embryonic transcription

RNA splicing Alive with DEAD proteins

Nuclear and cytoskeletal translocation and localization of heterotrimeric G-proteins

Immunity to avian pneumovirus infection in turkeys following in ovo vaccination with an attenuated vaccine

Regulation of transcription factor YY1 by acetylation and deacetylation

Gene expression associated with interferon alfa antiviral activity in an HCV replicon cell line

This work was supported by the Functional Pathogenomics of Mucosal Immunity (FPMI) project funded by Genome Prairie, Genome BC and their industry partners, Inimex Pharmaceuticals and Pyxis Genomics Inc. We are thankful to Dr. Andrew Ficzycz for providing information about YY1 protein from his PH.D thesis. We owe thanks to the scientific staff of the FPMI research group, adenovirus research group and avian research group at VIDO for their valuable suggestions and co-operation in designing and conducting the experiments and the data analysis. We are thankful to Stacy Miskolczi and Barry Carrol from the Animal Care facilities at VIDO. Published with permission of the director of VIDO as VIDO journal series number 382.