Yanmin Zhu
Guanliang Li
Jugpreet Singh3†
Awais Khan
Gennaro Fazio
Rui Xia- 1Tree Fruit Research Laboratory, USDA-ARS, Wenatchee, WA, United States
- 2College of Horticulture, South China Agricultural University, Guangzhou, China
- 3Plant Pathology and Plant-Microbe Biology Section, Cornell University, Geneva, NY, United States
- 4Plant Genetic Resources Unit, USDA-ARS, Geneva, NY, United States
Apple replant disease (ARD), incited by a pathogen complex including Pythium ultimum, causes stunted growth or death of newly planted trees at replant sites. Development and deployment of resistant or tolerant rootstocks offers a cost-effective, ecologically friendly, and durable approach for ARD management. Maximized exploitation of natural resistance requires integrated efforts to identify key regulatory mechanisms underlying resistance traits in apple. In this study, miRNA profiling and degradome sequencing identified major miRNA pathways and candidate genes using six apple rootstock genotypes with contrasting phenotypes to P. ultimum infection. The comprehensive RNA-seq dataset offered an expansive view of post-transcriptional regulation of apple root defense activation in response to infection from P. ultimum. Several pairs of miRNA families and their corresponding targets were identified for their roles in defense response in apple roots, including miR397-laccase, miR398-superoxide dismutase, miR10986-polyphenol oxidase, miR482-resistance genes, and miR160-auxin response factor. Of these families, the genotype-specific expression patterns of miR397 indicated its fundamental role in developing defense response patterns to P. ultimum infection. Combined with other identified copper proteins, the importance of cellular fortification, such as lignification of root tissues by the action of laccase, may critically contribute to genotype-specific resistance traits. Our findings suggest that quick and enhanced lignification of apple roots may significantly impede pathogen penetration and minimize the disruption of effective defense activation in roots of resistant genotypes. The identified target miRNA species and target genes consist of a valuable resource for subsequent functional analysis of their roles during interaction between apple roots and P. ultimum.
Introduction
Apple replant disease (ARD) is defined by the stunted growth or possible death of newly planted trees in replant sites, where continuous cultivation of apples or closely related species leads to pathogen inoculum accumulation. ARD form a soilborne pathogen complex that includes necrotrophic soilborne oomycetes (Phytophthora and Pythium) and fungi (Ilyonectria and Rhizoctonia) (Jaffee et al., 1982; Mazzola, 1998), has been a serious threat for the establishment of economically viable orchards. Pythium ultimum is a major component of this pathogen complex in orchard soils worldwide (Mazzola, 1997; Tewoldemedhin et al., 2011). Management of ARD depends almost exclusively on pre-plant fumigation of orchard soils to eradicate ARD pathogens, which also inadvertently eliminates other beneficial soil microbiota. Use of these chemical fumigants is under increasing regulatory restriction due to environmental and human health concerns. Maximized exploitation of host resistance through investigation of the molecular regulation controlling defense responses of apple roots, and thereafter development and deployment of resistant apple rootstocks can offer a cost-effective, environment friendly and durable strategy for ARD management (Zhu et al., 2014; Zhu and Saltzgiver, 2020).
To uncover the molecular mechanisms regulating apple roots resistance responses toward P. ultimum infection, two previous transcriptome analyses have provided the first panoramic description of genome-wide transcriptional networks during defense activation (Shin et al., 2016b; Zhu et al., 2017, 2019). The peak transcriptomic changes in apple roots in response to P. ultimum infection were shown to be at 48 hpi (hour post inoculation) when about 2% of the protein-encoding genes were differentially regulated (Shin et al., 2016b). A recent comparative transcriptome analysis identified a large number of differentially expressed genes between a P. ultimum resistant G.935® and a susceptible B.9 rootstock genotype (Zhu et al., 2019). The candidate genes, which may crucially contribute to resistance traits, include those with annotated function of pathogen perception, hormone signaling, secondary metabolite production, and resistance proteins (Zhu et al., 2017, 2019). Overall, it appears that an earlier and stronger molecular-level defense activation occur in the roots of resistant genotypes as they are challenged by P. ultimum. In contrast, deficient detoxification, weakened effector-triggered immunity (ETI), or existence of a susceptibility (S) gene may contribute to a disrupted defense activation process in the roots of the susceptible genotype. Parallel to these transcriptome analyses, a systematic phenotyping effort has recently identified a panel of apple rootstock germplasm with contrasting and reliable resistance traits against P. ultimum infection (Shin et al., 2016a; Zhu et al., 2018a,b). Besides the expression patterns of the protein coding genes, elucidating post-transcriptional control of genotype-specific defense responses may facilitate the identification of key genes underpinning apple roots resistance traits.
Plant endogenous small RNAs (sRNAs) are a class of short (∼21–25 nucleotides) non-coding single-stranded RNAs, which are ubiquitous in most cellular processes including plant-microbe interactions (Bartel, 2004; He and Hannon, 2004; Eckardt, 2012; Borges and Martienssen, 2015; Weiberg and Jin, 2015). sRNAs can be divided into two major categories based on their distinct biogenesis pathways, i.e., microRNAs (miRNAs) and small interfering RNAs (siRNAs) (Axtell, 2013). sRNAs are generated by DICER or DICER-like (DCL) endoribonucleases, and the mature sRNAs are loaded into Argonaute (AGO) proteins to induce silencing of genes with complementary sequences (Rogers and Chen, 2013; Borges and Martienssen, 2015), which is generally referred to as RNA interference (RNAi). By base-pairing to the target mRNAs of protein coding genes, sRNAs modulate gene expression through either translational repression or post-transcriptional silencing. Accumulating evidence indicates that plant miRNAs play a pivotal role in controlling plant immune responses to pathogen infection through numerous cellular processes (Katiyar-Agarwal and Jin, 2010; Fei et al., 2013, 2016; Bologna and Voinnet, 2014; Chen et al., 2018). For example, the NBS-LRR-miR482/2118 regulatory network has been demonstrated to have a role in fine-tuning the timing and intensity of defense responses in tomato-Phytophthora infestans pathosystem (Jiang et al., 2018; Hong et al., 2019). On the other hand, pathogens can hijack the host RNA silencing system to reduce plant defense responses (Weiberg and Jin, 2015; de Vries et al., 2018b, 2019). Profiling miRNAs during interaction between apple roots and necrotrophic oomycete pathogen P. ultimum, as well as identifying their corresponding target genes by degradome sequencing (German et al., 2008), should provide critical insights into regulation of resistance in apple roots to this pathogen.
The objective of this study was to identify key miRNA families and their specific target genes in apple roots during interaction with P. ultimum using high-throughput RNA sequencing analysis. Six apple rootstock genotypes showing contrasting phenotypes against P. ultimum infection were included in this study to identify genotype-specific patterns in miRNA expression and the identities of miRNA-targeted genes.
Materials and Methods
Apple Rootstock Genotypes and Plant Maintenance
Based on the results from a systematic phenotyping study (Zhu et al., 2018b), six apple rootstock genotypes, i.e., three resistant (or R group: R-#58, R-#161, and R-#164) and three susceptible genotypes (or S group: S-#106, S-#115, and S-#132), were identified and included in this study. These genotypes originated from an “Ottawa 3” × “Robusta 5” (O3R5) apple rootstock F1 population developed in the mid-1970s. Individual plants for each apple rootstock genotype and mock-inoculation controls were produced using synchronized plant tissue culture or micro-propagation procedures, as described previously (Zhu et al., 2018a). Shoot propagation and root induction/elongation required 6 weeks and 4 weeks, respectively. Individual plants with enough root system were transferred to pots containing autoclaved artificial soil consisting of construction sand, vermiculite and perlite in 1:1:1 ratio for a 1-week “in-soil acclimation” before pathogen infection assays and root tissue collection. For tissue culture generated root system, the in-soil acclimation is a critical step to allow further differentiation of root tissues and to fully express the inherent resistance traits. To minimize the effects of transplanting shock on plants, especially in roots from culture medium to soils in pots, a transparent 7′ vented humidity dome (Greenhouse Megastore, Danville, IL, United States) was used to cover the flat tray holding the pots to retain humidity. The temperature in the growth room was approximately 22 ± 1°C at night and 25 ± 1°C during daytime with 12 h light/dark photoperiod regime.
Inoculation of Apple Roots With Pythium ultimum and Root Tissue Collection
The Pythium ultimum isolate (#1062) used in this study was originally isolated from the roots of “Gala”/M.26 apple grown at Moxee, WA, United States. Inoculum of P. ultimum was prepared as previously described (Shin et al., 2016b). Briefly, the inoculum of P. ultimum was prepared by cultivation in potato-carrot broth (20 g of carrots and 20 g of peeled potatoes in 1 L of water boiled for 30 min) with two drops of wheat germ oil added per L of medium. The P. ultimum cultures were grown in potato-carrot broth in 9-cm petri dishes at 22°C for 4–6 weeks. Oospores and mycelium from the resultant hyphae mat were collected and ground in 0.5% methyl cellulose solution using a household electric blender for 30 s. The oospores and hyphal fragments were resuspended in 0.5% methyl cellulose to give a final concentration of approximately 2,000 oospores per milliliter (mL). The inoculation of seedlings with P. ultimum was performed by dipping the root system in the inoculum solution for 5 s. Inoculated plants were immediately transplanted into pots with autoclaved soils and watered thoroughly. Control plants were mock inoculated with 0.5% methyl cellulose solution and maintained similarly to the P. ultimum inoculated plants. All plants were maintained in an environmental growth room as specified above. Root tissues were collected at 48 hpi (hour post inoculation) by excising the roots after excavating plant from pot and quick rinsing off attached soil medium. Three independent biological replicates were set for each sample (genotype/treatment), and each replicate included the pooled root tissues of five plants. Tissues were store in −80°C freezer until RNA isolation.
Plant Survival Rate and Microscopic Root Necrotic Patterns
A detailed description of the resistance phenotypes for included six genotypes has been reported previously (Zhu et al., 2018b). Briefly, inoculated plants were allowed to grow for 4 weeks in autoclaved potting mix. Plant survival rate for each genotype was recorded at 3, 7, 10, 14, and 28 dpi (days post inoculation), although overall plant survival rate for most of the tested genotypes was stabilized after 7 dpi. Infection assays were repeated at least three times. For microscopic observation, plants were carefully excavated from the soil to minimize mechanical damage to the roots from 2 to 7 dpi. Residual soil along the root branches was gently removed under running tap water. Roots for both mock-inoculated control and P. ultimum inoculated plants were kept separately in 100-mL beakers filled with water until microscopic examination within 2 h. Individual root branches were carefully separated from each other, and a glass slide was used to hold the roots in petri dish filled with autoclaved water. A minimum of six plants per genotype from control and P. ultimum treated were examined with the assistance of a dissecting microscope (Olympus SXZ12). Images of healthy or necrotic root tissues were obtained using a DP73 digital camera installed on Olympus SXZ12 and the associated software suite of Celsense (Olympus, Center Valley, PA, United States). Digital images were only processed for resizing, cropping and adjusting overall brightness using a publicly available software FastStone Image Viewer 5.5.1
RNA Isolation, Sequencing, Target Genes Identification and phasiRNA Analysis
Total RNA was isolated as previously described (Zhou et al., 2018). Thirty-six libraries were constructed to represent six genotypes, two treatments (mock-inoculation and P. ultimum inoculation), and three independent biological replicates for each sample (genotype/treatment). Each biological replicate included the pooled root tissues from five plants. Frozen root tissue samples were ground to a fine powder in liquid nitrogen, and RNA quantity was determined using a NanoDrop spectrophotometer (ND-1000, Thermo Fisher Scientific). The quality of total RNA used for RNA-seq passed the required standard of RNA integrity number (RIN) x ≥ 8, which was evaluated using an Agilent 2100 Bioanalyzer. Total RNAs (∼2 μg for each sample) were gel purified and fragments within the 18–30 nucleotide length range were selected to construct a library. Thirty-six sRNA libraries were prepared and sequenced at the Center for Genome Research and Biocomputing in Washington State University using an Illumina HiSeq2500TM (Illumina Inc., San Diego, CA, United States). Clean sequence reads were obtained after removal of the low-quality tags, 3′ adapter null, insert null, 5′ adapter contaminants, poly A reads and reads shorter than 18 nt. The apple miRNAs were annotated as described previously (Xia et al., 2012, 2015). The total number of reads in a given library that perfectly matched the Malus × domestica genome was used for the normalization of read abundance. Each miRNA was normalized to 10 million reads. Malus × domestica genome sequences (Malus × domestica HFTH1 Whole Genome v1.0) were downloaded from the GDR.2 The microRNA sequence data was deposited in SRA (Sequence Read Archive) at the NCBI website under the accession number SRP295189.3 Statistical analysis of measured expression values between genotypes or treatments was performed using Mann–Whitney U test.
To identify which specific apple genes are targeted by miRNA-directed degradation after P. ultimum infection, three degradome libraries were sequenced and resulted sequence tags were annotated and quantified (German et al., 2008). First library (L1) was constructed using pooled total RNAs from roots of all mock-inoculation controls, second library (L2) from pooled total RNAs from infected root tissues of all three resistant genotypes, and third library (L3) from pooled total RNAs from root tissues of all three susceptible genotypes. Identification of miRNA target genes was conducted by analyzing three degradome libraries with CleaveLand (version 4) pipeline for degradome analysis (Addo-Quaye et al., 2009). Identification of PHAS loci is based on P-value calculation, which was developed and modified in previous studies (De Paoli et al., 2009; Xia et al., 2013). PHASMinerCLI from an in-house software “sRNAminer” was used to predict PHAS loci and TargetSoPipe from sRNAminer was used to identify the trigger miRNA of PHAS loci. The degradome p-adjusted value < 0.05 is adopted in assigning the identity of target genes.
The expression patterns of three laccase genes were examined using RT-qPCR as previously described (Zhu et al., 2019). The same set of total RNA samples for small RNA-seq were used for RT-qPCR analysis. The total RNA was treated with DNase I (Qiagen, Valencia, CA, United States) and then purified with RNeasy cleanup columns (Qiagen, Valencia, CA, United States). Two micrograms of DNase-treated RNA were used to synthesize first-strand cDNA using SuperScript II reverse transcriptase (Invitrogen, Grand Island, NY, United States) and poly dT (Operon, Huntsville, AL, United States) as the primer. The expression values of target genes were normalized to that of a previously validated internal reference gene (MDP0000095375) specific for gene expression analysis in apple roots 2–ΔΔCT method (the comparative Ct method) (Livak and Schmittgen, 2001). The gene sequences were retrieved from GDR.4 Forward and reverse primers were designed using web-based Primer3Plus5 and IDT OligoAnalyzer.6 Where possible, an optimum annealing temperature of 60°C, a GC content of 40–60%, an amplicon length of 150–180 bp, and a primer length of 20 bp were applied. The primer sequences of three laccase genes and a reference gene were listed in Table 1.
Table 1. Primer sequences of laccase and reference genes used for qRT-PCR analysis to study their expression patterns in apple roots in response to Pythium ultimum infection.
Genomic Location of miR397a and Variant Analysis
A de novo genome of “Robusta 5” parent was assembled using long-read PacBio sequences to identify the genomic location of miR397a precursors. The raw sequence quality cleaning and genome assembly was performed as per Canu assembly pipeline parameters with default settings (Koren et al., 2017). The miR397 sequences were aligned to the assembled “Robusta 5” genome to identify corresponding matching sequences using “Blastn” in the National Center of Biotechnology Information (NCBI) Blast software (Johnson et al., 2008). A 20 kb sequence around miR397a genomic locus on matching “Robusta 5” contigs was extracted using a custom script. Coding regions across these sequences were annotated by comparing them with the Golden Delicious double haploid (GDDH13 v1.1) genome (Daccord et al., 2017). In addition, a nucleotide blast search was conducted using online NCBI Blast tool (Johnson et al., 2008) to further annotate the miR397 matching sequences of “Robusta 5.” To identify nucleotide variants across the miR397a matching regions, sequence reads from “Ottawa 3” parent, and six P. ultimum susceptible and resistant genotypes were used for SNP (Single Nucleotide Polymorphism) variant calling as described earlier (Singh and Khan, 2019). First, sequence reads from individual genotypes were quality cleaned using the Trimmomatic software (Bolger et al., 2014) with parameters: leading:20, trailing:20, slidingwindow:4:15, avgqual:20, minlen:25. Read sequences below threshold quality of 20 were also removed. Remaining high-quality reads were aligned against the 20 kb miR397 matching sequences of “Robusta 5” as a reference using burrows-wheeler aligner (BWA) with default parameters (Li and Durbin, 2009). The alignment files were used as input to call SNP variants using Genome Analysis Toolkit (GATK version 3.8.0) (McKenna et al., 2010). SNP variant files for each sample were recorded in genotype variant call format (gVCF) with the HaplotypeCaller plugin in GATK. The sample-specific gVCF files were combined to obtain total SNP variants across all samples using the Genotype GVCF plugin in GATK. A final set of SNP variants supported by minimum 10 read sequences were retained using VCFtools software (Danecek et al., 2011).
Results
Resistance Phenotypes and Overall Experimental Design
Six apple rootstock genotypes included in this study were previously identified to have contrasting resistance traits based on repeated root infection assays with P. ultimum, as shown in Figure 1A (Zhu et al., 2018a, b). In addition to remarkable differences of plant survivability, microscopic observation revealed necrosis progression patterns and the intensity of pathogen hyphae growth along infected root that were distinguishable between resistant and susceptible apple rootstock genotypes. Compared to the white and intact tissue of mock-inoculated apple roots (Figure 1B), the sweeping spread of necrotic tissues and profuse growth of pathogen hyphae along the infected roots was frequently associated with the root system of inoculated susceptible genotypes (Figure 1C). In contrast, within the infected root system of the resistant genotypes, the defined lines separating healthy and necrotic sections (Figure 1D, red arrows) were commonly observed, suggesting an effective impediment of pathogen progression. For genotype-specific analysis of miRNA expression profiles and identification of miRNA target genes, three resistant (or R group: R-#58, R-#161, and R-#164) and three susceptible genotypes (or S groups: S-#106, S-#115, and S-#132) were directly compared at the critical stage of defense activation at 48 hpi (hour post inoculation) (Figure 1E).
Figure 1. Resistance phenotypes of apple rootstock in response to infection from Pythium ultimum and the experimental design for miRNA profiling and target gene identification. (A) Plant survival rates between a resistant (top, R-#58) and a susceptible (bottom, S-#106) genotype which were simultaneously inoculated with the same preparation of pathogen inoculum. Plants in the left column were the mock-inoculated control for each genotype. (B) Representing image of a mock-inoculated apple root system at 48 hpi; (C) Representing image of apple root branches of a susceptible rootstock genotype (S-#132) inoculated with P. ultimum at 48 hpi, showing profuse hyphae growth and widespread tissue necrosis and collapses. (D) A representing image of apple root branches for a susceptible rootstock genotype (R-#161) inoculated with P. ultimum at 48 hpi, showing abrupt necrosis arrest, with distinct lines separating healthy and necrotic sections along infected roots. Bars represent 300 μm. (E) Experimental design for microRNA profiling analysis and degradome sequencing. Root tissues at 48 hpi from three resistant (R-#58, -#161, and -#164) and susceptible genotypes (S-#106, -#115, and -#132) were used for small RNA library construction and subsequent high-throughput sequencing. Normalized read counts were compared within and between genotypes. Three vertical lines indicate three biological replicates for each sample (specific genotype/treatment combination). Three colored oval shapes represent three pooled samples for degradome libraries and sequencing analysis, where L1 stands for combined mock-inoculated root tissues (from both resistant and susceptible genotypes), L2 for combined root tissues consisting of infected resistant genotypes, and L3 for combined root tissues consisting of infected susceptible genotypes. Red arrows indicate the defined lines separating healthy and necrotic sections.
Overview of miRNA Sequencing Statistics in Apple Roots in Response to P. ultimum Infection
A total of ∼2.8 billion reads were generated by Illumina HiSeq2500 from 36 sequenced small RNA libraries. An average of 47% of clean reads among the thirty-six libraries were mapped to the apple genome (Zhang et al., 2019)7 according to an established protocol (Xia et al., 2012; Supplementary Table 1). The relative similarity of total read numbers from each sample indicated that neither sampling procedure, library construction, nor sequencing processes created significant bias or errors (Figure 2A). A similar distribution pattern of the size of identified sRNA from thirty-six libraries was observed, with 24 nt (nucleotide) sRNAs being the most abundant, followed by 23, 21, and 22 nt (Figure 2B). A total of 233 (including 50 novel) candidate miRNA species, belonging to 48 known and 39 novel miRNA families were mapped to the apple genome (Zhang et al., 2019; Figure 2C). The complementary miRNA∗ sequences for each novel candidate miRNA were also detected, even though most were present at lower copy levels than their corresponding mature miRNAs. The family size ranged from 1 to 22 species among identified miRNA families. Unsurprisingly, the expression levels, as expressed by reads per 10 million, showed vast differences from single digits to 106 of normalized read counts among miRNA species (Supplementary Table 1).
Figure 2. Overall statistics of identified miRNAs in apple roots. (A) The average numbers of sequence reads from three replicates for each sample (genotype/treatment). (B) Size distribution of small RNAs based on number of unique sequence reads. (C) Size of family of major miRNA in apple roots.
Expression Patterns of miRNA Families in Response to P. ultimum Infection
Among 43 major known miRNA (miR) families, most demonstrated downregulated expression in response to P. ultimum infection for all six genotypes, with rare exceptions, while only five miRNA families were upregulated in both resistant and susceptible genotypes, with a few more families specifically in infected root tissue of resistant genotypes (Table 2). Three upregulated miRNA families, miR160, miR167, and miR393 directly participate in auxin signaling, which may represent one of the significant transitions of cellular pathways from normal growth to defense responses in response to P. ultimum infection. Seven miRNA families (Table 2) showed differential regulation patterns between resistant and susceptible genotype groups, and five of them were upregulated in resistant genotypes; only the miR3627 family showed the opposite. Majority of the 43 known miRNA families exhibited low to medium levels of expression across the six genotypes, as judged arbitrarily by the normalized read counts at the level of 2–4 digits (Table 2). The remaining six families were highly expressed at the level of 105 to 106 normalized read counts.
Table 2. Genotype-specificity of miRNA regulation in apple roots in response to Pythium ultimum infection.
miRNA Families Regulating Apple Root Immune Responses to P. ultimum Infection
Several miRNA families appeared to be directly linked to plant immunity based on their target genes. Although these miRNA families all showed downregulated expression in response to P. ultimum infection, subtle variations in the level of downregulation existed between resistant and susceptible groups. Laccase-targeting miR397 showed the consistent and substantial levels of reduced expression. For example, compared to values of mock-inoculation controls, infection from P. ultimum resulted in significant reduction of detected read counts for miR397a across all genotypes (Figure 3A). More interestingly, even between the mock-inoculated controls, the resistant genotypes (except R#58) demonstrated strikingly lower levels of miR397a expression compared to those of susceptible genotypes. This observation suggests that a high level of laccase activity (from lower miRNA level) may also function as a part of the pre-formed defense system even before pathogen challenge. The targets of miR398 family include the genes encoding superoxide dismutase (SOD). Together with laccase and other groups, SODs belong to copper proteins, and their roles in early immune responses have been well elucidated (Pilon, 2017; Hong et al., 2019). The expression of miR398 in apple roots were uniformly downregulated upon pathogen inoculation (Figure 3B). However, miR398 had a slightly higher expression level in mock-inoculated tissues of resistant genotypes and a slightly stronger downregulation in infected root tissues of the susceptible genotypes. The regulatory mechanisms of miR398-SOD and their contribution to resistance traits of apple roots deserve further investigation.
Figure 3. Genotype-specific expression patterns for representing members in selected miRNA families potentially participating in apple roots resistance responses. (A) miR397a; (B) miR398b; (C) miR482d; (D) miR2118a. The numbers on Y axis represent the average values of the normalized read counts from three biological replicates. Genotype/treatment are aligned along X axis; with CK representing mock-inoculation control and “Pu” representing root samples inoculated using P. ultimum. Error bars are the standard deviation (SD) of count values within three biological replicates. Error bars represent SE of three independent biological replicates. Asterisks denote statistically significant differences in a one-tailed Mann–Whitney U test (*P < 0.05).
Members of the miR482/miR2118 super family are well-documented for their role in targeting specific NBS-LRR encoding genes (Fei et al., 2016; Islam et al., 2018). The entire miR482/miR2118 super family showed uniform downregulation in response to P. ultimum infection (Figures 3C,D and Supplementary Table 1), but a closer examination, both miR482a/b as the examples, indicated a subtle variation at the levels of downregulation between resistant and susceptible genotypes. A more substantial (>1x) downregulation was observed for all resistant genotypes, but in comparison, only one out of three susceptible genotypes showed a comparable level in response to pathogen infection. The minor variation, observed at miRNA level, could lead to amplified effect, through phasiRNA generation, on genotype-specific resistance trait toward P. ultimum infection.
miRNA Families Targeting Plant Hormone Signaling and TFs Involved in Defense Activation
Three miRNA families (miR160, miR167, and miR393) that target auxin signaling showed mostly upregulated expression upon P. ultimum infection (Table 3). Some members of miR160 family demonstrated differential expression patterns between R and S genotypes. Compared to the consistent upregulation of miR160b in all three susceptible genotypes, its downregulation was observed for two out of three resistant genotypes (R-#58 and R-#164) due to pathogen infection. miR172a, which targets AP2 and therefore potentially affects ethylene signaling, also demonstrated differential expression between resistant and susceptible genotypes, though its overall expression level was low. While miR172a was downregulated in the roots of all three susceptible genotypes, two out of three resistant genotypes (R-#58 and R-#161) showed upregulated expression. Five miRNA families, i.e., miR396, miR319, miR159, miR858, and miR164, are known to target various TF families including WRKY, TCP, MYB, and NAC that are implicated in JA biosynthesis and secondary metabolism under biotic or abiotic stress conditions (Table 3; Curaba et al., 2014; Liu and Axtell, 2015; Zhang, 2015). However, the overall comparability of their regulation patterns between R and S genotype groups suggested that their post-transcriptional regulations are unlikely the critical factor differentiating resistance and susceptibility, though they are the essential components of defense responses in apple roots. Notably, the targets of these miRNA families, such as those for TFs and hormone signaling, were the commonly identified genes from previous transcriptome analyses (Shin et al., 2016b; Zhu et al., 2019).
Table 3. Genotype-specific expression patterns for miRNA families regulating hormone signaling and TFs in secondary metabolisms in apple roots in response to Pythium ultimum infection.
Degradome Sequencing Identified Key Target Genes Due to P. ultimum Infection
Degradome sequencing identified a relatively wide, but generally expected, range of categories of target genes specific to this pathosystem between apple roots and P. ultimum (Table 4 and Supplementary Table 2). Most of the identified target genes fell into the typical spectrum of reported gene targets for a known miRNA family, although novel targets for several miRNA families were also identified. For example, “betaine aldehyde dehydrogenase” and “cation/calcium exchanger” were identified as the targets of miR156, in addition to its primary known targets of SPLs for this miRNA family. Similarly, “glucan endo-1,3-beta-glucosidase” was among the targets beyond the commonly predicated targets of “homeobox-leucine zipper protein” encoding genes for miR166, and “probable pectate lyase 22” as the target for miR395 (Table 3 and Supplementary Table 2). More practically, this experimental approach allowed the pinpointing of a small number of specific genes from a large gene family, which are specific to this pathosystem. For example, only five NBS-LRR encoding genes out of potentially a few hundred in the apple genome were identified as targets for miR482 from this pathosystem. On the other hand, only two ARF genes were shown to be the targets for miR160 out of several dozens of likely ARF family members. Great variation of the detected “tag abundance” existed per gene/miRNA species/library, from a single digit to over a thousand (HF12099, polyphenol oxidase), although this value can be dependent on the expression levels of both the target gene and miRNA. For some genes, such as those encoding laccases, it appeared the cleavage activity was reduced in P. ultimum infected tissues for both resistant and susceptible genotypes (L2 and L3, respectively), as compared with the mock-inoculation control tissue (L1). For other genes such as those encoding MYBs targeted by miR858, an elevated cleavage activity was observed for most target genes. In many cases, more than one member of the same miRNA family were shown to target the same gene ID, such as those genes encoding “growth regulating factors” that were targeted by miR396 species and genes encoding “auxin response factor” targeted by miR167 species. In contrast, many other genes showed a strict one (miRNA species) to one (gene ID) interaction pattern, such as those genes encoding laccase and targeted by miR397 species and R genes targeted by miR482 species. These identified target genes consist of valuable candidates for subsequent functional analysis to define their specific contribution to apple roots resistance to P. ultimum.
Table 4. Identified target genes, their functional annotations and detected tag abundance in apple roots in response to Pythium ultimum infection.
Three laccase encoding genes (laccase −3, −5, and −7), the target genes of miR397b based on degradome sequencing analysis, showed mostly upregulated expression patterns in response to P. ultimum infection (Figure 4). Among them, laccase-7 appeared to be more responsive to P. ultimum infection especially for resistance genotypes. For example, the detected transcript level of laccase-7 was increased three-fold in the infected root tissues of the resistant genotype #161, as compared to the mock-inoculated control tissue. The upregulated expression of these target genes corresponded to the downregulated expression of miR397b in response to P. ultimum infection. The observation of upregulated expression of these apple laccase genes in P. ultimum infected apple root tissues is consistent with the findings from a previous comparative transcriptome analysis between a resistant genotypes G.935 and a susceptible genotype B.9, which showed a peak response at 48 hpi and with stronger induction in the roots of the resistant genotype (Zhu et al., 2019). Therefore, the expression profiles of these laccase encoding genes are consistent with the findings from degradome sequencing, i.e., the reduced expression of miR397b and attenuated cleavage events likely contributed to the upregulated expression of these laccase genes in the infected apple root tissues.
Figure 4. Genotype-specific expression patterns of three laccase encoding genes in apple roots in response to Pythium ultimum infection. (A) Laccase-3; (B) Laccase-5 and (C) Laccase-7. The expression levels of these genes were measured by qRT-PCR. Genotypes are listed along X axis, with susceptible genotypes grouped at the left panel, and resistant genotypes grouped at right panel. Blue bar represents value of mock-inoculated root tissue and green bar represents the value of P. ultimum inoculated root tissue for each genotype. The values on Y axis represent the fold change of detected transcript level after normalizing to a reference gene (MD02G1221400) and calibrating by values of mock inoculated root tissue of the susceptible genotype #106. The values are the average of three biological replicates, and error bars are the standard deviation (SD) of fold changes. Error bars represent SE of three independent biological replicates. Asterisks denote statistically significant differences in a one-tailed Mann–Whitney U test (∗P < 0.05).
Genotype-Specific Expression of Selected miRNAs and Their Target Cleavage Activity
MiR397b-known-5p-mature is one of the four members from this family abundantly expressed in apple roots, and its cleavage activity on target (laccase) genes was only detected for miR397b (Tables 4, 5). One of the notable features regarding the genotype-specific miR397b expression, similar to miR397a (Figure 2), was that it expressed at a significantly lower level in mock-inoculated root tissues of resistant genotypes (except R-#58) compared to that of the susceptible genotypes. In response to P. ultimum infection, its expression was downregulated for all genotypes. Nevertheless, the lowest average value of read counts was detected in the infected root tissues of the resistant genotypes. The data from degradome sequencing indicated that three laccase encoding genes that are homologous to Arabidopsis thaliana laccase −3, −5, and −7, were the targets of miR397b. Furthermore, a substantially higher cleavage activity was observed in the infected root tissues of susceptible genotypes (L3), compared to that in resistant genotypes (L2). Taken together, the reduction of miR397b may result in less cleavage of laccase target mRNA, which could lead to higher laccase enzyme activity, possibly contributing to the observed resistance to P. ultimum infection.
Table 5. The genotype-specific expression of miR397b and target genes in response to Pythium ultimum infection of apple roots.
Resistance protein encoding genes, or R genes, are evidently the key players during plant-pathogen interactions. The differential expression patterns of miR482 between R and S genotypes and specifically identified targeted R genes by degradome sequencing were among the key findings from the current study (Table 1 and Supplementary Table 1). For example, miR482c-known-3p-mature is one of the three members in miR482 family which showed intriguing expression patterns among six genotypes and in response to P. ultimum infection (Table 6). In response to P. ultimum inoculation, the levels of miR482c were downregulated in both R and S genotypes, but with slightly larger reduction (>1x) in R group than that of S groups (<1x). The degradome sequencing identified the target of miR482c-known-3p-mature is Arabidopsis RPM1 homolog encoding gene, with two cleavage sites. Based on the detected values of tag abundance in three different libraries, it appeared that a more active cleavage activity occurred in susceptible genotypes than in resistant genotypes. Notably, the same RPM homolog gene was identified as one of the downregulated R genes specifically in susceptible B.9 plants from a previous comparative transcriptome analysis (Zhu et al., 2019). Therefore, at transcriptional or post-transcriptional levels, a consistent regulation scheme seemed to corroborate each other.
Table 6. The genotype-specific expression of miR482c and target genes in response to Pythium ultimum infection of apple roots.
The expression of miR10986 was generally detected at a low to moderate level in apple roots (Table 7), but a relatively large variation at read count values (2–3x) were observed among genotypes within the same R or S groups, or even within replicates for a given sample (genotype/treatment). Such uncommon variation among genotypes (within the same genotype group) probably indicated that its expression is also prone to certain abiotic stress conditions. In response to P. ultimum infection, downregulation was observed in all genotypes, but the degree of downregulation was more substantial (>1x) in resistant genotypes compared to that of susceptible genotypes. Similar to the regulation features for miR397, miR10986 also showed the lower basal expression in mock-inoculated tissues of R group as compared to that in susceptible genotypes. Degradome analysis indicated a polyphenol oxidase (PPO) is the cleavage target of miR10986. The overall cleavage activity on these PPO genes (for example, HF12100) was one of the strongest among the observed miRNA-target pairs from this dataset. However, a clear trend was missing, as the stronger cleavage activity was on either library L2 or L3 depending on the choices of cleavage sites. Therefore, the role of miR10986-PPO regulation during interaction between apple root and P. ultimum deserves future investigation.
Table 7. The genotype-specific expression of miR10986 and target genes in response to Pythium ultimum infection of apple roots.
Both the mature and star forms of miR7122a were detected with comparable abundance in apple roots (Table 8). In susceptible genotypes, their expression levels were slightly downregulated due to P. ultimum infection (<25% for both mature and star form). In contrast, variable regulation patterns were exhibited among resistant genotypes, i.e., slight upregulated for its mature form (∼5%) and substantially downregulated (47.8%) for its star form. Analysis of degradome sequencing data demonstrated that a homolog to the “putative pentatricopeptide repeat (PPR)-containing protein At1g12700” was the target of miR7122a. Slightly elevated cleavage activities on this target gene occurred in the resistant genotypes (L2). The roles of PPR proteins in plant immunity have been reported (Schmitz-Linneweber and Small, 2008; Barkan and Small, 2014) in other pathosystems, the definitive contribution of this miR-target pair in shaping up the resistance traits in apple roots to P. ultimum infection deserve subsequent research.
Table 8. Genotype-specific expression of miR7122a-known-5p-mature, 3P-star, and identified target genes in apple roots in response to Pythium ultimum infection.
PhasiRNA Analysis for miRNA390 and miRNA482/2118
PhasiRNAs are another main class of small RNAs in plants. Interestingly, we found that overall production of phasiRNAs was significantly lower in the P. ultimum-inoculated root tissues for both susceptible and resistant materials, suggesting that phasiRNA pathway contributes to the defense to P. ultimum. miR390 triggers the production of trans-acting siRNA3 (TAS3)-derived tasiRNAs to repress auxin responsive factor 2/3/4 (ARF2/3/4) genes, critical for auxin signaling; these tasiRNAs are known as tasiARFs. In plants, 5′ proximal miR390 target site on TAS3 is non-cleavable while 3′ proximal target site is sufficient for miR390-directed slicing, leading to a “two-hit, one-cleavage” model (Xia et al., 2017). In this study, the expression levels of miR390, tasiARF and all siRNAs from TAS3 transcripts were all significantly downregulated in P. ultimum-inoculated root tissues (Figure 5A), indicating that auxin signaling module regulated by the miR390-TAS3-ARF pathway were turned down by the P. ultimum infection. Detailed examination of the mapping profile of a TAS3 gene revealed that the miR390 cleavage site set the phase of the phasiRNA production, and indeed the abundance of sRNAs in mock-inoculated sample was much higher than that in P. ultimum-inoculated sample (Figure 5B).
Figure 5. miRNA triggers phasiRNA biogenesis in response to P. ultimum inoculation. (A) The expression levels of miR390, tasiARF, and siRNA from TAS3 transcripts in mock-inoculated (blue bars) and P. ultimum-inoculated (green bars) root issues. (B) Read abundance, sRNA abundance and phasing score distribution of a TAS3 locus in mock-inoculated and P. ultimum-inoculated root issues were viewed in IGV-sRNA. The yellow dashed lines and rectangle indicate the target site and slicing site of miR390a. (C) The expression levels of miR482, miR2118, and siRNA from NBS-LRR transcripts (HF37161) in mock-inoculated (blue bars) and P. ultimum-inoculated (green bars) root tissues. (D) Read abundance, sRNA abundance, and phasing score distribution of a representative PHAS locus in mock-inoculated and P. ultimum-inoculated root tissues were viewed in IGV-sRNA. The yellow dashed lines and rectangle indicate the target site of miR482c. Phasing score was calculated based on the mapping results of 21-nucleotide siRNAs (Xia et al., 2013). Error bars represent SE of three independent biological replicates. Asterisks denote statistically significant differences in a one-tailed Mann–Whitney U test (*P < 0.05).
MiR482/2118 is a well-known miRNA family important for disease resistance. All miRNA members of this family, including both miR482 and miR2118 variants are 22-nt long, and are capable to target NBS-LRR genes to trigger phasiRNAs biogenesis (Cuperus et al., 2010; Shivaprasad et al., 2012). In this study, a significant reduction on the expression levels of miR482, miR2118, and siRNAs from NBS-LRR transcripts were found in the P. ultimum-inoculated root tissues (Figure 5C), elucidating that miR482-mediated silencing of NBS-LRR genes was released in root tissues upon P. ultimum inoculation to activate defense. Detailed examination of the mapping profile of a NBS-LRR gene (R protein RPM1 homolog, HF37161) revealed that the miR482 cleavage site set the phase of the phasiRNA production and the abundance of phasiRNAs was greatly decreased in P. ultimum-inoculated roots (Figure 5D). These results regarding phasiRNAs suggested that apple trees turned down the endogenous sRNA silencing pathway to promote quick resistance response to the infection of P. ultimum. This kind of responses is likely widely present in plants.
Genomic Location of miR397a in Apple Genome
A comparison of miR397a precursor sequences against Robusta 5 de novo genome assembly identified four contig sequences (Figure 6A); two of which aligned to the same regions of chromosome 5 and chromosome 10, respectively, of the Golden Delicious Double Haploid genome (Daccord et al., 2017). The two overlapping contigs on same chromosomes represented the distinct haplotypes of the heterozygous apple genome. In addition, the regions of chromosomes 5 and 10 represent paralogs that were differentiated from a single ancestral chromosome through ancient whole genome duplication of the apple genome (Velasco et al., 2010). Annotation of 20 kb sequences around miR397a genomic loci identified a zinc-finger domain protein, an unknown protein, and a major facilitator superfamily protein. The unknown protein sequence appeared to match the promoter region of a miR397a gene from Pyrus pyrifolia (NCBI blast identifier: KY438936.1). Further comparison of Robusta 5 and Ottawa 3, across the miR397a targeted 20 kb genomic regions, identified 1,278 variants between them (Figure 6B). Co-linear contigs on chromosome 5 had relatively more variants than on chromosome 10 between Robusta 5 and Ottawa 3 on these regions. An alignment of P. ultimum susceptible (S) and resistant (R) genotype groups against Robusta 5 contigs resulted in 27 and 33 unique variants specific to each group, respectively (Figure 6B). A higher coverage of these regions with more sequencing reads can likely identify more variants in these regions.
Figure 6. Genomic location of miR397a, expression variant analysis. (A) Illustration of physical map around MIR397a. The two regions on chromosome 5 and 10 that matches MIR397a containing contigs are shown. The relative upstream and downstream sequences from MIR397a position was used to identify the annotated genes including zinc finger domain protein and a major facilitator superfamily protein on the two apple chromosomes. (B) Distribution of variants across the four MIR397a matching “Robusta 5” contigs after its comparison with the “O3” parent (gray dots), resistance (red dots), and susceptible (blue dots) genotypes.
Discussion
Upon pathogen infection, plants activate a sophisticated defense system that initiates massive reprogramming of global gene expression. As an integral part of the host transcriptome themselves, sRNAs are versatile and important post-transcriptional regulators of gene expression in almost all cellular processes including plant-pathogen interactions (Staiger et al., 2013; Weiberg et al., 2014; Yang and Huang, 2014; Chaloner et al., 2016). These sRNAs negatively regulate gene expression by binding and cleavage of their specific target mRNAs, or repressing the translation process (Bartel, 2004; Axtell and Bowman, 2008). The current study is part of an integrated effort to define the resistance mechanisms in apple roots to infection from P. ultimum, through the identification of key miRNA families and their corresponding targets.
The majority of identified miRNA families exhibited downregulated expression patterns in response to P. ultimum infection, although variable degrees of downregulation existed between resistant and susceptible genotypes. It is apparent that most of these target genes positively regulate defense activation in apple roots, and upon pathogen infection the reduced or attenuated cleavage activities from corresponding miRNAs lead to elevated defense activation. On the other hand, there is a possibility that the pathogen-derived effectors, toxins, or even mobile sRNAs, effectively and non-selectively sabotaged host miRNA pathways and many other cellular processes (Qiao et al., 2013; de Vries et al., 2017; Hou et al., 2019). Therefore, existence of non-host (or passive) resistance, or certain forms of preformed defense systems, such as the chemo/physical barrier impeding the initial pathogen penetration, may be a critical contributor to resistance to this necrotrophic pathogen. The microscopic features of necrosis progression or arrested necrotic progression could partially owe to the enhanced tissue fortification which impede pathogen hyphae penetration at the initial stage of infection and enhance the plant survivability for those resistant genotypes (Zhu et al., 2018b).
In response to pathogenic pressure, an increased level of tissue fortification such as lignin deposition has been hypothesized as a fundamental resistance mechanism (Vance et al., 1980; Bhuiyan et al., 2009; Miedes et al., 2014; Yang et al., 2018). In this dataset, four miRNA families which belong to copper (Cu) microRNAs (Pilon, 2017), i.e., laccase-targeting miR397, SOD-targeting miR398, nudix hydrolase-targeting miR408, and PPO-targeting miR10986, were markedly downregulated in infected root tissue. The targets of these miRNAs are known to directly participate in abiotic and biotic stresses through increased H2O2 generation, enhanced lignin content and modified cell wall composition (Lu et al., 2013). The levels of downregulation for both miR397 and miR10986 were more obvious in resistant genotypes, than in susceptible genotypes. Interestingly, even in mock-inoculated root tissues of resistant apple rootstock genotypes, much lower levels of miR397a/b and miR10986 were observed. The elevated activities of laccase and PPO appear to exist prior to pathogen exposure in resistant genotypes, which likely function as a preformed defense mechanism. Upon pathogen challenge, the stronger downregulation of laccase-targeting miRNA was observed in resistant genotypes, possibly leading to a quicker and stronger tissue fortification for an effective impediment of pathogen progression at the initial stages of infection. As an important part of the inducible defense mechanisms, production of secondary metabolites such as phenolics from phenylpropanoid pathway and the formation of lignin are believed to be responsible for delay or stop of the establishment of pathogen infection in many pathosystems (Danielsson et al., 2011; Oliva et al., 2015). Gene expression regulation such as that from MYB TFs and microRNAs as well as biochemical feature such as secondary metabolites from phenylpropanoid pathway and defense phytohormones appear to be evolutionarily conserved from non-flowering land plants and even streptophyte algae to angiosperms (Overdijk et al., 2016; de Vries et al., 2018a; Carella et al., 2019). Future studies are needed to define the detailed functional roles of Cu-miRNAs in shaping up resistance traits in apple roots to P. ultimum infection.
Post-transcriptional regulation of R genes is known to be a crucial aspect of plant immune response (Staiger et al., 2013; Weiberg et al., 2014; Liu et al., 2017; Islam et al., 2018). The roles of miR482/2118 on the expression of NBS-LRRs in tomato (Solanum lycopersicum) during infection by Phytophthora infestans infection are well investigated; and overexpression of mir482c was shown to induce enhanced susceptibility to late blight in tomato (de Vries et al., 2018b; Jiang et al., 2018; Canto-Pastor et al., 2019; Hong et al., 2019). Interestingly, in addition to a decreased expression of target NBS-LRR genes, the transgenic tomato plants overexpressing miR482c also showed lower peroxidase (POD), superoxide dismutase (SOD), and phenylalanine ammonia-lyase (PAL) activities, indicating the interconnected immune responses (Hong et al., 2019). Our data revealed a consistent downregulation of miR482/miR2118 superfamily across genotypes, which indicated that regulation of R gene expression, through miRNA function, is essential for the resistance traits of apple roots to P. ultimum infection. A slightly higher level of downregulation appeared to be the case in infected roots of resistant genotypes. In many plant species, including perennial woody plants, the miRNA–NB-LRR interactions result in the production of phased, secondary small interfering (phasi)RNAs, which function in both cis and trans, to amplify the regulatory effects across additional members of the target gene family (Katiyar-Agarwal and Jin, 2010). Therefore, it is likely that the subtle variation of expression of miR482/miR2118 leads to amplified effect on phasiRNA biogenesis, which may become an important factor differentiating apple roots resistance vs. susceptibility.
Hormone signaling and their interactions with corresponding transcription factors (TFs) are the integrated modules regulating plant defense response, and both elements are the preferred targets of miRNA regulation (Curaba et al., 2014). At the same time, both miRNAs and TFs are the primary regulators of gene expression, often for the same target genes at transcriptional and/or post-transcriptional levels (Li et al., 2014; Samad et al., 2017). The complicated interactions among these regulatory components (miRNA-TF-hormone signaling) are the common theme in many pathosystems (Balmer and Mauch-Mani, 2013; Weiberg and Jin, 2015; Chaloner et al., 2016; Liu et al., 2017). From our dataset, miR393, miR160, and miR167, all of which target auxin signaling, were among a few of the identified miRNA families with mostly upregulated expression upon P. ultimum infection. Interestingly, miR160 showed partially differential regulation between R and S groups, where two out of three resistant genotypes exhibited downregulated expression in response to P. ultimum infection. Results from previous transcriptome analyses also identified regulatory components related to JA, auxin, WRKY, MYB, and NAC, which are largely consistent with the findings in the current study; (Shin et al., 2016b; Zhu et al., 2017, 2019). The role of these regulatory components in secondary metabolites including monolignols or the precursors of lignin deserves further investigation.
The identification of bona fide targets of a given miRNA species represents one of the fundamental aspects of small RNA research. Plant miRNAs can target multiple non-paralogous but functionally related genes (Djami-Tchatchou et al., 2017; You et al., 2017), which are often from a large gene family. Utilizing powerful degradome sequencing, a few target genes specific to this pathosystem were experimentally identified (rather than by in silico prediction). For example, three miRNA families (miR159, miR319, and miR858) are known to target a large number of MYB TFs in the apple genome, this approach pinpointed less than a dozen of MYB TF encoding genes which are potentially essential in this pathosystem. As an extreme example, only a single SPL encoding gene was identified as the target of miR156. In addition, novel targets were identified beyond the “commonly” established spectrum of genes functioning in this pathosystem such as those genes encoding PPRs (pentatricopeptide repeat-containing proteins). Most plants contain several hundred PPRs, the RNA-binding proteins with proposed functions of RNA editing, RNA splicing, RNA cleavage and translation within mitochondria and chloroplasts (Barkan and Small, 2014; Park et al., 2014) which has been linked to plant resistance response (Schmitz-Linneweber and Small, 2008; Barkan and Small, 2014). Moreover, the quantified “tag abundance” also provides the estimated miRNA cleavage activities between resistant or susceptible genotype groups. For example, the detected tag abundance for two target gene IDs “HF12099-RA and HF12100-RA” (polyphenol oxidase) were about 50% in infected susceptible root tissues (L3), compared to the values of mock-inoculated control tissues (L1), or infected resistant root tissues (L2). These identified target genes constitute a focused candidate gene list, which may significantly contribute to apple root resistance traits through post-transcriptional regulations.
As part of the integrated effort to uncover the underpinned molecular mechanism controlling apple roots resistance to P. ultimum, the current dataset of miRNA profiling and degradome sequencing offers a unique perspective on its post-transcriptional regulation. The identified “miRNA-target gene” pairs likely represent the crucial components functioning in apple root defense activation to P. ultimum infection. As examples, genotype-specific expression patterns for miR397-laccse, miR10986-PPO, and miR398-SOD appeared to link with apple root resistance vs. susceptibility. Cell wall fortification through the function of copper protein encoding genes could be a key strategy to defend apple roots against this necrotrophic pathogen. Young apple roots, as a primary organ for taking up water and nutrients, lack effective protection layers such as cuticle or wax in aerial parts of a plant to deter or impede penetration from pathogen like P. ultimum. Winning the chemical war at the early stage of infection is crucial to thwart disruptive arsenals of effectors and toxins from this fast-growing necrotrophic pathogen. Therefore, quick and effective cell wall modification through the function of these Cu-miRNAs (Pilon, 2017) can significantly determine the outcome between the apple roots and P. ultimum. Such initial defense reactions are likely “buying the time” for a steady and robust defense activation to be developed. Inclusion of multiple genotypes with contrasting resistance phenotypes enhanced the reliability of this comprehensive dataset from this study. Those miRNA families targeting R genes, hormonal signaling, and various TFs represent the core regulation circuits for this pathosystem. Noticeably, many of the identified target genes in the current study are consistent with the results from previous transcriptome analyses (Shin et al., 2016b; Zhu et al., 2017, 2019). As an example, the same RPM1 homolog gene was also identified from a susceptible genotype B.9 with downregulated expression (Zhu et al., 2019). In summary, the identified miRNA species and their corresponding target genes constitute a valuable resource for subsequent functional analysis to define their detailed roles in shaping up the resistance traits in apple roots against P. ultimum infection. The findings from this study significantly contribute to the effort for deciphering the molecular mechanisms controlling defense responses for this less explored pathosystem between apple roots and a necrotrophic pathogen.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
YZ, RX, AK, and GF participated in the experimental design, data analysis, interpretations, and manuscript writing. MS performed most of the experiments of phenotyping and RNA preparation. GL and JS analyzed the dataset and contributed to the manuscript preparation. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the USDA-ARS base fund.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank India Amy Cain for language editing and reviewing the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.629776/full#supplementary-material
Footnotes
- ^ www.faststone.org
- ^ https://www.rosaceae.org/species/malus_x_domestica_HFTH1/genome_v1.0
- ^ ncbi.nlm.nih.gov/sra/review/SRP295189
- ^ http://www.rosaceae.org
- ^ www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi
- ^ www.idtdna.com/analyzer/Applications/OligoAnalyzer/
- ^ https://www.rosaceae.org/analysis/162
References
Addo-Quaye, C., Miller, W., and Axtell, M. J. (2009). CleaveLand: a pipeline for using degradome data to find cleaved small RNA targets. Bioinformatics 25, 130–131. doi: 10.1093/bioinformatics/btn604
Axtell, M. J. (2013). Classification and comparison of small RNAs from plants. Annu. Rev. Plant Biol. 64, 137–159.
Axtell, M. J., and Bowman, J. L. (2008). Evolution of plant microRNAs and their targets. Trends Plant Sci. 13, 343–349.
Balmer, D., and Mauch-Mani, B. (2013). Small yet mighty–microRNAs in plant-microbe interactions. MicroRNA 2, 73–80.
Barkan, A., and Small, I. (2014). Pentatricopeptide repeat proteins in plants. Annu. Rev. Plant Biol. 65, 415–442.
Bhuiyan, N. H., Selvaraj, G., Wei, Y., and King, J. (2009). Role of lignification in plant defense. Plant Signal. Behav. 4, 158–159.
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120.
Bologna, N. G., and Voinnet, O. (2014). The diversity, biogenesis, and activities of endogenous silencing small RNAs in Arabidopsis. Annu. Rev. Plant Biol. 65, 473–503. doi: 10.1146/annurev-arplant-050213-035728
Borges, F., and Martienssen, R. A. (2015). The expanding world of small RNAs in plants. Nat. Rev. Mol. Cell Biol. 16:727.
Canto-Pastor, A., Santos, B. A., Valli, A. A., Summers, W., Schornack, S., and Baulcombe, D. C. (2019). Enhanced resistance to bacterial and oomycete pathogens by short tandem target mimic RNAs in tomato. Proc. Natl. Acad. Sci. U.S.A. 116, 2755–2760. doi: 10.1073/pnas.1814380116
Carella, P., Gogleva, A., Hoey, D. J., Bridgen, A. J., Stolze, S. C., Nakagami, H., et al. (2019). Conserved biochemical defenses underpin host responses to oomycete infection in an early-divergent land plant lineage. Curr. Biol. 29, 2282.–2294. doi: 10.1016/j.cub.2019.05.078
Chaloner, T., van Kan, J. A., and Grant-Downton, R. T. (2016). RNA ‘Information Warfare’in pathogenic and mutualistic interactions. Trends Plant Sci. 21, 738–748. doi: 10.1016/j.tplants.2016.05.008
Chen, C., Zeng, Z., Liu, Z., and Xia, R. (2018). Small RNAs, emerging regulators critical for the development of horticultural traits. Horticult. Res. 5, 1–14. doi: 10.1038/s41438-018-0072-8
Cuperus, J. T., Carbonell, A., Fahlgren, N., Garcia-Ruiz, H., Burke, R.T., Takeda, A., et al. (2010). Unique functionality of 22-nt miRNAs in triggering RDR6-dependent siRNA biogenesis from target transcripts in Arabidopsis. Nat. Struct. Mol. Biol. 17, 997–1003. doi: 10.1038/s41438-018-0072-8
Curaba, J., Singh, M. B., and Bhalla, P. L. (2014). miRNAs in the crosstalk between phytohormone signalling pathways. J. Exp. Bot. 65, 1425–1438.
Daccord, N., Celton, J.-M., Linsmith, G., Becker, C., Choisne, N., Schijlen, E., et al. (2017). High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat. Genet. 49:1099. doi: 10.1038/ng.3886
Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., DePristo, M. A., et al. (2011). The variant call format and VCFtools. Bioinformatics 27, 2156–2158.
Danielsson, M., Lundén, K., Elfstrand, M., Hu, J., Zhao, T., Arnerup, J., et al. (2011). Chemical and transcriptional responses of Norway spruce genotypes with different susceptibility to Heterobasidion spp. infection. BMC Plant Biol. 11:154. doi: 10.1186/1471-2229-11-154
De Paoli, E., Dorantes-Acosta, A., Zhai, J., Accerbi, M., Jeong, D.-H., Park, S., et al. (2009). Distinct extremely abundant siRNAs associated with cosuppression in petunia. RNA 15, 1965–1970. doi: 10.1261/rna.1706109
de Vries, S., de Vries, J., and Rose, L. E. (2019). The elaboration of miRNA regulation and gene regulatory networks in plant–microbe interactions. Genes 10:310. doi: 10.3390/genes10040310
de Vries, S., de Vries, J., von Dahlen, J. K., Gould, S. B., Archibald, J. M., Rose, L. E., et al. (2018a). On plant defense signaling networks and early land plant evolution. Commun. Integr. Biol. 11, 1–14.
de Vries, S., Kukuk, A., von Dahlen, J. K., Schnake, A., Kloesges, T., and Rose, L. E. (2018b). Expression profiling across wild and cultivated tomatoes supports the relevance of early miR482/2118 suppression for Phytophthora resistance. Proc. R. Soc. B 285:20172560. doi: 10.1098/rspb.2017.2560
de Vries, S., von Dahlen, J. K., Uhlmann, C., Schnake, A., Kloesges, T., and Rose, L. E. (2017). Signatures of selection and host-adapted gene expression of the Phytophthora infestans RNA silencing suppressor PSR2. Mol. Plant Pathol. 18, 110–124. doi: 10.1111/mpp.12465
Djami-Tchatchou, A. T., Sanan-Mishra, N., Ntushelo, K., and Dubery, I. A. (2017). Functional roles of microRNAs in agronomically important plants—potential as targets for crop improvement and protection. Front. Plant Sci. 8:378. doi: 10.3389/fpls.2017.00378
Eckardt, N. A. (2012). A microRNA cascade in plant defense. Plant Cell 24:840. doi: 10.1105/tpc.112.240311
Fei, Q., Xia, R., and Meyers, B. C. (2013). Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks. Plant Cell 25, 2400–2415.
Fei, Q., Zhang, Y., Xia, R., and Meyers, B. C. (2016). Small RNAs add zing to the zig-zag-zig model of plant defenses. Mol. Plant Microbe Interact. 29, 165–169. doi: 10.1094/MPMI-09-15-0212-FI
German, M. A., Pillay, M., Jeong, D.-H., Hetawal, A., Luo, S., Janardhanan, P., et al. (2008). Global identification of microRNA–target RNA pairs by parallel analysis of RNA ends. Nat. Biotechnol. 26:941. doi: 10.1038/nbt1417
He, L., and Hannon, G. J. (2004). MicroRNAs: small RNAs with a big role in gene regulation. Nat. Rev. Genet. 5:522. doi: 10.1038/nrg1379
Hong, Y.-H., Meng, J., He, X.-L., Zhang, Y.-Y., and Luan, Y.-S. (2019). Overexpression of MiR482c in tomato induces enhanced susceptibility to late blight. Cells 8:822. doi: 10.3390/cells8080822
Hou, Y., Zhai, Y., Feng, L., Karimi, H. Z., Rutter, B. D., Zeng, L., et al. (2019). A Phytophthora effector suppresses trans-kingdom RNAi to promote disease susceptibility. Cell Host Microbe 25, 153–165 e155. doi: 10.1016/j.chom.2018.11.007
Islam, W., Qasim, M., Noman, A., Adnan, M., Tayyab, M., Farooq, T. H., et al. (2018). Plant microRNAs: front line players against invading pathogens. Microb. Pathog. 118, 9–17. doi: 10.1016/j.micpath.2018.03.008
Jaffee, B., Abawi, G., and Mai, W. (1982). Fungi associated with roots of apple seedlings grown in soil from an apple replant site. Plant Dis. 66, 942–944.
Jiang, N., Meng, J., Cui, J., Sun, G., and Luan, Y. (2018). Function identification of miR482b, a negative regulator during tomato resistance to Phytophthora infestans. Horticult. Res. 5, 1–11. doi: 10.1038/s41438-018-0017-2
Johnson, M., Zaretskaya, I., Raytselis, Y., Merezhuk, Y., McGinnis, S., and Madden, T. L. (2008). NCBI BLAST: a better web interface. Nucleic Acids Res. 36(suppl_2), W5–W9. doi: 10.1093/nar/gkn201
Katiyar-Agarwal, S., and Jin, H. (2010). Role of small RNAs in host-microbe interactions. Annu. Rev. Phytopathol. 48, 225–246.
Koren, S., Walenz, B. P., Berlin, K., Miller, J. R., Bergman, N. H., and Phillippy, A. M. (2017). Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736. doi: 10.1101/gr.215087.116
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760.
Li, J., Reichel, M., Li, Y., and Millar, A. A. (2014). The functional scope of plant microRNA-mediated silencing. Trends Plant Sci. 19, 750–756. doi: 10.1016/j.tplants.2014.08.006
Liu, Q., and Axtell, M. J. (2015). Quantitating plant microRNA-mediated target repression using a dual-luciferase transient expression system. Methods Mol. Biol. 1284, 287–303. doi: 10.1007/978-1-4939-2444-8_14
Liu, S. R., Zhou, J. J., Hu, C. G., Wei, C. L., and Zhang, J. Z. (2017). MicroRNA-mediated gene silencing in plant defense and viral counter-defense. Front. Microbiol. 8:1801. doi: 10.3389/fmicb.2017.01801
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 25, 402–408.
Lu, S., Li, Q., Wei, H., Chang, M.-J., Tunlaya-Anukit, S., Kim, H., et al. (2013). Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa. Proc. Natl. Acad. Sci. U.S.A. 110, 10848–10853. doi: 10.1073/pnas.1308936110
Mazzola, M. (1997). Identification and pathogenicity of Rhizoctonia spp. isolated from apple roots and orchard soils. Phytopathology 87, 582–587. doi: 10.1094/PHYTO.1997.87.6.582
Mazzola, M. (1998). Elucidation of the microbial complex having a causal role in the development of apple replant disease in washington. Phytopathology 88, 930–938. doi: 10.1094/PHYTO.1998.88.9.930
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The Genome Analysis Toolkit: a MAPREDUCE framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi: 10.1101/gr.107524.110
Miedes, E., Vanholme, R., Boerjan, W., and Molina, A. (2014). The role of the secondary cell wall in plant resistance to pathogens. Front. Plant Sci. 5:358. doi: 10.3389/fpls.2014.00358
Oliva, J., Rommel, S., Fossdal, C., Hietala, A., Nemesio-Gorriz, M., Solheim, H., et al. (2015). Transcriptional responses of Norway spruce (Picea abies) inner sapwood against Heterobasidion parviporum. Tree Physiol. 35, 1007–1015. doi: 10.1093/treephys/tpv063
Overdijk, E. J., De Keijzer, J., De Groot, D., Schoina, C., Bouwmeester, K., Ketelaar, T., et al. (2016). Interaction between the moss Physcomitrella patens and Phytophthora: a novel pathosystem for live-cell imaging of subcellular defence. J. Microsc. 263, 171–180. doi: 10.1111/jmi.12395
Park, Y. J., Lee, H. J., Kwak, K. J., Lee, K., Hong, S. W., and Kang, H. (2014). MicroRNA400-guided cleavage of pentatricopeptide repeat protein mRNAs renders Arabidopsis thaliana more susceptible to pathogenic bacteria and fungi. Plant Cell Physiol. 55, 1660–1668. doi: 10.1093/pcp/pcu096
Qiao, Y., Liu, L., Xiong, Q., Flores, C., Wong, J., Shi, J., et al. (2013). Oomycete pathogens encode RNA silencing suppressors. Nature Genet. 45, 330–333.
Rogers, K., and Chen, X. (2013). Biogenesis, turnover, and mode of action of plant microRNAs. Plant Cell 25, 2383–2399.
Samad, A. F., Sajad, M., Nazaruddin, N., Fauzi, I. A., Murad, A., Zainal, Z., et al. (2017). MicroRNA and transcription factor: key players in plant regulatory network. Front. Plant Sci. 8:565. doi: 10.3389/fpls.2017.00565
Schmitz-Linneweber, C., and Small, I. (2008). Pentatricopeptide repeat proteins: a socket set for organelle gene expression. Trends Plant Sci. 13, 663–670. doi: 10.1016/j.tplants.2008.10.001
Shivaprasad, P. V., Chen, H.-M., Patel, K., Bond, D. M., Santos, B. A., and Baulcombe, D. C. (2012). A microRNA superfamily regulates nucleotide binding site–leucine-rich repeats and other mRNAs. Plant Cell 24, 859–874. doi: 10.1016/j.tplants.2008.10.001
Shin, S., Lee, J., Rudell, D., Evans, K., and Zhu, Y. (2016a). Transcriptional regulation of auxin metabolism and ethylene biosynthesis activation during apple (Malus× domestica) fruit maturation. J. Plant Growth Regul. 35, 655–666.
Shin, S., Zheng, P., Fazio, G., Mazzola, M., Main, D., and Zhu, Y. (2016b). Transcriptome changes specifically associated with apple (Malus domestica) root defense response during Pythium ultimum infection. Physiol. Mol. Plant Pathol. 94, 16–26.
Singh, J., and Khan, A. (2019). Distinct patterns of natural selection determine sub-population structure in the fire blight pathogen, Erwinia amylovora. Sci. Rep. 9, 1–13. doi: 10.1038/s41598-019-50589-z
Staiger, D., Korneli, C., Lummer, M., and Navarro, L. (2013). Emerging role for RNA-based regulation in plant immunity. New Phytol. 197, 394–404. doi: 10.1111/nph.12022
Tewoldemedhin, Y. T., Mazzola, M., Botha, W. J., Spies, C. F., and McLeod, A. (2011). Characterization of fungi (Fusarium and Rhizoctonia) and oomycetes (Phytophthora and Pythium) associated with apple orchards in South Africa. Eur. J. Plant Pathol. 130, 215–229.
Vance, C., Kirk, T., and Sherwood, R. (1980). Lignification as a mechanism of disease resistance. Annu. Rev. Phytopathol. 18, 259–288.
Velasco, R., Zharkikh, A., Affourtit, J., Dhingra, A., Cestaro, A., Kalyanaraman, A., et al. (2010). The genome of the domesticated apple (Malus× domestica Borkh.). Nat. Genet. 42:833.
Weiberg, A., and Jin, H. (2015). Small RNAs—the secret agents in the plant–pathogen interactions. Curr. Opin. Plant Biol. 26, 87–94. doi: 10.1016/j.pbi.2015.05.033
Weiberg, A., Wang, M., Bellinger, M., and Jin, H. (2014). Small RNAs: a new paradigm in plant-microbe interactions. Annu. Rev. Phytopathol. 52, 495–516. doi: 10.1146/annurev-phyto-102313-045933
Xia, R., Meyers, B. C., Liu, Z., Beers, E. P., Ye, S., and Liu, Z. (2013). MicroRNA superfamilies descended from miR390 and their roles in secondary small interfering RNA biogenesis in eudicots. Plant Cell 25, 1555–1572. doi: 10.1105/tpc.113.110957
Xia, R., Xu, J., Arikit, S., and Meyers, B. C. (2015). Extensive families of miRNAs and PHAS loci in Norway spruce demonstrate the origins of complex phasiRNA networks in seed plants. Mol. Biol. Evol. 32, 2905–2918. doi: 10.1093/molbev/msv164
Xia, R., Zhu, H., An, Y.-Q., Beers, E. P., and Liu, Z. (2012). Apple miRNAs and tasiRNAs with novel regulatory networks. Genome Biol. 13:R47. doi: 10.1186/gb-2012-13-6-r47
Xia, R., Xu, J., and Meyers, B. C. (2017). The emergence, evolution, and diversification of the miR390-TAS3-ARF pathway in land plants. Plant Cell 29, 1232–1247. doi: 10.1186/gb-2012-13-6-r47
Yang, C., Liang, Y., Qiu, D., Zeng, H., Yuan, J., and Yang, X. (2018). Lignin metabolism involves Botrytis cinerea BcGs1-induced defense response in tomato. BMC Plant Biol. 18:103. doi: 10.1186/s12870-018-1319-0
Yang, L., and Huang, H. (2014). Roles of small RNAs in plant disease resistance. J. Integr. Plant Biol. 56, 962–970. doi: 10.1111/jipb.12200
You, C., Cui, J., Wang, H., Qi, X., Kuo, L.-Y., Ma, H., et al. (2017). Conservation and divergence of small RNA pathways and microRNAs in land plants. Genome Biol. 18:158.
Zhang, B. (2015). MicroRNA: a new target for improving plant tolerance to abiotic stress. J. Exp. Bot. 66, 1749–1761. doi: 10.1093/jxb/erv013
Zhang, L., Hu, J., Han, X., Li, J., Gao, Y., Richards, C. M., et al. (2019). A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour. Nat. commun. 10, 1–13. doi: 10.1038/s41467-019-09518-x
Zhou, Z., Tian, Y., Cong, P., and Zhu, Y. (2018). Functional characterization of an apple (Malus x domestica) LysM domain receptor encoding gene for its role in defense response. Plant Sci. 269, 56–65. doi: 10.1016/j.plantsci.2018.01.006
Zhu, Y., and Saltzgiver, M. (2020). A systematic analysis of apple root resistance traits to Pythium ultimum infection and the underpinned molecular regulations of defense activation. Horticult. Res. 7:62. doi: 10.1038/s41438-020-0286-4
Zhu, Y., Fazio, G., and Mazzola, M. (2014). Elucidating the molecular responses of apple rootstock resistant to ARD pathogens: challenges and opportunities for development of genomics-assisted breeding tools. Horticult. Res. 1:14043. doi: 10.1038/hortres.2014.43
Zhu, Y., Saltzgiver, M., and Zhao, J. (2018a). A phenotyping protocol for detailed evaluation of apple root resistance responses utilizing tissue culture micropropagated apple plants. Am. J. Plant Sci. 9:2183.
Zhu, Y., Shao, J., Zhou, Z., and Davis, R. E. (2017). Comparative transcriptome analysis reveals a preformed defense system in apple root of a resistant genotype of G.935 in the absence of pathogen. Int. J. Plant Genomics 2017:8950746. doi: 10.1155/2017/8950746
Zhu, Y., Shao, J., Zhou, Z., and Davis, R. E. (2019). Genotype-specific suppression of multiple defense pathways in apple root during infection by Pythium ultimum. Horticult. Res. 6:10. doi: 10.1038/s41438-018-0087-1
Keywords: apple root, soilborne pathogens, defense activation, plant resistance, small RNA profiling, degradome sequencing, post-transcriptional regulation
Citation: Zhu Y, Li G, Singh J, Khan A, Fazio G, Saltzgiver M and Xia R (2021) Laccase Directed Lignification Is One of the Major Processes Associated With the Defense Response Against Pythium ultimum Infection in Apple Roots. Front. Plant Sci. 12:629776. doi: 10.3389/fpls.2021.629776
Received: 15 November 2020; Accepted: 09 August 2021;
Published: 07 September 2021.
Edited by:
Magdalena Arasimowicz-Jelonek, Adam Mickiewicz University, PolandReviewed by:
Mickael Malnoy, Fondazione Edmund Mach, ItalyLaura Ellen Rose, Heinrich Heine University Düsseldorf, Germany
Sophie De Vries, Heinrich Heine University Düsseldorf, Germany, in collaboration with reviewer LR
Copyright © 2021 Zhu, Li, Singh, Khan, Fazio, Saltzgiver and Xia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yanmin Zhu, yanmin.zhu@ars.usda.gov; Rui Xia, rxia@scau.edu.cn
†These authors have contributed equally to this work