J Gen Virol Try IJSEM Online
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Originally published as JGV in Press, 10.1099/vir.0.008664-0 on March 4, 2009 J Gen Virol 90 (2009), 1505-1514; DOI 10.1099/vir.0.008664-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Table
Right arrow All Versions of this Article:
vir.0.008664-0v1
90/6/1505    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Google Scholar
Right arrow Articles by Rasoolizadeh, A.
Right arrow Articles by Cusson, M.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Rasoolizadeh, A.
Right arrow Articles by Cusson, M.
Agricola
Right arrow Articles by Rasoolizadeh, A.
Right arrow Articles by Cusson, M.

Tranosema rostrale ichnovirus repeat element genes display distinct transcriptional patterns in caterpillar and wasp hosts

Asieh Rasoolizadeh1,2, Catherine Béliveau1, Don Stewart1, Conrad Cloutier2 and Michel Cusson1,2

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
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
The endoparasitic wasp Tranosema rostrale transmits an ichnovirus to its lepidopteran host, Choristoneura fumiferana, during parasitization. As shown for other ichnoviruses, the segmented dsDNA genome of the T. rostrale ichnovirus (TrIV) features several multi-gene families, including the repeat element (rep) family, whose products display no known similarity to non-ichnovirus proteins, except for a homologue encoded by the genome of the Helicoverpa armigera granulovirus; their functions remain unknown. This study applied linear regression of efficiency analysis to real-time PCR quantification of transcript abundance for all 17 TrIV rep open reading frames (ORFs) in parasitized and virus-injected C. fumiferana larvae, as well as in T. rostrale ovaries and head–thorax complexes. Although transcripts were detected for most rep ORFs in infected caterpillars, two of them clearly outnumbered the others in whole larvae, with a tendency for levels to drop over time after infection. The genome segments bearing the three most highly expressed rep genes in parasitized caterpillars were present in higher proportions than other rep-bearing genome segments in TrIV DNA, suggesting a possible role for gene dosage in the regulation of transcription level. TrIV rep genes also showed important differences in the relative abundance of their transcripts in specific tissues (cuticular epithelium, the fat body, haemocytes and the midgut), implying tissue-specific roles for individual members of this gene family. Significantly, no rep transcripts were detected in T. rostrale head–thorax complexes, whereas some were abundant in ovaries. There, the transcription pattern was completely different from that observed in infected caterpillars, suggesting that some rep genes have wasp-specific functions.

