Gene Expression Analysis of Plum pox virus (Sharka) Susceptibility/Resistance in Apricot (Prunus armeniaca L.)

sRNA Analysis Evidenced the Involvement of Different Plant Viruses in the Activation of RNA Silencing-Related Genes and the Defensive Response Against Plum pox virus of 'GF305' Peach Grafted with 'Garrigues' Almond

Plum pox virus (PPV) causes sharka disease in Prunus trees. Peach (P. persica) trees are severely affected by PPV, and no definitive source of genetic resistance has been identified. However, previous results showed that PPV-resistant 'Garrigues' almond (P. dulcis) was able to transfer its resistance to 'GF305' peach through grafting, reducing symptoms and viral load in PPV-infected plants. A recent study tried to identify genes responsible for this effect by studying messenger RNA expression through RNA sequencing in peach and almond plants, before and after grafting and before and after PPV infection. In this work, we used the same peach and almond samples but focused the high-throughput analyses on small RNA (sRNA) expression. We studied massive sequencing data and found an interesting pattern of sRNA overexpression linked to antiviral defense genes that suggested activation of these genes followed by downregulation to basal levels. We also discovered that 'Garrigues' almond plants were infected by different plant viruses that were transferred to peach plants. The large amounts of viral sRNA found in grafted peaches indicated a strong RNA silencing antiviral response and led us to postulate that these plant viruses could be collaborating in the observed "Garrigues effect."

TRAF proteins as key regulators of plant development and stress responses

Tumor necrosis factor receptor-associated factor (TRAF) proteins are conserved in higher eukaryotes and play key roles in transducing cellular signals across different organelles. They are characterized by their C-terminal region (TRAF-C domain) containing seven to eight anti-parallel β-sheets, also known as the meprin and TRAF-C homology (MATH) domain. Over the past few decades, significant progress has been made toward understanding the diverse roles of TRAF proteins in mammals and plants. Compared to other eukaryotic species, the Arabidopsis thaliana and rice (Oryza sativa) genomes encode many more TRAF/MATH domain-containing proteins; these plant proteins cluster into five classes: TRAF/MATH-only, MATH-BPM, MATH-UBP (ubiquitin protease), Seven in absentia (SINA), and MATH-Filament and MATH-PEARLI-4 proteins, suggesting parallel evolution of TRAF proteins in plants. Increasing evidence now indicates that plant TRAF proteins form central signaling networks essential for multiple biological processes, such as vegetative and reproductive development, autophagosome formation, plant immunity, symbiosis, phytohormone signaling, and abiotic stress responses. Here, we summarize recent advances and highlight future prospects for understanding on the molecular mechanisms by which TRAF proteins act in plant development and stress responses.

Deciphering Prunus Responses to PPV Infection: A Way toward the Use of Metabolomics Approach for the Diagnostic of Sharka Disease

Sharka disease, caused by Plum pox virus (PPV), induces several changes in Prunus. In leaf tissues, the infection may cause oxidative stress and disrupt the photosynthetic process. Moreover, several defense responses can be activated after PPV infection and have been detected at the phytohormonal, transcriptomic, proteomic, and even translatome levels. As proposed in this review, some responses may be systemic and earlier to the onset of symptoms. Nevertheless, these changes are highly dependent among species, variety, sensitivity, and tissue type. In the case of fruit tissues, PPV infection can modify the ripening process, induced by an alteration of the primary metabolism, including sugars and organic acids, and secondary metabolism, including phenolic compounds. Interestingly, metabolomics is an emerging tool to better understand Prunus–PPV interactions mainly in primary and secondary metabolisms. Moreover, through untargeted metabolomics analyses, specific and early candidate biomarkers of PPV infection can be detected. Nevertheless, these candidate biomarkers need to be validated before being selected for a diagnostic or prognosis by targeted analyses. The development of a new method for early detection of PPV-infected trees would be crucial for better management of the outbreak, especially since there is no curative treatment.

