Genome-wide analysis of Dof transcription factors reveals functional characteristics during development and response to biotic stresses in pepper

The DNA-binding with one zinc finger proteins (Dofs) are a plant-specific family of transcription factors. The Dofs are involved in a variety of biological processes such as phytohormone production, seed development, and environmental adaptation. Dofs have been previously identified in several plants, but not in pepper. We identified 33 putative Dof genes in pepper (CaDofs). To gain an overview of the CaDofs, we analyzed phylogenetic relationships, protein motifs, and evolutionary history. We divided the 33 CaDofs, containing 25 motifs, into four major groups distributed on eight chromosomes. We discovered an expansion of the CaDofs dated to a recent duplication event. Segmental duplication that occurred before the speciation of the Solanaceae lineages was predominant among the CaDofs. The global gene-expression profiling of the CaDofs by RNA-seq analysis showed distinct temporal and pathogen-specific variation during development and response to biotic stresses (two TMV strains, PepMoV, and Phytophthora capsici), suggesting functional diversity among the CaDofs. These results will provide the useful clues into the responses of Dofs in biotic stresses and promote a better understanding of their multiple function in pepper and other species.


Results
Identification of the Dof gene family in pepper. We identified 33 putative CaDofs from the Capsicum annuum 'CM334' (hereafter 'CM334') genome. Additionally, we obtained 21 putative CaDofs from the Plant Transcription Factor Database (Plant TFDB). We aligned the 33 putative CaDofs from the 'CM334' genome with 21 corresponding genes from the Plant TFDB. All of the candidate CaDofs from the Plant TFDB were matched to corresponding genes among the 33 putative CaDofs. Therefore, we used only the 33 CaDofs from C. annuum 'CM334' for further investigation. To verify the sequences of the CaDofs, we aligned all of the CaDofs with corresponding transcripts or genes from the 'CM334' RNA-seq 40 , pepper EST 41 , and C. annuum 'Zunla-1' CDS 42 databases and manually corrected. We designated the generated CaDofs as CaDof01-33 according to their locations on the chromosomes. The number of Dofs in the pepper genome was similar to those in Arabidopsis (36 Dofs), rice (30 Dofs) and tomato (34 Dofs) 16,21 . All of the CaDofs contained the highly conserved Dof domain including the zinc finger motif ( Supplementary Fig. S1). The conserved Dof domain consists of 52 basic residues located in the N-terminal region. Twenty-six of the 52 residues were perfectly conserved (100% conserved) in all 33 CaDofs. The Dof domains were also highly conserved among the Dofs of other plants, such as Arabidopsis, rice, and tomato ( Supplementary Fig. S2). Information about the CaDofs; such as the gene ID, position on the chromosomes, number of nucleotides and amino acids, isoelectric point (pI), and molecular weight (MW); is presented in Table 1. The size of the coding sequence (CDS) of the CaDofs ranged from 525 bp to 1,512 bp, with deduced polypeptide sizes ranging from 173 amino acids to 503 amino acids. The pIs of the deduced CaDofs ranged from 3.61 to 9.79, and the MWs ranged from 15.8 kD to 54.27 kD.

Phylogenetic relationships of the Dof family in pepper and tomato plants.
To examine the phylogenetic relationships between the Dof families in pepper and tomato, we performed multiple sequence alignment of the 33 CaDofs and 34 tomato Dofs (SlDofs) using ClustalW. We constructed an unrooted phylogenetic tree using the neighbor-joining method with alignment of the full-length Dofs from tomato and pepper. Based on the resulting phylogenetic tree, we categorized the CaDofs and SlDofs into four major groups (A, B, C, and D) and seven subgroups ( Fig. 1 and Supplementary Fig. S3). We named the pepper groups based on the corresponding groups in tomato 21 . Groups B and C were the largest major groups, each comprising 17 genes (accounting for 25.3%) in pepper and tomato. Groups B and C were divided into two subgroups: B1 and B2 and C1 and C2, respectively. Subgroup B1, which consisted of nine CaDofs and eight SlDofs, was the largest subgroup in both pepper and tomato. Group A and subgroups C1, C2, and D1 contained 10,9,8, and 13 members, respectively, from the two plants. Subgroups B2 and D2 were not clearly separated and clustered into the same subgroup with four CaDofs and four SlDofs. Most of the CaDofs were clustered with SlDofs, and each subgroup had a high bootstrap value, suggesting that the CaDofs have a close phylogenetic relationship with that of the SlDofs.