Primer sequences are available with the online version of this paper.


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Hymenopteran endoparasitoids deposit their eggs within the haemocoele of arthropods, most of which are insects (Eggleton & Belshaw, 1993Down). To protect their eggs from detection by the host immune system and to provide an appropriate developmental and physiological milieu for survival of their immature progeny within the host, female wasps typically inject their eggs along with various materials capable of disguising the egg surface and/or altering host physiology. For example, some members of the families Ichneumonidae and Braconidae transmit, to their caterpillar hosts, a virus that is essential for survival of the immature wasp within the parasitized insect (reviewed by Stoltz, 1993Down). These viruses, known as polydnaviruses (PDVs), feature a segmented, circular dsDNA genome, with individual genome segments varying in size and genetic content. A copy of the viral genome is present as a provirus within the wasp's chromosomes, thus providing a mechanism for the vertical transmission of PDVs within parasitoid populations. Virus replication is restricted to the calyx region of the wasp ovary, from which virions are released into the lumen of lateral oviducts. There, they form the particulate fraction of the calyx fluid (CF). During oviposition, a female wasp injects one or more eggs, along with CF and other secreted proteins and venom, into the lepidopteran host. Although no virus replication occurs in the parasitized caterpillar, expression of PDV genes causes developmental and immune dysfunctions that protect the egg and wasp larvae from encapsulation by host haemocytes and/or lead to retardation or arrest of host metamorphosis, thus providing more time for the wasp larva to complete its development in advance of host pupation (reviewed by Kroemer & Webb, 2004Down; Stoltz, 1993Down).

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, 1988Down). 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., 2007Down; Kroemer & Webb, 2005Down; Thoetkiattikul et al., 2005Down).

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, 1987Down). Members of the rep gene family encode non-secreted proteins that are conserved among several ichnovirus species (Tanaka et al., 2007Down; Volkoff et al., 2002Down; Webb et al., 2006Down). Expression of rep genes has been detected in both parasitoids and their parasitized hosts (Galibert et al., 2006Down; Theilmann & Summers, 1988Down). 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., 2007Down). In an earlier study, two TrIV rep genes (TrFrep1 and TrFrep2) were shown to be expressed from TrIV genome segment F (Volkoff et al., 2002Down); 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., 2007Down). 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., 2006Down). Using a recently developed and powerful quantitative real-time PCR (qPCR) approach (Rutledge & Stewart, 2008aDown, bDown), 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
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
RNA and DNA extraction.
Within 24 h of the moult to the last (sixth) instar, C. fumiferana larvae were either parasitized once by T. rostrale or injected with 0.5 female equivalents (FE) of T. rostrale CF, as described previously (Doucet & Cusson, 1996aDown, bDown). For each sampling point [1, 3 and 5 days post-parasitization (p.p.) or post-injection (p.i.)], total RNA was extracted and pooled from three to five whole C. fumiferana larvae, using TRIzol reagent (Invitrogen), according to the manufacturer's instructions (Béliveau et al., 2000Down). RNA was also extracted from the fat body, cuticular epithelium, midgut and haemocytes obtained from a pool of three to five larvae 48 h p.p. In addition, total RNA was extracted and pooled from five ovary pairs dissected from post-emergence 5–10-day-old T. rostrale females, using a QIAshredder and RNeasy Mini kit (Qiagen), according to the manufacturer's instructions. The head–thorax complexes of the same five females were also subjected to total RNA extraction using TRIzol reagent.

TrIV DNA was extracted from the CF of 16 T. rostrale female wasps as described previously (Stoltz et al., 1986Down). 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 {lambda} genomic DNA as a quantitative standard (Rutledge & Stewart, 2008aDown, bDown). 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, 1998Down) 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:


Formula

where nil is the number of amplification reactions failing to produce an amplification profile and total is the total number of reactions [see Rutledge & Stewart (2008b)Down for additional details about LDA]. Multiplication of Nmean by the dilution factor provides the LDA estimate.

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., 1997Down) using default settings.


   RESULTS AND DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Critical assessment of the LRE methodology
To assess the reliability of the qPCR estimates made in this study using the LRE approach, target quantities in three of our samples were determined using both LRE and LDA analysis. The LDA method generated estimates that were congruent with those obtained by LRE analysis, for one DNA and two RNA samples (Table 1Down), confirming the accuracy of the LRE methodology. In addition, amplification efficiencies were high and uniform across all 17 rep genes, in all treatment groups, and across all ten TrIV genome segments, with maximal amplification efficiencies (see Rutledge & Stewart 2008aDown, bDown) of 101.3±1.7 % and 101.5±1.5 % (mean±SD) for transcript and genome segment abundance, respectively. Thus, in assessing transcript levels for large multi-gene families such as those found in PDVs, or for measuring the relative abundance of many PDV genome segments, application of the LRE approach to qPCR determinations provides unprecedented accuracy and substantially improves analytical throughput over methods requiring the production of standard curves for each DNA examined.


View this table:
[in this window]
[in a new window]

 
Table 1. Critical assessment of the accuracy of LRE-based qPCR determinations (Rutledge & Stewart, 2008aDown, bDown) by comparison with estimates obtained by application of the LDA method (Wang & Spadoro, 1998Down)

