Previous Article | Next Article ![]()
Journal of Virology, April 2003, p. 4899-4910, Vol. 77, No. 8
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.8.4899-4910.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Departments of Microbiology and Veterinary Pathobiology, and Biomedical Genomics Center, University of Minnesota, St. Paul, Minnesota 55108
Received 14 October 2002/ Accepted 16 January 2003
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Avian species are highly susceptible to respiratory tract infections due to the unusual anatomical and immunological features of their respiratory tracts (45). Despite the great morbidity and mortality associated with these infections, very little is known about the molecular basis of disease pathogenesis and the host's response to infection. Among the major pathogens associated with respiratory disease in birds is the avian pneumovirus (APV), a member of the Paramyxoviridae family of enveloped viruses and the causal agent of acute and highly contagious respiratory tract infections in poultry species and wild bird populations worldwide (2). Infected birds typically exhibit severe inflammatory changes such as swollen sinuses, nasal and ocular discharge, and respiratory distress. The disease is characterized histopathologically by acute respiratory tract inflammatory changes including diffuse cellular infiltration, hemorrhage, and dense submucosal infiltration with heterophils (29).
We describe here the application of suppression subtractive hybridization (SSH) in conjunction with cDNA microarray analysis to explore the alterations in host cell transcriptional program in an APV model of infection. The present study provides intriguing information on the changes in host gene expression upon infection with a pneumovirus, which could serve as a foundation for future investigations on the roles of various genes and pathways in the pathogenesis of this important disease in avian species.
| MATERIALS AND METHODS |
|---|
|
|
|---|
95% confluent; at this time they were infected with the Colorado strain of APV subtype C (5) at a multiplicity of infection of 0.25 PFU per cell in a volume of 10 ml of medium and were allowed to adsorb at 37°C for 1 h. Following adsorption, 25 ml of Leibovitz McCoy's medium was added, and flasks were incubated at 37°C. Cells used as uninfected controls were mock treated (with sterile media) and processed in a fashion similar to that of the infected cells. Cells were harvested at 2.5, 24, 48, and 96 h postinfection (p.i.) from APV- and mock-infected cultures, respectively. Purification of poly(A)+ mRNA. Adherent APV-infected and control cells were washed once with phosphate-buffered saline, followed by lysis and homogenization in Trizol reagent (Invitrogen Life Technologies, Carlsbad, Calif.). Total RNA was purified according to the manufacturer's instructions. Poly(A)+ mRNA was extracted from the total RNA by using Oligotex mRNA spin columns (Qiagen, Valencia, Calif.) per the manufacturer's instructions.
Construction of normalized subtracted cDNA libraries. A forward subtracted library was constructed in which poly(A)+ mRNA from APV-infected cells served as the "tester," or the population whose upregulated transcripts were to be identified, and poly(A)+ mRNA from mock-infected control cells was designated the "driver," or the population whose transcripts served as a reference for cDNA subtraction. APV-infected and mock-treated cell transcripts were also used conversely as driver and tester samples, respectively, to construct a reverse subtracted library. The differentially expressed cDNAs (targets) are present in the tester cDNA but are absent (or present at lower levels) in driver cDNAs. The SSH procedure (9) was followed to construct forward and reverse subtracted libraries by using the PCR-Select cDNA subtraction kit (Clontech, Palo Alto, Calif.). Briefly, 2 µg of the poly(A)+ mRNA from each of the tester and driver, collected at 96 h p.i., was converted to double-stranded cDNA by reverse transcription (RT). Both tester and driver cDNAs were digested with RsaI to produce shorter blunt-ended fragments. Only the digested tester cDNA was subdivided into two portions, each of which was ligated with a different adapter sequence, resulting in two populations of tester. After ligation, a series of two hybridization steps was performed. During the first hybridization, an excess of driver was added to each tester, denatured, and allowed to anneal. Due to the second-order kinetics of hybridization and the suppression PCR effect, the concentrations of high- and low-abundance target single-stranded molecules in the tester are normalized (43) and the target sequences in the tester become significantly enriched for differentially expressed genes. In the second hybridization, the two first-hybridization reaction products were mixed with each other and with fresh denatured driver cDNA. The populations of normalized and subtracted single-stranded target cDNAs anneal with each other, forming double-stranded hybrids with different adapter sequences at their 5' ends. The adapter ends were filled in with DNA polymerase, and the subtracted hybrid molecules were specifically amplified by nested PCR using adapter-specific primer pairs.
Cloning of the SSH libraries into the TA vector.
Amplified products from the secondary nested PCR, constituting the subtracted target cDNAs, were ligated with the pGEM-T plasmid vector (Promega, Madison, Wis.) and transformed into maximum-efficiency Escherichia coli DH5
cells (Life Technologies). Subsequently, transformed bacteria were plated onto Luria-Bertani agar plates containing ampicillin (100 µg/ml), 5-bromo-4-chloro-3-indolyl-ß-D-galactoside (X-Gal; 50 µg/ml), and isopropyl-ß-D-thiogalactoside (IPTG; 100 µM), followed by overnight incubation at 37°C. Random recombinant white colonies were selected and cultured in Luria-Bertani broth containing ampicillin (100 µg/ml), and plasmid extraction was performed with the QIAwell 96-well plasmid purification system (Qiagen). A total of 960 forward and reverse subtracted clones were sequenced, and the sequences were identified based on homology searches with public genetic databases (Swiss-Prot, TrEMBL, PIR, NRL-3D, and GenPept). Sequence data were analyzed and edited for quality and vector sequences by using Phred-Phrap/Consed analysis software.
cDNA microarray construction. A 2,950-element cDNA microarray chip containing 960 SSH clones was constructed as described elsewhere (34, 35). In brief, 576 forward SSH and 384 reverse SSH clone inserts were amplified by PCR using a primer pair corresponding to the flanking adapter sequences (Clontech). All PCR products were visualized on 1% agarose gels to ensure quality and adequate amplification, followed by purification with MultiScreen PCR plates (Millipore, Bedford, Mass.). PCR products were printed in triplicate onto poly-L-lysine-coated glass slides by employing a Microgrid II robot (BioRobotics, Boston, Mass.). Three bacterial genes (from Pasteurella multocida) and one plant gene (the Arabidopsis thaliana gene encoding the chlorophyll a/b-binding protein [Cab]) were included as negative controls, and total cellular cDNA was included as a positive control. All control elements were spotted 14 times each on the array.
Fluorescent probe labeling and hybridization. Poly(A)+ RNA (300 ng) purified from APV-infected and mock-infected cells at 2.5, 24, 48, and 96 h p.i. was reverse transcribed by using an oligo(dT)12-18 primer, deoxynucleoside triphosphates, aminoallyl dUTP, and Superscript II reverse transcriptase (Invitrogen Life Technologies). cDNAs from infected and control samples were labeled with the monofunctional dyes Cy5 and Cy3 (Amersham, Piscataway, N.J.), respectively, followed by hybridization with the spotted array at 67°C for 5 h. For all hybridizations, A. thaliana Cab gene mRNA (Stratagene, La Jolla, Calif.) was spiked into the cDNA synthesis reactions of both samples, and information from this mRNA was later used for data normalization. Prior to final hybridizations, hybridization conditions, such as the optimum amount of labeled probe, the hybridization buffer composition, temperature, and washing stringencies, were standardized to ensure maximum specificity and sensitivity, and all hybridizations were performed with the optimized conditions, at which there was no appreciable cross-hybridization with the control spots. The complete experiment, including primary culture of cells, virus infection, and RNA isolation, was performed twice, with probe synthesis and hybridizations performed at least twice per independent experiment. This plan resulted in 12 independent Cy5/Cy3 intensity ratio data points for each spotted cDNA at each time point. Images of the hybridized arrays were acquired by laser confocal scanning (Scanarray 5000; GSI Lumonics, Watertown, Mass.) and analyzed with Quantarray, version 3.0 (GSI Lumonics), and Spotfire Decision Site, version 6.5 (http://www.spotfire.com), software.
Data analysis. Prior to the import of gene identifiers and descriptions into the data set, the raw data were analyzed as follows for quality control: (i) individual spots were flagged for quality and discarded during analysis and quantification of spot intensities, (ii) local background fluorescence was subtracted from the fluorescence intensity of each spot for both the Cy5 and Cy3 channels, (iii) the entire data set was normalized for both channels based on the A. thaliana (Cab gene) spike control, (iv) spots with high background intensity for either dye were eliminated, (v) replicate spots that had Cy5/Cy3 intensity ratios 2 or more standard deviations higher than the mean intensity ratio were discarded, and (vi) replicate spot ratios were averaged, and differential expression ratios were determined. Further data analysis and visualizations of expression profiles were performed using Spotfire Decision Site software, version 6.5.
Differential expression analysis by real-time RT-PCR. The differential expression of a group of 20 selected genes was further validated by real-time quantitative RT-PCR using SYBR green-based detection (ABI) that was run on an ABI 7700 fluorescent sequence detection system (Perkin-Elmer, Foster City, Calif.). Selection of genes for further validation was based on (i) previously known critical functions of the encoded proteins in the context of the host-virus relationship and (ii) their representation of the whole spectrum of expression patterns from very strong induction to weak induction and significant repression by APV infection. The quality and specificity of amplified products were confirmed by visualization on a 2% agarose gel.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
cDNA microarray analysis revealed that the expression levels of 352 genes were significantly altered over the 96-h period. Of these genes, 268 transcripts showed an increase and 84 showed a decrease in expression levels; a partial list is given in Table 1. Transcripts representing several functional classes were perturbed by APV infection; these included many genes with known roles in blocking virus replication and others with a likely role in facilitating virus propagation. The major functional classes of genes that were transcriptionally altered during infection are summarized below.
|
The gene encoding the Mx protein was one of the most strongly induced (more than 60-fold) during infection. Mx proteins constitute a group of large antiviral GTPases that are key components of the IFN-induced antiviral defenses (39) and mediate their effect by sequestering viral nucleocapsids, thus rendering them inaccessible for replication (21). Owing to their significance, Mx proteins are evolutionarily conserved in all vertebrates analyzed so far, including chickens, mammals, and fish. Activity against orthomyxoviruses, paramyxoviruses, rhabdoviruses, bunyaviruses, and togaviruses has previously been demonstrated for human MxA, and activity against orthomyxoviruses has been demonstrated for mouse Mx1 (33). Although up-regulation of Mx expression has been observed during infection with a wide range of viruses such as hepatitis C virus, herpes simplex virus, and human cytomegalovirus (1, 26, 53), its antiviral potency in these cases remains to be determined. Similarly, it is not known if the avian Mx protein has an antiviral effect toward APV or other avian respiratory viruses. Sequence analysis suggests that the avian Mx gene contains an IFN-stimulated response element (ISRE), common to all ISG promoters, that is strongly activated upon IFN stimulation (39). The massive induction of Mx gene expression during the course of APV infection suggests its possible activity as an antiviral ISG, and thereby it could have an important role in APV pathogenesis.
The gene encoding the IFI-56K protein was up-regulated early (2.5 h) in infection and continued to accumulate over time. It is known that IFN creates an antiviral state, in part, by halting mRNA translation. The IFI-56K protein contains eight tetratricopeptide repeat (TPR) motifs that mediate protein-protein interaction by binding the P48 subunit of translation initiation factor eIF-3, thus inhibiting protein synthesis (14). By using the gene array approach, enhanced expression of IFI-56K was also found during infections with hepatitis C virus, herpes simplex virus, and cytomegalovirus (1, 26, 53). The considerable (6- to 32-fold) induction of the IFI-56K gene suggests a role for this protein in inhibiting the translation complex during APV infection. Interestingly, a gene encoding one of the eIF-3 subunits in the translational machinery predicted to be inhibited by IFI-56K showed increased (two- to sixfold) expression after the onset of the increase in IFI-56K transcription; this gene may represent a virus-induced countermeasure to overcome the IFI-56K-mediated block of eIF-3, an intriguing hypothesis that requires further investigation.
The analysis also showed that the genes encoding RI58 and IFN-induced protein with TPR 4 were up-regulated. Both of these are ISGs whose roles in the IFN pathway remain to be elucidated. Interestingly, however, both possess TPR motifs, have high homology with IFI-56K (50), and in our studies were found to be up-regulated with profiles similar to that of IFI-56K (Fig. 1B). It is therefore tempting to speculate that they might have a collaborative regulatory function similar to that of IFI-56K during infection.
|
/ß signaling pathway (6), was strongly induced by APV (>50-fold). STAT1 not only serves as an IFN signaling molecule, but its own expression is also induced by IFN (7). Although the IFN transcript was not isolated by SSH, up-regulation of STAT1, a hallmark of the IFN signaling cascade, suggests that the induction of ISGs was perhaps due to IFN stimulation. However, this does not exclude the possibility that induction of a few or all ISGs identified here may be partially attributed to APV infection itself and its double-stranded RNA replication intermediates by their own IFN-independent signaling pathways. Overall, the data suggest that the increased expression of STAT1 may account for the early and strong up-regulation of the ISGs observed soon after APV infection. It is interesting that PKR and 2-5A were not isolated here by SSH as differentially expressed genes. Either they were not induced during infection, were activated to a level below the assay detection limit, or were indeed isolated but are among the unknown differentially expressed cDNAs that have no homologue in the database, since the chicken sequences for these genes are not known.
Together, the data reveal massive induction of ISGs during APV infection. While the exact role of this enhanced ISG expression during APV infection is unclear, it is noteworthy that paramyxoviruses in general are known to block the IFN signaling cascade so as to circumvent IFN antiviral responses in a virus-specific fashion (10). Similarly, the human and bovine respiratory syncytial viruses (HRSV and BRSV, respectively), both close relatives of APV, although they do not inhibit IFN-
/ß or IFN-
production or signaling, do antagonize IFN-induced antiviral proteins such as Mx through the action of the viral NS1 and NS2 proteins and thus are resistant to the effects of IFN (37). Our analyses suggest that, like HRSV and BRSV, APV does not block IFN production or signaling, as indicated by up-regulation of several ISGs. However, unlike these other viruses, APV lacks the genes encoding NS1 and NS2 proteins and thus may be sensitive to the antiviral effects of Mx and other ISGs. Alternatively, if APV does possess any IFN resistance mechanism, it could conceivably operate as an event downstream of the IFN signaling cascade and at the level of action of ISG-encoded proteins. Future investigations to elucidate the effects of IFN-triggered proteins identified in this study on APV replication will be of immense significance for understanding their role in APV biology.
Proinflammatory chemokines, complement, and adhesion molecules. Chemokines are a superfamily of chemoattractant cytokines that recruit immune cells to the site of infection with viruses, bacteria, or fungi. Based on the number and relative positions of the conserved cysteine residues, they are classified into the CXC, CC, C, and CX3C subfamilies (32). Although well known for their immunomodulatory and chemoattractive activities, chemokines have been implicated in pathological injury caused by severe inflammation. Here we observed significant induction by APV of genes that are well-known mediators of proinflammatory responses. The genes encoding the proinflammatory CXC chemokine mRNAs, including 9E3 (chicken homologue of human interleukin 8 [IL-8]) (19) and K60 (42), two of the four chemokines known in chickens, were considerably up-regulated (53- and 18-fold, respectively). While these proteins have not been previously implicated in APV pathogenesis, their chemoattractant properties are in agreement with the type of cellular infiltrate present in the avian respiratory tract during infection. 9E3 and K60 are both CXC chemokines and strong chemoattractants for heterophils (the avian equivalents of neutrophils) (42), which are the predominant infiltrating leukocytes in the APV-infected respiratory tract. This observation is consistent with the HRSV-induced up-regulation of IL-8 expression (22, 52) and neutrophil infiltration seen during human infections; along with other chemokines, IL-8 is believed to be the cause of airway inflammation and immunopathogenic injury caused by neutrophils (47, 52). While the exact mechanisms governing the inflammation and recruitment of circulating heterophils into the respiratory tracts of APV-infected birds remain unknown, induction of 9E3 and K60 provides a plausible explanation for the extensive inflammation and swelling of the face, inflammatory heterophil infiltration, exudation, and edema of sinuses of birds during APV infection.
There was strong induction of the gene encoding the complement pathway protein C3 (2- to 12-fold induction over the time course). The complement system is a critical component of the host innate defense mechanisms against a variety of pathogens. It is composed of multiple factors that act in concert to kill virus-infected cells, bacteria, and parasites (reviewed in reference 11). The C3 complement is the most abundant protein, has leukocyte chemoattractant properties, and plays a pivotal role in initiating the complement cascade. The transcriptional activation and role of C3 complement in APV infection have not been reported previously. In human infections with HRSV, it has been shown that HRSV stimulates C3 production in infected cells, which mediates cell lysis through neutrophil-dependent cytotoxicity (20); this, in conjunction with chemokines, is considered to be the basis of infection-related pathology.
The results also show elevated expression (>2-fold) of the gene encoding the ornithokinin receptor, the only known chicken receptor for proinflammatory peptides known as kinins, which are responsible for changes associated with inflammation such as vasodilation, increased vascular permeability, and pain (38). Although further investigation is required to confirm the role of activated 9E3, K60, C3, and ornithokinin receptor genes in APV pathogenesis, the results strongly suggest that the inflammation and migration of heterophils into the respiratory tract triggered by these factors are likely to play an important role in APV-related pathology.
Ub-proteasome pathway.
Several genes belonging to the ubiquitin (Ub)-proteasome pathway, including those encoding the 26S proteasome subunit p42, COP9 complex subunit 4, and E3 ubiquitin ligase CUL-1, the gene similar to sequestosome 1, and the gene encoding deubiquitinating enzyme 18, were induced during APV infection. The functions of many of these genes have been well characterized in simple and higher eukaryotes. The Ub-proteasome pathway comprises a major proteolytic system that regulates the dynamics of proteins (17). Multi-Ub chains are added to proteins to tag them for degradation by the 26S proteasome, a multiproteolytic enzyme complex. Ubiquitination is accomplished by a set of enzymes including E1 (Ub activating), E2 (Ub conjugating), and E3 (Ub-substrate ligase). A gene similar to sequestosome 1 encodes a nonproteasomal Ub binding protein, p62, that binds preferentially with multi-Ub and forms a novel cytoplasmic structure called a sequestosome, which contains Ub-protein complexes and is implicated in Ub-mediated signaling and/or protein degradation by virtue of its role as a temporary storage place for Ub-protein complexes (41, 46). Induction of several Ub-proteasome pathway proteins by APV indicates enhanced activity of this proteolytic machinery during infection. In conjunction with the IFI-56K inhibition of protein translation described above, this may be a parallel or alternative host antiviral response to enhance proteolysis. This hypothesis is supported by the fact that several of the proteasome system genes are induced by IFN-
/ß stimulation (8). Alternatively, this response may serve to enhance the production of viral antigenic peptides for display in the context of major histocompatibility complex class I molecules to virus-specific T lymphocytes for triggering of an antiviral immune response.
Given the induction of the Ub system genes that mediate proteolysis, an intriguing observation is the up-regulation of the Ub protease (Ubp) or deubiquitinating enzyme during APV infection. Ubp is known to cleave Ub tags from proteins and render them resistant to proteasome degradation (17). Induction of this enzyme together with the gene similar to sequestosome 1 occurred by 2.5 h p.i., and expression remained elevated (maximum, 11.9-fold) throughout the experiment. While Ubps are substrate specific and the functional relevance of the observed up-regulation remains to be determined, the enhanced expression of Ubp can be interpreted in several ways. (i) It may be a virus-directed virulence mechanism to prevent Ub-mediated degradation of nascent viral proteins in order to ensure successful production of virus progeny. In this context, it is noteworthy that Ubp and IFI-56K have similar induction profiles (Fig. 1B). Both genes are up-regulated very soon after infection, at 2.5 h p.i., and remain activated throughout the experiment, suggesting that activation of Ubp may be a virus-induced response to stabilize proteins in the presence of an IFI-56K-mediated translation block. (ii) Since proteasome-generated peptides are presented in the context of major histocompatibility complex class I molecules on the surfaces of infected cells to T lymphocytes, Ubp activation may be an exquisite virulence mechanism to prevent antigen presentation and escape host immune recognition right from the beginning of host invasion. This could possibly be the basis of the poor adaptive immunity of poultry flocks, which repeatedly experience APV outbreaks with the same virus strain even during their short life spans and are more susceptible to secondary infections when exposed to APV. (iii) The increased activity of Ubp leads to higher levels of free intracellular Ub, which along with other vesicular-system proteins, plays a key role in the maturation and budding of a variety of viruses and could have similar functions in APV replication.
Interestingly, a variety of other virus types exploit the Ub-proteasome pathway to facilitate their replication and evade host immunity. For instance, porcine reproductive and respiratory syndrome (PRRS) virus (51) and human cytomegalovirus (53) infections cause up-regulation of Ubp expression, which in the case of the PRRS virus is thought to be a mechanism for escaping host immunity (51).
Vesicular protein trafficking and fusion. Enveloped viruses are packaged, mature, and exit from the host cell by budding (pinching off) from the cell membrane. Viral proteins and nucleocapsids utilize the host cell machinery for their transport to the plasma membrane. An interesting group of vacuolar protein sorting (Vps) genes showed altered expression upon APV infection. These include genes involved in protein sorting, trafficking, vesicle formation, membrane tethering, and fusion such as those encoding Vps45, Rab5c-like protein, syntaxin 7, clathrin heavy chain, signal recognition particle receptor, and the endocytic receptor Endo180. The genes encoding the Vps45 and signal recognition particle receptor were induced >5- and >2-fold, respectively, and are thought to play critical roles in vacuolar protein sorting from the Golgi apparatus through the endosomal pathway to the lysosome (vacuole in yeast) (4). The gene encoding the Rab5c-like protein was also substantially induced during APV infection. It functions in the transport of cargo in clathrin-coated vesicles to the early endosomes (36). Rab5c-like proteins are small GTPases that complex with other proteins, perform functions such as transport vesicle formation and tethering of vesicles at their target membranes, and are also known to recruit SNARE proteins (such as syntaxins) that are responsible for fusion of membranes at the Rab5-tethered sites (24). Endo180 functions as a recycling molecule and mediates the transport of glycosylated proteins to the endosomal compartment (40). The membrane-bound proteins are sorted by the Vps pathway for ultimate degradation in the lysosome, which serves as a mechanism for regulating cellular proteins. Syntaxin 7, in addition to fusing membranes, is involved along with Vps45 in the formation of multivesicular bodies, which are endosomal vesicles containing cargos brought in by the Vps machine for lysosomal degradation (27). It appears, therefore, that upon infection, several genes that are all known to serve various functions in the Vps pathway are up-regulated.
The role of the Vps and Ub pathways in the replicative biology of viruses is being increasingly recognized. Over recent years, fairly diverse viruses such as those in the retrovirus, filovirus, and rhabdovirus families have been witnessed to hijack and exploit the host Vps and Ub pathways for their packaging and egress from the cell (31). It is believed that the virus budding phenomenon is topologically similar to multivesicular-body formation in the Vps pathway (in both mechanisms, the membrane invaginates away from the cytoplasm) and that viruses therefore can usurp this machinery in an analogous fashion for viral bud formation and exocytosis (12). Though far from being completely understood, most of the information in this arena comes from studies on human immunodeficiency virus type 1 (HIV-1) and other quite divergent viruses (reviewed in reference 31). Exploitation of the Vps pathway for budding, however, requires a conserved functional consensus sequence known as the late assembly (L) domain. First identified in Rous sarcoma virus (48), L domains have been identified in many RNA viruses, and thus use of the Vps pathway proteins is believed to be a broadly applicable mechanism of virus budding. Although all L domains have the same budding function, the motif sequence is variable across viruses. The sequence motif PPxY in the Rous sarcoma virus Gag and vesicular stomatitis virus matrix (M) proteins and the sequence motifs PTAPP and YxxL in Gag proteins of HIV-1 and equine infectious anemia virus, respectively, have been identified as L domain consensus sequences (16, 31). The current model suggests that the L domain mediates ubiquitination of viral Gag or M protein by Ub ligases. The Ub-tagged protein then recruits Vps pathway proteins (such as Vps4 and Tsg101 in HIV-1) that ultimately mediate the budding process (31). The budding machinery is also known to require active proteasomes, free Ub (15, 16, 30), and the recruitment of the Vps-endosomal protein machinery to target viral proteins for final budding (12, 31). Given the use of an analogous mechanism for budding by fairly divergent virus families, it could conceivably be a pathway for virus egress from the host cell for other viruses as well that have not yet been studied from this perspective.
Importantly, in addition to Vps pathway protein transcripts, many genes with functions in the Ub pathway were also transcriptionally induced upon infection with APV. The molecular mechanism of budding in paramyxoviruses in general is virtually unknown, and the existence of the L domain also has not been reported for any of the family members. Aside from the knowledge that the paramyxovirus matrix (M) protein has a role in virus assembly, the molecular mechanisms that drive M protein-directed virus assembly at the limiting membrane remain undescribed. However, the M proteins of Sendai virus (a paramyxovirus) and human parainfluenza virus (a close relative of APV) (both classified under the same order, Mononegavirales) and that of influenza virus (an orthomyxovirus) have been shown to be capable of budding out of cells in the form of virus-like particles when expressed alone in the absence of other viral proteins (3, 13, 44), a feature reminiscent of the Gag and M proteins of retroviruses, rhabdoviruses, and filoviruses that is attributed to the functional L domain in their sequences (reference 15 and references therein). Significant evidence, therefore, suggests that M proteins of paramyxoviruses (3, 44) complete the budding process perhaps in a fashion similar to that of M and Gag proteins of rhabdoviruses, filoviruses, and retroviruses, which employ their L domains (16, 31) in recruiting host machinery for late stages of budding. Hence, paramyxovirus M proteins might also have a functional domain(s) that could be utilized in a similar manner.
In order to investigate if APV and other pneumoviruses possess any of the known L domain consensus sequences, we aligned the M protein sequences of all known pneumoviruses that infect a variety of host species. Our analysis shows that the M proteins of all pneumoviruses possess the L domain consensus sequence YxxL, located in the C-terminal region (Fig. 2). The putative L domain was found to span residues 197 to 200 with a consensus sequence of YAGL for APV-A, -B, and -C, HRSV B, human metapneumovirus, BRSV, ovine RSV, and pneumonia virus of mice. The sequence YSGL was present instead at the same location in HRSV A and had a single residue substitution of A to S at position 198 but still conformed to the YxxL L domain consensus. Therefore, based on (i) the fact that the M protein is sufficient for mediating budding in a variety of paramyxoviruses, (ii) the presence of a highly conserved putative L domain consensus in M protein sequences of all known pneumoviruses, and (iii) the enhanced expression of several Vps and Ub pathway proteins during APV infection, it is reasonable to hypothesize that APV employs the M protein to usurp the Vps and Ub pathways during late assembly and egress from the host cell. Based on this, we have developed a model for the recruitment of Vps proteins by the APV M protein (Fig. 3) for virus budding. Further studies testing this model and examining the role of the putative L domain of the APV M protein in recruiting Vps proteins are likely to provide important insights into pneumovirus biology and also identify the mechanisms involved in the final stages of paramyxovirus budding and maturation.
|
|
stimulation (49). Hence, up-regulation of B7-1 transcription in APV infection could be a sequel of IFN induction and may have a role in antigen presentation by the infected cells.
|
The differentially expressed genes were also analyzed by principal-component analysis to identify gene clusters based on the similarity of their expression patterns over time. Twelve distinct clusters of genes whose expression changed in a coordinated manner were obtained. Results for the cluster of IFN response genes and the cluster of chemokine and adhesion molecule genes are shown in Fig. 1B and C, respectively. There appears to be a strong correlation between functional groups and clusters, i.e., genes that perform similar functions share comparable expression profiles. For instance, the dynamics for all IFN response genes were very similar and clustered together as shown (Fig. 1B). This is not altogether surprising, since all IFN response genes are known to be coregulated by transcription factors activated by IFN-
/ß or virus infection. Of great interest is the finding that several differentially regulated genes that have no sequence homology in the database and were thus designated "unknowns" followed expression profiles very similar to those of the IFN-responsive genes and clustered with them (Fig. 1B). Based on this observation, it is likely that these unknown transcripts could be (i) novel IFN response factors or (ii) novel chicken ISGs with only limited homology to their known orthologues in other species. Future detailed analysis of these genes in the avian system, particularly examination of their upstream promoter sequences for the presence of ISREs and the ability to be activated by IFN stimulation, would help determine if they belong to the IFN-activated group of genes. Similarly, the chemokine (9E3 and K60), complement C3, and adhesion molecule (CD80 and ornithokinin receptor) mRNAs also had very similar expression profiles and clustered together into a distinct group. Interestingly, 9E3, K60, and C3, all three of which are chemoattractants, have profiles similar to each other, whereas the surface adhesion molecules (CD80-like protein and ornithokinin receptor), though they have patterns of induction similar to those of 9E3, K60, and C3, show dynamics more closely clustered with each other than with 9E3, K60, and C3. Again, a common transcription coregulatory mechanism, possibly mediated by IFN and nuclear factor kappa B (NF-
B), may be responsible for a similar temporal pattern of these genes.
Real-time RT-PCR analysis. RT-PCR analysis was performed on a set of 20 genes that were found to be either up-regulated, down-regulated, or unchanged in expression by microarray analysis of APV-infected cells. Results obtained by RT-PCR (Table 2) for the genes examined were in accordance with the microarray findings for 90% of the genes. The two genes showing discordance were ring finger protein 6 and CUL-1. They showed elevated expression by microarray analysis at 4.38- and 2.13-fold in APV-infected cells, respectively. Although both genes also showed up-regulation in expression at 1.66- and 1.23-fold, respectively, by RT-PCR, the change was not significant compared to that found by microarray analysis. The difference in the results of these two methods for 10% of the genes tested may possibly be due to differences in the kinetics and sensitivities of these two techniques for various target sequences.
A "synthesis" model of the host response to APV infection. The availability of expression profiles of host cell genes during the course of infection with APV enables the development of a "synthesis" model depicting how the cell responds to viral infection (Fig. 3). What emerges is a picture of the changes unfolding at the gene expression level in cells infected with APV and a revelation of key genes and pathways with a role in the antiviral and pathophysiological response to APV infection. The model presents a powerful means of visualizing the global changes that occur in an APV-infected cell and provides a facile means of developing testable hypotheses and designing specific experiments to test the roles of individual genes and pathways in an infected cell. For instance, the model predicts that IFN and proinflammatory genes are a predominant immediate host response to infection and that proinflammatory molecules could be the basis of disease-related pathology and likely targets for developing disease intervention strategies. The model also predicts the involvement of the Vps pathway in pneumovirus maturation and budding and provides a mechanism for this process. Further investigation of this mechanism is expected to advance our knowledge of the biology of pneumovirus budding, which could also be extended to other paramyxoviruses.
In sum, the results of our investigations have shed considerable light on host cell molecular responses to APV infection and have provided insights into host defensive measures as well as putative virus-directed mechanisms at the cellular level with implications for disease pathogenesis. The data reported here provide a myriad of mRNAs whose levels change upon host-APV interaction and that suggest an array of testable hypotheses should serve as a foundation for future research to unravel the roles of these genes in the pathogenic mechanisms of APV. Our investigations also demonstrate that study of the molecular details of host-virus interaction is crucial to our understanding of the complexities of the host-virus relationship and is likely to facilitate the development of antiviral agents, immunological interventions, and strategies to alleviate respiratory diseases in avian species.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
B and NF-IL-6. J. Infect. Dis. 174:262-267.[Medline]
. Chin. Med. J. 112:61-66.[Medline]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Mol. Cell. Biol. | Microbiol. Mol. Biol. Rev. |
|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|