Genetic dissection of soybean partial resistance to sclerotinia stem rot through genome wide association study and high throughout single nucleotide polymorphisms

Sclerotinia stem rot (SSR) is a disease of soybean [Glycine max (L.) Merr] that causes severe yield losses. We studied 185 representative soybean accessions to evaluate partial SSR resistance and sequenced these by the specific-locus amplified fragment sequencing method. In total, 22,048 single-nucleotide polymorphisms (SNPs), with minor allele frequencies (MAF) ≥5% and missing data <3%, were developed and applied to genome-wide association study of SSR responsiveness and assess linkage disequilibrium (LD) level for candidate gene selection. We identified 18 association signals related to SSR partial resistance. Among them, six overlapped the regions of previous quantitative trait loci, and twelve were novel. We identified 243 candidate genes located in the 200 kb genomic region of these peak SNPs. Based on quantitative real-time polymerase chain reaction and haplotype analysis, Glyma.03G196000 and Glyma.20G095100, encoding pentatricopeptide repeat proteins, might be important factors in the resistance response of soybean to SSR.

Gene Expression Analysis of Induced Plum pox virus (Sharka) Resistance in Peach (Prunus persica) by Almond (P. dulcis) Grafting

No natural sources of resistance to Plum pox virus (PPV, sharka disease) have been identified in peach. However, previous studies have demonstrated that grafting a “Garrigues” almond scion onto “GF305” peach rootstock seedlings heavily infected with PPV can progressively reduce disease symptoms and virus accumulation. Furthermore, grafting a “Garrigues” scion onto the “GF305” rootstock has been shown to completely prevent virus infection. This study aims to analyse the rewiring of gene expression associated with this resistance to PPV transmitted by grafting through the phloem using RNA-Seq and RT-qPCR analysis. A total of 18 candidate genes were differentially expressed after grafting “Garrigues” almond onto healthy “GF305” peach. Among the up-regulated genes, a HEN1 homolog stands out, which, together with the differential expression of RDR- and DCL2-homologs, suggests that the RNA silencing machinery is activated by PPV infection and can contribute to the resistance induced by “Garrigues” almond. Glucan endo-1,3-beta D-glucosidase could be also relevant for the “Garrigues”-induced response, since its expression is much higher in “Garrigues” than in “GF305”. We also discuss the potential relevance of the following in PPV infection and “Garrigues”-induced resistance: several pathogenesis-related proteins; no apical meristem proteins; the transcription initiation factor, TFIIB; the speckle-type POZ protein; in addition to a number of proteins involved in phytohormone signalling.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

BACKGROUND

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

RESULTS

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

CONCLUSIONS

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

Background

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

Results

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

Conclusions

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Translatome Profiling of Plum Pox Virus-Infected Leaves in European Plum Reveals Temporal and Spatial Coordination of Defense Responses in Phloem Tissues

Plum pox virus (PPV) is the causative agent of sharka, a devastating disease of stone fruits including peaches, apricots, and plums. PPV infection levels and associated disease symptoms can vary greatly, depending upon the virus strain, host species, or cultivar as well as developmental age of the infected tissues. For example, peaches often exhibit mild symptoms in leaves and fruit while European plums typically display severe chlorotic rings. Systemic virus spread into all host tissues occurs via the phloem, a process that is poorly understood in perennial plant species that undergo a period of dormancy and must annually renew phloem tissues. Currently, little is known about how phloem tissues respond to virus infection. Here, we used translating ribosome affinity purification followed by RNA sequencing to identify phloem- and nonphloem-specific gene responses to PPV infection during leaf development in European plum (Prunus domestica L.). Results showed that, during secondary leaf morphogenesis (4- and 6-week-old leaves), the phloem had a disproportionate response to PPV infection with two- to sixfold more differentially expressed genes (DEGs) in phloem than nonphloem tissues, despite similar levels of viral transcripts. In contrast, in mature 12-week-old leaves, virus transcript levels dropped significantly in phloem tissues but not in nonphloem tissues. This drop in virus transcripts correlated with an 18-fold drop in phloem-specific DEGs. Furthermore, genes associated with defense responses including RNA silencing were spatially coordinated in response to PPV accumulation and were specifically induced in phloem tissues at 4 to 6 weeks. Combined, these findings highlight the temporal and spatial dynamics of leaf tissue responses to virus infection and reveal the importance of phloem responses within a perennial host.