Three example runs are shown, two for transcript levels and one for viral DNA (sample id). Nil, number of amplification reactions failing to produce an amplification profile; Nmean, mean number of molecules per aliquot (see Methods for details). The LRE values reported here for the transcripts are the copy number per 2.5 ng total RNA, whilst the value for C166 genomic DNA is the number of genome segments per 0.01 ng DNA (i.e. the concentrations at which the LRE measurements were made).

 
Transcript abundance in parasitized larvae
TrIV rep genes displayed important differences in gene-specific, time-dependent and tissue-specific levels of transcripts in naturally parasitized last-stadium C. fumiferana larvae. In whole caterpillars, transcripts were detected for almost all genes examined, but transcript abundance was generally low [<550 transcripts (ng total RNA)–1] except for F1-1 (~16 000 at 3 days p.p.) and, to a lesser extent, F1-2 (~2200 at 1 day p.p.; Fig. 1Down). Whether these differences in transcript abundance among rep genes are indicative of their relative importance in the subjugation of C. fumiferana hosts is not clear, but the strong predominance of F1-1 transcripts suggests that the product of this gene plays a vital role in the success of parasitism.


Figure 1
View larger version (12K):
[in this window]
[in a new window]

 
Fig. 1. Transcript levels of 17 TrIV rep genes in naturally parasitized C. fumiferana sixth-instar larvae, as determined by qRT-PCR using total RNA extracted from whole caterpillars at 1, 3 and 5 days p.p. Larvae were parasitized within 24 h of the moult to the sixth (last) stadium. Actual transcript numbers are provided above each bar for values ≤50 transcripts (ng total RNA)–1. Each value presented here is the mean±SD of four technical replicates carried out on an RNA extract obtained from a pool of three to five parasitized larvae.

 
Levels of rep gene transcripts were not stable during the course of parasitism and tended to decrease by >50 % between the first (day 1) and last (day 5) sampling points p.p., although three genes (most notably F1-1) displayed higher levels of transcripts at 3 days p.p. than at the other two sampling times (Fig. 1Up). A temporal pattern of expression similar to that observed here for F1-1 has been reported previously for another TrIV gene, TrV1, in parasitized C. fumiferana larvae (Béliveau et al., 2000Down). Differences in temporal patterns of expression among viral genes, in a given host, have been observed for other PDVs, including examples where maximal transcript levels were seen several days after oviposition (e.g. Ibrahim et al., 2007Down). Although such differences suggest that individual PDV gene products may target specific phases of parasitism, the observed transcriptional patterns may be dictated, at least in part, by the stability of the viral genome segments from which transcripts are generated, a variable that could differ considerably according to whether or not the developing wasp larva feeds on infected tissues supporting viral gene expression (Beck et al., 2007Down).

In our investigation of tissue-specific transcription at 2 days p.p., the overall gene-specific pattern of transcript abundance (Fig. 2Down) was similar to that observed in whole larvae (Fig. 1Up), 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., 2002Down), 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. 2Down). These results are in contrast to those obtained by Galibert et al. (2006)Down, 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.


Figure 2
View larger version (12K):
[in this window]
[in a new window]

 
Fig. 2. Transcript levels of 17 TrIV rep genes in naturally parasitized sixth-instar larvae, as determined by qRT-PCR using total RNA extracted from four different tissues: FB, fat body; CE, cuticular epithelium; HL, haemocytes; MG, midgut. The larvae were parasitized within 24 h of the moult to the sixth (last) stadium, and RNA was extracted from individual tissues 2 days p.p. Each value is the mean±SD of four technical replicates carried out on an RNA extract obtained from a pool of three to five parasitized larvae.

 
Tissue-specific differences in PDV gene transcript abundance in parasitized hosts have also been observed for ichnovirus ank genes (Kroemer & Webb, 2005Down) and bracovirus PTP genes (Gundersen-Rindal & Pedroni, 2006Down; Provost et al., 2004Down). Such tissue-specific expression suggests that the diversity of genes within a given PDV gene family may be associated with the existence of tissue-specific roles for different family members in the caterpillar hosts, or that some of these related gene products, whilst having the same function, are more effective in one tissue than in another. Irrespective of its functional significance, tissue-specific variation in transcript levels implies that there exist tissue-specific host factors modulating the transcription of specific rep genes.

