|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, 1055 du PEPS, PO Box 10380, Stn Sainte-Foy, QC G1V 4C7, Canada
2 Département de Biologie, Université Laval, QC G1V 0A6, Canada
Correspondence
Michel Cusson
michel.cusson{at}nrcan.gc.ca
| ABSTRACT |
|---|
|
|
|---|
Primer sequences are available with the online version of this paper.
| INTRODUCTION |
|---|
|
|
|---|
PDV genes are divided into three broad categories based on whether they are expressed in the carrier wasp (class I), the infected caterpillar (class II) or both (class III) (Theilmann & Summers, 1988
). Because of their potentially major significance in the success of parasitism, class II genes have been studied more extensively than those of the other two groups. A number of these genes encode proteins displaying motifs or structural and sequence features observed in previously characterized eukaryotic proteins. Based on these similarities, it has been possible to generate and test hypotheses about their probable functions. Such an approach has led to proposed functions for various PDV genes (e.g. the vankyrins) (Falabella et al., 2007
; Kroemer & Webb, 2005
; Thoetkiattikul et al., 2005
).
Other PDV genes, however, display no known similarity to other eukaryotic or viral (non-PDV) genes, rendering their functional analysis more difficult. Such is the case for the repeat element (rep) gene family, the largest gene family identified to date in the genus Ichnovirus (PDVs associated with ichneumonid wasps). These genes consist of imperfectly conserved repeats of
540 bp, arranged either singly or in direct tandem arrays (Theilmann & Summers, 1987
). Members of the rep gene family encode non-secreted proteins that are conserved among several ichnovirus species (Tanaka et al., 2007
; Volkoff et al., 2002
; Webb et al., 2006
). Expression of rep genes has been detected in both parasitoids and their parasitized hosts (Galibert et al., 2006
; Theilmann & Summers, 1988
). The Tranosema rostrale ichnovirus (TrIV) genome contains at least 17 different open reading frames (ORFs) identified as belonging to the rep gene family; they are located on ten different genome segments (Tanaka et al., 2007
). In an earlier study, two TrIV rep genes (TrFrep1 and TrFrep2) were shown to be expressed from TrIV genome segment F (Volkoff et al., 2002
); this genome segment has been renamed F1 and the two rep genes it contains are now referred to as F1-1 and F1-2 (Tanaka et al., 2007
). As a first step towards elucidating the function(s) of these gene products, we initiated a study of the temporal and tissue-specific transcription of all known and putative TrIV rep genes. A similar study of ten rep genes from the ichnovirus of Hyposoter didymator (HdIV) has revealed important differences in gene-specific transcript abundance but minor differences in host and tissue specificity (Galibert et al., 2006
). Using a recently developed and powerful quantitative real-time PCR (qPCR) approach (Rutledge & Stewart, 2008a
, b
), the present study examined transcriptional patterns in the host Choristoneura fumiferana, either naturally parasitized by T. rostrale or injected with its CF, as well as in the wasp's ovaries and head–thorax complexes. We also examined the possible effect of gene dosage (i.e. relative genome segment abundance) on rep gene transcription. Finally, we have presented the results of new bioinformatics analyses conducted with the aim of detecting rep homologues among more recent entries in the NCBI protein database.
| METHODS |
|---|
|
|
|---|
TrIV DNA was extracted from the CF of 16 T. rostrale female wasps as described previously (Stoltz et al., 1986
). The DNA was first ethanol precipitated and then resuspended in 100 µl Tris/EDTA buffer (pH 7.6).
Reverse transcription and qPCR.
To remove DNA contaminants from RNA extracts, 500 ng total RNA was treated with 2 U amplification-grade DNase I (Invitrogen) for 15 min at 25 °C. We ran controls with no reverse transcriptase for the four most highly transcribed ORFs and detected no significant amplification, indicating the absence of genomic DNA contamination in the extracts. RNA (500 ng) from parasitized and CF-injected C. fumiferana larvae, as well as from T. rostrale head–thorax complexes, and 200 ng RNA from ovarian tissue was reverse transcribed using 0.5 µg of an oligo(dT) primer and 200 U Superscript II RNase H– reverse transcriptase (Invitrogen). The reaction was carried out in 1x PCR buffer, with 0.5 mM each dNTP and 40 U RNAguard RNase inhibitor (Amersham Biosciences), at 42 °C for 50 min.
For qPCR analysis, four primers were initially designed for each rep gene, using diverse regions among aligned rep nucleotide sequences. These four primer pairs were used to assess primer performance and quantitative precision. Initial amplification tests were conducted on reverse-transcribed RNA obtained from parasitized C. fumiferana larvae. A single primer pair was then selected for each rep gene (see Supplementary Table S1, available in JGV Online), based on high amplification efficiency and lack of non-specific amplification products, and used for analysis of the remaining samples.
PCR amplifications were carried out on aliquots of individual reverse transcription reactions containing cDNA in amounts equivalent to 2.5 ng RNA, except for ovarian samples, which contained amounts of cDNA equivalent to 1 ng RNA. Four replicate amplification reactions containing 500 nM each primer were conducted for each sample, using an MX3000P spectrofluorometric thermal cycler (Stratagene) and a QuantiTect SYBR Green PCR kit (Qiagen), initiated with a 15 min incubation at 95 °C, followed by a cycling regime (50 cycles) of 95 °C for 10 s and 65 °C for 2 min. Each run was completed with a melting-curve analysis to confirm the specificity of amplification and the absence of primer dimers. Amplification efficiency was determined for each amplification reaction using linear regression of efficiency (LRE) analysis, and the number of target molecules was calculated using
genomic DNA as a quantitative standard (Rutledge & Stewart, 2008a
, b
). LRE is a powerful methodology recently developed for modelling real-time qPCR amplification. It provides absolute target amounts without the need to produce standard curves and can generate absolute accuracies of less than ±25 %, whilst displaying single-molecule sensitivity.
To assess the proportion of each rep gene-bearing genome segment within a TrIV DNA extract, the same qPCR approach was applied directly to 0.01 ng TrIV DNA, using one of the primer pairs designed for transcript quantification for each genome segment (see Supplementary Table S1).
To evaluate the accuracy of the measurements made here using LRE analysis, we applied a limiting-dilution assay (LDA; Wang & Spadoro, 1998
) approach to three of our samples, and compared the estimates obtained with each method. Briefly, based on values determined by LRE, samples were diluted so that each of 20 replicate aliquots would contain approximately one copy of cDNA or genomic DNA. As dictated by Poisson distribution, a large proportion of aliquots will not contain a target molecule and will fail to produce an amplification profile. The mean number of molecules per aliquot (Nmean) can be calculated using the equation:
|
|
Bioinformatics.
To explore the possibility that sequences recently deposited in GenBank may be homologous to ichnoviral rep genes, all TrIV rep proteins were submitted to a BLASTP analysis. Alignments of amino acid sequences were performed with CLUSTAL_X (Thompson et al., 1997
) using default settings.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
16 000 at 3 days p.p.) and, to a lesser extent, F1-2 (
2200 at 1 day p.p.; Fig. 1
|
In our investigation of tissue-specific transcription at 2 days p.p., the overall gene-specific pattern of transcript abundance (Fig. 2
) was similar to that observed in whole larvae (Fig. 1
), but with some notable exceptions. For example, in the four tissues examined, F1-2 displayed lower proportions of transcripts relative to F1-1 than in whole larvae, whilst the opposite trend was observed for F3-2. This suggests that the tissues supporting high levels of F1-2 transcription were not sampled in the present study, whereas some of the sampled tissues were enriched for F3-2 transcripts. More significantly, TrIV rep genes exhibited important differences in their tissue specificity: whereas F1-1 transcripts were most abundant in C. fumiferana cuticular epithelium and the fat body, corroborating previous assessments made by Northern blot analysis (Volkoff et al., 2002
), the transcripts of several other genes were at higher levels in haemocytes (B2-2, C7-2 and F3-2) or the midgut (C4-2, D5-2, D6-1 and F1-2) than in the other three tissues (Fig. 2
). These results are in contrast to those obtained by Galibert et al. (2006)
, who found that the fat body and cuticular epithelium of parasitized Spodoptera littoralis hosts had the highest levels of HdIV rep transcripts for all ten rep ORFs examined, followed by nervous tissue, which was not investigated in the present study. It remains to be seen whether the observed trend in HdIV rep gene expression was influenced by the choice of rep ORFs that were studied, as we now know that the HdIV genome contains many additional rep genes (A. N. Volkoff, personal communication). Thus, this apparent difference between the two biological systems could be due to a gene-sampling bias.
|
Transcript abundance in CF-injected larvae
The TrIV rep transcript levels observed in CF-injected larvae (Fig. 3
) displayed gene-specific and time-dependent differences similar to those observed for parasitized whole larvae (Fig. 1
), with the exception that absolute transcript levels were generally much higher than those observed at equivalent sampling times in parasitized larvae, particularly 1 day after treatment (>85 times higher in the case of F1-1), indicating that the virus dose contained in 0.5 FE of CF is much higher than that injected by a female wasp during natural parasitization. As a point of comparison, the dose of virus injected by the wasp Microplitis demolitor into its host has been estimated to be between 0.04 and 0.005 FE of CF per ovipositional event (Beck et al., 2007
).
|
Transcript abundance in wasp ovary and head–thorax complexes
The pattern of TrIV rep gene transcription in T. rostrale ovaries was markedly different from that seen in naturally parasitized or CF-injected C. fumiferana larvae. Whereas F1-1 and F1-2 were the most highly expressed rep genes in infected caterpillars (Figs 1
, 2
and 3
), transcripts generated from these two genes displayed low abundance in wasp ovaries compared with other genes such as C166-1, the transcript levels of which were by far the highest (Fig. 4
). Interestingly, the C3-1 gene, whose transcription was barely detectable in infected caterpillars (Figs 1
, 2
and 3
), was the second most highly transcribed gene in wasp ovaries. In addition, the transcript levels of C3-2, C7-2, D5-2 and F3-2, which were modest in infected C. fumiferana larvae (Figs 1
, 2
and 3
), varied between
5000 and
10 000 transcripts (ng total RNA)–1 in wasp ovaries (Fig. 4
). In comparison, all TrIV rep genes had undetectable or very low transcript levels in wasp head–thorax complexes (Fig. 4
).
|
Gene dosage
In earlier work examining the relationship between the abundance of PDV gene transcripts and the proportion of the genome segments bearing these genes within the packaged viral genome, no clear correlation between the two variables was observed (Beck et al., 2007
; Galibert et al., 2006
). Here, the three most highly expressed TrIV rep genes in parasitized caterpillars, F1-1, F1-2 and C166 (Fig. 1
), were found to be borne by the two most abundant TrIV genome segments (Fig. 5
), suggesting that gene dosage, in this particular instance, may have some impact on transcript abundance. Yet, when all TrIV rep genes were considered, we observed no significant correlation between transcript levels and the proportion of the originating genome segments. Clearly, factors other than, or in addition to, gene dosage affect transcript levels, including possible differences in promoter strength, the presence or absence of host factors that may affect the transcription of individual rep genes and/or differences in mRNA stability. For example, there were important differences in the abundance of F1-1 and F1-2 transcripts, which are generated from genes present on the same genome segment. Integration of genome segment F1 into the lepidopteran host genomic DNA could also be a factor resulting in the enhancement of F1-1 and F1-2 transcription. Although the integration of genome segment F1 has not been demonstrated in the parasitized host, it clearly occurs in infected C. fumiferana CF-124T cells in culture (Doucet et al., 2007
). Such an integration event would permit sustained expression of the integrated genes when titres of episomal DNA decrease. The question of whether other rep-containing genome segments undergo integration into C. fumiferana genomic DNA remains to be examined.
|
18 aa that also appears well conserved among rep proteins from other ichnoviruses, including Hyposoter fugitivus ichnovirus (HfIV) and HdIV (Fig. 7
|
|
Although the new bioinformatics analyses performed here provided few new insights into the function(s) of rep genes, the presence of rep homologues in the recently sequenced genome of a granulovirus could eventually provide an indirect means of assessing their role through the production of a HearGV rep knockout, followed by an assessment of the effect of this genetic alteration on viral replication or other aspects of the infection cycle. Deployment of this strategy, however, would require the prior development of an efficient in vitro system for HearGV.
In summary, the present study suggests that the very high level of diversification seen within the ichnoviral rep gene family may have evolved in response to the necessity to fine-tune the function(s) and/or effectiveness of rep proteins for expression in different hosts and tissues. Given that rep genes encode proteins that are not secreted and that some of them are expressed at relatively high levels in wasp ovaries without any overt pathological consequence, the possibility exists that their function has more to do with cell homeostasis (in ichnovirus- or granulovirus-infected lepidopteran cells or in ovarian wasp cells supporting virus replication) than virulence. Some PDV-encoded proteins are secreted and display deleterious effects on other cells (e.g. Béliveau et al., 2003
); because PDVs do not replicate in the lepidopteran host, sustained viral gene expression for the duration of immature parasitoid development is predicted to require a mechanism preventing infected cells from being negatively affected by secreted PDV proteins and/or suppressing breakdown of viral DNA and transcripts by host cells. Some CsIV ank gene products appear to have such a function, given that they have been shown to delay lysis of baculovirus-infected cells (Fath-Goodin et al., 2006
). We are currently examining the effect of TrIV rep gene expression on C. fumiferana host-cell gene expression, with the aim of identifying the pathway(s) targeted by rep proteins.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Béliveau, C., Laforge, M., Cusson, M. & Bellemare, G. (2000). Expression of a Tranosema rostrale polydnavirus gene in the spruce budworm, Choristoneura fumiferana. J Gen Virol 81, 1871–1880.
Béliveau, C., Levasseur, A., Stoltz, D. & Cusson, M. (2003). Three related TrIV genes: comparative sequence analysis and expression in host larvae and Cf-124T cells. J Insect Physiol 49, 501–511.[CrossRef][Medline]
Cusson, M., Barron, J. R., Goulet, H., Régnière, J. & Doucet, D. (1998). Biology and status of Tranosema rostrale rostrale (Hymenoptera: Ichneumonidae), a parasitoid of the eastern spruce budworm (Lepidoptera: Tortricidae). Ann Entomol Soc Am 91, 87–93.
Doucet, D. & Cusson, M. (1996a). Alteration of developmental rate and growth of Choristoneura fumiferana parasitized by Tranosema rostrale: role of the calyx fluid. Entomol Exp Appl 81, 21–30.[CrossRef]
Doucet, D. & Cusson, M. (1996b). Role of calyx fluid in alterations of immunity in Choristoneura fumiferana larvae parasitized by Tranosema rostrale. Comp Biochem Physiol 114, 311–317.
Doucet, D., Levasseur, A., Béliveau, C., Lapointe, R., Stoltz, D. & Cusson, M. (2007). In vitro integration of an ichnovirus genome segment into the genomic DNA of lepidopteran cells. J Gen Virol 88, 105–113.
Eggleton, P. & Belshaw, B. (1993). Comparisons of dipteran, hymenopteran and coleopteran parasitoids: provisional phylogenetic explanations. Biol J Linn Soc Lond 48, 213–226.[CrossRef]
Falabella, P., Varricchio, P., Provost, B., Espagne, E., Ferrarese, R., Grimaldi, A., Eguileor, M., Fimiani, G., Ursini, M. V. & other authors (2007). Characterization of the I
B-like gene family in polydnaviruses associated with wasps belonging to different braconid subfamilies. J Gen Virol 88, 92–104.
Fath-Goodin, A., Kroemer, J., Martin, S., Reeves, K. & Webb, B. A. (2006). Polydnavirus genes that enhance the baculovirus expression vector system. Adv Virus Res 68, 75–90.[CrossRef][Medline]
Galibert, L., Devauchelle, G., Cousserans, F., Rocher, J., Cérutti, P., Barat-Houari, M., Fournier, P. & Volkoff, A. N. (2006). Members of the Hyposoter didymator ichnovirus repeat element gene family are differentially expressed in Spodoptera frugiperda. Virol J 3, 48[CrossRef][Medline]
Gundersen-Rindal, D. E. & Pedroni, M. J. (2006). Characterization and transcriptional analysis of protein tyrosine phosphatase genes and an ankyrin repeat gene of the parasitoid Glyptapanteles indiensis polydnavirus in the parasitized host. J Gen Virol 87, 311–322.
Harrison, R. L. & Popham, H. J. R. (2008). Genomic sequence analysis of a granulovirus isolated from the Old World bollworm, Helicoverpa armigera. Virus Genes 36, 565–581.[CrossRef][Medline]
Ibrahim, A. M., Choi, J. Y., Je, Y. H. & Kim, Y. (2007). Protein tyrosine phosphatases encoded in Cotesia plutellae bracovirus: sequence analysis, expression profile, and a possible biological role in host immunosuppression. Dev Comp Immunol 31, 978–990.[CrossRef][Medline]
Kroemer, J. A. & Webb, B. A. (2004). Polydnavirus genes and genomes: emerging gene families and new insights into polydnavirus replication. Annu Rev Entomol 49, 431–456.[CrossRef][Medline]
Kroemer, J. A. & Webb, B. A. (2005). I
B-related vankyrin genes in the Campoletis sonorensis ichnovirus: temporal and tissue-specific patterns of expression in parasitized Heliothis virescens lepidopteran hosts. J Virol 79, 7617–7628.
Marti, D., Grossniklaus-Bürgin, C., Wyder, S., Wyler, T. & Lanzrein, B. (2003). Ovary development and polydnavirus morphogenesis in the parasitic wasp Chelonus inanitus. I. Ovary morphogenesis, amplification of viral DNA and ecdysteroid titres. J Gen Virol 84, 1141–1150.
Provost, B., Varricchio, P., Arana, E., Espagne, E., Falabella, P., Huguet, E., La Scaleia, R., Cattolico, L., Poirié, M. & other authors (2004). Bracoviruses contain a large multigene family coding for protein tyrosine phosphatases. J Virol 78, 13090–13103.
Rutledge, R. G. & Stewart, D. (2008a). A kinetic-based sigmoidal model for the polymerase chain reaction and its application to high-capacity absolute quantitative real-time PCR. BMC Biotechnol 8, 47[CrossRef][Medline]
Rutledge, R. G. & Stewart, D. (2008b). Critical evaluation of methods used to determine amplification efficiency refutes the exponential character of real-time PCR. BMC Mol Biol 9, 96[CrossRef][Medline]
Stoltz, D. B. (1993). The polydnavirus life cycle. In Parasites and pathogens of insects, vol. 1, pp. 80–101. Edited by N. E. Beckage, S. N. Thompson & B. A. Federici. San Diego, CA: Academic Press.
Stoltz, D. B., Guzo, D. & Cook, D. (1986). Studies on polydnavirus transmission. Virology 155, 120–131.[CrossRef][Medline]
Tanaka, K., Lapointe, R., Barney, W. E., Makkay, A. M., Stoltz, D., Cusson, M. & Webb, B. A. (2007). Shared and species-specific features among ichnovirus genomes. Virology 363, 26–35.[CrossRef][Medline]
Theilmann, D. A. & Summers, M. D. (1987). Physical analysis of the Campoletis sonorensis virus multipartite genome and identification of a family of tandemly repeated elements. J Virol 61, 2589–2598.
Theilmann, D. A. & Summers, M. D. (1988). Identification and comparison of Campoletis sonorensis virus transcripts expressed from four genomic segments in the insect hosts Campoletis sonorensis and Heliothis virescens. Virology 167, 329–341.[Medline]
Thoetkiattikul, H., Beck, M. H. & Strand, M. R. (2005). Inhibitor
B-like proteins from a polydnavirus inhibit NF-
B activation and suppress the insect immune response. Proc Natl Acad Sci U S A 102, 11426–11431.
Thompson, J. D., Gibson, T. J., Plewniak, F., Jeanmougin, F. & Higgins, D. G. (1997). The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25, 4876–4882.
Volkoff, A. N., Béliveau, C., Rocher, J., Hilgarth, R., Levasseur, A., Duonor-Cérutti, M., Cusson, M. & Webb, B. A. (2002). Evidence for a conserved polydnavirus gene family: ichnovirus homologs of the CsIV repeat element genes. Virology 300, 316–331.[CrossRef][Medline]
Wang, Z. & Spadoro, J. (1998). Determination of target copy number of quantitative standards used in PCR-based diagnostic assays. In Gene quantification, pp. 31–43. Edited by F. Ferré. Boston: Birkhäuser.
Webb, B. A. & Summers, M. D. (1992). Stimulation of polydnavirus replication by 20-hydroxyecdysone. Experientia 48, 1018–1022.[CrossRef][Medline]
Webb, B. A., Strand, M. R., Dickey, S. E., Beck, M. H., Hilgarth, R. S., Barney, W. E., Kadash, K., Kroemer, J. A., Lindstrom, K. G. & other authors (2006). Polydnavirus genomes reflect their dual roles as mutualists and pathogens. Virology 347, 160–174.[CrossRef][Medline]
Received 7 November 2008;
accepted 2 February 2009.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |