Previous Article | Next Article ![]()
Journal of Virology, April 2004, p. 4145-4155, Vol. 78, No. 8
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.8.4145-4155.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Center for Biosystems Research, University of Maryland Biotechnology Institute, College Park, Maryland 20742-4450
Received 16 September 2003/ Accepted 21 November 2003
|
|
|---|
|
|
|---|
Fungal strains infected with the two hypoviruses also exhibit a number of distinguishing phenotypic traits when cultured under defined laboratory conditions. CHV1-EP713-infected C. parasitica strains exhibit reduced growth rates on solid media and produce colonies characterized by reduced production of aerial hyphae, irregular margins, and the general absence of asexual spores (4). C. parasitica strains infected with the mild hypovirus CHV1-Euro7 actually grow faster than the corresponding virus-free strain, and produce colonies with abundant aerial hyphae and regular margins and that produce conidia at levels intermediate between those of virus-free and CHV1-EP713-infected strains (4). Both hypoviruses suppress production of orange pigment, resulting in white colonies, and cause loss of female fertility.
Efforts to understand the mechanisms underlying the differences in symptom expression by the mild and severe hypovirus isolates revealed differential modulation of cyclic AMP-regulated signaling and laccase production by the mild and severe isolates (25). These results suggested the possibility that the two hypoviruses, even though closely related at the nucleotide and amino acid sequence levels, may have quite different effects on host gene expression. However, caution must be exercised in making such predictions since these studies relied on single gene readouts or differences in the activity of a specific family of enzymes.
Allen et al. (1) recently used a spotted cDNA microarray containing ca. 2,200 C. parastica genes to monitor the change in host transcriptional profile after infection with severe hypovirus isolate CHV1-EP713. That analysis confirmed that CHV1-EP713 infection results in a persistent reprogramming of a significant portion (estimated at 13.4%) of the C. parasitica transcriptome and expanded the number of identified C. parasitica genes that respond to hypovirus infection from less than 20 (9, 17) to nearly 300. We have taken advantage of the C. parasitica cDNA microarray platform to determine here whether the mild and severe hypovirus isolates cause similar or dissimilar changes in the host transcriptional profile and to identify virus-common and virus-specific host-responsive genes.
|
|
|---|
Total RNA isolation. Cultures used for RNA isolation were grown on PDA-cellophane for 6 days and harvested by freezing the mycelia in liquid nitrogen with immediate grinding of the mycelia into a fine powder by using a mortar and pestle. RNA isolation was performed as described by Allen et al. (1).
Microarray analysis. Microarray slide printing, fluorescent probe generation, microarray hybridization, scanning, and data analyses were performed as described by Allen et al. (1). Integrated pixel intensity values for each spot were calculated by using TIGR Spotfinder software and saved in tab-delimited format for use by MIDAS software (http://www.tigr.org/software) from The Institute for Genomic Research (Rockville, Md.). All hybridization data among three sets of dye-swap experiments were normalized simultaneously in MIDAS to correct for experimental error within a specific hybridization and between repeated hybridizations.
Selection of differentially expressed clones in each hybridization was performed in the Functional Genomics module of Spotfire DecisionSite 7.0 (Spotfire, Inc., Somerville, Mass.) by calculating the log2 (cy3/cy5) ratio for each clone. The clones were divided into "groups" based on the number of standard deviations that the corresponding log2 ratio varied from the data sets' average log2 ratio. Clones with log2 ratios of
2 standard deviations in a minimum of four of six hybridizations were considered differentially expressed. For experiments comparing RNA isolated from EP155/CHV1-EP713 to that from EP155/CHV1-Euro7, a clone was considered differentially expressed if its log2 ratio was
2 standard deviations in a minimum of three of five hybridizations. Genes identified as differentially expressed were subjected to hierarchical cluster analysis in Spotfire DecisionSite 7.0.
Microarray data management.
The arrayed EST (AEST) library used to construct the cDNA microarray chip consisted of 3,864 EST clones representing ca. 2,200 unique C. parasitica genes, of which roughly 1,600 were estimated to be present as a single clone (1, 8). Thus, ca. 27% of the genes present on the microarray are represented by multiple independent EST clones. Consequently, a subset of the genes on the primary list of differentially expressed genes generated in a microarray experiment will also be represented by more than one related AEST clone (8). To reduce the redundancy of the primary lists for publication purposes (1), the scored AEST clones are sorted according to NCBI gene index (gi) numbers assigned by basic local alignment search tool (BLAST) analysis. The AEST clone with the strongest BLAST hit is subsequently chosen to represent all of the scored AEST clones with the same gi number. This level of redundancy can present some minor cross-referencing problems when two published lists are compared. That is, if the representative AEST of a gi family for some technical reason fails to be scored, it will be replaced by the related EST with the next strongest BLAST hit. This can result in the same gene being represented by different related AEST clones on two different published, nonredundant lists of differentially expressed genes. A list of related AEST clones sorted according to gi numbers is available to facilitate cross-referencing (http://www.umbi.umd.edu/
cbr/AESTRedundant.pdf). Nonredundant lists for hybridizations involving EP155/CHV1-EP713 versus EP155, EP155/CHV1-Euro7 versus EP155, and EP155/CHV1-EP713 versus EP155/CHV1-Euro7 are available at http://www.umbi.umd.edu/
cbr/155-713NonRedundant.pdf, http://www.umbi.umd.edu /
cbr/155-Euro7NonRedundant.pdf, and http://www.umbi.umd.edu/
cbr/713-Euro7NonRedundant.pdf, respectively.
It is important that primary lists of responsive genes and not the published, nonredundant lists are used in all list comparison analyses. For example, host genes responsive to both mild and severe hypoviruses were identified by comparing the primary CHV1-EP713- and CHV1-Euro7-responsive lists in Microsoft Excel by using the MATCH command. The resulting common list was reduced in redundancy as described above for presentation in Table 1. Occasionally, members of a gi family are found to be present on two primary lists, but the representative of the family with the strongest BLAST hit is absent on one or both of the lists. In such cases, the family members common to both lists are sorted according to BLAST scores, and the AEST clone with the strongest BLAST hit is included in the published, nonredundant common list to represent the gene corresponding to that gi family.
|
View this table: [in a new window] |
TABLE 1. C. parasitica genes responsive to infection by both CHV1-EP713 and CHV1-Euro7a
|
Ct method (16) relative to the amount of 18S rRNA in the sample with primers and conditions as described by Parsley et al. (25). Differential expression based on real-time RT-PCR measurements was defined as a change in transcript accumulation of
2-fold. |
|
|---|
Of the 2,200 C. parasitica genes represented on the microarray, only 166 were scored as being differentially expressed in the CHV1-Euro7-infected strain (90 upregulated and 76 downregulated [available at http://www.umbi.umd.edu/
cbr/155-Euro7NonRedundant.pdf]). This is roughly one-half the number of genes found to be altered in transcript accumulation after infection by the severe isolate CHV1-EP713 (295 genes, with 132 upregulated and 163 downregulated [1]). That is, CHV1-Euro7 caused only half as many changes in the transcriptional profile as CHV1-EP713. Comparison of the 295 genes differentially expressed after CHV1-EP713 infection (EP155/CHV1-EP713 versus EP155) with the 166 genes differentially expressed after CHV1-Euro7 infection (EP155/CHV1-Euro7 versus EP155) indicated an overlap of 80 genes that were altered in transcript accumulation after infection of EP155 independently by CHV1-EP713 and CHV1-Euro7 (Fig. 1). The list of virus-common genes is presented in Table 1 under headings of putative biological processes as assigned by Dawe et al. (8), according to the classification guidelines outlined by the Gene Ontology Consortium (http://www.genontology.org) and as previously reported for CHV1-EP713-responsive genes by Allen et al. (1).
![]() View larger version (40K): [in a new window] |
FIG. 1. Venn diagram illustrating the total number of differentially expressed genes identified in hybridizations between EP155/CHV1-EP713 versus EP155 (1) and EP155/CHV1-Euro7 versus EP155 (the present study). A total of 80 genes were found on both lists of differentially expressed clones, and these genes are described in Table 1.
|
cbr/713-Euro7NonRedundant.pdf). Of these, 195 were increased in transcript abundance in EP155/CHV1-EP713 relative to EP155/CHV1-Euro7; 216 were decreased.
![]() View larger version (11K): [in a new window] |
FIG. 2. Illustration of the differences predicted in the list of differentially regulated genes generated from EP155/CHV1-EP713 versus EP155/CHV1-Euro7 hybridizations from those lists generated from EP155/CHV1-EP713 versus EP155 and EP155/CHV1-Euro7 versus EP155 hybridizations. The length of each colored arrow represents transcript abundance measurements of various magnitude relative to virus-free EP155. The red and green arrows represent transcript changes too small to be reliably detected by microarray profiling for EP155/CHV1-Euro7 versus EP155 and EP155/CHV1-EP713 versus EP155 hybridizations, respectively. The blue arrow represents the same transcript that can be reliably detected by microarray profiling for EP155/CHV1-EP713 versus EP155/CHV1-Euro7 hybridizations. Clones identified in this manner include clones [35], [23], and [12] highlighted in Fig. 3 and Table 3.
|
![]() View larger version (80K): [in a new window] |
FIG. 3. Visual representation of hierarchically clustered hybridization data sorted according to similarities in gene expression patterns. Column A represents the average log2 (cy3/cy5) ratio for each cDNA clone measured in six hybridizations of EP155/CHV1-EP713 versus EP155. Column B represents the average log2 (cy3/cy5) ratio for each cDNA clone measured in five hybridizations of EP155/CHV1-EP713 versus EP155/CHV1-Euro7. Column C represents the average log2 (cy3/cy5) ratio for each cDNA clone measured in six hybridizations of EP155/CHV1-Euro7 versus EP155. In columns A and C, red lines indicate an increase in transcript abundance in hypovirus-infected strains relative to virus-free EP155. Green lines indicate a decrease in transcript accumulation. In column B, red lines indicate an increase in transcript abundance in EP155/CHV1-EP713 relative to EP155/CHV1-Euro7. Green lines indicate a decrease. In all columns, black lines indicate no significant change in transcript accumulation between biological samples. Clones of interest are highlighted to the right of the cluster diagram. Each clone is preceded by a number in brackets, which refers the reader to the real-time RT-PCR data in Table 3.
|
![]() View larger version (38K): [in a new window] |
FIG. 4. Graphical representation of the relative magnitudes of microarray-predicted changes in transcript accumulation for the 80 identified virus-common C. parasitica-responsive genes (from Table 1). Blue-shaded bars indicate the magnitude of transcript accumulation change following CHV1-EP713 infection (fold change [y axis]), whereas the magnitude of change for the same genes after CHV1-Euro7 infection is indicated by the yellow-shaded bars. Specific genes tested by real-time RT-PCR are indicated by bracketed numbers that refer the reader to data listed in Table 3.
|
|
View this table: [in a new window] |
TABLE 2. Performance of control spots under three different hybridization conditionsa
|
|
View this table: [in a new window] |
TABLE 3. Validation of microarray measurements using real-time RT-PCRa
|
|
|
|---|
There is considerable evidence that hypoviruses are major contributors to the diversity in virulence levels and phenotypic traits observed for C. parasitica field isolates (6, 11-14, 22, 26). The constellation of altered traits that can accompany hypovirus-mediated attenuation of C. parasitica virulence is thought to significantly influence the effectiveness of hypovirulent strains as biological control agents for chestnut blight (9). For example, the severely reduced levels of asexual sporulation and loss of female fertility observed for many hypovirus-infected strains limit hypovirus transmission through and persistence within a C. parasitica population (22). The ability to monitor global changes in transcriptional responses to different hypoviruses provides unprecedented opportunities to gain insights into molecular mechanisms underlying differences in hypovirus-mediated symptom expression and virulence attenuation.
The high level of sequence similarity for CHV1-EP713 and CHV1-Euro7 coupled with the quite different phenotypic changes caused by the two viruses suggested several possible outcomes for microarray analysis. These ranged from (i) altered expression of a similar set of genes by both viruses but with a much greater magnitude of change by the severe virus to (ii) each virus causing alterations in transcript accumulation for an entirely distinct set of host genes. As indicated in Fig. 1, 3, and 4, the observed outcome fell between these two extremes. The number of host genes differentially expressed as a result of CHV1-Euro7 infection equaled approximately half the number that were responsive to CHV1-EP713 infection (166 versus 295, respectively). Of these, 80 genes were found to be responsive to both viruses.
Considered from another perspective, of the ca. 2,200 C. parasitica genes represented on the microarray chip, a total of 301 (13.6%) were scored as responsive to one virus or the other, and only 80 (3.6%) were responsive to both viruses, with 73 of those responding in the same direction. Extrapolations based on the estimate that the EST clones contained on the microarray represent ca. 22% of the predicted C. parasitica gene coding capacity (1) suggest that a total of only 363 host genes are altered in transcript accumulation as a general response to hypovirus infection. One could imagine that these virus-common responsive genes might include genes that are (i) altered in expression as part of a cellular defense response, (ii) involved in the control of hypovirus copy number, or (iii) altered in expression levels in order for viral RNA replication to proceed.
Differences were also evident in the magnitudes by which transcript accumulation changed in response to infection by the mild and severe hypovirus isolates. CHV1-Euro7 infection generally caused smaller changes in transcript accumulation than were observed for CHV1-EP713 infection. This difference in magnitude is clearly illustrated for the virus-common responsive genes in Fig. 4 and was supported by real-time RT-PCR analysis for a selected number of genes (Table 3). Additional experimental evidence for this general trend was generated by comparative hybridizations with cDNAs derived from RNA isolated from CHV1-EP713- and CHV1-Euro7-infected strain EP155. The expression change profile for the EP155/CHV1-EP713 versus EP155/CHV1-Euro7 hybridization was clearly more similar to that of EP155/CHV1-EP713 versus EP155 than to the profile derived from the EP155/CHV1-Euro7 versus EP155 hybridization (Fig. 3).
Allen et al. (1) previously discussed the potential relevance of a number of CHV1-EP713-responsive genes with high database sequence matches to hypovirus-mediated symptom expression and virus replication. Interestingly, several of these genes also appeared in lists of hypovirus-responsive genes generated in the present study. Transcripts for the glutathione S-transferase (GST) homologue AEST-12-G-04 and the HSP70 homologue AEST-10-H-10 were found to be constitutively increased (6- to 13-fold and 2.7- to 2.8-fold, respectively) in CHV1-EP713-infected mycelia (1). GSTs are a superfamily of isoenzymes responsible for detoxifying the cellular environment by removing reactive oxygen through conjugation of thiol reduced glutathione to various harmful ligands, including plant phenols and aflatoxins (29). Microarray profiling studies (using Affymetrix GeneChips) conducted in Arabidopsis thaliana leaves independently infected with five different RNA viruses (31) identified four different GST homologues that are each induced in the presence of all five viruses. Cellular HSP70 proteins have been reported to be recruited during infection by a number of viruses to facilitate virion assembly or genome replication (7, 15, 30). Since the increased transcript accumulation of GST and HSP70 homologues was not accompanied by altered transcript accumulation for many additional heat shock and classical stress response genes represented on the C. parasitica EST microarray chip, Allen et al. (1) advanced the possibility that these homologues belong to a subset of stress-related genes induced by hypovirus infection to facilitate viral functions. As indicated in Table 1, the HSP70 homologue was increased in transcript accumulation after infection by both viruses. However, the GST homologue was downregulated by 4.5-fold after CHV1-Euro7 infection (Table 3 and Fig. 3); one of just seven genes that were found to be regulated in opposite direction by CHV1-EP713 and CHV1-Euro7. The observation that GST transcript accumulation is strongly increased in EP155/CHV1-EP713 and strongly decreased in EP155/CHV1-Euro7 justifies an examination of relative effects of the mild and severe hypovirus isolates on the cellular redox state and possible ramifications for differences in gene expression. The upregulation of HSP70 by both viruses adds some additional support for the speculation that HSP70 may play a role in virus replication.
S-Adenosyl-L-methionine synthetase (SAMS; AEST-08-F-10 and S-adenosyl-L-homocysteine hydrolase (SAHH; represented by AEST-22-B-11 on the CHV1-EP713-responsive list and by AEST-02-A-07 on the CHV1-Euro7-responsive list, see legend to Table 3) are responsible for generation of the primary methyl donor S-adenosyl-L-methionine (SAMe) and removal of S-adenosyl-L-homocysteine (SAH), the by-product of transmethylation reactions, respectively. Because of the central role played by SAMe in cellular metabolism, Allen et al. (1) considered the possibility that the constitutive increase in transcript accumulation for genes encoding both SAMS and SAHH in response to CHV1-EP713 infection could have significant metabolic or physiologic consequences for the host. The association between abnormal intracellular levels of SAMe and altered genome stability (18-20, 27), coupled with the reported influence of intracellular SAH levels on senescence and cell growth (21, 32), also prompted the suggestion that hypovirus infection may provide a useful model for examining the consequences of chronic RNA virus infection on host genome stability. These considerations are reinforced by the observation that homologues of both of these enzymes are upregulated by both the mild and severe isolates but to a greater extent by CHV1-EP713 (Tables 1 and 3).
Allen et al. (1) made the interesting observation that only 3 of 26 genes on the EST chip that fell under the molecular function category "transcription regulation/transcription factors" in Dawe et al. (8) were responsive to CHV1-EP713; each was reduced by at least threefold. These included the homologues of Mst12 (AEST-30-C-09) from Magnaporthe grisea, shown to be important for regulating infectious hyphae growth (24); Pro1 (AEST-27-F-10), involved in controlling sexual sporulation in several filamentous fungi (23); and HoxX (AEST-05-C-02), part of a bacterial two-component regulatory system (10). HoxX appears on the virus-common list (Table 1), and its downregulation by both viruses was confirmed by real-time RT-PCR (Table 3). Pro1 and Mst12 homologues were present only on the CHV1-EP713-responsive list (Table 3). However, a slight reduction in transcript accumulation (below the twofold cutoff) was observed for these genes by real-time RT-PCR measurements after CHV1-Euro7 infection. Finally, the homologue of the general transcription factor AP1 (AEST-12-F-11) was found by microarray analysis to be upregulated by CHV1-Euro7. This result was confirmed by real-time RT-PCR (Table 3), which also revealed a modest level of upregulation in CHV1-EP713-infected mycelia. This collection of hypovirus-regulated putative transcription factors represent prime candidates for future functional studies related to hypovirus-mediated alterations of host gene expression.
The microarray analyses described here have resulted in the initial assignment of >20% of the C. parasitica-expressed genes to one of the following categories for the experimental conditions used in the present study: (i) nonresponsive to hypovirus infection, (ii) upregulated by both mild and severe hypoviruses, and (iii) downregulated by both mild and severe hypoviruses, and (iv) regulated in opposite directions by the two viruses. As indicated in Table 1 and online supplementary lists, the genes assigned an "unknown" biological function form the largest category in both the virus-specific and virus-common responsive gene lists reported in the present study. Consequently, a comprehensive understanding of the significance of the transcriptional responses to hypovirus infection will require additional functional studies of the differentially expressed genes identified through microarray profiling. However, it is anticipated that further refinement of the virus-specific and virus-common lists through the addition of differential expression data for other experimental conditions and hypoviruses, e.g., chimeric viruses, will reveal trends and patterns that will drive the direction of future mechanistic studies.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»