Transcript abundance in CF-injected larvae
The TrIV rep transcript levels observed in CF-injected larvae (Fig. 3Down) displayed gene-specific and time-dependent differences similar to those observed for parasitized whole larvae (Fig. 1Up), 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., 2007Down).


Figure 3
View larger version (13K):
[in this window]
[in a new window]

 
Fig. 3. Transcript levels of 17 TrIV rep genes in sixth-instar larvae injected with 0.5 FE of T. rostrale calyx fluid, as determined by qRT-PCR using total RNA extracted from whole caterpillars at 1, 3 and 5 days p.i. Larvae were injected within 24 h of the moult to the sixth (last) stadium. Actual transcript numbers are provided above each bar for values <2000 transcripts (ng total RNA)–1. Each value is the mean±SD of four technical replicates carried out on an RNA extract obtained from a pool of three to five injected larvae.

 
Another difference between patterns found for parasitized and CF-injected larvae was the rise seen in F1-1 transcript abundance on day 3 p.p., an increase that was not observed in injected caterpillars, although absolute levels of F1-1 transcripts were higher in the latter than in the former group, at all three sampling points. With few exceptions, transcript levels decreased substantially from day 1 to days 3 and 5 p.i., suggesting that the unusually high inoculum injected in larvae may have triggered faster clearance or breakdown of some viral DNA in the host than in parasitized caterpillars. The present qPCR findings for F1-1 (TrFrep1) are in agreement with an earlier assessment made by Northern blot analysis, which showed F1-1 to be transcribed at much higher levels in CF-injected larvae than in parasitized caterpillars (Volkoff et al., 2002Down).

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 1Up, 2Up and 3Up), 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. 4Down). Interestingly, the C3-1 gene, whose transcription was barely detectable in infected caterpillars (Figs 1Up, 2Up and 3Up), 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 1Up, 2Up and 3Up), varied between ~5000 and ~10 000 transcripts (ng total RNA)–1 in wasp ovaries (Fig. 4Down). In comparison, all TrIV rep genes had undetectable or very low transcript levels in wasp head–thorax complexes (Fig. 4Down).


Figure 4
View larger version (12K):
[in this window]
[in a new window]

 
Fig. 4. Transcript levels of 17 TrIV rep genes in T. rostrale ovaries and head–thorax complexes, as determined by qRT-PCR. Total RNA was extracted from five ovary pairs dissected from post-emergence 5–10-day-old T. rostrale females and from the head–thorax complexes of the same females. Each value is the mean±SD of four technical replicates carried out on each RNA extract. Actual transcript numbers are provided above each bar for values <500 transcripts (ng total RNA)–1.

 
Using Northern blot analysis, Theilmann & Summers (1988)Down provided the first report on the transcription of Campoletis sonorensis ichnovirus (CsIV) rep genes in Cam. sonorensis female reproductive tissues. These authors observed that some rep genes were transcribed exclusively in the parasitized host, whilst others produced transcripts only in wasp ovaries or in both hosts. The quantitative transcriptional data provided here for 17 TrIV rep genes in both parasitized hosts and wasp ovaries are in agreement with this earlier finding. The distinct transcriptional patterns of rep genes in T. rostrale ovaries (Fig. 4Up) and parasitized or CF-injected larvae (Figs 1Up, 2Up and 3Up) suggest that individual rep genes may play either wasp- or caterpillar-specific roles. In contrast, HdIV rep1 was the most highly expressed rep gene in both infected caterpillar hosts and wasp ovaries (Galibert et al., 2006Down), suggesting that the host-specific expression reported here may not apply to all ichnoviruses. With respect to TrIV, the observation that some rep genes may be expressed only in the wasp (e.g. C3-1) raises the questions as to (i) why such genes are found in a packaged virus meant to be delivered to the lepidopteran host and (ii) whether there are additional, unpackaged rep genes in the T. rostrale genome, expressed only in wasp ovaries. Of course, the possibility exists that some of the TrIV rep genes that were found to be weakly expressed in parasitized C. fumiferana larvae would be expressed strongly in other lepidopteran hosts (e.g. Choristoneura rosaceana; Cusson et al., 1998Down) or tissues not sampled yet, including genes that were found here to be expressed only in the ovary. Additionally, it is not quite clear whether rep genes that are expressed in the wasp ovary are transcribed from episomal or chromosomal DNA, or both. Interestingly, none of the rep genes that were found to be expressed in T. rostrale ovaries was transcribed at significant levels in the other wasp tissues examined (Fig. 4Up), thus suggesting an ovary-specific role for those that are transcribed in that tissue. Given that rep gene products are not predicted to be secreted, rep proteins expressed in wasp ovaries are not expected to be released in the lumen of the oviduct for subsequent injection into the caterpillar during parasitization. For this reason, their expression in the ovary suggests that they could play a role in virus replication, a hypothesis that could be tested by following developmental changes in ovarian rep transcript abundance in pupae, the stage at which virus replication begins (Marti et al., 2003Down; Webb & Summers, 1992Down).

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., 2007Down; Galibert et al., 2006Down). Here, the three most highly expressed TrIV rep genes in parasitized caterpillars, F1-1, F1-2 and C166 (Fig. 1Up), were found to be borne by the two most abundant TrIV genome segments (Fig. 5Down), 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., 2007Down). 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.


Figure 5
View larger version (11K):
[in this window]
[in a new window]

 
Fig. 5. Assessment of genome segment abundance within the TrIV packaged genome, as determined by qPCR using viral DNA as template. The same primer pairs used for transcript quantification were used to quantify genome segments. C166 and C289 are contigs associated with genome segments that have not been cloned and that remain partially sequenced (Tanaka et al., 2007Down). Error bars indicate SD.

 
Comparison of TrIV rep proteins and identification of non-PDV rep homologues
A CLUSTAL_X alignment of all 17 deduced TrIV rep proteins revealed regions that are well conserved across all members of this family, including five cysteine residues that are present in all proteins except C7-2; the latter lacks the second and third cysteines, and its N terminus is substantially truncated relative to the other TrIV rep proteins (Fig. 6Down). The most conserved region is observed in the vicinity of the fifth cysteine residue (Fig. 6Down), comprising a segment of ~18 aa that also appears well conserved among rep proteins from other ichnoviruses, including Hyposoter fugitivus ichnovirus (HfIV) and HdIV (Fig. 7Down). A BLASTP search using all TrIV rep proteins as query sequences revealed the existence of two putative rep homologues in Helicoverpa armigera granulovirus (HearGV), the genome of which has recently been sequenced and annotated (Harrison & Popham, 2008Down). One of these two proteins, hear76, has only 70 aa and displays a modest level of similarity to ichnoviral rep proteins (e.g. BLASTP expect value of 0.36 for similarity to TrIV F3-1); however, the other predicted protein, hear75, has 171 aa, contains four of the five conserved cysteine residues referred to above and shows significant similarity to many ichnoviral rep proteins, most notably within the aforementioned highly conserved region (Fig. 7Down). BLASTP expect values for similarity between hear75 and ichnoviral rep proteins varied between 6x10–09 and 3x10–05 for HfIV D3-2 and TrIV F3-1, respectively. No rep homologues have been detected in the other baculovirus genomes sequenced to date; thus, their presence in HearGV may well be the result of lateral gene transfer from an ichnovirus genome (Harrison & Popham, 2008Down).