Molecular Plant-Plum Pox Virus Interactions

Plum pox virus, the agent that causes sharka disease, is among the most important plant viral pathogens, affecting Prunus trees across the globe. The fabric of interactions that the virus is able to establish with the plant regulates its life cycle, including RNA uncoating, translation, replication, virion assembly, and movement. In addition, plant-virus interactions are strongly conditioned by host specificities, which determine infection outcomes, including resistance. This review attempts to summarize the latest knowledge regarding Plum pox virus-host interactions, giving a comprehensive overview of their relevance for viral infection and plant survival, including the latest advances in genetic engineering of resistant species.

Resistance to Sharka in Apricot: Comparison of Phase-Reconstructed Resistant and Susceptible Haplotypes of 'Lito' Chromosome 1 and Analysis of Candidate Genes

Sharka, a common disease among most stone fruit crops, is caused by the Plum Pox Virus (PPV). Resistant genotypes have been found in apricot (Prunus armeniaca L.), one of which-the cultivar 'Lito' heterozygous for the resistance-has been used to map a major quantitative trait locus (QTL) on linkage group 1, following a pseudo-test-cross mating design with 231 individuals. In addition, 19 SNP markers were selected from among the hundreds previously developed, which allowed the region to be limited to 236 kb on chromosome 1. A 'Lito' bacterial artificial chromosome (BAC) library was produced, screened with markers of the region, and positive BAC clones were sequenced. Resistant (R) and susceptible (S) haplotypes were assembled independently. To refine the assembly, the whole genome of 'Lito' was sequenced to high coverage (98×) using PacBio technology, enabling the development of a detailed assembly of the region that was able to predict and annotate the genes in the QTL region. The selected cultivar 'Lito' allowed not only to discriminate structural variants between the two haplotypic regions but also to distinguish specific allele expression, contributing towards mining the PPVres locus. In light of these findings, genes previously indicated (i.e., MATHd genes) to have a possible role in PPV resistance were further analyzed, and new candidates were discussed. Although the results are not conclusive, the accurate and independent assembly of R and S haplotypes of 'Lito' is a valuable resource to predict and test alternative transcription and regulation mechanisms underpinning PPV resistance.

The apricot (Prunus armeniaca L.) genome elucidates Rosaceae evolution and beta-carotenoid synthesis

Apricots, scientifically known as Prunus armeniaca L, are drupes that resemble and are closely related to peaches or plums. As one of the top consumed fruits, apricots are widely grown worldwide except in Antarctica. A high-quality reference genome for apricot is still unavailable, which has become a handicap that has dramatically limited the elucidation of the associations of phenotypes with the genetic background, evolutionary diversity, and population diversity in apricot. DNA from P. armeniaca was used to generate a standard, size-selected library with an average DNA fragment size of ~20 kb. The library was run on Sequel SMRT Cells, generating a total of 16.54 Gb of PacBio subreads (N50 = 13.55 kb). The high-quality P. armeniaca reference genome presented here was assembled using long-read single-molecule sequencing at approximately 70× coverage and 171× Illumina reads (40.46 Gb), combined with a genetic map for chromosome scaffolding. The assembled genome size was 221.9 Mb, with a contig NG50 size of 1.02 Mb. Scaffolds covering 92.88% of the assembled genome were anchored on eight chromosomes. Benchmarking Universal Single-Copy Orthologs analysis showed 98.0% complete genes. We predicted 30,436 protein-coding genes, and 38.28% of the genome was predicted to be repetitive. We found 981 contracted gene families, 1324 expanded gene families and 2300 apricot-specific genes. The differentially expressed gene (DEG) analysis indicated that a change in the expression of the 9-cis-epoxycarotenoid dioxygenase (NCED) gene but not lycopene beta-cyclase (LcyB) gene results in a low β-carotenoid content in the white cultivar "Dabaixing". This complete and highly contiguous P. armeniaca reference genome will be of help for future studies of resistance to plum pox virus (PPV) and the identification and characterization of important agronomic genes and breeding strategies in apricot.