Identification of conserved motifs and analysis of gene structures among the CaDofs.
To identify the diversity and conservation among the CaDofs, we investigated the conserved motifs using the MEME suite (http://meme-suite.org/tools/meme). A total of 25 distinct motifs were predicted among all 33 CaDofs Scientific RepoRts | 6:33332 | DOI: 10.1038/srep33332 (Fig. 2). The size of the identified motifs ranged from 8 amino acids to 113 amino acids. The consensus sequences of the motifs are presented in Supplementary Table S1. The motifs of closely related genes within each group in the phylogenetic tree shared common sequences and positions (Fig. 2). Most of the subgroups included several specific motifs in the C-terminal regions. Motif 1 was uniformly found in all 33 CaDofs and represented the pepper Dof domain ( Fig. 2 and Supplementary Table S1). Group A and subgroups B2 and D2 represent the simplest composition of motifs, with no specific motifs except for the Dof domain. Although subgroup B1 had several motifs (motifs 2, 7, 17, 18, and 19) at the C-terminal regions, it appeared to be characterized by motifs 11 and 12 at the N-terminal region. One CaDof (CaDof31) in subgroup B1 did not contain motifs 11 and 12. Subgroup C1 contained four specific motifs on both sides of the Dof domain: motifs 7 and 14 in the N-terminal region and motifs 24 and 23 in the C-terminal region. Subgroup C2 contained one specific motif (motif 25). Subgroup D1, which showed the most diverse composition of motifs among the subgroups, was characterized by five motifs (motifs 3, 4, 5, 8, and 10). Motif 8 was located in the N-terminal region in all of the D1 members except for CaDof19. Conserved motifs 3, 4, 5, and 10 at the C-terminus were detected in seven of the eight Dofs in subgroup D1. CaDof04 in subgroup D1 contained only two motifs (motif 8 and the Dof domain), possibly because it had a shorter gene sequence than the other D1 members.
In order to obtain further insight into the structural diversity of the CaDofs, we analyzed the exon-intron organization. As shown in Fig. 3, 22 of the CaDofs had no introns, whereas seven and four CaDofs had one and two introns, respectively. Those exon-intron structures were similar to those of the Dofs in Arabidopsis, rice, and tomato 16,21 . Most of the introns were generally located upstream of the Dof domain; only three genes (CaDof04, 29, and32) contained introns downstream of the Dof domain. Most of the genes within a given group showed similar organizations of exons and introns. Most of the genes from subgroup D1 contained one intron, and the intron sizes and positions were comparable among the genes in D1. Subgroups B2 and D2 contained no introns. The structures of the CaDofs in group A and subgroups B1, C1, and C2 were more variable than those in subgroups D1, B2, and D2. The distribution of motifs together with the gene structure within each group indicated the phylogenetic relationships among the CaDofs. Those results indicate that the conservation and specificity of gene characters among the subgroups could be related to functional diversity among the Dofs.

Chromosomal location and gene duplication of CaDofs.
To determine the genomic localization of the CaDofs, we aligned each CaDof as query sequence against the pepper genome database using BLAST. The positions of the CaDofs on the chromosomes had an uneven distribution. Twenty-six genes were located on 8 of the 12 chromosomes, with no genes mapping to chromosomes 7, 8, 9, and 12 ( Fig. 4). Nine CaDofs were mapped to chromosome 2. Among them, CaDof05, 06, and 07 and CaDof08 and 09 were clustered together, respectively. Chromosome 6 included six CaDofs, which were clustered in two groups (CaDof18, 19, and 20 and CaDof21 and 22, respectively) on the end of the long arm. Chromosomes 3 and 10 each contained three genes, and CaDof23 and 24 were clustered together. Chromosome 5 contained two genes. Chromosomes 1, 4 and 11 each contained only one gene. Seven of the 33 CaDofs (CaDof27-33) were assigned to chromosome 0 that means unassigned pepper sequences to any of the 12 pepper chromosomes.
To detect potential duplications of the CaDofs, we identified 10 putative pairs of CaDof paralogs based on the chromosomal localization and the neighbor-joining tree (Figs 1 and 4). Each pair of paralogs shared high sequence similarity. The identified paralogs are connected by the lines shown in Fig. 4. Seven of the paralog pairs were putative segmental duplication events, and three pairs were distributed at the same chromosome level. Only one putative pair (CaDof23 and 24) was possibly the result of tandem duplication. Overall, those results suggest that segmental duplication predominated in the expansion of the CaDof family in peppers, although tandem duplication was also involved.
In order to estimate the dates of the duplication events, we estimated the nonsynonymous substitution rate (Ka), the synonymous substitution rate (Ks), and the Ka/Ks ratio. The approximate dates of the duplication events are presented in Table 2. The nine segmental duplications (Fig. 4) of the CaDofs took place from 194.48 million years ago (Mya) to 37.06 Mya, and one tandem duplication (Fig. 4) occurred approximately 1.45 Mya. With the speciation data between pepper and tomato around 19.2 Mya 30 , all the segmental duplications occurred before the pepper/tomato split, which is also supported by the phylogenic tree shown in Fig. 1. The Ka/Ks ratios of five segmental duplication pairs and the tandem duplication pair were > 0.3, suggesting that significant functional divergence among those CaDofs might have occurred after the duplication events. The Ka/Ks ratios of all the duplicate pairs were < 1.0, indicating that the CaDofs evolved under negative selection acting against protein-coding changes. Expression patterns of CaDofs in various developmental stages. To investigate the involvement of the CaDofs in pepper growth and development, we generated a heat map of the global transcription patterns of the CaDofs based on publicly available RNA-seq data for five tissues (root, stem, leaf, pericarp, and placenta), including seven developmental stages of the pericarp and placenta (Fig. 5). The expression profiles of each CaDof revealed various patterns in different organs and stages. CaDof13, 17,18,19,21,27, and 30 were expressed at relatively high levels, while CaDof14, 15, 16, 25, 31, and 32 were expressed at relatively low levels in all the tested tissues. The expression levels of CaDof11, 29, and 08 were high in the root and stem and higher in the placenta from 6 days post anthesis (PL_6DPA) to the mature green (MG) stage compared with those in the placenta at other stages, but they were low in the pericarp at all stages. CaDof22 was expressed at higher levels in the root and stem than in the other tissues. Although CaDof21 was expressed at relatively high levels in all the tested tissues, the expression levels in the placenta at several stages [PL_16 and 25DPA, MG, and breaker (B)] were higher than those in the other tissues. In comparison, CaDof30 was also highly expressed in the placenta, especially at the B stage. CaDof18 was constitutively expressed in most of the tested tissues, with especially high expression in some PR stages (PR 6, 16, and 25DPA and MG). CaDof27 showed the highest expression level in the roots among all the CaDofs. CaDof04 was expressed at a high level only in roots and was undetectable in the other organs. Taken together, those results indicate that the CaDofs are possibly involved in biological functions during plant development.  Table S1.   Responses of CaDofs to TMV, PepMoV, and P. capsici pathogen challenges. 'CM334' is well known as a powerful resource for the development of resistant lines against several viruses and P. capsici [37][38][39] . To examine whether CaDofs are involved in responses to biotic stresses in pepper, we analyzed the expression profiles of the CaDofs using RNA-seq data from 'CM334' challenged with TMV Pathotype 0 (TMV-P0), TMV-P2, PepMoV, and P. capsici (Supplementary Table S2). 'CM334' is resistant to TMV-P0, PepMoV, and P. capsici but susceptible to TMV-P2. We inoculated 'CM334' leaves with TMV-P0, TMV-P2, PepMoV, and P. capsici and then sampled at 0, 1, 2, and 3 days post inoculation. We used the leaf samples to construct RNA-seq libraries, which we subsequently sequenced. We then generated heat maps of the transcript profiles of the CaDofs with each gene expression value normalized by each control sample. The expression patterns of the CaDofs were differentially regulated in the resistant and susceptible responses against TMV-P0 and TMV-P2, respectively (Fig. 6). CaDof32 was highly up-regulated only in the resistant response against TMV-P0, while the expression levels of CaDof07, 12, 14, and 15 were increased in the susceptible response against TMV-P2. Some of the CaDofs, such as CaDof9, 10, 11, 17, and 20, were up-regulated in both responses but showed temporally different expression patterns. Among those, CaDof10 showed earlier and higher levels of expression in the resistant response. After PepMoV inoculation, some CaDofs were up-regulated, as in the responses to TMV-P0 and TMV-P2, but showed different patterns, such as those of CaDof09, 10, 11, 20, and 32 (Fig. 6). CaDof04 showed earlier and higher levels of expression after PepMoV inoculation but did not show significant changes after other pathogen inoculations. Other CaDofs showed no detectable changes compared with controls and among treatments. In response to P. capsici inoculation, nine CaDofs were expressed at higher levels compared with the controls (Fig. 6). Among those, CaDof17 showed higher expression levels compared with those in control and the other viruses treatments. CaDof04 was significantly down-regulated. In addition, two genes (CaDof05 and 31) showed particularly strong expression levels after P. capsici inoculation compared with those after the other treatments. Overall, the expression profiles in response to the biotic stresses suggest diversity and conservation of the biological functions of CaDofs during plant-microbe interactions.

Discussion
The Dofs are a major group of plant-specific TFs involved in diverse functions 4,8 . Several studies have addressed the genome-wide identification, function, and evolution of Dofs in Arabidopsis, rice 16 , poplar 17 , algae,  19 , soybean 20 , Chinese cabbage 22 , and tomato 21 . No previous study has reported the genome-wide identification and transcriptional profiles in response to biotic stresses of Dofs in pepper, an agriculturally important member of the Solanaceae. We identified a total of 33 CaDofs based on two pepper genomes, C. annuum 'CM334' and C. annuum 'Zunla-1' . The number of CaDofs identified was similar to those in other plants, such as the 36 AtDofs in Arabidopsis, the 30 OsDofs in rice 16 , and the 34 SlDofs in tomato 21 . In addition, the number of reported Dofs in many plant species ranges from 20 to 50 genes according to information from the Plant TFDB. Taken together, the numbers of Dofs among plant species appear to be regular, regardless of genome size. Consistent with previous results from other plants 16,21 , the CaDofs were divided into four major groups and seven subgroups based on a neighbor-joining phylogenetic tree. Phylogenetic analysis of the Dofs in pepper and tomato revealed particular clusters of paralogous and orthologous gene pairs within each subgroup supported by high bootstrap values ( Fig. 1 and Supplementary Fig. S3). Although the clade patterns are consistent with previous results from Arabidopsis, rice, and tomato, the number of duplicated genes within some groups of Dofs in pepper were different from those in the other three species. Group B contained the most genes in pepper, whereas group C contains the most genes in Arabidopsis and tomato. The numbers of genes in group A and subgroup D1 in pepper were the same as those in tomato. Previous studies indicated that orthologous Dofs in the same clades between monocots and dicots might be derived from common ancestors that evolved before the monocot/dicot split 16,21 . Taken together, the results suggest that there were unequal gene duplication events in each clade of the Dof family among species, leading to expansion and diversity of the Dof repertoire during the course of evolution.
Structural analysis showed that the CaDofs contained few introns, with intron numbers ranging from zero to two per gene (Fig. 3). Most of the genes had no introns, and those genes all clustered into the same groups. The lengths of the introns within the genes showed similar patterns within the same groups. Similar characteristics of intron/exon structure were observed in Arabidopsis, rice, and tomato 16,21 , suggesting evolutionary conservation of the Dofs. The chromosomal locations of the CaDofs showed that the CaDofs were located on 8 of the 12 pepper 'CM334' is resistant to TMV-P0, PepMoV, and P. capsici but susceptible to TMV-P2.The 29 CaDof genes were used to construct heat map. The 4 CaDof genes (CaDof03, 23, 24, and 33), which are absent data, were excluded. For expression pattern after virus inoculation, mean of RPKM value in triplicate data sets were used. This data was normalized by control samples (Mock_virus and Mock_P.c). Blue and red colors represent relatively low and high expression (log 2 RPKM value), respectively. P.c, Phytophthora capsici; D, days post inoculation. chromosomes, with chromosomes 7, 8, 9, and 12 containing no CaDofs. Chromosome 2 contained nine CaDofs, which was the highest number among the chromosomes. Likewise, chromosomes 7 and 12 in tomato lack Dofs, and the highest number of Dofs (nine SlDofs) in tomato is found on chromosome 2. The chromosomal locations of the orthologous gene pairs between pepper and tomato were compared based on a phylogenetic tree. That comparative analysis revealed that most of the orthologous gene pairs were syntenic between the two plant species. Those results are consistent with previous reports of conserved locations of homologous genes on chromosomes between species 40,44 . Hence, we assumed that the locations of the seven unassigned CaDofs (CaDof27, 28, 29, 30, 31, 32, and 33) might correspond to syntenic positions of orthologous SlDofs. For instance, CaDof31, which is orthologous to SlDof14, might be located on chromosome 3 in pepper.
The motif analysis of CaDofs revealed that the CaDofs shared conserved motifs in the same subgroup. For example, C1 has specific motifs 7, and 14, whereas C2 contains motif 25. We also found specific motifs in other subgroups: motif 11 and 12 in B1, and motif 3, 4, 5, 8, and 10 in D1 (Fig. 2). The distribution of motifs also showed conserved position of each gene in the same subgroup. These conserved motifs of Dof genes indicated that such structures have been preserved by evolution suggesting that these motifs might play important roles in subgroup functions.
Tandem duplication, segmental duplication, and transposition are major evolutionary mechanisms of gene-family expansion 45 . In plants, segmental duplication occurs more frequently than tandem duplication and transposition, because most plants are diploidized polyploids and retain numerous duplicated chromosomal blocks within their genomes 45 . We obtained 10 pairs of paralogs in the pepper Dof family based on chromosomal distribution, phylogenetic analysis, and sequence similarity. One pair of paralogs (CaDof23 and 24) was the result of a putative tandem duplication event. Among the remaining nine pairs, two pairs (CaDof02 and 10 and  (Table 2). Therefore, the duplication of most of the CaDofs occurred before the pepper/tomato divergence, and only one paralogous pair was duplicated recently, after the pepper/tomato split. Similar cases have been reported in tomato. The duplication of one SlDof was estimated to have occurred before the pepper/tomato split (19.1-61.1 Mya) 21 . The Ka/Ks ratio provides a measure of selective pressure. If the gene pair were under purifying selection, it would have been purged by natural selection, presumably because of deleterious effects. Conversely, if the gene pair was under positive selection, it would have been advantageous during the evolution of the two duplicates 47 . We calculated the Ka/Ks ratio of the duplicated pairs in the CaDof family. The Ka/Ks ratios of all 10 duplicated pairs in the CaDof family were < 1.0, suggesting that the duplicated pairs evolved under purifying selection. Duplicate pairs of Dofs in tomato and soybeans were also estimated to have evolved under purifying selection 20,21 . Overall, the duplication analysis suggests that the segmental expansion of the CaDofs might have taken place before the speciation event separating pepper and tomato and that the functional divergence of the CaDofs was retained after the duplications.
Transcriptional control of gene expression is an important part of plant growth, development, and response to biotic and abiotic stresses 48 . RNA-seq has been a useful and powerful tool for the analysis of gene expression profiles since next-generation sequencing has become very straightforward. The RNA-seq data showed that most of the CaDofs had distinct tissue-specific expression patterns among the five different tissues tested. The CaDof expression profiles revealed various patterns among different organs and stages (Fig. 5). CaDof17, 18, 21, 27, and 30 were expressed at relatively high levels in various tissues and at several developmental stages. Seventeen of 33 genes showed relatively high expression levels in the root. Similarly, 22 SlDofs in tomato and 22 GmDofs in soybeans were expressed at relatively high levels in the root 20,21 . Ten (CaDof08, 11, 13, 17, 18, 21, 22, 27, 29, and 30) and seven (CaDof08, 11, 13, 18, 22, 27, and 30) genes showed relatively high expression levels in the stem and leaf, respectively. The expression levels of CaDof11, 21, and 30 in the placenta and that of CaDof18 in the pericarp were higher than those in other organs, suggesting that those genes could play key roles in fruit development, such as capsaicinoid biosynthesis. Taken together, the results indicate that the CaDofs are possibly involved in biological functions during plant development.
Recently, Dof family has also been reported to be possibly involved in biotic stresses 49,50 . However, still their role associated with biotic stress were largely unknown. Additional clues and works should be needed to understand and confirm their newer roles in biotic stress. In our study, the CaDof family showed distinct expression patterns in response to different biotic stresses (Fig. 6). The transcript levels of the CaDofs showed temporal and pathogen-specific variation. CaDof04, 05, and 32 showed changes in expression levels in the resistant response to TMV-P0 and/or PepMoV and P. capsici, while CaDof07 was up-regulated only in the susceptible response to TMV-P2. In particular, CaDof10 and 11 showed relatively high expression levels in response to each of the viruses, suggesting that they could play a role in defense response against viruses. In addition, some CaDofs such as CaDof09 and 17, showed similar expression patterns and temporal variation following various pathogen inoculations, indicating a shared function in response to biotic stresses. Thus, the differentially regulated CaDofs might play a crucial role in defense against biotic stresses. Those results will provide useful information for the functional characterization of Dofs to understand biological processes in pepper and other plants.
Scientific RepoRts | 6:33332 | DOI: 10.1038/srep33332 Methods Identification and annotation of CaDofs. We obtained the consensus amino-acid sequence of the Dof protein as a seed from the Pfam (PF02701) database (http://pfam.sanger.ac.uk/). We then used that sequence to search all of the Dof proteins against the 'CM334' genome 40 using the HMMER software package (Version 3.0, http://hmmer.org/). We confirmed the collected CaDofs by a BLAST search against NCBI. We then collected the CaDofs manipulated by alignment of the 'CM334' genome sequence, the 'Zunla-1' genome sequence, the pepper ESTs, and the Dof sequences from the Plant TFDB. If the alignment result showed the same sequence among the Dofs from several databases, we chose the longest sequence as a representative. We computed the pIs and MWs of the predicted protein sequences using the ExPASy Compute pI/MW tool (http://web.expasy.org/compute_pi/). We obtained the full amino-acid sequences of the tomato Dofs, designated SlDofs, from the Plant TFDB.
Motif analysis of the CaDofs and exon/intron structural analysis. We used the MEME suite program (http://meme-suite.org/tools/meme) to identify conserved motifs among the CaDofs. The analysis settings of the MEME suite were a motif length of 6-200 amino acids and 2-120 motif sites. We set the maximum number of motifs to 25. We set the other conditions to the default values. We performed the alignment of the identified CaDof domains using ClustalW (http://www.ebi.ac.uk/Tools/clustalw/). We predicted the zinc finger structure of the Dof domain based on information from the tomato and Arabidopsis Dof domains 16,21 . We searched and confirmed the exon/intron structures of the CaDofs by aligning the CDSs with the corresponding genomic regions. We presented the exon/intron structures using the online display tool Gene Structure Display Server (http://gsds. cbi.pku.edu.cn/).

Phylogenetic tree analysis.
We performed a multiple alignment of the 33 full-length CaDofs and 34 full-length SlDofs using ClustalW. We obtained the Dof sequences for tomato from previous research 21 and the Plant TFDB. We used the alignment result to construct a phylogenetic tree using the neighbor-joining method of MEGA6 software (http://www.megasoftware.net/). We created the phylogenetic tree using the following parameters: Poisson correction, pairwise deletion, and a bootstrap value of 1,000.
Chromosomal location and calculation of duplication events. We mapped the CDSs of the Dofs to the 'CM334' chromosomes in the pepper genome database 40 using BLASTn. Information about the chromosomal locations of the CaDofs was based on the identified physical positions on the chromosomes. We used MEGA6 to make the pairwise alignments of paralogs with the removal of gaps. We estimated the Ka, Ks, and Ka/Ks using K-Estimator 6.1v 51 . We then used the Ks value of each pair to estimate the duplication time by following formula: duplication time = Ks/2λ , with the clock-like rate (λ ) = 6.1 × 10 −9 based on previous research 52 .

Plant materials and pathogen inoculation.
We used 'CM334' to analyze the responses to TMV-P0, TMV-P2, PepMoV, and P. capsici. We conducted the virus inoculations as described by Kang et al. 53 . We inoculated Nicotiana benthamiana leaves infected with TMV-P0, TMV-P2, and PepMoV mechanically with 0.1 M phosphate buffer (pH 7.0) into true leaves of 'CM334' seedlings. For virus inoculation, three independent experiments of each virus were each performed. We performed the preparation and inoculation of P. capsici as described by Yeom et al. 37 . We inoculated 'CM334' seedlings with 1.0 × 10 6 zoospores/ml P. capsici. The pepper plants inoculated with viruses and P. capsici were kept in a growth chamber at 23-25 °C under a 16 h/8 h light/ dark photoperiod.

Transcriptome analysis.
To obtain RNA-seq data from pathogen-infected pepper tissues, we harvested the leaves from various time points after pathogen inoculation to prepare the total RNA. We used RNA samples extracted from three biological replicates of TMV-P0, TMV-P2, and PepMoV inoculated tissues at various time points for library construction using a modified protocol for strand-specific library construction 54 (Illumina Inc., San Diego, USA). We sequenced the constructed libraries (insert size: 150-200 bp) using Illumina HiSeq 2000 (Illumina Inc., San Diego, USA). We first aligned the RNA-seq reads to rRNA and tRNA sequences to remove possible contamination. Low-quality reads were filtered by in-house trimming scripts. We then aligned the resulting reads to the reference pepper genome using CLC Assembly Cell (CLC bio, Aarhus, Denmark). Following alignment to each gene model, we normalized the number of mapped reads to the reads per kilobase per million mapped reads (RPKM).
To analyze the expression patterns of the CaDofs in specific tissues and at different developmental stages, we used the RNA-seq data obtained from 'CM334' 40 . We generated transcript data for five tissues: root, stem, leaf, pericarp, and placenta. The pericarp and placenta data included seven developmental stages: 6, 16, and 25 days post inoculation; mature green (MG); breaker (B); and 5 and 10 days post B. We illustrated the log 2 RPKM values of the expression data for the tissues and developmental stages using the heat-map R package (http://bioconductor.org/). To evaluate the expression patterns in response to the virus and P. capsici inoculations, we compared the normalized transcripts obtained under each biotic stress with each negative control (Supplementary Table S2). For data analysis of virus inoculation, three RNA-seq data sets corresponding three biological replicates of each virus were analyzed. Then, mean of values in triplicate datasets were used for analysis of expression pattern and construction of heat maps. We created the heat map using the heat-map R package (http://bioconductor.org/).