Resistance to Plum Pox Virus (PPV) in apricot (Prunus armeniaca L.) is associated with down-regulation of two MATHd genes
BMC Plant Biology volume 18, Article number: 25 (2018)
Plum pox virus (PPV), causing Sharka disease, is one of the main limiting factors for Prunus production worldwide. In apricot (Prunus armeniaca L.) the major PPV resistance locus (PPVres), comprising ~ 196 kb, has been mapped to the upper part of linkage group 1. Within the PPVres, 68 genomic variants linked in coupling to PPV resistance were identified within 23 predicted transcripts according to peach genome annotation. Taking into account the predicted functions inferred from sequence homology, some members of a cluster of meprin and TRAF-C homology domain (MATHd)-containing genes were pointed as PPV resistance candidate genes.
Here, we have characterized the global apricot transcriptome response to PPV-D infection identifying six PPVres locus genes (ParP-1 to ParP-6) differentially expressed in resistant/susceptible cultivars. Two of them (ParP-3 and ParP-4), that encode MATHd proteins, appear clearly down-regulated in resistant cultivars, as confirmed by qRT-PCR. Concurrently, variant calling was performed using whole-genome sequencing data of 24 apricot cultivars (10 PPV-resistant and 14 PPV-susceptible) and 2 wild relatives (PPV-susceptible). ParP-3 and ParP-4, named as Prunus armeniaca PPVres MATHd-containing genes (ParPMC), are the only 2 genes having allelic variants linked in coupling to PPV resistance. ParPMC1 has 1 nsSNP, while ParPMC2 has 15 variants, including a 5-bp deletion within the second exon that produces a frameshift mutation. ParPMC1 and ParPMC2 are adjacent and highly homologous (87.5% identity) suggesting they are paralogs originated from a tandem duplication. Cultivars carrying the ParPMC2 resistant (mutated) allele show lack of expression in both ParPMC2 and especially ParPMC1.
Accordingly, we hypothesize that ParPMC2 is a pseudogene that mediates down-regulation of its functional paralog ParPMC1 by silencing. As a whole, results strongly support ParPMC1 and/or ParPMC2 as host susceptibility genes required for PPV infection which silencing may confer PPV resistance trait. This finding may facilitate resistance breeding by marker-assisted selection and pave the way for gene edition approaches in Prunus.
Sharka disease, caused by Plum pox virus (PPV), is currently the most important viral disease affecting Prunus species . PPV is a member of the Potyvirus genus in the Potyviridae, one of the largest families of plant viruses, and has been included in the ‘Top 10’ ranking of scientific/economically relevant plant viruses . Described for the first time infecting plums (Prunus domestica L.) in Bulgaria around 1917 , PPV spread into most temperate fruit crop-growing areas since then . The growth of PPV-resistant Prunus cultivars is pointed out as the ideal long-term solution, especially in endemic areas where fruit trees cannot be efficiently protected from Sharka infection . However, resistant sources are scarce. Germplasm screenings have just identified a handful of North American apricot (Prunus armeniaca L.) PPV resistant cultivars  currently used as donors in the breeding programs, and a few plum genotypes showing tolerance or hypersensitive response to PPV infection .
Genetic control of PPV resistance in apricot has long been a source of controversy. Differences in the phenotyping methods and the use of distinct PPV strains hampered to reach firm conclusions . However, most studies currently support the involvement of one major dominant locus (PPVres locus) located in the upper part of the apricot linkage group 1, including genetic mapping [15, 25, 26, 34, 41, 50] and genome-wide association approaches . In a previous work, we narrowed down the PPVres locus to ~ 196 kb according to the peach genome syntenic region . This study also pointed out a cluster of meprin and TRAF-C homology domain (MATHd)-containing genes as responsible for PPV resistance. Genetic evidence primarily supported the 5-bp insertion mutation disrupting one of these genes as the candidate variant causing resistance. Thus, since PPVres heterozygous genotypes confer resistance through the mutated PPVres allele, a gain-of-function or a dominant negative mutation was hypothesized . In Arabidopsis thaliana, another MATHd-only protein-encoding gene, RTM3, is one of the dominant RTM genes involved in the restriction of PPV long distance movement . In addition, the non-functionality of one or more RTM alleles is sufficient to abolish the resistance phenotype [10, 40]. On the contrary, PPV resistance in apricot is suggested to be associated with a MATHd mutated allele encoding a truncated non-functional protein . Similarly, loss-of-function of the host eukaryotic translation initiation factor 4E isoform (eIF(iso)4E) has been shown to confer PPV resistance in A. thaliana and plum [14, 56]. However, resistance in these cases is due to ‘recessive homozygosity’ while, in apricot, natural PPV resistance is present in heterozygosis .
Genetic engineering technologies have been explored to obtain PPV resistant Prunus cultivars and rootstocks overcoming breeding limitations such as incompatibility barriers and long generation periods. Much of the work has involved pathogen derived resistance using constructs of the PPV coat protein (CP) gene (reviewed by ). This strategy was successfully implemented to develop the PPV resistant ‘HoneySweet’ C5 plum by RNA silencing activation against the PPV CP coding sequence (; Patent No. US PP15154 P2). Host-derived resistance also has been exploited by silencing the host susceptibility (S) gene eIF(iso)4E to obtain PPV resistant transgenic plums . However, transformation and regeneration in Prunus still are limiting factors for gene transfer technologies highly genotype-dependent . No genome-edited Prunus species has been reported to date, but this technology will presumably have a great impact in PPV resistance breeding.
In this work we combined transcriptomic and genomic data to validate candidate genes previously identified by Zuriaga et al.  and to characterize the global transcriptome response to PPV infection in apricot. As a whole, results allowed us to support that PPV resistance relies on two down-regulated PPVres locus MATHd genes and to propose a dominant gene action via gene-silencing, both findings being relevant for future breeding in Prunus.
Transcriptome profiles differ between resistant/susceptible apricot cultivars
To analyze the apricot response to PPV infection, PPV resistant cultivars ‘Goldrich’ and ‘Stella’ and the susceptible ‘Canino’ were grafted onto susceptible GF-305 peach rootstocks, and half of the plants were inoculated by chip budding with PPV-D. Gene expression was analyzed in a total of 16 leaf RNA samples using pair-ends (PE) HiSeq2000 Illumina sequencing. More than 490 M of 101 bp raw sequence reads were obtained, averaging 30.63 M per sample (Additional file 1: Table S1). After quality trimming and adapter clipping, a total of ~ 48,885 M of high quality sequenced bases (99.75% of raw sequenced bases) belonging to 70–100 bp cleaned reads were obtained. As suggested by Haas et al. , normalized cleaned sequences were assembled by Trinity software. After refinement (see Experimental procedures section for details), 91,735 transcripts were generated and grouped into 61,096 genes, with 98,391,234 bases and a contig mean length of 1072.56 bases. Regarding other Prunus species, comparable results were obtained studying the response to PPV infection in plum , the response of reproductive tissues to frost stress in almond (Prunus dulcis (Mill.) D.A. Webb) , and the dynamics of fruit development  and the anthocyanin biosynthesis in sweet cherry (Prunus avium L.) . Prunus persica reference transcriptome (peach v.1.0, rosaceae.org), obtained from different tissues (i.e. fruits, roots, leaves, embryos and cotyledons) has 28,689 transcripts . Up to 10,894 peach transcripts are represented in this leaf RNA only-based apricot transcriptome with a length coverage over 80% (Additional file 2: Table S2). Putative orthologs were detected in peach for 34% of the total apricot assembled transcripts by using the reciprocal best Blast hit (RBH) criterion (Additional file 3: Data S1).
Plum pox virus genome sequence (NCBI Reference Sequence: NC_001445.1) was blasted (e-value >1e-07) against the apricot assembled contigs and just one contig (c34934_g0_i1) showed similarity with the virus sequence (Additional file 4: Table S3). PPV genome was almost completely assembled into this single contig, with 9741 bases length and more than 98% identity, except for a small portion of both ends and a 45 bases insertion. Reads mapping against this ‘PPV contig’ were used to check the presence of the virus (Additional file 5: Table S4). As expected, PPV was significantly present in the 3 replicates of the PPV-inoculated susceptible ‘Canino’ but almost absent in the rest. PPV symptoms were clearly observed in all inoculated GF305 rootstock indicators.
Transcript abundance estimates for each sample were obtained by aligning trimmed PE reads to the assembled transcriptome. In order to identify trends or detect putative biases in the data set, relationships between samples were checked using multidimensional scaling (MDS) plots (Additional file 6: Figure S1). As a whole, transcriptome profile variability was higher between cultivars than it was between infection conditions (I: inoculated; NI: non-inoculated) (Additional file 6: Figure S1a). PPV-inoculated ‘Goldrich’ replicate ‘Go_I_rep3’ was clearly separated from the rest of the ‘Goldrich’ samples and therefore it was eliminated for subsequent analyses to prevent background noise. Within cultivar, ‘Canino’ (Additional file 6: Figure S1b) and ‘Stella’ (Additional file 6: Figure S1d) samples appear separated according to the infection conditions but not in ‘Goldrich’ (Additional file 6: Figure S1c). Technical replicate pairs (‘Ca_NI_rep2a’/'Ca_NI_rep2b’; ‘Go_I_rep1a’/'Go_I_rep1b’; ‘Go_I_rep2a’/'Go_I_rep2b’) cluster together respectively and were considered as single biological replicates (‘Ca_NI_rep2’, ‘Go_I_rep1’, and ‘Go_I_rep2’) for subsequent analyses.
Number of differentially expressed genes (DEGs) between I and NI infection conditions differs for each cultivar: 793 in ‘Canino’ (homozygous for PPV susceptibility), 194 in ‘Stella’ (homozygous for PPV resistance) and 23 in ‘Goldrich’ (heterozygous) (Fig. 1a). ‘Canino’ and ‘Stella’ share in common 102 DEGs while 688 appear exclusively in ‘Canino’ and 92 in ‘Stella’. Two common DE cell-wall related genes exhibit opposed behaviors in both cultivars, c29444_g0 appears down-regulated in inoculated ‘Canino’ and over-expressed in inoculated ‘Stella’, while the opposite occurs with c33041_g0 (Fig. 1b). To gain insights into the biological roles of DEGs in both cultivars, a Gene Ontology (GO) categories enrichment analysis (Fisher exact test, with a FDR < 0.05) was performed using Blast2GO (Table 1). In ‘Canino’, we found GO terms for molecular functions significantly enriched in binding (iron ion, chromatin and DNA) and catalytic activities (oxidoreductase, hydrolase, lyase, transferase and transporter). For biological processes, enriched GO terms were related to abiotic stimulus response. Finally, for cellular component the enriched GO terms were related to chloroplasts (thylakoids and stroma), nucleus and nucleosome. In ‘Stella’ we only found enriched GO terms for biological processes related to the response to abiotic and endogenous stimulus, circadian rhythm and positive development regulation.
Two PPVres locus MATHd genes are down-regulated in apricot resistant cultivars
The major dominant PPVres locus in apricot comprises ~ 196 kb according to the peach genomic syntenic region . Nineteen assembled apricot genes were RBH and other ten showed highly similarity with some peach annotated genes within the PPVres locus (Additional file 7: Table S5). None of them was DE within cultivar comparing PPV infection conditions but 6 (ParP-1 to ParP-6) appear significantly DE between cultivars (Fig. 2, Additional file 8: Table S6). ParP-3 and -4 were highly down-regulated in ‘Stella’ (absolute logFC values between 5.19–7.4), while ParP-1, − 2, − 5 and − 6 were up-regulated (absolute logFC values between 0.78–1.45). ParP-3, ParP-4 and ParP-5 are homologous to 3 peach MATHd genes (ppa022254m, ppb0221 95 m and ppa008951m, respectively) previously pointed as PPV resistance candidate genes . RNAseq results were further confirmed by qRT-PCR analyses that revealed similar gene-expression patterns (Fig. 3). For instance, no statistically significant differences (p < 0.05) were detected within cultivars under different infection conditions. ParP-3 and ParP-4 were clearly down-regulated in ‘Stella’ with regard to ‘Canino’, while ‘Goldrich’ occupied an intermediate position between them (Fig. 3a-b). ParP-3 showed the highest relative expression differences, ~ 300-fold increase in non-inoculated ‘Canino’ compared with ‘Stella’ (Fig. 3a, Additional file 9: Table S7). On the contrary, ParP-5 showed no significant differences neither between nor within samples (Fig. 3c). Correlation between these gene-expression patterns and the PPV-resistance phenotype was confirmed by testing two additional resistant cultivars (‘Harlayne’ and ‘Orange Red’, both PPVres heterozygous) and four susceptible cultivars (‘Currot’, ‘Ginesta’, ‘Katy’ and ‘Mitger’) (Fig. 3 ; Additional file 9: Table S7). ParP-3 and ParP-4 were found clearly down-regulated in all resistant cultivars (Fig. 3a-b), and ParP-3 again showed more striking differences ranging between ~ 300 to ~ 4267-fold higher gene-expression in susceptible cultivars with regard to ‘Stella’ (Additional file 9: Table S7). Moreover, ParP-5 did not show consistent differences between susceptible and resistant cultivars (Fig. 3c).
ParP-3 and ParP-4 accumulate genomic variants linked to PPV resistance
Whole genome sequences (WGS) of 24 apricot cultivars (10 PPV-resistant and 14 PPV-susceptible) and two apricot relatives (PPV-susceptible), as well as reads derived from six BAC clones corresponding to the ‘Goldrich’ PPVres locus R-haplotype, were aligned against the peach genome (peach v1.0; http://www.rosaceae.org) (Additional file 10: Table S8). Variant calling to detect SNPs or small insertions and deletions (INDELs) within the PPVres locus identified 2424 to 3928 putative SNPs and 425 to 711 putative INDELs per genotype (Additional file 10: Table S8). Afterwards, three filters were sequentially applied to discriminate SNPs/INDELs associated with PPV resistance from all 7459 detected variants: i) variants should be linked in coupling with PPV-resistance, as confirmed by their presence in ‘Goldrich’ R-haplotype BACs, being heterozygous in diploid ‘Goldrich’ WGS; ii) they should be homozygous in ‘Stella’ and heterozygous in the rest of the resistant cultivars; iii) they had to be absent in all 16 susceptible cultivars (Fig. 4a, Additional file 11: Table S9). A total of 44 SNP/INDELs fulfilled these three conditions. Twenty-eight of these variants were found in intergenic regions, being 14 and 11 in the putative promoter regions of ParP-3 and ParP-4, respectively (Fig. 4b). In addition, 1 filtered variant was present within ParP-3 (nsSNP: Ile109Leu) and 15 within ParP-4 (7 in intronic regions, 5 sSNPs, 2 nsSNPs: Glu194Lys and Thr266Ala, and 1 5-bp deletion) (Fig. 4b). Interestingly, this latter 5-bp deletion is located in the second exon and produces a frameshift mutation that creates a premature stop codon (Fig. 4c). The 5-bp deletion could be consistently screened on agarose gel electrophoresis by allele-specific PCR in susceptible and resistant cultivars providing a useful tool for marker-assisted selection (MAS) into apricot breeding programs (Fig. 4c-d). Following the gene naming guidelines of the Genome Database for Rosaceae (GDR) , ParP-3 and -4 were respectively named ParPMC1 and ParPMC2 (Prunus armeniaca PPVres locus MATH-domain containing (PMC) genes).
Five apricot genes were identified as putative orthologs or highly similar to the 9 peach PPVres locus MATHd genes according to reciprocal Blast results (Additional file 7: Table S5). Maximum likelihood based phylogeny of these 14 genes revealed 3 highly supported sub-clusters (Fig. 5). ParPMC1, ParPMC2 and ParP-5, and their putative peach orthologs, ppa022254m, ppb0221 95 m and ppa008951m, grouped together in the same sub-cluster. ParPMC1 and ParPMC2 showed the shortest genetic distance among all apricot pairs (0,148) having 87,5% of sequence identity (Additional file 12: Table S10).
Apricot response against PPV infection
Differences in transcriptome profiles were more striking between apricot cultivars than between viral infection conditions. Different factors related or not with PPV infection may explain this. For instance, ‘Goldrich’ and ‘Stella’ are North American apricot cultivars adapted to cold-growing conditions while the Spanish ‘Canino’ is mainly grown through the temperate Mediterranean Basin . PPV inoculation procedure entails a cold treatment to break dormancy . Accordingly, distinct cultivar-dependent responses to this treatment could be expected as reflected by their global gene-expression profiles. As expected, deeper differences in transcriptome profiles between PPV inoculated and non-inoculated tissues were observed in the susceptible materials. In accordance with previous works, this expression-pattern might be due to changes suffered by cells experiencing pathogenic stress. Taking into account the broad differences among previous studies, similar biological processes and molecular functions were found to be affected in this work. In apricot, Rubio et al.  identified DEGs involved in biological processes associated with responses to different stimuli. Wang et al.  observed up-regulation of genes involved in defense, cellular transport, development, protein synthesis and binding functions in peach infected leaves. Studies on the response to PPV infection in Arabidopsis leaves identified altered genes belonging to major groups of metabolism, transcription/splicing/RNA processing proteins, defense and development/storage proteins . Enriched GO terms for cellular components are in agreement with papers describing the effect of PPV on the photosynthetic processes producing physical and biochemical changes in the chloroplasts [11, 22, 44, 49]. As a whole, transcriptomic data provided in this work have sharpened our knowledge on the response to PPV at gene-expression level and might be helpful to search for additional genes involved in PPV resistance.
PPV resistance in apricot is mediated by MATHd gene/s
The PPVres locus syntenic region in peach contains 29 genes (one of them encoding for 3 alternative transcripts) as annotated by the International Peach Genome Initiative. Gene-expression analysis in this work showed that none of the corresponding apricot genes was DE between infected/non-infected tissues, suggesting that PPV-D presence does not modulate their expression. However, six DEGs were identified between the susceptible ‘Canino’ and the resistant ‘Stella’. Two of them, ParPMC1 and ParPMC2, were significantly DE emerging as the best candidates from expression data. Previously, using genomic data of three PPV-resistant and four PPV-susceptible cultivars, a total of 68 variants linked in coupling with PPV resistance were identified in 23 apricot genes located at the ~ 196 kb PPVres locus . In this work, a much deeper variant calling analysis based on 26 apricot WGS, allowed us to confirm 44 variants matching with the imposed criteria by the genetics of PPV resistance in apricot [50, 53]. ParPMC1 and ParPMC2 genes contain 1 and 15 of these variants, respectively, while another 14 and 11 are located in their putative promoter regions. Remarkably, both genes are homologous to members of the peach MATHd genes cluster previously suggested as the most likely PPV resistance candidate genes . ParPMC1 has only one nsSNP (Ile109Leu) within the first MATH domain but ParPMC2 accumulates 15 variants including a 5 nt deletion that results in a premature stop codon. ParPMC1 and ParPMC2 are highly similar (87.5% identity), close to ParP-5, and all 3 occupy adjacent positions in the MATHd genes cluster. Therefore, we hypothesize that they are paralogs originated by tandem duplications from a single ancestral gene. However, their expression patterns are different, ParPMC1 and ParPMC2 are down-regulated in PPV resistant cultivars but ParP-5 expression seems to be similar in both PPV resistant and susceptible cultivars. ParPMC1 remains basically unaltered while the ParPMC2 resistant allele has accumulated deleterious mutations probably as a consequence of a pseudogenization process . It could be speculated that this, in turn, affects the expression of ParPMC1 and ParPMC2 susceptible alleles but not ParP-5 expression due to sequence differences. Thus, it may also be suggested that ParP-5 function is not redundant to ParPMC1 and ParPMC2. Altogether, these results support that ParPMC1 and/or ParPMC2 genes are required for PPV susceptibility in apricot. Other loci contributing to PPV resistance have been proposed in studies using PPV-Rec  and PPV-M strains , but consistent candidate genes remain to be identified. Moreover, Decroocq et al. [12, 13] observed that some apricot genotypes carrying the ParPMC2 resistant allele are susceptible to PPV, mainly to the M strain. These reports suggest that ParPMC2 resistant allele does not suffice to confer PPV resistance being necessary additional genes. This point cannot be discarded on the basis of the evidence provided in this work. However, all available studies (including these latter) undoubtedly show a tight linkage between ParPMC2 resistant allele and PPV-D resistance. Moreover, misfits may be due to other genes but also to the PPV strain (different genes may underlie PPV-M resistance).
Down-regulation of MATHd gene/s expression causes PPV resistance in apricot
Down-regulation of ParPMC1 and ParPMC2 genes expression is the differential factor between PPV resistant and susceptible apricot cultivars regarding the PPVres locus. Mutations accumulated in the promoter regions of the PPV resistant alleles may be directly affecting their expression, but half-gene dosage can not account for the observed differences between resistant and susceptible cultivars. Nevertheless, ParPMC1 and ParPMC2 gene expression may also be prevented by gene-silencing. For instance, nonsense-mediated mRNA decay (NMD) is a conserved mechanism that targets aberrant mRNAs carrying premature termination codons for destruction, preventing the accumulation of potentially harmful truncated proteins . According to these authors, a typical target has a termination codon positioned more than 50–55 nt upstream of the last exon-exon junction or has a long 3′ untranslated region. Interestingly, ParPMC2 has a 5-bp loss-of-function deletion in the second exon, 73 nt upstream of an exon–exon junction, being a potential target for NMD. Pseudogenes are suggested to potentially regulate their protein-coding cousins , as exemplified by the silencing of a gene required for sexual seed formation in apomictic Paspalum simplex associated with its homolog pseudogene . Moreover, overlaps between RNA silencing and NMD pathways have been reported in plants. Christie et al.,  pointed out that when aberrant mRNA formation is frequent (i.e. due to nonsense mutations) there is a threshold-dependent induction of RNA silencing. According to gene-expression data, apricot cultivars carrying only one ParPMC2 resistant (mutated) allele showed lack of expression of both ParPMC2 and ParPMC1 genes, but especially ParPMC1. This scenario is consistent with a dominant mutation in ParPMC2 conferring PPV resistance by silencing functional homologs such as the susceptible (non-mutated) ParPMC2 allele and/or ParPMC1. Therefore, ParPMC1 and ParPMC2 can be considered host susceptibility (recessive) genes which silencing may confer PPV resistance trait.
Together with the eIF(iso)4E–like factor , ParPMC1 and ParPMC2 are, to our knowledge, the first factors required for PPV susceptibility identified in Prunus, the natural hosts of PPV. Functional analyses are currently in progress to elucidate the exact role of ParPMC1 and ParPMC2 genes in PPV resistance. Meanwhile, silencing of these and other recessive genes, such as the eIF(iso)4E–like factor, by genome editing technologies seems a promising strategy to exploit host-derived resistance in Prunus pyramiding durable resistance to PPV. Furthermore, closer analysis of other DEGs in the susceptible cultivar ‘Canino’ may allow the identification of additional host S-genes involved in viral infection.
Plant material and PPV inoculation
The PPV resistant cultivars ‘Goldrich’ and ‘Stella’ and the susceptible ‘Canino’ were used in the RNAseq experiment. These cultivars are maintained as part of the germplasm collection at IVIA (Valencia, Spain). Each genotype was grafted onto PPV susceptible ‘GF305’ peach seedlings growing in pots under controlled greenhouse conditions as described by Moustafa et al. . Half of the plants were inoculated by chip budding on the GF305 rootstock using the PPV Dideron strain 3.3 RB . Four weeks after grafting, plants were subjected to an artificial period of dormancy in darkness at 5 °C for 8 weeks. Subsequently, the plants were transferred again to the greenhouse. After 8 weeks, young and fully developed leaves, randomly distributed along the sprout, were harvested and directly frozen in liquid nitrogen in the greenhouse before RNA extraction. The presence of virus was scored by visual inspection of symptoms in susceptible ‘GF-305’ leaves and by detecting PPV transcripts through RNAseq analysis. Three replicates for each condition (inoculated/non-inoculated) were sampled from ‘Goldrich’ and ‘Canino’ and two from ‘Stella’.
Additionally, 2 resistant (‘Harlayne’ and ‘Orange Red’) and 5 susceptible (‘Canino’, ‘Currot’, ‘Ginesta’, ‘Katy’ and ‘Mitger’) cultivars were employed for gene expression analyses. All of them also maintained as part of the germplasm collection at IVIA (Valencia, Spain).
Total RNA isolation and sequencing
Total RNA was isolated from 100 mg of powdered leaves with the RNeasy Plant Mini Kit (Qiagen, Valencia, CA, USA), adding 1% (w:v) polyvinylpyrrolidone (PVP-40) to the kit extraction buffer before use, followed by a DNase treatment with the RNase-Free DNase Set (Qiagen, Valencia, CA, USA). RNA quality and quantity were checked by agarose gel electrophoresis and Nanodrop ND-1000 spectrophotometer (Nanodrop Technologies, Wilminton, DE, USA). Library construction and sequencing were performed by the Beijing Genomics Institute (BGI Tech Solutions (Hong Kong) Co., LTD). RNA quality and quantity were rechecked before sequencing by BGI using both Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and Nanodrop. ‘Short-insert’ libraries were sequenced using an Illumina HiSeq2000 instrument for PE sequencing of 101-bp. Some samples were sequenced twice to obtain the amount of clean data needed and treated as technical replicates in subsequent analyses. Cleaned paired-end sequence dataset was deposited in the NCBI Short Read Archive (SRA) under the accession numbers SRR5591366 to SRR5591375, associated with the BioProject PRJNA387702.
De novo transcriptome assembly and quality control
FastQC v.0.10.1 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) software was used to assess the quality of raw and clean read sets. Reads were quality trimmed using FASTX-toolkit (http://hannonlab.cshl.edu/fastx_toolkit) with a minimum quality score of 25 and a minimum length of 40. Adaptor sequences were trimmed using the ‘trim_blast_short’ script available as part of seq_crumbs (http://bioinf.comav.upv.es/seq_crumbs/). High quality reads of all samples were combined in order to perform the de novo assembly using the Trinity v.20140413p1 software  (see Additional file 13 : Methods S1 for details).
In order to check the quality of the assembly, obtained transcripts were blasted against the GDR peach transcripts database v.1.0 (ftp://ftp.bioinfo.wsu.edu/species/Prunus_persica/Prunus_persica-genome.v1.0/genes/) for similarity searches (top 25 hits at e-value >1e-07). Results allowed us to estimate the number and length of the genes recovered using the ‘analyze_blastPlus_topHit_coverage.pl’ script from the Trinity software. Assembled transcriptome was deposited in the NCBI Transcriptome Shotgun Assembly Sequence (TSA) Database associated to the BioProject PRJNA387702.
De novo transcriptome annotation
Gene ontology (GO) annotation was performed using the Blast2GO software  and the Blastx results against the NCBI non-redundant protein (nr) database (http://www.ncbi.nlm.nih.gov/protein). Blastx Reciprocal Best Hits (RBH) analysis was performed to obtain a set of putative orthologs between apricot and peach using peptides contained in the GDR peach v1.0 database. Plum pox virus genome sequence (NCBI Reference Sequence: NC_001445.1) was blasted (e-value >1e-07) against the apricot assembled transcripts in order to identify virus sequences. The ngs_backbone software  was employed for annotations from Blast results in all cases, except for the GO term annotation. Blast analyses were conducted using Picasso supercomputer from the Supercomputing and Bioinnovation Center at the University of Málaga (http://www.scbi.uma.es). Rest of the analyses were made using the Bioinformatics Department server (http://bioinf.comav.upv.es/) of the Instituto de Conservación y Mejora de la Agrodiversidad Valenciana (COMAV) at the Polytechnic University of Valencia.
Differential expression analysis
Cleaned RNA-seq reads were aligned to the assembled transcriptome using Bowtie  through the Trinity software . Transcript quantification was performed with RSEM  and the edgeR package  was used to call differentially expressed genes (DEGs). Samples considered as technical replicates were analyzed both independently and combined as a single sample. False discovery rate (FDR) ≤0.05 was used to determine the threshold of the p-value in multiple tests. Relations between samples were observed using multidimensional scaling (MDS) plots, where distance between each pair of samples can be interpreted as the leading log-fold change between the samples for the genes that best distinguish that pair of samples. By default, leading fold-change is defined as the root-mean-square of the largest 500 log2-fold changes between that pair of samples. In our case, MDS plots were obtained using those 500 genes and also using all genes, with no significant differences between both approaches. GO enrichment analysis of DEGs was performed using Blast2GO software with a cutoff value of FDR ≤ 0.05. Venn diagrams were obtained using the tool disposable in http://bioinfogp.cnb.csic.es/tools/venny/index.html. Heat-map was performed using a custom R script.
Four PPV resistant (‘Goldrich’, ‘Stella’, ‘Harlayne’ and ‘Orange Red’) and 5 PPV susceptible (‘Canino’, ‘Currot’, ‘Ginesta’, ‘Katy’ and ‘Mitger’) cultivars were analyzed by qRT-PCR. Total RNA (500 ng) was reverse transcribed with the PrimeScript RT reagent kit using an Oligo-d(T) primer (Takara Bio, Otsu, Japan) in a total volume of 10 μl. Two microliters of 10X diluted first-strand cDNA were used for PCR reactions in a final volume of 20 μl. qRT-PCR was performed on a StepOnePlus Real-Time PCR System (Life Technologies, Carlsbad, CA, USA), using SYBR premix Ex Taq (Tli RNaseH plus) (Takara Bio). Primer pairs are listed in Additional file 14: Table S11. Cycling protocol consisted of 10 min at 95 °C, followed by 40 cycles of 15 s at 95 °C for denaturation and 1 min at 60 °C for annealing and extension. PCR reaction specificity was assessed by the presence of a single peak in the dissociation curve after amplification and through size estimation of the amplified products by agarose electrophoresis. Normalized gene expression levels were measured by the relative standard curve procedure using the geometric mean of two reference genes, Actin and Sand-like . Results were the average of 1–3 independent biological replicates with 3 technical replicates each one. Comparisons of multiple samples were evaluated by the non-parametric Kruskal-Wallis test, with a confidence level of 95%, using the Statgraphics Centurion XVII v. 17.2.00 software (Statpoint Technologies, Warrenton, VA, USA). Significantly different samples were labelled with different letters.
WGS mapping, variant calling and filtering
WGS of 10 PPV resistant and 14 PPV susceptible cultivars and 2 PPV susceptible apricot relatives were used in this study (Additional file 10: Table S8). Ten of these WGS, and the 454 sequenced BAC clones belonging to the ‘Goldrich’ PPVres locus R-haplotype, were available from our previous works [39, 58]. Other 16 WGS were downloaded from the SRA repository (https://www.ncbi.nlm.nih.gov/sra). All raw reads were processed using the ‘run_trimmomatic_qual_trimming.pl’ script from the Trinity software. After removing the low-quality regions as well as vector and adaptor contaminants, cleaned reads were aligned to the peach genome v1.0 (ftp://ftp.bioinfo.wsu.edu/species/Prunus_persica/Prunus_persica-genome.v1.0/) using Bowtie2 v2.2.4 software . Variant calling to detect SNPs or small INDELs was performed using HaplotypeCaller tool from the Genome Analysis Toolkit (GATK) v3.5–0-g36282e4 software  setting the minimum phred-scaled quality score of 10 in emission confidence and of 30 in calling confidence. Following the GATK Best Practices (https://software.broadinstitute.org/gatk/best-practices/), variant discovery analysis was performed as cohort of samples. In order to discriminate variants linked in coupling with PPV resistance from all the detected variants, 3 filters were sequentially applied using homemade python scripts: (i) variants must be present in the PPV resistant ‘Goldrich’ haplotype (454 BAC contig sequences) and be heterozygous in the diploid ‘Goldrich’ WGS; (ii) variants should be homozygous in ‘Stella’ but heterozygous in the other PPV resistant cultivars; (iii) variants must not be present in any PPV susceptible cultivar. Peach genome annotation available from the GDR was used as a reference to identify polymorphisms associated with the predicted transcripts.
ParPMC2 cDNA amplification and sequencing
Complete coding DNA sequences (CDS) of ParPMC2 susceptible (S) and resistant (R) alleles were PCR-amplified with the primers pair cDNA_EcoRI_F/cDNA_BamHI_R (Additional file 14: Table S11) using ‘Goldrich’ leaf cDNA as template. This cDNA was synthesized from 500 ng of total RNA using the SuperScript III First-Strand Synthesis System kit (Invitrogen). Cycling protocol consisted of 2 min at 95 °C; 10 cycles of 30 s at 95 °C, 30 s at 49 °C (+ 0.5 °C every cycle) and 1 min at 60 °C; 25 cycles of 30 s at 95 °C, 30 s at 57 °C and 1 min (+ 10 s every cycle) at 72 °C; and finally 72° for 10 min. PCRs were performed in a final volume of 25 μL containing 10 × buffer, 1.8 mM MgCl2, 0.2 mM of each dNTP, 400 mM of each primer, 1 U of FastStart Taq DNA Polymerase (Roche) and 2 μL of cDNA template. PCR products were electrophoresed in 1% (w/v) agarose gels. PCR fragments were cloned into pGEM-T vector (Promega) according to the manufacturer’s instruction. DNA sequencing was performed in an ABI Prism 3130XL genetic analyzer, following manufacturer instructions for the BigDye terminator v.3.1 cycle sequencing kit (Applied Biosystems, Foster City, CA), using T7 and SP6 universal primers. Forward and reverse sequences were assembled and edited with the Staden package v.1.6-r (http://staden.sourceforge.net/). Pairwise alignment was made using the ClustalW function in BioEdit version 7.0.9 . Complete CDS of ParPMC2 resistant (R) and susceptible (S) alleles are deposited under GenBank accession numbers MF346726 and MF346727, respectively.
ParPMC2 allele-specific PCR assay
ParPMC2 R- and S-alleles were specifically PCR-amplified using two forward specific primers (ParP-4_F_alleleS or ParP-4_F_alleleR), differing at the 3′ end, the reverse primer ParP-4_R (Fig. 4, Additional file 14: Table S11). PCRs were performed in a final volume of 20 μL containing 1 × DreamTaq buffer, 0.2 mM of each dNTP, 5 μM of each primer, 1 U of DreamTaq DNA polymerase (Thermo Fisher) and 100 ng of DNA. Cycling conditions were as follows: an initial denaturing of 95 °C for 5 min; 35 cycles of 95 °C for 30 s, 55 °C for 45 s and 72 °C for 45 s; and a final extension of 72 °C for 10 min. PCR products were electrophoresed in 1% (w/v) agarose gels.
CDS of the 9 MATHd peach genes clustered within the PPVres locus  were downloaded from the GDR database. Multiple sequence alignment using peach and apricot MATHd genes was performed using the ClustalW function in BioEdit version 7.0.9 . Poorly aligned positions and divergent regions of the alignment were eliminated using Gblocks v.0.91b . Model of nucleotide substitution comparison was performed by using the Akaike information criterion. The best-fitting evolutionary model (The Tamura 3-parameter model (T92) + G) was implemented in the Maximum Likelihood phylogenetic analysis using 500 bootstrap replications. Evolutionary divergences between sequences were estimated using the same evolutionary model and all codon positions, but removing ambiguous positions for each sequence pair. Pairwise identity matrix was also obtained. Evolutionary analyses were conducted in MEGA6 .
Basic local alignment search tool
Complementary deoxyribonucleic acid
Coding DNA sequences
Diferentially Expressed Gene
False discovery rate
Genome Database for Rosaceae
Molecular Assisted Selection
Meprin and TRAF-C homology domain
Multidimensional scaling plot
NCBI non-redundant protein sequence
Prunus armeniaca PPVres
Prunus armeniaca PPVres MATHd-containing
Plum Pox virus
PPV resistance locus
Quantitative real-time polymerase chain reaction
Reciprocal best Blast Hit;
Single Nucleotide Polymorphisms
Short Read Archive
Transcriptome Shotgun Assembly Sequence Database
Whole genome sequences
Alkio M, Jonas U, Declercq M, Van Nocker S, Knoche M. Transcriptional dynamics of the developing sweet cherry (Prunus avium L.) fruit: sequencing, annotation and expression profiling of exocarp-associated genes. Hortic Res. 2014;1:11.
Asensio M. El virus de la sharka (plum pox virus): caracterización, diagnóstico y detección mediante anticuerpos monoclonales específicos. PhD dissertation. Spain: University of Valencia; 1996. p. 193.
Atanasoff D. Plum Pox. A New Virus Disease. In: Yearbook Faculty Agricultural University 1932/1933, Sofia, Bulgaria. 1933;11:49–69.
Babu M, Griffiths JS, Huang T, Wang A. Altered gene expression changes in Arabidopsis leaf tissues and protoplasts in response to Plum pox virus infection. BMC Genomics. 2008;9:325.
Blanca JM, Pascual L, Ziarsolo P, Nuez F, Cañizares J. ngs_backbone: a pipeline for read cleaning, mapping and SNP calling using next generation sequence. BMC Genomics. 2011;12:285.
Capote N, Cambra M, Llácer G, Petter F, Platts L, Roy A, Smith I. Current status of Plum pox virus and sharka disease worldwide. EPPO Bull. 2006;36:205–18.
Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17:540–52.
Christie M, Brosnan CA, Rothnagel JA, Carroll BJ. RNA decay and RNA silencing in plants: competition or collaboration? Front Plant Sci. 2011;2:99.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.
Cosson P, Sofer L, Le QH, Leger V, Schurdi-Levraud V, Whitham SA, Yamamoto ML, Gopalan S, Le Gall O, Candresse T, Carrington JC, Revers F. RTM3, which controls long-distance movement of Potyviruses, is a member of a new plant gene family encoding a Meprin and TRAF homology domain-containing protein. Plant Physiol. 2010;154:222–32.
Dardick C. Comparative expression profiling of Nicotiana benthamiana leaves systemically infected with three fruit tree viruses. Mol Plant-Microbe Interact. 2007;20:1004–17.
Decroocq S, Chague A, Lambert P, Roch G, Audergon J, Geuna F, Chiozzotto R, Bassi D, Dondini L, Tartarini S, Salava J, Krska B, Palmisano F, Karayiannis I, Decroocq V. Selecting with markers linked to the PPVres major QTL is not sufficient to predict resistance to Plum Pox virus (PPV) in apricot. Tree Genet Genomes. 2014;10:1161–70.
Decroocq S, Cornille A, Tricon D, Babayeva S, Chague A, Eyquard JP, Karychev R, Dolgikh S, Kostritsyna T, Liu S, Liu W, Geng W, Liao K, Asma BM, Akparov Z, Giraud T, Decroocq V. New insights into the history of domesticated and wild apricots and its contribution to Plum pox virus resistance. Mol Ecol. 2016;25:4712–29.
Decroocq V, Sicard O, Alamillo JM, Lansac M, Eyquard JP, García JA, Candresse T, Le Gall O, Revers F. Multiple resistance traits control Plum pox virus infection in Arabidopsis thaliana. Mol Plant-Microbe Interact. 2006;19:541–9.
Dondini L, Lain O, Vendramin V, Rizzo M, Vivoli D, Adami M, Guidarelli M, Gaiotti F, Palmisano F, Bazzoni A, Boscia D, Geuna F, Tartarini S, Negri P, Castellano M, Savino V, Bassi D, Testolin R. Identification of QTL for resistance to plum pox virus strains M and D in Lito and Harcot apricot cultivars. Mol Breed. 2011;27:289–99.
Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J. Preservation of duplicate genes by complementary, degenerative mutations. Genetics. 1999;151:1531–45.
García JA, Cambra M. Plum pox virus and sharka disease. Plant Viruses. 2007;1:69–79.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2011;29:644–52.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M, MacManes MD, Ott M, Orvis J, Pochet N, Strozzi F, Weeks N, Westerman R, William T, Dewey CN, Henschel R, LeDuc RD, Friedman N, Regev A. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protocols. 2013;8:1494–512.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucl Acids Symp Ser. 1999;41:95–8.
Hartmann W, Neumüller M. Breeding for resistance: breeding for plum pox virus resistant plums (Prunus domestica L.) in Germany. EPPO Bull. 2006;36:332–6.
Hernández JA, Díaz-Vivancos P, Rubio M, Olmos E, Ros-Barceló A, Martínez-Gómez P. Long-term plum pox virus infection produces an oxidative stress in a susceptible apricot, Prunus armeniaca, cultivar but not in a resistant cultivar. Physiol Plantarum. 2006;126:140–52.
Ilardi V, Tavazza M. Biotechnological strategies and tools for plum pox virus resistance: trans-, intra-, cis-genesis, and beyond. Front Plant Sci. 2015;6:379.
Jung S, Bassett C, Bielenberg DG, Cheng C, Dardick C, Main D, Meisel L, Slovin J, Troggio M, Schaffer RJ. A standard nomenclature for gene designation in the Rosaceae. Tree Genet Genomes. 2015;11:108.
Lalli DA, Abbott AG, Zhebentyayeva TN, Badenes ML, Damsteegt V, Polák J, Krška B, Salava J. A genetic linkage map for an apricot (Prunus armeniaca L.) BC1 population mapping plum pox virus resistance. Tree genet. Genomes. 2008;4:481–93.
Lambert P, Dicenta F, Rubio M, Audergon JM. QTL analysis of resistance to sharka disease in the apricot (Prunus armeniaca L.) ‘Polonais’ x ‘stark early Orange’ F1 progeny. Tree genet. Genomes. 2007;3:299–309.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Meth. 2012;9:357–9.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Layne REC, Biley CH, Hough LF. Apricots. In: Janick J, Moore JN, editors. Fruit breeding: tree and tropical fruits, vol. 2. NY: John Wiley and Sons; 1996. p. 79–11.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
Llácer G, Badenes ML, Romero C. Problems in the determination of inheritance of Plum Pox virus resistance in apricot. Acta Hort. 2006;781:263–7.
Lloret A, Conejero A, Leida C, Petri C, Gil-Muñoz F, Burgos L, Badenes ML, Ríos G. Dual regulation of water retention and cell growth by a stress-associated protein (SAP) gene in Prunus. Sci Rep. 2017;23:332.
Lykke-Andersen S, Jensen TH. Nonsense-mediated mRNA decay: an intricate machinery that shapes transcriptomes. Nat Rev Mol Cell Biol. 2015;16:665–77.
Marandel G, Pascal T, Candresse T, Decroocq V. Quantitative resistance to Plum pox virus in Prunus davidiana P1908 linked to components of the eukaryotic translation initiation complex. Plant Pathol. 2009;58:425–35.
Mariette S. Wong Jun tai F, Roch G, Barre a, Chague a, Decroocq S, Groppi a, Laizet Y, Lambert P, Tricon D, Nikolski M, Audergon J, Abbott AG, Decroocq V. Genome-wide association links candidate genes to resistance to Plum Pox Virus in apricot (Prunus armeniaca). New Phytol. 2016;209:773–84.
Martínez-Gómez P, Dicenta F, Audergon JM. Behaviour of apricot (Prunus armeniaca L.) cultivars in the presence of Sharka (Plum pox potyvirus): a review. Agronomie. 2000;20:407–22.
Mousavi S, Alisoltani A, Shiran B, Fallahi H, Ebrahimie E, Imani A, Houshmand S. De novo transcriptome assembly and comparative analysis of differentially expressed genes in Prunus Dulcis mill. In response to freezing stress. PLoS One. 2014;14:e104541.
Moustafa TA, Badenes ML, Martínez-Calvo J, Llácer G. Determination of resistance to sharka (plum pox) virus in apricot. Sci Hort. 2001;91:59–70.
Muñoz-Sanz JV, Zuriaga E, Badenes ML, Romero C. A disulfide bond A-like oxidoreductase is a strong candidate gene for self-incompatibility in apricot (Prunus armeniaca) pollen. J Exp Bot. 2017;68:5069–78.
Pagny G, Paulstephenraj PS, Poque S, Sicard O, Cosson P, Eyquard J, Caballero M, Chague A, Gourdon G, Negrel L, Candresse T, Mariette S, Decroocq V. Family-based linkage and association mapping reveals novel genes affecting Plum pox virus infection in Arabidopsis thaliana. New Phytol. 2012;196:873–86.
Pilarova P, Marandel G, Decroocq V, Salava J, Krka B, Abbott AG. Quantitative trait analysis of resistance to Plum pox virus in the apricot F1 progeny ‘Harlayne’ x ‘Vestar’. Tree Genet Genomes. 2010;6:467–75.
Pink RC, Wicks K, Caley DP, Punch EK, Jacobs L, Carter DRF. Pseudogenes: pseudo-functional or key regulators in health and disease? RNA. 2011;17:792–8.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Rodamilans B, San León D, Mühlberger L, Candresse T, Neumüller M, Oliveros JC, García JA. Transcriptomic analysis of Prunus domestica undergoing hypersensitive response to Plum Pox Virus infection. PLoS One. 2014;9:e100477.
Rubio M, Ballester AR, Olivares PM. Castro dM, Dicenta F, Martínez-Gómez P. Gene expression analysis of Plum pox virus (Sharka) susceptibility/resistance in apricot (Prunus armeniaca L.). PLoS One. 2015;10:e014e4670.
Scholthof KG, Adkins S, Czosnek H, Palukaitis P, Jacquot E, Hohn T, Hohn B, Saunders K, Candresse T, Ahlquist P, Hemenway C, Foster GD. Top 10 plant viruses in molecular plant pathology. Mol Plant Pathol. 2011;12:938–54.
Scorza R, Callahan A, Dardick C, Ravelonandro M, Polak J, Malinowski T, Zagrai I, Cambra M, Kamenova I. Genetic engineering of plum pox virus resistance: ‘HoneySweet’ plum: from concept to product. Plant Cell Tissue Organ Cult. 2013;115:1–12.
Siena LA, Ortiz JPA, Calderini O, Paolocci F, Cáceres ME, Kaushal P, Grisan S, Pessino SC, Pupilli F. An apomixis-linked ORC3-like pseudogene is associated with silencing of its functional homolog in apomictic Paspalum simplex. J Exp Bot. 2016;67:1965–78.
Sochor J, Babula P, Adam V, Krska B, Kizek R. Sharka: the past, the present and the future. Viruses. 2012;4:2853–901.
Soriano JM, Vera-Ruiz E, Vilanova S, Martínez-Calvo J, Llácer G, Badenes ML, Romero C. Identification and mapping of a locus conferring Plum Pox Virus resistance in two apricot-improved linkage maps. Tree Genet Genomes. 2008;4:391–402.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evolution. 2013;30:2725–9.
Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J, Banks E, Garimella KV, Altshuler D, Gabriel S, DePristo MA. From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013;43:11.10.1–33.
Vera Ruiz EM, Soriano JM, Romero C, Zhebentyayeva T, Terol J, Zuriaga E, Llácer G, Abbott AG, Badenes ML. Narrowing down the apricot Plum pox virus resistance locus and comparative analysis with the peach genome syntenic region. Mol Plant Path. 2011;12:535–47.
Verde I, Abbott AG, Scalabrin S, Jung S, Shu S, Marroni F, Zhebentyayeva T, Dettori MT, Grimwood J, Cattonaro F, Zuccolo A, Rossini L, Jenkins J, Vendramin E, Meisel LA, Decroocq V, Sosinski B, Prochnik S, Mitros T, Policriti A, Cipriani G, Dondini L, Ficklin S, Goodstein DM, Xuan P, Fabbro CD, Aramini V, Copetti D, Gonzalez S, Horner DS, Falchi R, Lucas S, Mica E, Maldonado J, Lazzari B, Bielenberg D, Pirona R, Miculan M, Barakat A, Testolin R, Stella A, Tartarini S, Tonutti P, Arus P, Orellana A, Wells C, Main D, Vizzotto G, Silva H, Salamini F, Schmutz J, Morgante M, Rokhsar DS. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nat Genet. 2013;45:487–94.
Wang A, Chapman P, Chen L, Stobbs LW, Brown DCW, Brandle JE. A comparative survey, by expressed sequence tag analysis, of genes expressed in peach leaves infected with Plum pox virus (PPV) and free from PPV. Can J Plant Pathol 2005;27:410-419.
Wang X, Kohalmi SE, Svircev A, Wang A, Sanfaçon H, Tian L. Silencing of the host factor eIF(iso)4E gene confers Plum Pox Virus resistance in plum. PLoS One. 2013;8:e50627.
Wei H, Chen X, Zong X, Shu H, Gao D, Liu Q. Comparative Transcriptome analysis of genes involved in Anthocyanin biosynthesis in the red and yellow fruits of sweet cherry (Prunus avium L.). PLoS One. 2015;10:e0121164.
Zuriaga E, Soriano JM, Zhebentyayeva T, Romero C, Dardick C, Cañizares J, Badenes ML. Genomic analysis reveals MATH gene(s) as candidate(s) for plum pox virus (PPV) resistance in apricot (Prunus armeniaca L.). Mol. Plant Path. 2013;14:663–77.
We thank University of Málaga (Spain) for data processing through the Picasso supercomputer. We acknowledge Dr. Gabino Rios and Alba Lloret for advice on qPCR, Dr. Chris Dardick, Dr. Tetyana Zhebentyayeva and Dr. Albert Abbot for providing some WGS data and Dr. Joaquin Cañizares for his helpful comments on the manuscript.
This research was supported by the Spanish Ministry of Economy, Industry and Competitiveness (Research Project RTA2013–00026-C03–01).
Availability of data and materials
Cleaned paired-end sequence dataset was deposited in the NCBI Short Read Archive (SRA) under the accession numbers SRR5591366 to SRR5591375, associated with the BioProject accession number PRJNA387702. Sanger sequences were deposited into the GenBank database with the accession numbers MF346726 and MF346727.
EZ designed and conducted the experiments, analyzed the data and drafted the manuscript. CR participated in the analyses of the data and revised the manuscript. JMB contributed to RNAseq computational analyses and provided technical expertise in bioinformatics. MLB contributed materials, reagents and analysis tools, supervised the data analysis and revised the manuscript. All authors read and approved the final manuscript
Ethics approval and consent to participate
Consent for publication
The authors declare that they do not have competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of RNA-seq data. (PDF 73 kb)
Distribution of percent length coverage of the assembled apricot transcripts against peach annotated transcripts (Peach v.1.0). (PDF 49 kb)
Data S1. Putative peach orthologs of the apricot assembled transcriptome. (TXT 2064 kb)
Blastn analysis of the Plum pox virus genome sequence (NCBI Reference Sequence: NC_001445.1) against the apricot assembled transcripts. (PDF 56 kb)
PPV abundance estimation based on RNA-seq sequences mapped against the apricot PPV assembled contig (c34934_g0_i1). (PDF 72 kb)
Multidimensional scaling (MDS) plots of RNA-seq expression profiles showing sample and replicate relationships. Distance between each pair of samples is the leading log-fold change between them, defined as the root-mean-square of the largest 500 log2-fold changes between that pair of samples. (JPEG 1101 kb)
Similarity analysis between apricot and peach PPVres locus genes. (XLS 17 kb)
Apricot PPVres locus DEGs identified by RNA-seq data. (PDF 107 kb)
qRT-PCR analysis of PPVres locus MATHd genes showing differential expression according to RNA-seq data. (PDF 102 kb)
Summary of genome sequences used for variant calling analysis. (PDF 90 kb)
Filtered variants putatively associated with Plum pox virus (PPV) resistance. (XLS 121 kb)
Estimates of evolutionary divergence among the apricot and peach MATHd genes clustered in the PPVres locus. (PDF 87 kb)
Methods S1. De Novo Transcriptome Assembly. (PDF 76 kb)
Primers used in this study for qRT-PCR and PCR-genotyping. (PDF 46 kb)
About this article
Cite this article
Zuriaga, E., Romero, C., Blanca, J.M. et al. Resistance to Plum Pox Virus (PPV) in apricot (Prunus armeniaca L.) is associated with down-regulation of two MATHd genes. BMC Plant Biol 18, 25 (2018). https://doi.org/10.1186/s12870-018-1237-1