The significance of methionine cycle enzymes in plant virus infections

Both biotic and abiotic stresses cause changes in the activities of plant methionine cycle (MTC) enzymes. These changes contribute to the ability of the plant to manage stress. On the other hand, viruses utilize MTC enzymes to promote infection. Here, we review the growing but still limited knowledge of the interactions between plant viral proteins and MTC enzymes. Virus-induced changes in S-adenosyl methionine synthetase and S-adenosyl homocysteine hydrolase activities debilitate transcriptional and post-transcriptional RNA silencing and affect antiviral defense reactions connected to ethylene and polyamine biosynthesis pathways. Viral perturbations of host methionine homeostasis couple trans-sulfuration and gluthathione biosynthesis pathways to MTC functions. Large multiprotein complexes, which contain viral proteins and MTC enzymes, may represent metabolons assembled for specific viral functions or host defense responses. Proper understanding of the MTC-associated metabolic and regulatory interactions will reveal those with potential to create virus resistance in plants.

Identification of phloem-associated translatome alterations during leaf development in Prunus domestica L

Phloem plays a fundamental role in plants by transporting hormones, nutrients, proteins, RNAs, and carbohydrates essential for plant growth and development. However, the identity of the underlying phloem genes and pathways remain enigmatic especially in agriculturally important perennial crops, in part, due to the technical difficulty of phloem sampling. Here, we used two phloem-specific promoters and a translating ribosome affinity purification (TRAP) strategy to characterize the phloem translatome during leaf development at 2, 4, and 6 weeks post vernalization in plum (Prunus domestica L.). Results provide insight into the changing phloem processes that occur during leaf development. These processes included the early activation of DNA replication genes that are likely involved in phloem cell division during leaf expansion, as well as the upregulation of phloem genes associated with sink to source conversion, induction of defense processes, and signaling for reproduction. Combined these results reveal the dynamics of phloem gene expression during leaf development and establish the TRAP system as a powerful tool for studying phloem-specific functions and responses in trees.

Phytohormone Signaling of the Resistance to Plum pox virus (PPV, Sharka Disease) Induced by Almond (Prunus dulcis (Miller) Webb) Grafting to Peach (P. persica L. Batsch)

Plum pox virus (PPV, sharka) is a limiting factor for peach production, and no natural sources of resistance have been described. Recent studies, however, have demonstrated that grafting the almond cultivar "Garrigues" onto the "GF305" peach infected with Dideron-type (PPV-D) isolates progressively reduces disease symptoms and virus accumulation. Furthermore, grafting "Garrigues" onto "GF305" prior to PPV-D inoculation has been found to completely prevent virus infection, showing that resistance is constitutive and not induced by the virus. To unravel the phytohormone signaling of this mechanism, we analyzed the following phytohormones belonging to the principal hormone classes: the growth-related phytohormones cytokinin trans-zeatin (tZ) and the gibberellins GA₃ and GA₄; and the stress-related phytohormones ethylene acid precursor 1-aminocyclopropane-1-carboxylic acid (ACC), abscisic acid (ABA), salicylic acid (SA), and jasmonic acid (JA). PPV inoculation produced a significant increase in GA₃ and ABA in peach, and these imbalances were related to the presence of chlorosis symptoms. However, grafting "Garrigues" almond onto the PPV-inoculated "GF305" peach produced the opposite effect, reducing GA₃ and ABA contents in parallel to the elimination of symptoms. Our results showed the significant implication of SA in this induced resistance in peach with an additional effect on tZ and JA concentrations. This SA-induced resistance based in the decrease in symptoms seems to be different from Systemic Acquired Resistance (SAR) and Induced Systemic Resistance (ISR), which are based in other reactions producing necrosis. Further studies are necessary, however, to validate these results against PPV-D isolates in the more aggressive Marcus-type (PPV-M) isolates.

