Geminiviruses are a large and important family of plant viruses that infect a wide range of crops throughout the world. The Begomovirus genus contains species that are transmitted by whiteflies and are distributed worldwide causing disease on an array of horticultural crops. Symptom remission, in which newly developed leaves of systemically infected plants exhibit a reduction in symptom severity (recovery), has been observed on pepper (Capsicum annuum) plants infected with Pepper golden mosaic virus (PepGMV). Previous studies have shown that transcriptional and post-transcriptional gene silencing mechanisms are involved in the reduction of viral nucleic acid concentration in recovered tissue. In this study, we employed deep transcriptome sequencing methods to assess transcriptional variation in healthy (mock), symptomatic, and recovered pepper leaves following PepGMV infection.
Differential expression analyses of the pepper leaf transcriptome from symptomatic and recovered stages revealed a total of 309 differentially expressed genes between healthy (mock) and symptomatic or recovered tissues. Computational prediction of differential expression was validated using quantitative reverse-transcription PCR confirming the robustness of our bioinformatic methods. Within the set of differentially expressed genes associated with the recovery process were genes involved in defense responses including pathogenesis-related proteins, reactive oxygen species, systemic acquired resistance, jasmonic acid biosynthesis, and ethylene signaling. No major differences were found when compared the differentially expressed genes in symptomatic and recovered tissues. On the other hand, a set of genes with novel roles in defense responses was identified including genes involved in histone modification. This latter result suggested that post-transcriptional and transcriptional gene silencing may be one of the major mechanisms involved in the recovery process. Genes orthologous to the C. annuum proteins involved in the pepper-PepGMV recovery response were identified in both Solanum lycopersicum and Solanum tuberosum suggesting conservation of components of the viral recovery response in the Solanaceae.
These data provide a valuable source of information for improving our understanding of the underlying molecular mechanisms by which pepper leaves become symptomless following infection with geminiviruses. The identification of orthologs for the majority of genes differentially expressed in recovered tissues in two major solanaceous crop species provides the basis for future comparative analyses of the viral recovery process across related taxa.
Geminiviruses are a large and important family of plant viruses that infect a wide variety of crops around the world. The family Geminiviridae is divided into four genera (Mastrevirus, Curtovirus, Begomovirus, Topocuvirus) based on genome organization (mono- or bipartite), insect vector (whiteflies, leafhoppers, treehoppers), and host range (monocotyledonous or dicotyledonous plant species). Geminiviral genomes are composed of circular, single-stranded DNA (ssDNA) molecules encapsidated in twin icosahedral virions for which the family is named [1, 2]. The ssDNA viral genomes are transcribed, replicated, and encapsidated in the nuclei of infected cells. Geminiviruses also traffic within and between host cells moving systemically throughout the infected plant and are dependent on host machinery for both replication and movement [3, 4].
Species of the genus Begomovirus are transmitted by whiteflies (Bemisia tabaci Genn.) and distributed worldwide causing diseases in horticultural crops such as tomato and pepper [5–7]. Crop losses of up to 100% have been reported for geminivirus diseases . Pepper golden mosaic virus (PepGMV) is one of the most important and widely distributed Begomovirus throughout Mexico and infects several Solanaceae crops including pepper (Capsicum annuum), tomato (Solanum lycopersicum), and tomatillo (Physalis ixocarpa) [8–10]. The bipartite genome of PepGMV encodes six proteins. The DNA-A component encodes proteins involved in replication (Rep and REn), trans-activation (TrAP), and the capsid protein (CP) whereas viral DNA-B encodes proteins related to movement (NSP and MP) . PepGMV infection results in bright yellow mosaic symptoms on leaves that is associated with twisting and distortion of leaves and fruits, stunted plants, and reduced yield .
In studies under controlled conditions, symptom remission has been observed on PepGMV-infected pepper plants . Thus, in the pepper-PepGMV recovery system, pepper plants show severe symptoms around 9 to 10 days post-inoculation (dpi) (Figure 1A,B, “Symptoms”). Following this initial expression of disease symptoms, newly emerged leaves show a reduction in symptom severity (Figure 1C, “Pre-Recovery”). By 20 dpi, the third set of emerging leaves is nearly symptomless (Figure 1D, “Recovery”). This process has been termed symptom remission or recovery [12, 13]. Similar observations with several plant viruses were reported as early as 1928  and more recently, in plants infected with geminiviruses [15–17]. In PepGMV-infected pepper plants, the concentrations of viral ssDNA and viral mRNA in recovered leaves are dramatically reduced compared to the viral nucleic acid concentrations in severely symptomatic leaves [12, 13]. Based on the presence and characterization of small RNA of viral origin (svRNA), as well as the levels of methylated viral DNA by host methylation machinery, the recovery process has been associated with transcriptional and post-transcriptional gene silencing mechanisms . However, little is known in the pepper-PepGMV recovery system with respect to changes in host gene expression that occur during infection.
For organisms with no publicly available genome sequence, such as pepper, transcriptome analyses can provide major insights into genes involved in important biological processes. Deep transcriptome sequencing technologies, such as pyrosequencing provided through the Roche/454 sequencing platform , is a powerful tool for the identification of transcripts and transcript variation in plant-pathogen interactions [19–21]. In this study, transcriptional variation in pepper during PepGMV infection was analyzed in healthy (mock), symptomatic, and recovered pepper leaves. To provide insight into the geminivirus-host molecular interactions in the pepper-PepGMV recovery system, we employed next generation transcriptome sequencing of viral-infected pepper leaves and through differential gene expression, identified key transcripts (genes) involved in this phenomenon. We report a total of 309 differentially expressed (DE) genes in the pepper-PepGMV recovery system with major differences in up- and down-regulated genes observed between healthy (mock) and symptomatic or recovered tissues. Of these, 246 have a known function including genes that are associated with defense responses such as pathogenesis-related (PR) proteins, reactive oxygen species, jasmonic acid, and ethylene signaling pathways. A set of genes with a novel role in defense responses were also identified thereby expanding our understanding of the molecular interactions that underlie the PepGMV-pepper recovery system.
Results & discussion
Sampling & transcriptome sequencing
Four-leaf stage Chili pepper plants (C. annuum L cv. Sonora Anaheim) were inoculated by bombardment with dimeric PepGMV clones (Figure 1A) [12, 13]. Symptom appearance was observed at 9 dpi in newly developed leaves (Figure 1B, “Symptoms”) and a decreased severity of symptoms was observed at 15 dpi in the next set of newly emerged leaves (Figure 1C, “Pre-Recovery”). The third pair of new leaves after inoculation emerged at 20 dpi and was symptomless (or showed an important reduction in the severity of the symptoms) (Figure 1D, “Recovery”). Symptomatic and recovered leaf tissues were sampled at 9 and 20 dpi, respectively, RNA isolated, and cDNA libraries constructed. As a control, total RNA from a pool of tissues (9, 15, and 20 dpi) was isolated from mock-inoculated healthy pepper leaves and a single cDNA library was constructed.
All of the libraries were sequenced using pyrosequencing on the Roche 454 GS20 platform . A total of 1,838,567 reads were obtained from nine pyrosequencing runs. The number of runs and reads for each condition were mock-inoculated (1 run, 222,558 reads), symptomatic tissue at 9 dpi (5 runs, 865,103 reads), and recovered tissue at 20 dpi (3 runs, 750,906 reads) (Table 1).
Differential expression analyses
To measure transcript abundances for each condition (Mock (M), Symptomatic (S), and Recovered (R)), pyrosequencing reads were aligned using the BLASTN algorithm  to a Capsicum annuum Reference Transcriptome (CaRT) that contains 32,220 total transcripts representing 12.5 Mb . Pyrosequencing reads were mapped against the CaRT dataset and alignments with an identity equal to or greater than 96.6% and alignment length equal to or greater than 30 bp were retained . The number of pyrosequencing reads mapped to a specific CaRT transcript was used to estimate transcript levels in each condition (Mock, Symptomatic and Recovered) . Read counts were normalized by using a relative frequency of reads, i.e., number of reads mapped for a given contig relative to the total number of reads for a specific library, a method used previously in efficient detection of DE genes .
To assess the technical reproducibility of the pyrosequencing runs, we compared different runs from a single condition (technical replicates) by calculating the transcript abundances from each independent run (p ≤ 2.26e-16) . Correlation within the technical replicates from the recovered leaf library (e.g., R-run1 vs. R-run2; R-run2 vs. R-run3; R-run1 vs. R-run3) revealed no differential gene expression between the 454 runs indicating a high degree of technical replication (Additional file 1). Examination of the number of reads, average length, and average quality on the technical replicates of the symptomatic library revealed a single run with low quality; this single run was discarded leaving a total of 4 runs with a total of 770,694 reads from the symptomatic library that were used in downstream analyses. The high degree of technical replication between the different 454 runs of the symptomatic leaf library is shown in Additional file 2.
The Fisher’s exact test was used to estimate the probability of significant differences among the transcript abundances in pair-wise comparisons (Mock vs. Symptomatic; Mock vs. Recovered; Symptomatic vs. Recovered) and significant differences were estimated using p ≤ 1.5e-6 (α = 0.1; probability of Error Type I) and the Bonferroni correction for multiple testing (Additional file 3). Comparisons of the symptomatic and recovered tissue relative to control leaves (mock-inoculated), revealed a total of 309 differentially expressed (DE) genes with 254 and 264 DE genes in symptomatic and recovered leaves relative to mock-inoculated leaves, respectively (Figure 2A). In symptomatic tissue, 147 and 107 genes out of 254 were up- and down-regulated relative to mock-inoculated leaves, respectively; whereas the recovered leaves, of 264 genes, 145 were up-regulated and 119 were down-regulated relative to mock-inoculated leaves (Figure 2A). Two genes were found up-regulated in recovered tissue relative to symptomatic tissue; one encodes a phosphodiesterase and the second gene has no known function. Expression patterns of the 309 DE genes in mock-inoculated, symptomatic and recovered tissues (Figure 2B), as well as hierarchical clustering of fold expression of the 309 DE genes in symptomatic and recovered tissue vs. mock-inoculated tissue (S/M, R/M) and recovered vs. symptomatic (R/S), revealed that expression profiles in symptomatic and recovered tissues are similar to each other (Additional file 4). Principal Component Analysis (PCA), showed that comparisons of fold-change of the 309 DE genes between the two stages (R/S) were significantly different than comparisons of the two stages relative to mock-inoculated tissue (S/M, R/M) (Figure 2C).
A comparative analysis of DE genes (up- and down-regulated) in both comparisons, symptomatic vs. mock and recovered vs. mock, showed that nearly 80% of the genes were in common to both groups; 124 up-regulated genes and 85 down-regulated genes were DE in both symptomatic vs. mock-inoculated and recovered tissue vs. mock-inoculated (Figure 2D). Some genes were differentially expressed only in symptomatic or recovered vs. mock-inoculated. Of these, 23 were up-regulated and 22 were down-regulated in symptomatic tissue, and in the recovered tissue, 21 up- and 34 down-regulated genes were identified. However, less stringent parameters revealed that most of these genes were also differentially regulated in both conditions (data not shown). Thus, the Fisher’s exact test for sensitivity failed to detect these genes as DE and examination of the expression differences indicate they were near the cutoff for annotating genes as differentially expressed (data not shown). Due to restrictive criteria used to define differential expression in this study, these data should be considered an under-estimation of differential expression in the response of pepper to PepGMV.
A BLAST search (E ≤ 1e-06) with the 309 (168 up- and 141 down-regulated) DE genes revealed that 246 (79.6%) matched an entry in the National Center for Biotechnology Information (NCBI) non-redundant database. Gene Ontology (GO, [26, 27]) associations were assigned by performing a BLASTX  search against the predicted Arabidopsis thaliana proteome [27, 28] (E ≤ 1e-06) and transitively assigning GO terms for the biological process category from A. thaliana to the corresponding pepper genes. GO terms were further reduced to GO Slim terms by using GOTermMapper (http://go.princeton.edu/cgi-bin/GOTermMapper). Analysis of the DE genes revealed that 219 (70.8%) of the 309 transcripts had a significant alignment to the A. thaliana proteome. Of these, 136 were up-regulated and 83 were down-regulated. A total of 32 and 58 pepper up- and down-regulated genes, respectively, did not have a significant match to A. thaliana proteome and as a consequence, were not assigned a GO annotation.
GO Slim terms associated with up-regulated transcripts (168 total genes) in symptomatic and recovered relative to the mock-inoculated tissues revealed 48 and 50 genes with GO associations classified as “responses to stress” in symptomatic and recovered tissues, respectively (Figure 3A). GO Slim terms assigned to down-regulated genes in symptomatic and recovered leaves (141 total genes) included 28 and 33 genes associated with “biosynthetic process”, respectively. In both recovered and symptomatic tissue, 18 down-regulated genes were associated with “catabolic process” and 18 down-regulated genes were associated with “photosynthesis” (Figure 3B). The GO results are consistent with the hypothesis that biotic stress marks a transition from growth and reproduction to physiology and metabolism tailored for defense responses  including a programmed down-regulation of primary metabolism.
Signaling and pathogen response genes
To further characterize the PepGMV-pepper recovery process or phenomenon, genes encoding proteins implicated in defense functions, such as PR proteins, reactive oxygen species, jasmonic acid biosynthesis, and ethylene signal transduction were examined in depth in our set of differentially expressed genes. Although not all of the components for each of defensive pathway or response were present in our datasets, several well-known markers for these pathways were identified suggesting that these defense responses are activated in the pepper-PepGMV interaction. Transcripts encoding the pathogenesis-related proteins, PR5 (Pepper00302) and two members of PR1 gene family (Pepper27140, Pepper31625), were highly represented during PepGMV infection (Table 2). Genes involved in jasmonic and ethylene-mediated defense responses, b-CHI (basic chitinase; Pepper32368) and HEL (hevein-like; Pepper32274), were up-regulated in infected leaves. Interestingly, the mRNA encoding the final step of ethylene biosynthesis , ACC oxidase (Pepper01276), was up-regulated in both conditions (Table 2). In contrast, the EIN3 (ethylene insensitive 3; Pepper00009) transcription factor, which is part of the ethylene signaling pathway  was down-regulated in infected tissues (Table 2). The transcripts for components of the jasmonic acid pathway, OPR1 (12-oxophytodienaote reductase 1; Pepper26071) and LOX1 (lipoxygenase; Pepper31749), were elevated in both conditions (Table 2).
Transcripts involved in oxidative stress and redox signaling were up-regulated in both conditions relative to mock-inoculated tissue including CAT (catalase; Pepper28784) which is involved in scavenging H2O2, consistent with results reported by  in which CAT is induced in PepGMV-infected Capsicum chinense plants. Recent data suggests that annexins may play an important role in conferring oxidative protection  and elevated levels of mRNAs that encode for annexin (ANN) (Pepper25924 and Pepper30410) were observed during PepGMV infection. GST1 (glutathione S-transferase) is hypothesized to enhance oxidative stress tolerance  and is induced by salicylic acid, jasmonic acid, ethylene, and hydrogen peroxide . In this study, GST1 (Pepper28222) was induced in both symptomatic and recovered tissues relative to mock-inoculated tissues (Table 2), consistent with previous studies with common bean (Phaseolus vulgaris) and the geminivirus Bean dwarf mosaic virus (BDMV) .
One of the most up-regulated genes in this study encodes a kiwelling ripening-related protein 1 (Pepper05849; from here on referred to as RRP1). The mRNA for this gene was up-regulated 69.5 and 66.9 fold in symptomatic and recovered tissues relative to mock-inoculated tissues, respectively. RRP1 has not been previously reported as related as a plant defense-related gene, although, based on its characterized biochemical features, it may be a new type of PR protein .
Overall, these results are consistent with observations reported by  in which expression profiles of A. thaliana infected with the geminivirus Cabbage leaf curl virus (CaLCuV) were analyzed utilizing microarrays. The suite of genes identified in CaLCuV-infected A. thaliana leaves mirrored the DE genes identified in the pepper-PepGMV recovery system in this study . One exception is EIN3, which was down-regulated at 9 and 20 dpi in the pepper-PepGMV recovery system (Table 2) yet in the CaLCuV-Arabidopsis interaction the Arabidopsis ortholog was up-regulated at 12 dpi  suggesting differences in the response of these two dicotyledonous plants to viral infection with respect to the ethylene signaling pathway.
Experimental validation of differentially expressed genes by quantitative real-time PCR
Computationally determined differential expression was validated using quantitative reverse-transcription PCR assays (qRT-PCR). For this purpose, seven up-regulated genes from the Mock vs. Symptomatic and Mock vs. Recovered tissue comparisons were selected for qRT-PCR validation. These genes spanned a number of metabolic pathways as well as levels of computationally predicted expression levels (high, medium, and low; Table 3). As shown in Figure 4, the qRT-PCR assay confirmed the up-regulation of the seven selected genes computationally-predicted to be differentially expressed. An interesting difference, however, between the qRT-PCR results and the computationally-predicted results was noticed in the recovered tissue. In the qRT-PCR, the RNA levels detected for the seven genes analyzed were substantially reduced in the recovered tissue. On the other hand, the expression levels observed in the computationally-prediction analysis were basically similar in both symptomatic and recovered tissues.
RRP1 (Pepper05849) was computationally predicted to be up-regulated 69.5-fold and 66.9-fold in symptomatic and recovered tissue, respectively, relative to mock-inoculated leaves. qRT-PCR showed up to 138.6-fold expression difference in symptomatic tissue and 38.3-fold in recovered tissue (Table 3). Computational predictions and qRT-PCR analysis indicated that this transcript is highly abundant in symptomatic leaves (9 dpi), however the abundance decreases in recovered tissue (20 dpi) (Table 3, Figure 4). The role of PR genes in plant defense has been widely documented [39–43] and PR5 (thaumatin-like protein) has been used as a marker of systemic acquired resistance, SAR [42, 44]. Our computational analyses showed that Pepper00302 which encodes a PR5 was up-regulated 25-fold in symptomatic and 24.4-fold in recovered tissues. qRT-PCR results showed 6-fold and 2.3-fold up-regulation in symptomatic and recovered tissues, respectively, relative to mock-inoculated tissues (Figure 4, Table 3). Interestingly, similar results have been shown in the BG-3821 accession of C. chinense Jacq., that displays resistance to PepGMV infection .
Pepper25924, which encodes a predicted annexin (ANN4), was computationally predicted to be 7-fold up-regulated in symptomatic and recovered compared to mock-inoculated tissue. qRT-PCR assays showed 6.3-fold and 2.4-fold up-regulation in symptomatic and recovered tissues, respectively, compared to the mock-inoculated tissue (Figure 4, Table 3). Annexin genes have been reported to have peroxidase activity, and it has been hypothesized that annexins can sense reactive oxygen species and modulate endogenous reactive oxygen species responses . Studies in the PepGMV-resistant accession (BG-3821-R) of C. chinense have shown that PepGMV infection is able to trigger a reactive oxygen species-mediated response . Additionally, it is hypothesized that a failure to control reactive oxygen species can lead to cell death [42, 45] as the cellular damage resulting from high reactive oxygen species levels show hallmarks of necrosis . It is important to note that the recovery phenomenon in the pepper-PepGMV system does not involve programmed cell death or the hypersensitive response. Thus, the annexin peroxidase activity could enhance oxidative tolerance by regulating reactive oxygen species levels in the PepGMV-infected leaves thereby resulting in the absence of a hypersensitive phenotype.
The computational prediction for Pepper27731, which encodes a serine carboxypeptidase (SCP), is 6-fold up-regulated in symptomatic and recovered tissues relative to mock-inoculated tissues. qRT-PCR results showed this gene 6.7- and 2.7-times up-regulated in symptomatic and recovered tissues, respectively, relative to mock-inoculated leaves (Figure 4, Table 3). Serine carboxypeptidases have been identified in many plant species . The physiological role of serine carboxypeptidases in plant defense is unclear; however, characterization of this protein in different plant-pathogen systems suggests that it is required for the synthesis of defense compounds [46, 47].
Pepper28222 encodes a GST1 (glutathione-S-transferase 1) which has been demonstrated to be transcriptionally activated by reactive oxygen species . Computational predictions revealed that Pepper28222 was up-regulated 4.8-fold in symptomatic and 5.4-fold in recovered leaves relative to mock-inoculated leaves. qRT-PCR results were similar, 4.2- and 2.2-fold up-regulation in symptomatic and recovered tissues, respectively, relative to mock-inoculated leaves (Figure 4, Table 3). Studies in A. thaliana infected with Cauliflower mosaic virus, a DNA virus, have shown that the activity of this gene is associated with both local and systemic accumulation of H2O2.
Jasmonic acid is a key signal transducer in the production of phytoalexins , which are important compounds for plant defense. Pepper26071 encodes a 12-oxophytodienoic acid reductase (OPR1). OPR catalyzes the NADPH-dependent reduction of 12-oxophytodienoic acid (OPDA) into 3-oxo-2[(Z)-2’-pentyl]cyclopentane-1-octanoic acid (OPC-8:0), and this reaction is part of the biosynthetic pathway leading to jasmonic acid. Computationally, Pepper26071 was 2.7-fold up-regulated in symptomatic tissue and 2.5-fold up-regulated in recovered tissue relative to mock; qRT-PCR results showed this gene up-regulated 3.4-fold in symptomatic and 2.2-fold in recovered tissue relative to mock-inoculated leaves (Figure 4, Table 3).
Heat-shock proteins (HSP) are a central component of the cellular chaperone network and play a crucial role in maintaining protein homeostasis by re-establishing functional native conformations under environmental stress conditions . Interestingly, yeast two-hybrid studies in the geminivirus-plant interaction with the bipartite Abutilon mosaic virus (AbMV) suggested that these proteins may interact with the viral movement protein, and therefore, have a role during geminiviral cell-to-cell transport . Pepper31770, which encodes a HSP70 chaperone, was 2.6-fold up-regulated in symptomatic and 2.8-fold up-regulated in recovered tissues relative to the mock-inoculated leaves. qRT-PCR revealed 2.9- and 1.3-fold up-regulation in symptomatic and recovered tissues, respectively, relative to the mock-inoculated leaves (Figure 4, Table 3).
Overall, a good correlation (Pearson’s coefficient of determination; r2 = 0.7) between differential expression values obtained in our computational predictions and those obtained by qRT-PCR was observed (Additional file 5). Discrepancies between the methods to determine DE genes may be related to methodological differences including the normalization methods used in our RNA-seq dataset as it has been reported that different normalization procedures impact differential expression detection .
Temporal expression patterns of differentially expressed genes during PepGMV infection
The analysis of a subset of genes was further extended to a wider time course to determine the expression profiles during viral infection. Total RNA was isolated from pre-symptomatic (6 dpi), symptomatic (9 dpi), pre-recovery (15 dpi) and recovered (20 dpi) leaves (Figure 1) and pepper transcripts 05849, 00302 and 28222, encoding RRP1, PR5, and GST1 genes, respectively, were examined for expression through the time course. qRT-PCR assay results showed that transcript levels of these genes increased significantly prior to the symptomatic stage (6 dpi) suggesting that a major differential transcriptional activity occurs in infected plants prior to symptom development (Figure 5A-C). As the PepGMV infection progressed, the levels of RRP1, GST1 and PR5 transcripts levels decreased (Figure 5A-C). The expression patterns of RRP1, GST1 and PR5 were correlated with levels of two PepGMV transcripts. Rep (Replication-associated protein) and NSP (Nuclear shuttle protein) from viral component A and B, respectively (Figure 5D). The results showed high levels of the Rep and NSP transcripts at 6 and 9 dpi, which gradually decreased during the recovery process (15–20 dpi). These results are consistent with previous studies [12, 13, 52] in which viral ssDNA and mRNA concentrations peaked at 9–10 dpi when symptoms develop in the first set of leaves following inoculation (Figure 1). By 20 dpi (typically the third set of leaves after inoculation), both viral DNA and RNA concentrations showed a reduction of 40 to 60% of the highest peak [12, 13, 52].
Orthologous and paralogous clusters
To determine if orthologs of our DE genes were present in the related Solanaceae crop species, tomato (S. lycopersicum) and potato (S. tuberosum), we generated orthologous clusters of the predicted proteomes of C. annuum, S. lycopersicum and S. tuberosum using the OrthoMCL algorithm  (Figure 6A). A total of 7,470 clusters containing 24,635 proteins were identified (Additional file 6) that represents the core proteome for all three species. Lineage-specific and two-species specific clusters were also identified. A total of 271 clusters containing 660 proteins were unique to C. annuum and S. tuberosum whereas 690 clusters with 1,522 proteins were unique between C. annuum and S. lycopersicum (Figure 6A). The OrthoMCL algorithm also identifies close paralogs within a species and we were able to identify 36 clusters grouping 80 proteins restricted to C. annuum (Figure 6A).
A total of 210 out of the 309 DE genes were identified in 191 clusters (Figure 6B; Additional file 6); of these, 187 of the DE genes were within 175 clusters that were shared by all three species (Figure 6B). Not all the transcripts could be clustered and 99 pepper transcripts remained as singleton transcripts (data not shown). A total of 11 DE pepper genes lacked a S. tuberosum ortholog and were restricted to a C. annuumS. lycopersicum orthologous group whereas seven DE pepper genes lacked a S. lycopersicum ortholog and were restricted to a C. annuumS. tuberosum orthologous cluster. The lineage-specific C. annuum DE genes could be clustered into 5 clusters containing 10 genes (Figure 6B). Pepper transcripts that encode proteins related to the oxidative response (CAT, ANN4, ANN1 and GST1), pathogenesis-related protein 5 (PR-5), ethylene and jasmonic acid signaling (ACC, EIN3, OPR1 and LOX1) and the novel gene PRR1 were present within the three species orthologous clusters (Table 2; Additional file 6). Interestingly, transcripts encoding HEL (hevein-like protein) grouped in a lineage-specific cluster (Additional file 6). These data are consistent with previous reports that show a high degree of conservation within the Solanaceae family  and provide candidate genes for further investigation of the recovery process in solanaceous species.
Identification of new components in the PepGMV-pepper recovery system
It has been reported that next generation sequencing methods are an excellent tool for the discovery of novel genes [19–21, 57–59]. The results of this study suggest that the RRP1 gene (Pepper05849) may have a role in plant defense. Interestingly, RRP1 is not expressed in mock-inoculated (healthy) plants but in PepGMV-inoculated plants, expression is high in pre-symptomatic, symptomatic and recovered leaves during the establishment of viral infection (Figure 5A). To assess whether the RRP1 gene is exclusively expressed in response to viral infection, qRT-PCR assays were performed with pepper plants challenged with different pathogens. Pepper plants were infected with two DNA viruses (PepGMV and Pepper huasteco yellow vein virus (PHYYV)), two RNA viruses ((Tobacco etch virus (TEV) and Tobacco rattle virus, (TRV)), a bacterium (Xanthomonas campestris pv. vesicatoria), and infested with whiteflies (Bemisia tabaci), the insect vector for PepGMV (Figure 7). Total RNA from leaves infected with the different pathogens was isolated at 10 dpi and from leaves exposed to whiteflies for 20 days. Pepper-RRP1 expression levels were measured in all samples. Infection with the DNA viruses (PepGMV and PHYVV) and TEV, as well as whitefly infestation, resulted in induction of RRP1 expression (Figure 7). This result is consistent with transcriptome analysis of A. thaliana challenged with whiteflies in which genes involved in systemic acquired resistance, such as PR1 and PR5, were induced . The highest level of induction was observed in TEV-infected pepper plants suggesting that RRP1 may have a generalized role in the response to a wide range of biotic stress defense mechanisms.
The core histones (H3, H4, H2A, H2B) are the components of the nucleosome complex, the basic unit of chromatin . Geminiviral genomes are replicated and transcribed in infected plant cells through double-stranded DNA intermediates, which are then assembled into mini-chromosomes [3, 4, 62]. Recently, it has been shown that the host histone H3 interacts with both viral movement proteins, NSP and MP, creating a complex composed of H3, NSP, MP and viral DNA, suggesting that H3 may play a role in geminivirus cell-to-cell movement through the formation of a movement-competent complex . Interestingly, genes encoding all four core histones were identified as differentially expressed during PepGMV infection. Pepper transcripts 01647, 30474, 31867 and 31973 are predicted to encode histone H3, whereas pepper contigs 26934, 28934, 30839, 31472 and 31544 are predicted to encode H4. Similarly, seven pepper contigs are predicted to encode histone H2A-related proteins (25732, 31056, 31911, 32303, 32359, 32453 and 32284) and three contigs were predicted to encode histone H2B (32370, 32371 and 32456). All of these contigs were up-regulated in symptomatic and recovered tissues (Additional file 3). In the OrthoMCL analysis with the 309 DE genes, one of the large gene families within the C. annuum-S. lycopersicum-S. tuberosum orthologous clusters contains 24 histone H4-encoding proteins, of these, three were from C. annuum. Two C. annuum-S. lycopersicum-S. tuberosum orthologous clusters contained histone H3 while histone H2A was contained within a single C. annuum-S. lycopersicum-S. tuberosum cluster. Interestingly, two genes encoding histone H4 protein were present within a C. annuum-specific paralogous cluster (Additional file 6). Detection of DE genes encoding core histones in our dataset suggests that these might have an important role in the response of pepper to virus, including the recovery process.
It has been also reported that host plants methylate viral chromatin as a defense against geminiviruses [13, 64–66]. Thus, viral chromatin is a target for transcriptional gene silencing and post-transcriptional gene silencing  and previous studies have suggested that the geminivirus recovery system requires the host RNA-directed DNA methylation pathway [13, 64]. Indeed, methylation-deficient mutants of A. thaliana are hyper-susceptible to geminivirus infection and histone H3 methylated at lysine 9 (H3K9) is highly represented in viral chromatin . In addition, studies on chromatin structure and gene regulation have characterized a remodeling process called “histone replacement”, in which canonical histones (i.e., H2A) are substituted by histone variants . One of the variants of H2A is H2A.Z and recently it was shown that histone H2A.Z may function to maintain the repressed or active transcriptional states of a number of genes related to the systemic acquired resistance response in A. thaliana.
Recent studies suggest that histone H2A.Z may be involved in the induction of ERF1 (ETHYLENE RESPONSIVE FACTOR1) and b-Chi, genes that act downstream in the ethylene and jasmonic acid signaling pathways [61, 68, 69]. Pepper32368 encodes b-Chi and our computational predictions revealed that this gene is up-regulated 39.6 and 41.3 fold in symptomatic and recovered tissue, respectively (Additional file 3). Collectively, these results have yielded candidate genes to further investigate the potential role of host chromatin modification and post-transcriptional gene silencing in the pepper-PepGMV recovery system.
Virus-induced gene expression in plants has been studied in many host-virus models, including several geminivirus [32, 36, 38]. In most cases, however, the models permit examination of changes between healthy or non-inoculated plants and the infected, symptomatic tissues . The pepper-PepGMV recovery system affords the opportunity to examine susceptibility and recovery in the same system. Quantification and comparison of transcript abundances using deep transcriptome sequencing and the C. annuum reference transcriptome  allowed us to analyze the transcriptional status of PepGMV-infected plants during the initial symptom stage and subsequent recovered condition. Modification of transcript levels for many genes occurs prior to the appearance of symptoms with the highest peak of expressed genes observed around 6 dpi. Newly emerged leaves are nearly symptomless (recovery stage), thus analysis of DE genes suggests that several elements related to the defense machinery of the plant (PR proteins, reactive oxygen species, as well as jasmonic acid and ethylene signaling pathways) may contribute in pepper-PepGMV recovery system along with the previously reported PTGS and TGS mechanisms. Interestingly, novel genes, such as Pepper-RRP1 and histone proteins, were identified which may have a role in plant defense. The results presented in this study provide valuable information for our understanding of the underlying molecular mechanisms by which PepGMV-infected pepper plants recover from geminiviral infection.
Material & methods
Biological material & sampling
C. annuum L cv. Sonora Anaheim seeds, a susceptible cultivar [10, 12], were germinated in plant growth chambers under an 18 hr light/6 hr dark photoperiod at 26 to 28°C for two weeks. At the four true leaf stage, plants were inoculated with PepGMV (A + B) dimeric clones using biolistics as previously reported . Mock-inoculated plants were bombarded with gold particles alone. Leaf tissues were sampled at 6, 9, 15 and 20 dpi, frozen immediately in liquid nitrogen, and stored at −80°C until use.
Inoculation of PHYVV and TEV
Pepper plants at the four true leaf stage were inoculated with dimeric clones using biolistics as previously reported . Leaf tissue was sampled at 6 and 10 dpi, frozen immediately in liquid nitrogen, and processed. TEV was inoculated using sap from a infected tobacco plant and carborundum.
Agroinoculation of TRV vectors
Agrobacterium tumefaciens carrying the TRV vector was grown in LB medium containing kanamycin 50 mg/L (Bristol-Myers Squibb de México), carbenicillin 100 mg/L (Pfizer, New York, NY, USA) and rifampicin 50 mg/L (Bristol-Myers Squibb de México). The bacterial pellet was diluted in buffer (10 mM MgCl2, 10 mM MES pH 5.6, 150 μM acetosyringone) to achieve an optical density of 1. The bacterial suspension was infiltrated into the main vein of pepper leaves using a 1 cc syringe. Leaf tissues were sampled at 6 and 10 dpi, frozen immediately in liquid nitrogen, and processed.
Inoculation of X. campestris
X. campestris pv. vesicatoria was grown for 48 hrs in 5 ml of NYB medium (Casein Peptone 10 g/L, yeast extract 5 g/L, NaCl 5 g/L pH 7) at 27°C. The bacterial suspension was infiltrated into the main vein of pepper leaves using a 1 cc syringe. Inoculated plants were incubated at 27°C and leaf tissues were sampled at 6 and 10 dpi, frozen in liquid nitrogen, and processed.
Infestation with B. tabaci(whitefly)
Pepper plants at the four true leaf stage were exposed to a colony of B. tabaci (approximately 30 whiteflies). Leaf tissues were sampled at 20 days post exposure, frozen in liquid nitrogen, and processed.
RNA extraction & cDNA library preparation
Total RNA was extracted from frozen tissue using the combined method of Trizol and PureLink Micro-to-Midi Total RNA Purification System kit (Invitrogen, Carlsbad, CA). RNA purity was checked using the Agilent 2100 Bioanalyzer RNA 6000 Nano Assay chip (Agilent Technologies, Stockport, U.K). cDNA synthesis was performed from 3.5 μg of total RNA using the Message Amp-II kit (Ambion, Foster City, CA) following the manufacturer’s protocol as described earlier . For RRP1 expression, total RNA was extract from frozen tissues using Trizol (Invitrogen, Carlsbad, CA, U.S.A.).
cDNA samples were prepared for GS20-454 pyrosequencing as described previously [18, 70]. Nine runs of three cDNA libraries were performed resulting in 1,838,567 total reads. Sequences of the 454-G20 reads are publicly available for download from http://www.bioingenios.ira.cinvestav.mx:81/Joomla/ and in NCBI Sequence Read Archive (accession number SRA052606).
In silico differential expression analysis
Stand-alone BLAST software  was obtained from NCBI (http://www.ncbi.nih.gov). The 454-pyrosequencing reads were aligned using BLAST to the C. annuum Reference Transcriptome . A custom local MySQL database was constructed to store and query information from BLAST alignments. Using criteria described in , an alignment was considered significant if ≥ 30 bp aligned at ≥ 96.6% identity. One impact on differential expression detection is the normalization method, e.g., FPKM values (fragments per kilobase of exon model per million mapped reads). FPKM values are heavily affected by a relatively small proportion of highly-expressed genes and, as such, can introduce biased estimates of differential expression if these genes are differentially expressed across the conditions under comparison . Therefore, we elected to normalize transcript levels using a relative frequency of reads, i.e., number of reads mapped for a given contig relative to the total number of reads for a specific library. This method has proven to be efficient in the identification of DE genes . The fold change of DE genes was estimated by obtaining the ratio between the relative frequencies for the two conditions. The probability (p-value) and significant differences between the samples (symptomatic vs. mock; recovered vs. mock; symptomatic vs. recovered) were estimated using the Fisher’s exact test. To reduce false positives, the QVALUE software was used to adjust p-values obtained from the Fisher’s exact test . Changes in signal intensity of ± 1.45 or higher/lower between treatments were considered highly significant (p-value 1.5e-6; 90% confidence). However, we focused on DE genes using traditional criteria from microarray experiments in which the cut-off threshold for up-regulated is ≥ 2 and down-regulated genes is ≤ 0.5.
Functional annotation and OrthoMCL analysis
Transctripts were annotated using BLASTX  searches against non-redundant polypeptides database from NCBI and the A. thaliana proteome (TAIR10; arabidopsis.org) as described previously . GO associations [26, 27] were made by GOTermMapper (http://go.princeton.edu/cgi-bin/GOTermMapper).
Orthologs and close paralogs were identified in the three predicted proteomes using OrthoMCL (v 1.4)  using the default parameters with an E-value cutoff of 1e-10. Transposable elements were filtered out to avoid clusters comprised entirely of transposable elements.
Evaluation of genes expression by quantitative reverse-transcription PCR (qRT-PCR)
The same total RNA samples used for 454-pyrosequencing sequencing were used in the qRT-PCR validation experiments. Biological replicates confirmed the qRT-PCR results. DNA contamination was determined by running a PCR under the same conditions for the RNA samples. DNA-free RNA (1 μg) was used for cDNA synthesis using Superscript II Reverse Transcriptase (Invitrogen, Carlsbad, CA). qRT-PCR was carried out as previously described . The primers used in this study are described in Additional file 7.
Capsicum annuum Reference Transcriptome
Pepper golden mosaic virus
Pepper huasteco yellow vein virus
Post-transcriptional gene silencing
Quantitative reverse transcription PCR
Reactive oxygen species
Single stranded DNA
Small RNA of viral origin
Tobacco etch virus
Transcriptional gene silencing
Tobacco rattle virus.
Bisaro DM: Geminivirus DNA Replication. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press; 1996.
Rojas MR, Hagen C, Lucas WJ, Gilbertson RL: Exploiting chinks in the plant's armor: evolution and emergence of geminiviruses.Annu Rev Phytopathol 2005, 43: 361-394. 10.1146/annurev.phyto.43.040204.135939
Torres-Pacheco I, Garzon-Tiznado JA, Brown JK, Becerra-Flora A, Rivera-Bustamante RF: Detection and distribution of geminiviruses in Mexico and Southern United States.Phytopathology 1996, 86: 1186-1192. 10.1094/Phyto-86-1186
Mendez-Lozano J, Torres-Pacheco I, Fauquet CM, Rivera-Bustamante RF: Interactions between geminiviruses in a naturally occurring mixture:Pepper huasteco virusandPepper golden mosaic virus.Phytopathology 2003, 93: 270-277. 10.1094/PHYTO.2003.93.3.270
Carrillo-Tripp J, Lozoya-Gloria E, Rivera-Bustamante RF: Symptom remission and specific resistance of pepper plants after infection byPepper golden mosaic virus.Phytopathology 2007, 97: 51-59. 10.1094/PHYTO-97-0051
Góngora-Castillo E, Fajardo-Jaime R, Fernández-Cortés A, Jofre-Garfias AE, Lozoya-Gloria E, Martinez O, Ochoa-Alejo N, Rivera-Bustamante R: The capsicum transcriptome DB: a "hot" tool for genomic research.Bioinformation 2012, 8: 43-47. 10.6026/97320630008043
Heil M, Ibarra-Laclette E, Adame-Alvarez RM, Martinez O, Ramirez-Chavez E, Molina-Torres J, Herrera-Estrella L: How plants sense wounds: damaged-self recognition is based on plant-derived elicitors and induces octadecanoid signaling.PLoS One 2012, 7: e30537. 10.1371/journal.pone.0030537
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium.Nat Genet 2000, 25: 25-29. 10.1038/75556
Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, et al.: The Gene Ontology (GO) database and informatics resource.Nucleic Acids Res 2004, 32: D258-D261. 10.1093/nar/gkh036
Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander DL, Garcia-Hernandez M, et al.: The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools.Nucleic Acids Res 2012, 40: D1202-D1210. 10.1093/nar/gkr1090
Konopka-Postupolska D, Clark G, Goch G, Debski J, Floras K, Cantero A, Fijolek B, Roux S, Hennig J: The role of annexin 1 in drought stress in Arabidopsis.Plant Physiol 2009, 150: 1394-1410. 10.1104/pp.109.135228
Lieberherr D, Wagner U, Dubuis PH, Metraux JP, Mauch F: The rapid induction ofglutathione S-transferasesAtGSTF2 and AtGSTF6 by avirulentPseudomonas syringaeis the result of combined salicylic acid and ethylene signaling.Plant Cell Physiol 2003, 44: 750-757. 10.1093/pcp/pcg093
Seo YS, Jeon JS, Rojas MR, Gilbertson RL: Characterization of a novel Toll/interleukin-1 receptor (TIR)-TIR gene differentially expressed in common bean (Phaseolus vulgariscv. Othello) undergoing a defence response to the geminivirusBean dwarf mosaic virus.Mol Plant Pathol 2007, 8: 151-162. 10.1111/j.1364-3703.2007.00379.x
Tuppo L, Giangrieco I, Palazzo P, Bernardi ML, Scala E, Carratore V, Tamburrini M, Mari A, Ciardiello MA: Kiwellin, a modular protein from green and gold kiwi fruits: evidence of in vivo and in vitro processing and IgE binding.J Agric Food Chem 2008, 56: 3812-3817. 10.1021/jf703620m
Bienert MD, Delannoy M, Navarre C, Boutry M: NtSCP1 from tobacco is an extracellular serine carboxypeptidase III that has an impact on cell elongation.Plant Physiol 2012, 158: 1220-1229. 10.1104/pp.111.192088
Mugford ST, Qi X, Bakht S, Hill L, Wegel E, Hughes RK, Papadopoulou K, Melton R, Philo M, Sainsbury F, et al.: A serine carboxypeptidase-like acyltransferase is required for synthesis of antimicrobial compounds and disease resistance in oats.Plant Cell 2009, 21: 2473-2484. 10.1105/tpc.109.065870
Sobajima H, Takeda M, Sugimori M, Kobashi N, Kiribuchi K, Cho EM, Akimoto C, Yamaguchi T, Minami E, Shibuya N, et al.: Cloning and characterization of a jasmonic acid-responsive gene encoding 12-oxophytodienoic acid reductase in suspension-cultured rice cells.Planta 2003, 216: 692-698.
Bullard JH, Purdom E, Hansen KD, Dudoit S: Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments.BMC Bioinformatics 2010, 11: 94. 10.1186/1471-2105-11-94
Rensink WA, Lee Y, Liu J, Iobst S, Ouyang S, Buell CR: Comparative analyses of six solanaceous transcriptomes reveal a high degree of sequence conservation and species-specific transcripts.BMC Genomics 2005, 6: 124. 10.1186/1471-2164-6-124
Wang B, Guo G, Wang C, Lin Y, Wang X, Zhao M, Guo Y, He M, Zhang Y, Pan L: Survey of the transcriptome ofAspergillus oryzaevia massively parallel mRNA sequencing.Nucleic Acids Res 2010, 38: 5075-5087. 10.1093/nar/gkq256
Fernandez D, Tisserant E, Talhinhas P, Azinheira H, Vieira A, Petitot AS, Loureiro A, Poulain J, Da Silva C, Silva Mdo C, Duplessis S: 454-pyrosequencing ofCoffea arabicaleaves infected by the rust fungusHemileia vastatrixreveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant-rust interaction.Mol Plant Pathol 2012, 13: 17-37. 10.1111/j.1364-3703.2011.00723.x
Kempema LA, Cui X, Holzer FM, Walling LL: Arabidopsis transcriptome changes in response to phloem-feeding silverleaf whitefly nymphs. Similarities and distinctions in responses to aphids.Plant Physiol 2007, 143: 849-865.
Zhou Y, Rojas MR, Park MR, Seo YS, Lucas WJ, Gilbertson RL: Histone H3 interacts and colocalizes with the nuclear shuttle protein and the movement protein of a geminivirus.J Virol 2011, 85: 11821-11832. 10.1128/JVI.00082-11
March-Diaz R, Garcia-Dominguez M, Lozano-Juste J, Leon J, Florencio FJ, Reyes JC: Histone H2A.Z and homologues of components of the SWR1 complex are required to control immunity in Arabidopsis.Plant J 2008, 53: 475-487.
Solano R, Stepanova A, Chao Q, Ecker JR: Nuclear events in ethylene signaling: a transcriptional cascade mediated by ETHYLENE-INSENSITIVE3 and ETHYLENE-RESPONSE-FACTOR1.Genes Dev 1998, 12: 3703-3714. 10.1101/gad.12.23.3703
Larsen PB, Cancel JD: A recessive mutation in the RUB1-conjugating enzyme, RCE1, reveals a requirement for RUB modification for control of ethylene biosynthesis and proper induction of basic chitinase and PDF1.2 in Arabidopsis.Plant J 2004, 38: 626-638. 10.1111/j.1365-313X.2004.02068.x
Vega-Arreguin JC, Ibarra-Laclette E, Jimenez-Moraila B, Martinez O, Vielle-Calzada JP, Herrera-Estrella L, Herrera-Estrella A: Deep sampling of the Palomero maize transcriptome by a high throughput strategy of pyrosequencing.BMC Genomics 2009, 10: 299. 10.1186/1471-2164-10-299
This work was supported by Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación (SAGARPA) and Consejo Nacional de Ciencia y Tecnologia (Conacyt) . We acknowledge support from Conacyt-Mexico to Elsa Góngora-Castillo (Ph. D. fellowship). We acknowledge and thank Dr. Robin Buell for critical reading of the manuscript and helpful comments. We would also like to acknowledge the anonymous reviewers that made important suggestions to improve this manuscript.
Authors and Affiliations
Departamento de Ingeniería Genética, Carretera Irapuato-León, Centro de Investigación y de Estudios Avanzados del I.P.N (Cinvestav)-Unidad Irapuato, Km 9.6 Libramiento Norte, Irapuato, Gto, 36821, México
Elsa Góngora-Castillo, Diana L Trejo-Saavedra & Rafael F Rivera-Bustamante
Laboratorio Nacional de Genómica para la Biodiversidad (Langebio), Cinvestav-Irapuato, Carretera Irapuato-León, Km 9.6 Libramiento Norte, Irapuato, Gto, 36821, México
The authors declare that they have no competing interests.
EGC designed methods and experiments; performed the plant inoculations, RNA isolations, qRT-PCR assays, and bioinformatics analysis; analyzed and interpreted the data; and wrote the manuscript. EIL designed methods and experiments, constructed the cDNA library for pyrosequencing, performed the statistical analysis and collaborated to draft the manuscript. DTS designed methods and experiments, carried out plant inoculations, RNA isolation and qRT-PCR assays. RRB coordinated the project, designed experiments and collaborated in analyzing data and writing the manuscript. All the authors have read and approved the final manuscript.
Additional file 1:Correlation of aligned 454-reads to theCapsicum annuumReference Transcriptome (CaRT) between different 454-runs from the recovered leaf cDNA library. The best fitting linear correlations for each pair-wise comparison was calculated (p < 2.26e-16): run1 vs. run2 (r2 = 0.9666), run1 vs. run3 (r2 = 0.9698) and run2 vs. run3 (r2 = 0.9703). (PDF 198 KB)
Additional file 2:Correlation of aligned 454-reads to theCapsicum annuumReference Transcriptome (CaRT) between different runs derived from the symptomatic leaf cDNA library. The best fitting linear correlations for each pair-wise comparison was calculated (p < 2.2e-16): run1 vs. run2 (r2 = 0.9657), run1 vs. run3 (r2 = 0.9644), run1 vs. run4 (r2 = 0.9705), run2 vs. run 3 (r2 = 0.9658), run2 vs. run5 (r2 = 0.9667), run3 vs. run5 (r2 = 0.9663). Low quality was observed on run4 and therefore it was discarded from the downstream analysis (data not shown). (PDF 270 KB)
Additional file 4:Hierarchical clustering of differentially expressed genes identified in PepGMV-infected pepper plants. A total of 309 genes were identified (fold change of at least ± 2 and p-value ≤ 1.57e-06) and the ratio for each comparison (Symptomatic (S) vs. Mock (M); Recovered (R) vs. Mock; Recovered vs. Symptomatic) was used for the analysis. Clustering was performed using the Smooth correlation and average linkage clustering in GeneSpring GX 7.3.1 software (Agilent Technologies®). Green indicates down-regulated, red up-regulated and black unchanged values, as shown on the color scale at the side of the figure. (PDF 120 KB)
Additional file 6:Table S1. Clusters of orthologous and paralogous genes families inS. lycopersicum, S. tuberosumandC. annuumspecies. Table S2. Clusters of orthologous and paralogous of differentially expressed genes in S. lycopersicum, S. tuberosum and C. annuum species. Predicted peptides from S. lycopersicum, S. tuberosum transcripts were clustered using OrthoMCL . (XLSX 1004 KB)
This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License (
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Góngora-Castillo, E., Ibarra-Laclette, E., Trejo-Saavedra, D.L. et al. Transcriptome analysis of symptomatic and recovered leaves of geminivirus-infected pepper (Capsicum annuum).
Virol J9, 295 (2012). https://doi.org/10.1186/1743-422X-9-295