Figure 6
View larger version (130K):
[in this window]
[in a new window]

 
Fig. 6. CLUSTAL_X alignment of all known and predicted TrIV rep family proteins. Arrows indicate the positions of conserved cysteine residues. Asterisks (*), double dots (:) and single dots (.) above letters in the alignments denote identical residues, and conserved and semi-conserved substitutions, respectively. NCBI protein database accession nos: C7-1, BAF45598 [GenBank] ; C7-2, BAF45599 [GenBank] ; D5-1, BAF45610 [GenBank] ; D5-2, BAF45611 [GenBank] ; C289-1, BAF45769 [GenBank] ; C3-2, BAF45588 [GenBank] ; F3-2, BAF45626 [GenBank] ; F3-3, BAF45627 [GenBank] ; F3-1, BAF73402 [GenBank] ; C4-1, BAF45589 [GenBank] ; C4-2, BAF45590 [GenBank] ; B2-2, BAF45579 [GenBank] ; C3-1, BAF45585 [GenBank] ; D6-1, BAF45614 [GenBank] ; F1-1, AAN32723 [GenBank] ; F1-2, ACJ72220 [GenBank] ; C166-1, BAF45767 [GenBank] .

 

Figure 7
View larger version (50K):
[in this window]
[in a new window]

 
Fig. 7. CLUSTAL_X alignment of selected ichnoviral rep proteins from HfIV, TrIV and HdIV, as well as a rep-like protein from the granulovirus HearGV. Black arrows indicate the positions of conserved cysteine residues; the grey arrow points to the third conserved cysteine residue, which is replaced by an alanine in the HearGV rep-like protein. Asterisks (*), double dots (:) and single dots (.) above letters in the alignments denote identical residues, and conserved and semi-conserved substitutions, respectively. NCBI protein database accession nos: HfIV D3-2, BAF45718 [GenBank] ; HfIV D10-2, BAF45741 [GenBank] ; TrIV F3-1, BAF73402 [GenBank] ; HdIV rep5, AAR89177 [GenBank] ; HearGV hear75, ABY47766 [GenBank] .

 
As observed in earlier analyses of rep proteins, no conserved domains were detected in any of the 17 TrIV representatives of this family, with the exception of F1-2, in which a PIWI-like domain was detected in the region between residues 60 and 150, although with a low (0.001) expect value. The same protein was also found to display a modest level of similarity to a bacterial transposase (NCBI protein database accession no. ABM04822 [GenBank] ) within its C terminus, an interesting observation given that TrIV genome segment F1 has been shown to integrate spontaneously into the genome of C. fumiferana cells in culture, through an unidentified mechanism (Doucet et al., 2007Down).

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., 2003Down); 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., 2006Down). 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
 
The authors thank A. N. Volkoff for fruitful discussions about the work of her group on HdIV rep genes and D. Stoltz for helpful comments on an earlier version of the manuscript. This research was supported by grants from the Canadian Forest Service (CFS) and a Discovery grant from the Natural Sciences and Engineering Research Council of Canada to M. C.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Beck, M. H., Inman, R. B. & Strand, M. R. (2007). Microplitis demolitor bracovirus genome segments vary in abundance and are individually packaged in virions. Virology 359, 179–189.[CrossRef][Medline]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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{kappa}B-like gene family in polydnaviruses associated with wasps belonging to different braconid subfamilies. J Gen Virol 88, 92–104.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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{kappa}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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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 {kappa}B-like proteins from a polydnavirus inhibit NF-{kappa}B activation and suppress the insect immune response. Proc Natl Acad Sci U S A 102, 11426–11431.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Table
Right arrow All Versions of this Article:
vir.0.008664-0v1
90/6/1505    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Google Scholar
Right arrow Articles by Rasoolizadeh, A.
Right arrow Articles by Cusson, M.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Rasoolizadeh, A.
Right arrow Articles by Cusson, M.
Agricola
Right arrow Articles by Rasoolizadeh, A.
Right arrow Articles by Cusson, M.


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