Molecular insights into Cassava brown streak virus susceptibility and resistance by profiling of the early host response

Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) are responsible for significant cassava yield losses in eastern sub-Saharan Africa. To study the possible mechanisms of plant resistance to CBSVs, we inoculated CBSV-susceptible and CBSV-resistant cassava varieties with a mixed infection of CBSVs using top-cleft grafting. Transcriptome profiling of the two cassava varieties was performed at the earliest time point of full infection (28 days after grafting) in the susceptible scions. The expression of genes encoding proteins in RNA silencing, salicylic acid pathways and callose deposition was altered in the susceptible cassava variety, but transcriptional changes were limited in the resistant variety. In total, the expression of 585 genes was altered in the resistant variety and 1292 in the susceptible variety. Transcriptional changes led to the activation of β-1,3-glucanase enzymatic activity and a reduction in callose deposition in the susceptible cassava variety. Time course analysis also showed that CBSV replication in susceptible cassava induced a strong up-regulation of RDR1, a gene previously reported to be a susceptibility factor in other potyvirus-host pathosystems. The differences in the transcriptional responses to CBSV infection indicated that susceptibility involves the restriction of callose deposition at plasmodesmata. Aniline blue staining of callose deposits also indicated that the resistant variety displays a moderate, but significant, increase in callose deposition at the plasmodesmata. Transcriptome data suggested that resistance does not involve typical antiviral defence responses (i.e. RNA silencing and salicylic acid). A meta-analysis of the current RNA-sequencing (RNA-seq) dataset and selected potyvirus-host and virus-cassava RNA-seq datasets revealed that the conservation of the host response across pathosystems is restricted to genes involved in developmental processes.

Resistance to Plum Pox Virus (PPV) in apricot (Prunus armeniaca L.) is associated with down-regulation of two MATHd genes

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

Quantitative proteome-level analysis of paulownia witches' broom disease with methyl methane sulfonate assistance reveals diverse metabolic changes during the infection and recovery processes

Paulownia witches' broom (PaWB) disease caused by phytoplasma is a fatal disease that leads to considerable economic losses. Although there are a few reports describing studies of PaWB pathogenesis, the molecular mechanisms underlying phytoplasma pathogenicity in Paulownia trees remain uncharacterized. In this study, after building a transcriptome database containing 67,177 sequences, we used isobaric tags for relative and absolute quantification (iTRAQ) to quantify and analyze the proteome-level changes among healthy P. fortunei (PF), PaWB-infected P. fortunei (PFI), and PaWB-infected P. fortunei treated with 20 mg L^-1 or 60 mg L^-1 methyl methane sulfonate (MMS) (PFI-20 and PFI-60, respectively). A total of 2,358 proteins were identified. We investigated the proteins profiles in PF vs. PFI (infected process) and PFI-20 vs. PFI-60 (recovered process), and further found that many of the MMS-response proteins mapped to "photosynthesis" and "ribosome" pathways. Based on our comparison scheme, 36 PaWB-related proteins were revealed. Among them, 32 proteins were classified into three functional groups: (1) carbohydrate and energy metabolism, (2) protein synthesis and degradation, and (3) stress resistance. We then investigated the PaWB-related proteins involved in the infected and recovered processes, and discovered that carbohydrate and energy metabolism was inhibited, and protein synthesis and degradation decreased, as the plant responded to PaWB. Our observations may be useful for characterizing the proteome-level changes that occur at different stages of PaWB disease. The data generated in this study may serve as a valuable resource for elucidating the pathogenesis of PaWB disease during phytoplasma infection and recovery stages.