1. Introduction
Capsicum or pepper, belonging to the genus Capsicum and family Solanaceae, is a nutritionally important vegetable that originated in Central and South America, the Caribbean, and Mexico, with at least 27 reported species [1,2]. Capsicum production has increased due to strong demand and culinary use of its fruit [3,4]; however, its growth is frequently threatened by different abiotic and biotic stresses and can be affected by those stresses simultaneously [3,4,5,6]. Therefore, it is of increasing importance to study the impact of combined stresses on capsicum plants, such as virus infection and high temperature.
Capsicum chlorosis virus (CaCV) is a serious pathogen that infects not only capsicum but also chili, tomato, pineapple, and peanut crops [7,8,9]. The virus has been reported in Australia, China, Greece, India, Iran, Taiwan, Thailand, and the USA [10]. In southern China, CaCV is considered a significant disease in peanuts, with incidence reaching up to 20% and causing noticeable yield losses, particularly when plants are infected at early growing stages [10,11]. In India, the virus affects the production of chili peppers, with disease incidences over 20% reported [10,12]. In Australia, the emergence of CaCV has been reported in large commercial capsicum production areas, including Bundaberg and the northern dry tropics of Queensland [9]. The typical CaCV symptoms on capsicum plants include stunting, marginal and interveinal chlorosis and leaf deformation on emerging leaves, and concentric chlorotic or necrotic lesions on mature leaves [9,13]. Recently, CaCV-resistant advanced breeding lines (PI 290972 × C. annuum cv. Mazurka and cv. Warlock inbred lines) were selected by the Queensland Department of Agriculture and Fisheries (DAF) breeding program [13,14]. CaCV is taxonomically classified in the species Capsicum chlorosis orthotospovirus in the genus Orthotospovirus, family Tospoviridae, order Bunyavirales [15]. The genome of this orthotospovirus consists of three segments of single-stranded RNA (ssRNA). Among them, the large (L) segment is of negative polarity and encodes the RNA-dependent RNA polymerase (RdRp), while both medium (M) and small (S) segments are of ambisense coding polarity with viral movement protein (NSm) and glycoproteins (Gn/Gc) encoded on the M segment, and a viral silencing suppresser (NSs) and nucleocapsid protein (N) encoded on the S segment [13,16,17].
The progression of virus diseases in plants is significantly affected by temperature since both virus propagation and plant development are temperature-dependent [18]. Given current concerns about global warming, elevated temperatures have been recognized as an important climate-changing variable impacting plant–virus interactions [19,20,21]. To date, evidence shows that exposure to elevated temperatures suppresses effector-triggered immunity (ETI)-mediated plant resistance [22,23,24]. ETI is a robust plant defense response triggered by recognizing pathogen effector molecules through nucleotide-binding leucine-rich repeat (NLR) proteins [25]. Typical ETI responses are associated with calcium production, salicylic acid (SA) accumulation, and a burst of oxidative reactive oxygen species (ROS). Furthermore, these plant responses result in two typical manifestations of disease resistance, hypersensitive response (HR) and systemic acquired resistance (SAR) [25]. Examples of temperature-mediated breaking of disease resistance include the potato virus Y (PVY)–potato and tomato spotted wilt virus (TSWV)–capsicum pathosystems [24,26,27,28,29]. In these pathosystems, virus-resistant plant varieties, carrying resistance (R) genes, are unable to mount effective ETI at elevated temperatures. The capsicum Tsw gene and potato R gene, Ny-DG—which confer resistance to TSWV and PVY, respectively—are compromised at 32 °C or above 28 °C [26,29].
MicroRNAs (miRNAs), which are post-transcriptional regulators, are found in most eukaryotes [30]. In plants, miRNAs are conserved regulators of plant developmental processes and stress responses [31,32]. The 20–24 nucleotide (nt) miRNAs originate from miRNA genes (MIR). These genes are transcribed by RNA polymerase II (Pol II) into a stem–loop structure, named primary (pri)-miRNAs. Pri-miRNAs are processed by DICER-LIKE 1 (DCL1) and several assisting proteins to form the precursor (pre)-miRNAs and mature miRNA duplex. At the last step of this pathway, the guide strand of the miRNA duplex is loaded onto specific ARGONAUTE (AGO) proteins for target mRNA cleavage or translational inhibition [33]. MiRNAs are essential elements affecting the plant–virus arms race [34,35]. Several miRNAs have been reported to act as regulators to support plant antiviral resistance [36,37]. For example, the expression of miR6019/miR6020, which downregulates expression of the NLR receptor N gene in tobacco plants, is reduced, allowing effective induction of NLR-mediated ETI during tobacco mosaic virus (TMV) infection [37]. During infection of tomato leaf curl New Delhi virus, Sly-miR159 appears to be downregulated, which increases the expression of its target SlMyB33 and the downstream R gene SlSw5a in resistant tomatoes. This modulation enhances the plant defense against the virus by triggering HR [38]. Conversely, some miRNAs act as regulators that are favorable for virus infections, like miR319 [39,40]. The expression of miR319 in rice plants is induced by rice ragged stunt virus infections. This leads to a suppression of jasmonic acid (JA)-mediated defenses, which results in an increased susceptibility of rice plants to virus disease [39].
Although the interplay between plant miRNAs and negative-sense RNA viruses has been explored [41], the role of miRNAs in capsicum resistance to CaCV at different temperatures remains unknown. Therefore, this study investigated the effects of high temperature on capsicum resistance to CaCV through phenotypic and miRNA analysis. Small RNA (sRNA) high-throughput sequencing (HTS) was used for a systemic comparative analysis of differentially expressed miRNAs in plants challenged with or without CaCV at ambient (25 °C) and elevated temperatures (35 °C/30 °C), elucidating miRNA-related underlying mechanisms involved in the capsicum–CaCV pathosystem.
2. Materials and Methods
2.1. Plants and Growth Conditions
CaCV-resistant capsicum plants from the third backcross generation of Capsicum chinense PI 290972 × commercial CaCV-susceptible C. annuum cultivars [14] were grown in a temperature-controlled glasshouse compartment at an ambient temperature (AT) of approximately 25 °C. Seedlings aged four weeks with two true leaves were mechanically inoculated with CaCV isolate QLD 3432 [17]. Inoculum was prepared by grinding fresh CaCV-infected symptomatic capsicum leaves of the susceptible cultivar Yolo Wonder in 10 mM phosphate buffer, pH 7.6, containing freshly added 20 mM sodium sulphite with a mortar and pestle. Mock treatment used buffer only. Inoculum and buffer were rub-inoculated onto the carborundum-dusted first two leaves and two cotyledons. Subsequently, half of each plant group was transferred to a growth cabinet at high temperature (HT) of 35 °C/30 °C (16 h day/8 h night), light intensity of 230 μmol·m−2·s−1, and relative humidity of 60%.
2.2. RNA Isolation, Small RNA Library Construction. and Sequencing
Capsicum plants, treated under four conditions, were used for sRNA library construction: mock-inoculated plants grown at AT (‘AM’ hereafter); CaCV-inoculated plants grown at AT (‘AV’ hereafter); mock-inoculated plants grown at HT (‘HM’ hereafter); and CaCV-inoculated plants grown at HT (‘HV’ hereafter). Systemic leaves from three individual plants, each grown in one of those four conditions, were collected at 10 dpi. For each biological replicate, three leaf disks from the third leaf from the top were sampled. Total RNA from sample lysates in lysis buffer (ISOLATE II RNA Plant Kit, Bioline, London, UK) was purified by Direct-zol RNA Miniprep kit (Zymo Research, Irvine, CA, USA) according to the manufacturer’s instructions. On-column DNase I treatment was performed using the same kit. Total RNA (2–3 µg per sample) was sent to the Australian Genome Research Facility (AGRF, Melbourne, VIC, Australia) for sRNA library construction and HTS. The sRNA libraries were prepared using NEBNext® Small RNA Library Prep Set (New England Biolabs, Ipswich, MA, USA) for Illumina® following the manufacturer’s instructions. Single-end 50 bp sequencing was performed on an Illumina NOVASEQ 6000, AGRF, Melbourne, VIC, Australia.
2.3. Small RNA Sequencing Data Pre-Processing and microRNA Identification
Raw reads generated by sRNA sequencing were pre-processed using the source code of iwa software (
The MiRNA Compilation and MiRNA Selection modules in iwa software were used to identify the putative known and novel miRNAs [42]. Raw reads were uploaded into iwa software and were pre-processed by the same parameters mentioned above. In the MiRNA Compilation module, miRDeep-P2 was selected for miRNA prediction [44]. The PmiREN (
2.4. Analysis of Virus-Derived siRNAs
The genome of CaCV-Qld-3432, containing 8913 nt in the L segment, 4846 nt in the M segment, and 3944 nt in the S segment (GenBank accession numbers KM589495, KM589494, KM589493) [17], was used for sRNA alignment. The profile of 21-nucleotide, 22-nucleotide, and 24-nucleotide vsiRNA normalized read counts was obtained using SCRAM aligner module (
2.5. Differential Expression Analysis of miRNAs
Differentially expressed (DE) miRNAs in the two treatments were analyzed using DEseq2 (Galaxy Version 2.11.40.6 + galaxy1) [53]. Here, four pairwise comparisons were performed: AV vs. AM; HV vs. HM; HM vs. AM; and HV vs. AV. Total read counts for each miRNA were normalized by median ratio normalization. The resulting p-values were corrected for multiple testing using Benjamini and Hochberg’s false discovery rate (FDR) [54]. MiRNAs were judged to be DE if the FDR-adjusted p-value (Padj) was ≤0.01 and log2 fold-change was >1.0 or <−1.0. Heatmaps were created using bioinfokit python package (Version 2.0.8) to visualize significantly DE miRNAs [55].
2.6. miRNA Target Prediction and Enrichment Analysis
PsRNAtarget (2017 release,
Single enrichment analysis (SEA) was performed using agriGO v2.0 (
2.7. Evaluation of the Expression of miRNAs and Their Potential Target Genes
The abundance of selected miRNAs was examined using linear specific (S)-poly (A)-tailed real-time RT-PCR, as described by Xie and collaborators [61,62,63]. Total RNA purification and DNase treatment were performed using the methods described in Section 2.2 above. DNase-treated RNAs were then polyadenylated using poly (A) polymerase (New England BioLabs, Ipswich, MA, USA) according to the manufacturer’s instructions. For each miRNA, cDNA was synthesized from polyadenylated RNA (1 µg) by using Superscript IV reverse transcriptase (Invitrogen, Carlsbad, CA, USA) with S-poly (A)-tailed RT primers (Table S1). These RT primers (SRT) that target specific miRNAs were designed using sRNAprimerDB software version 1.0 (
Northern blot hybridization was used for the detection and quantification of mature Can-miR164b/c since a specific product of mature miRNA164b/c was unable to be amplified by using linear S-poly (A)-tailed real-time RT-PCR. Total RNA was separated on a 17% polyacrylamide gel [67]. Subsequentially, gel-separated RNA was transferred onto a Hybond-N+ nylon membrane (Roche) using a Bio-Rad mini trans-blot system. For U6 and Can-miR164b/c detection, DIG-labeled DNA probes of U6 (5′-TCATCCTTGCGCAGGGGCCA) and Can-miR164b/c (5′-TGCACGTGCCCTGCTTCTCCA) were generated using the DIG Oligonucleotide 3′-End Labelling Kit (Roche, Basel, Switzerland). The relative accumulation level of Can-miR164b/c was then calculated by densitometry of the Northern blot chemiluminescent images using iBright analysis software (Version 5.0.0). The densitometry value of the AM condition was set as reference in each blot. The normalized U6 values served as an internal control for the normalized values of Can-miR164b/c.
Statistical tests were performed using GraphPadPrism software (Version 9.3.1). Data were analyzed with unpaired Student’s t-tests and considered significantly different if the two-tailed p-value was <0.05.
2.8. Cleavage Sites Validation of Target Genes of miRNAs
Modified 5′ RNA ligase-mediated amplification of cDNA ends (5′-RLM-RACE) was used to validate the cleavage sites of predicted miRNA target genes [68,69]. This modified procedure starts with ligating the 5′-RNA adapter (RLM_RNA_adapter, Table S1) to total RNAs using T4 RNA ligase (New England Biolabs, Ipswich, MA, USA). Then, the adapter-ligated RNAs were used to synthesize first-strand cDNA using Superscript IV reverse transcriptase (Invitrogen) with the oligo dT primer following the manufacturer’s instructions. The cDNA was used as the template for touchdown PCR using Phusion High-Fidelity DNA Polymerase (New England Biolabs, Ipswich, MA, USA) with the primers RLM, gene-specific primers (RLM_GSP)—including RLM_AGO1b_GSP, RLM_LAC4_GSP, and RLM_CLAVATA_GSP (Table S1)—and GeneRacer_5P. The cycling program of touchdown PCR was set as the initial annealing temperature of 68 °C with 3 cycles, followed by a reduced annealing temperature of 66 °C with 3 cycles, 63 °C with 7 cycles, and 60 °C with 22 cycles. The products of the first amplification were used as the templates for the second PCR with nested primer pairs RLM_nGSP/GeneRacer_n5P (Table S1) that annealed internally to the first primer pairs. The products of the second PCR were then cloned into pGEM®-T Easy vector (Promega, Madison, WI, USA) and sequenced to analyze the cleavage sites.
3. Results
3.1. CaCV Symptom Development in Resistant Capsicum Plants at Elevated Temperatures
To evaluate the effect of elevated temperature on capsicum resistance to CaCV, symptoms were observed over a 10-day period in CaCV-infected capsicum plants grown at HT and AT. Four-week-old CaCV-resistant capsicum plants grown at AT were inoculated with CaCV prior to half of them being transferred to HT. At 5 dpi, all resistant plants challenged with CaCV showed HR-mediated necrotic spots on virus-inoculated leaves at both AT and HT (Figure 1a). At 10 dpi, 8 out of 14 (57%) CaCV-infected plants grown at HT showed distinct necrotic spots on upper systemic leaves, while none of the CaCV-infected plants grown at AT had symptoms on the systemic leaves (Figure 1b). All mock-inoculated CaCV-resistant plants grown at HT or AT remained symptomless and had a similar appearance at the two temperatures. These results indicate that the HR was successfully triggered by CaCV infection in the inoculated leaves of CaCV-resistant capsicum, which was independent of growth temperature. However, the HR appears to be unable to restrict CaCV infection in some plants grown at HT, allowing CaCV to move to systemic leaves and trigger HR there.
3.2. Small RNAs in CaCV-Infected and Virus-Free Capsicum Grown at Different Temperatures
To identify miRNAs that may be involved in capsicum resistance to CaCV at HT and/or AT, systemic leaves from three CaCV-infected and three mock-inoculated plants grown at either temperature were collected. For HV treatment, leaf disks were sampled from three plants with systemic necrotic spot symptoms. A total of twelve sRNA libraries, containing three libraries in each of the four conditions, were generated for Illumina sequencing and the data are summarized in Table S2. In total, 123,888,579, 95,340,112, 160,227,368, and 127,753,253 raw reads were obtained from the libraries of the AV, AM, HV, and HM treatments, respectively. After trimming adaptor sequences and low-quality reads, 121,643,325, 93,441,595, 157,681,773, and 125,749,412 cleaned reads, and 58,305,917, 45,565,384, 102,316,388, and 65,410,603 length-limited reads remained for AV, AM, HV, and HM treatments, respectively (Table S2). In the analysis of length distribution, 24-nucleotide sRNAs were the largest sRNA group among all length-limited reads in the libraries of AV, AM, and HM treatments (Figure 2a,b,d), whereas 21-nucleotide and 22-nucleotide sRNAs accounted for the highest percentage of all sRNAs in the HV libraries in addition to the 24-nucleotide sRNA peak (Figure 2c). To further analyze vsiRNAs and miRNAs, length-limited reads were aligned to the capsicum genome. An average of 70.0%, 70.2%, 56.4%, and 73.3% of reads for the AV, AM, HV, and HM treatments, respectively, mapped to the capsicum reference genome. The reads that failed to map to the capsicum genome in the AV and HV libraries were retrieved for vsiRNA analysis and indicated that the expression of CaCV-derived vsiRNAs was significantly higher in CaCV-infected capsicum grown at HT than in those plants grown at AT (Figure 3, compare panels a and b). Visualization of vsiRNA coverage in CaCV-infected capsicum grown at HT showed that vsiRNAs were more evenly dispersed across the entire L segment than across M and S segments (Figure 3a). Viral siRNAs were less abundant across the intergenic region of both M and S segments, indicating vsiRNA hotspots in the coding regions of M and S segments (Figure 3a). Moreover, the abundance of 21-nucleotide, 22-nucleotide, and 24-nucleotide vsiRNAs at the higher temperature was highest in several hotspots in the CaCV S segment in all three HV libraries (Figure 3a).
3.3. Identification of Known and Novel miRNAs
The sRNA datasets were grouped according to the four treatments (AV, AM, HV, and HM), and the tool miRDeep-P2 in iwa-miRNA software was used to predict miRNA clusters. The capsicum miRNAs were annotated using the miRNA annotations retrieved from PmiREN, sRNAanno, and PsRNA databases. After aggregating annotated and predicted miRNAs into a candidate miRNAs list, miRNAs that passed both high-throughput-based and machine learning-based criteria were selected. From the four grouped datasets, 105 known and 83 novel miRNAs were identified. Among all identified miRNAs, 89 known and 47 novel miRNAs were present in libraries across more than one treatment. In addition, 69 known and 21 novel miRNAs were identified in all libraries (Figure S1). The mature, star, and precursor sequences of identified miRNAs present in libraries across more than one treatment are listed in Table S3. The abundance of miRNAs in individual sRNA libraries was determined by aligning reads against the mature, star, and precursor sequences of those identified miRNAs. They were given the prefix ‘Can-’ for Capsicum annuum. Overall, Can-MIR159 and Can-MIR166 were the most abundant known miRNA families and Can-MIRN19 was the most abundant novel miRNA family in all libraries across all treatments.
3.4. Differentially Expressed miRNAs
The differential expression of miRNAs was calculated using DEseq2 in four pairwise comparisons: AV vs. AM; HV vs. HM; HM vs. AM; and HV vs. AV. Five DE miRNAs were identified at AT when comparing CaCV-infected and mock-inoculated treatments, which included one upregulated and four downregulated miRNAs (Figure 4a). In HV vs. HM and HM vs. AM comparisons, 34 and 48 DE miRNAs, respectively, were detected. Of those DE miRNAs, 18 upregulated and 16 downregulated miRNAs were identified in the HT CaCV-infected treatment versus mock treatment; 8 upregulated and 40 downregulated miRNAs were identified in mock treatments comparing HT and AT conditions (Figure 4b,c).
To identify miRNAs that are potentially involved in temperature-sensitive resistance to CaCV, we focused on the comparative expression levels of miRNAs in HV and AV conditions. A total of 35 DE miRNAs, including 21 upregulated and 14 downregulated miRNAs, were identified at HT compared to AT in CaCV-infected plants (Figure 4d). Among those, DE miRNAs, Can-MIR408, Can-MIR398, Can-MIR397, Can-MIR393, Can-MIRN492, and Can-MIRN7 were the most highly upregulated miRNA families with log2-fold changes above 2.0 in HV compared to AV conditions (Table 1). Interestingly, Can-miR408a-3p, Can-miR398b-3p, Can-miR397-3p, and Can-miR397-5p were upregulated when comparing the HV and HM treatments but were downregulated when comparing the AV and AM treatments (Table 1). For those downregulated miRNAs in the CaCV-infected treatment comparing HT and AT growth conditions, some miRNAs in the Can-MIR169, Can-MIR477, Can-MIRN482, and Can-MIR399 families were responsive to HT and showed a reduced expression in HM versus AM treatment. Notably, Can-miR157b-5p, Can-miR157c-5p, Can-miR164b-5p, Can-miR164c-5p, and Can-miR_N80a-3p were downregulated only in HV compared to AV treatment but not in HM compared to AM treatment (Table 1). These data indicate that increased or reduced expression of some miRNAs in Can-MIR408, Can-MIR398, Can-MIR397, Can-MIR157, Can-MIR164, and Can-MIR_N80a families may be associated with CaCV-induced changes between HT and AT treatments.
3.5. Functional Analysis of Predicted Target Genes of DE miRNAs in Response to CaCV under Ambient and High-Temperature Conditions
For a better understanding of the potential regulatory roles of DE miRNAs during CaCV infection and at elevated temperature, target genes of all identified miRNAs were predicted by psRNAtarget software with stringent criteria (expectation score ≤ 3) and were functionally annotated with GO terms. An overview of all DE miRNAs and their predicted mRNA targets in the four pairwise comparisons AV vs. AM, HV vs. HM, HM vs. AM, and HV vs. AV is shown in the Venn diagrams depicted in Figure S2. A large proportion of the upregulated miRNAs and their predicted targets in HV compared to AV overlapped with those upregulated miRNAs and their predicted targets in HV compared to HM treatments (Figure S2a,c), suggesting a significant virus effect. However, downregulated miRNAs and their predicted targets in HV compared to AV overlapped with those downregulated miRNAs and their predicted targets in HM compared to AM (Figure S2b,d), suggesting a significant temperature effect.
To identify pathways that are relevant to the observed phenotypic differences between HV and AV treatments, GO enrichment analyses were conducted for the targets that were predicted from up- or downregulated miRNAs. Besides analyzing targets of DE miRNAs in HV vs. AV treatments, targets of DE miRNAs in other pairwise comparisons (AV vs. AM, HV vs. HM, and HM vs. AM) were also included. All significantly enriched (Padj < 0.05) GO terms describing biological process (BP), molecular function (MF), and cellular component (CC) for the predicted targets of DE miRNAs are listed in Tables S4–S9. A total of 156 targets were predicted from 21 miRNAs that were upregulated in CaCV-infected plants at HT compared to AT. Among them, 143 miRNA targets could be functionally annotated. Subsequently, 8 and 12 GO terms listed in MF and BP were enriched through the SEA of those miRNA targets. The hierarchical relationships of the significantly enriched terms for BP and MF are presented in Figure 5 and Figure 6, respectively. DE miRNAs and their corresponding target genes with enriched GO terms are shown in Table S4. Briefly, miRNA targets that were annotated with significantly enriched GO terms in BP and MF categories—including the phenylpropanoid metabolic process (GO:0009698), lignin catabolic process (GO:0046274), lignin metabolic process (GO:0009808), phenylpropanoid catabolic process (GO:0046271), hydroquinone: oxygen oxidoreductase activity (GO:0052716), and copper ion binding (GO:0005507)—were regulated by Can-miR408a-3p and Can-miR397-5p (Figure 5 and Figure 6 and Table S4). In agreement with the relevance of DE miRNAs among the pairwise comparisons HV vs. AV, HV vs. HM, and AV vs. AM, the most highly enriched GO terms were also enriched through SEA of targets of downregulated miRNAs in AV compared to AM as well as targets of upregulated miRNAs in HV compared to HM (Figures S3 and S4). In addition, the GO-enriched term, cysteine-type peptidase activity (GO:0008234) for MF, was associated with a group of targets that were affected by upregulated miRNAs (Can-miR_N7b-3p, Can-miR_N7c-3p, Can-miR_N7d-3p, miRN37-Novel-5p, and Can-miR408a-5p) in HV compared to both AV and HM treatments (Figure 6, Figure S4 and Tables S4 and S6).
For miRNAs that were downregulated at higher temperature in CaCV-infected plants, a total of 48 target mRNAs were functionally annotated. Among those targets, only two GO terms (GO:0003677 and GO:0003676) were enriched in the MF category. Moreover, these two enriched terms were associated with targets that were predicted from eight miRNAs, including Can-miR157b-5p, Can-miR157c-5p, Can-miR164b-5p, Can-miR164c-5p, Can-miR169h-3p, Can-miR477a-5p, Can-miR_N80a-3p, and Can-miRN482-5p (Table S5). Based on the functional prediction of miRNA targets and the DE miRNAs, Can-miR408a-3p, Can-miR397-5p, Can-MIR157, and Can-MIR164 families may represent important known miRNAs involved in temperature-dependent resistance-breaking.
3.6. Validation of the Expression Patterns of Selected miRNAs and Their Target Genes
To confirm the sRNA sequencing data, S- poly(A)-tailed real-time RT-PCR or Northern blot with signal intensity quantification was used to quantify the levels of selected miRNAs at 10 dpi. Based on the results of the differential expression analysis of miRNAs and functional analysis of targets, Can-miR408a-3p, Can-miR397-5p, and Can-miR164b/c-5p were selected based on their potential involvement in temperature-dependent resistance-breaking. In addition, Can-miR168-5p, which was upregulated with a log2-fold change of 1.0 and Padj of 0.017 in HV compared to the AV condition, was selected despite the cut-off threshold of Padj ≤ 0.01. This exception was based on the reported importance of miR168. Accumulation of miR168 compromises AGO1-mediated antiviral RNA silencing in several virus-infected plant species [70,71,72]. In addition, accumulation of miR168 may affect the expression of many target genes since AGO1 plays a central role in RNA silencing [73]. The results of the real-time RT-PCR and Northern blot (Figure S5) confirmed the increased expression levels of Can-miR408a-3p, Can-miR397-5p, and Can-miR168-5p as well as the decreased expression level of Can-miR164b/c-5p in HV compared to AV treatments, and showed that they were consistent with the expression profiles obtained by sRNA sequencing (Figure 7a). When HV was compared to HM, the expression trends of Can-miR408a-3p, Can-miR397-5p, and Can-miR164b/c-5p were similar between the results observed in real-time RT-PCR and sRNA sequencing (Figure 7b). However, expression trends of Can-miR408a-3p, Can-miR397-5p, and Can-miR168-5p in the AV and AM comparison were opposite between the results of real-time RT-PCR and sRNA sequencing (Figure 7c).
In addition to analyzing the expression of mature miRNAs, quantitative real-time RT-PCR was also used to investigate the effect of these miRNAs on their target genes in the HV condition. Targets were selected based on the target prediction data obtained from psRNAtarget and WPMIAS software [74] as well as the target validation data reported by Zhang and collaborators [75]. CAN.G1061.9 (CLAVATA1-related), CAN.G671.1 (Laccase 5, LAC5), and CAN.G1305.35 (Plantacyanin) were chosen for Can-miR408a-3p target analysis; CAN.G394.71 (LAC2), CAN.G351.2 (LAC4), and CAN.G355.8 (LAC4) were chosen for Can-miR397-5p target analysis; CAN.G637.6 (AGO1b) was chosen for Can-miR168-5p target analysis; NAM/ATAF1,2/CUC2 (NAC) transcription factors—CAN.G943.28 (NAC1) and CAN.G587.8 (NAC5)—were chosen for Can-miR164b/c-5p target analysis. Unexpectedly, expression patterns of the targets CAN.G671.1, CAN.G1305.35, CAN.G394.71, CAN.G351.2, CAN.G355.8, and CAN.G637.6 were not inverse to the expression of their corresponding miRNAs when comparing the HV and AV conditions (Figure 8c,d,f–h,m). Conversely, the expected inverse expression pattern was observed in Can-miR408a-3p/CAN.G1061.9, Can-miR164b/c-5p/CAN.G943.28, and Can-miR164b/c-5p/CAN.G587.8 pairs (Figure 8b,j,k). Notably, we found that the transcript levels of all three selected targets of Can-miR397-5p in mock-inoculated plants drastically increased at HT (Figure 8f–h). Moreover, the transcript level of one of the targets, CAN.G355.8, was downregulated, while that of the other two targets (Figure 8h), CAN.G394.71 and CAN.G351.2, remained unchanged in CaCV-infected plants as compared with mock-inoculated plants at HT (Figure 8f,g). These results suggest that the downregulation of CAN.G355.8 may be due to the negative regulation of upregulated Can-miR397-5p in HV compared to HM. Similarly, an inverse expression pattern between Can-miR168-5p and CAN.G637.6 was also observed in HV vs. HM (Figure 8m). Different from other targets that showed inverse or no correlation tendency with miRNAs, CAN.G1305.35 showed a positively correlated expression pattern with Can-miR408a-3p in three pairwise comparisons (AV vs. AM, HV vs. HM, and HV vs. AV) (Figure 8d). Overall, diverse regulatory effects of miRNAs on their targets were observed when plants were treated at different temperatures and biotic stress conditions.
3.7. Identification of miRNA-Mediated Cleavage of Target Genes
To confirm putative miRNA targets, the specific cleavage sites were identified experimentally using 5′-RACE. Three targets, including CAN.G1061.9, CAN.G355.8, and CAN.G637.6, were chosen for the 5′-RACE verification based on the observed inverse expression patterns of miRNA/target pairs. CAN.G1061.9 was selected for its inverse expression pattern at HV vs. AV, while CAN.G355.8 and CAN.G637.6 were chosen for their inverse patterns at HV vs. HM. The other targets, CAN.G943.28 and CAN.G587.8, which showed opposite correlation with Can-miR164b/c-5p, were verified by degradome-based miRNAs-target analysis using WPMIAS (Table S10). Based on the results of 5′-RACE, the cleavage sites were all located in the regions that were predicted as miRNA-binding sites (Figure 9). The majority cleavage site in CAN.G1061.9 was mapped to the 11th position of the complementary sequence at the 5′-end of Can-miR408a-3p (Figure 9a). On the other hand, the main cleavage site in CAN.G355.8 and CAN.G637.6 was mapped to the 10th position of the complementary sequences at the 5′-end of Can-miR397-5p and Can-miR168-5p, respectively (Figure 9b,c). These results suggest that CAN.G1061.9, CAN.G355.8, and CAN.G637.6 can be cleaved by Can-miR408a-3p, Can-miR397-5p, and Can-miR168-5p, respectively, in capsicum.
4. Discussion
The suppressive effect of elevated temperatures on ETI-mediated resistance or its associated HR response has been demonstrated in both non-viral [20,76,77] and viral pathosystems [20,78]. For example, the N gene-mediated HR and resistance to TMV in tobacco are compromised by exposing plants to temperatures above 28 °C [20,78]. In the present study, we found that ETI-mediated resistance to CaCV was overcome at 35 °C/30 °C (day/night) in 8 out of 14 capsicum plants. Interestingly, instead of abolishing HR, distinct local lesions developed in inoculated leaves, and systemic necrotic lesions were also observed in these plants in which resistance is compromised at 35 °C. This result agrees with the previous finding suggesting that HR may fail to restrict virus movement in local lesions at elevated temperatures [26,27,79]. At 32 °C, the Tsw-mediated resistance to TSWV in capsicum is compromised, leading to virus movement into uninoculated leaves. Nonetheless, local HR persists in the inoculated leaves of those plants [26,27,79]. Notably, there is limited knowledge regarding the mechanisms that underpin temperature-sensitive ETI-mediated resistance in viral pathosystems [20,29]. In the TMV–tobacco pathosystem, HR is suppressed at elevated temperatures by preventing NLR proteins from translocating into the nucleus to orchestrate plant immune signaling. Tobacco plants mount an effective HR response when the host N protein accumulates in the nucleus after recognition of the TMV coat protein at 22 °C, while plants are unable to induce an HR response at 28 °C due to decreased N protein nuclear localization [20]. In addition, changes in miRNA accumulation that affect the expression of their target genes appear to be linked to the temperature sensitivity of Ny-DG-mediated resistance in the PVY–potato pathosystem. Resistance to PVYNTN was compromised at 28 °C and an enhanced miRNA-mediated downregulation of a specific NLR transcript was observed in PVY-inoculated potato leaves [29].
In previous studies, several miRNAs have been identified in Capsicum spp. [41,80,81,82,83,84]. Some of these miRNAs were potentially involved in the interactions between capsicum and plant pathogens [41,84]. In the present study, Illumina sequencing of sRNAs was conducted to investigate the involvement of miRNAs in temperature-sensitive antiviral resistance. By predicting and retrieving miRNA annotations from three sRNA databases—PmiREN, sRNAanno, and PsRNA—105 known and 83 novel miRNAs were identified across four treatment groups: AV, AM, HV, and HM. Given the distinct phenotypic differences observed between CaCV-infected capsicum plants grown at HT and AT, the DE miRNAs revealed by comparing HV and AV treatments allowed us to identify miRNAs potentially involved in temperature-sensitive resistance. Among the DE miRNAs, a significant number expressed higher in the HV than AV were responsive to CaCV infection at HT. In addition, in the enrichment analysis of the upregulated miRNAs in HV compared to AV, targets of Can-miR408a-3p and Can-miR397-5p—which are associated with biological processes, including the metabolism and catabolism of phenylpropanoids and lignin—were identified as the most enriched GO terms.
Lignin deposition plays a crucial role in plant resistance to biotrophic pathogens by confining them to the infection site [85]. In A. thaliana, lignification is induced during the ETI response, triggered by infection of avirulent bacterial strains, Pseudomonas syringae pv. tomato (Pst) DC3000 (AvrRpm1) and Pst DC3000 (AvrRpt2). This lignification is associated with localized cell death, restricting the bacterial pathogens to the infection site [85]. Increased lignin accumulation has also been associated with plant resistance to viral pathogens, including TMV, tobacco necrosis virus, southern bean mosaic viruses, and sweetpotato virus disease (SPVD) [86,87]. We, therefore, speculate that the upregulation of Can-miR408a-3p and Can-miR397-5p in HV compared to AV may reduce lignin content by targeting genes in the LAC family, which in turn may interfere with restrictions to virus spread from the infection site. Although there is no direct evidence showing the involvement of miR408 and miR397 in lignin-related HR, reduced lignin deposition has been observed in plants overexpressing these two miRNAs (Lu et al. 2013; Song et al. 2018). Moreover, inhibition of miR397 has been shown to increase lignin deposition in the cell wall and reduce the accumulation of SPVD [87].
Based on computational results and their involvement in plant responses to various stresses [88], Can-miR408a-3p, Can-miR397-5p, Can-miR168-5p, and Can-miR164b/c-5p were selected for further analysis of sRNA and target genes expression profiles. Real-time RT-PCR revealed that Can-miR408a-3p, Can-miR397-5p, and Can-miR168-5p were upregulated, whereas Can-miR164b/c-5p was downregulated in HV compared to AV. Additionally, Can-miR408a-3p and Can-miR397-5p were upregulated in HV compared to HM, consistent with our sRNA sequencing data. However, discrepancies were observed between the results of sRNA sequencing and real-time RT-PCR in some pairwise comparisons. For instance, in AV and AM comparisons, real-time RT-PCR indicated an increased level of Can-miR408a-3p and an unchanged level of Can-miR397-5p. Conversely, sRNA sequencing showed decreased levels for both. This inconsistency between sequencing and RT-PCR data has also been documented in prior studies [89,90]. Several factors may contribute to discrepancies between computational sequence analysis and Northern blot or real-time RT-PCR. For instance, biases during sRNA library preparation and amplification, such as adapter ligation bias and varying GC content among miRNAs, can introduce distortions [91,92,93].
Genes that encode copper-containing proteins, including proteins in the phytocyanin family—cupredoxin, plantacyanin, and uclacyanin—or proteins involved in lignin polymerization such as LACs, have been validated as targets of miR408 in A. thaliana [94,95]. For miR397, genes in the LAC family are its major target genes in various plant species [96,97,98,99,100]. In the present study, GO enrichment analysis suggests that LAC family genes, predicted to be regulated by Can-miR408a-3p and Can-miR397-5p, may be crucial for resistance breakdown at HT. Unexpectedly, real-time RT-PCR analysis did not reveal significant differences in the expression levels of the selected LAC genes, including CAN.G671.1, CAN.G394.7, CAN.G351.2, and CAN.G355.8, between plants treated with HV and AV. Their expression was not inversely correlated with the increased pattern of Can-miR408a-3p and Can-miR397-5p. However, negative regulation of CAN.G355.8 by Can-miR397-5p was observed in HV compared to HM conditions, wherein the CAN.G355.8 level was reduced, while the Can-miR397-5p level was increased in HV-treated plants. Furthermore, transcript levels of selected Can-miR397-5p targets, including CAN.G394.7, CAN.G351.2, and CAN.G355.8, were increased by HT in mock-inoculated plants. Given that reduced lignin content is detrimental to Medicago truncatula growth at HT [101,102], an induction of LACs at HT may be crucial for heat tolerance in capsicum plants. Taken together, we speculate that the suppression of CAN.G355.8 expression by CaCV at HT is likely to reduce plant tolerance to stresses, which may lead to CaCV resistance breakdown at HT.
In addition to the genes in the LAC family, the expression of other targets predicted for Can-miR408a-3p was investigated. CAN.G1061.9 (leucine-rich receptor-like kinase CLV1-related) was shown to be cleaved by Can-miR408a-3p through 5′-RACE in the present study. The transcript level of CAN.G1061.9 was lower in CaCV-infected plants grown at HT than those grown at AT, which showed an opposite tendency to Can-miR408a-3p in the HV vs. AV comparison. CAN.G1061.9 encodes a CLV1-related receptor kinase-like protein, which appears to be similar to BARELY ANY MERISTEM 1 (BAM2) in the Plaza 4.0 database [57]. Interestingly, BAM2, which functions redundantly with BAM1, participates in the regulation of cell-to-cell RNAi movement [103]. Therefore, a downregulation of CAN.G1061.9 may affect systemic spread of antiviral RNAi, resulting in reduced resistance of capsicum plants to CaCV at HT. Unlike CAN.G1061.9, a negative correlation between miRNA and its target did not occur for Can-miR408a-3p and CAN.G1305.35 (Plantacyanin). The transcript level of CAN.G1305.35 increased during Can-miR408a-3p induction in CaCV-infected plants either grown at HT or AT. This agrees with the prior finding that the levels of miR408 and its targets, Plantacyanin and Uclacyanin, increased simultaneously in late stages of natural senescence in A. thaliana [104]. Emerging evidence suggests that a negative correlation between miRNAs and their targets is not strictly required for target validation [105]. Such exceptions may occur when miRNAs and their targets are expressed partially overlapped or in a cell-type-specific manner [105,106,107]. For example, miR395 and its target, SULTR2;1, are mainly expressed in phloem companion cells and xylem parenchyma cells in roots, respectively [107]. This, therefore, prevents miR395 from downregulating SULTR2;1, leading to an upregulation of both SULTR2;1 and miR395 in roots during sulfur starvation [107].
In plants, ROS function as central regulators in complex signaling networks [108,109]. The ROS burst can trigger HR and induce plant resistance to viruses [25]. However, excessive ROS accumulation occurs when the homeostasis of ROS is imbalanced, leading to irreversible oxidative damage and accelerated senescence [109]. In our study, miRNA profiles were investigated in systemic, uninoculated leaves that displayed necrotic spots. This necrotic spot phenotype may be associated with an increased accumulation of ROS in systemic leaves at elevated temperature [110]. Interestingly, not only an upregulation of miR408 but also a downregulation of miR164, which have been shown to be associated with accelerated leaf senescence [88,111], were observed in HV compared to AV. Repressing miR164 expression was found to induce the expression of its target gene, ORE1, leading to early-senescence phenotypes [111,112]. In addition, NAC transcription factors targeted by miR164, such as ANAC021/22 (NAC1), ANAC079/80 (NAC4), and ANAC100 (NAC5), were upregulated during senescence [113]. MiRNA168 is crucial in miRNA and siRNA pathways by targeting AGO1, the core component of the RNA-induced silencing complex [114]. A well-known auto-regulatory feedback loop of miR168 and AGO1 is involved in maintaining their homeostasis [114]. During several viral infections in A. thaliana and N. benthamiana, increased miR168 levels disrupt this homeostasis, resulting in an inhibition of AGO1-mediated antiviral RNA silencing [71,72]. In our study, the transcript levels of two selected targets, CAN.G943.28 (NAC1) and CAN.G587.8 (NAC5), were upregulated in HV compared to AV, which was negatively correlated with the downregulated Can-miR164b/c-5p. However, no significant expression differences of CAN.G637.6 (AGO1b) were observed in HV compared to AV. Overall, our findings suggest that the upregulation of Can-miR408a-3p and Can-miR397-5p, as well as the downregulation of Can-miR164b-5p and Can-miR164c-5p, may underpin the temperature-sensitive resistance-breaking phenotype of CaCV-infected capsicum at elevated temperature. The functions of these miRNAs, including regulating lignin deposition and leaf senescence, may affect important cellular mechanisms in CaCV-infected capsicum plants, which may lower plant fitness under stress conditions and further impact the capacity of plant resistance to viruses.
5. Conclusions
In this study, we investigated the effect of elevated temperature on a CaCV-resistant capsicum breeding line. Our results showed that ETI-mediated resistance to CaCV is compromised at elevated temperatures (35 °C/30 °C) in some capsicum plants. Through sRNA HTS and computational analysis, we identified 105 known and 83 novel miRNAs across different treatment groups. Of particular interest was the upregulation of Can-miR408a-3p and Can-miR397-5p in HV plants compared to AV among all DE miRNAs. GO enrichment analysis of the target genes (LACs) predicted from these miRNAs suggested the potential involvement of phenylpropanoid (GO:0009698) and lignin (GO:0009808) metabolic processes in the capsicum response to CaCV under different temperatures. MiRNAs, including Can-miR408a-3p, Can-miR397-5p, Can-miR164b/c-5p, and Can-miR168-5p, may underlie the temperature-sensitivity of ETI-mediated resistance in CaCV-infected capsicum plants.
Conceptualization, W.-A.T., N.M. and R.G.D.; methodology, W.-A.T., N.M., C.A.B. and R.G.D.; validation, W.-A.T., N.M., C.A.B. and R.G.D.; investigation, W.-A.T.; resources, N.M. and R.G.D.; data curation, W.-A.T.; writing—original draft preparation, W.-A.T.; writing—review and editing, N.M., C.A.B. and R.G.D.; visualization, W.-A.T.; supervision, N.M., C.A.B. and R.G.D.; project administration, W.-A.T. and R.G.D. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
All data can be accessed on the NCBI database under Bioproject PRJNA1128345.
The Authors acknowledge Ritesh Jain, Fernanda Borges Naito, and Stephen Fletcher for scientific discussion; Shirani Widana Gamage for data availability of transcriptome of capsicum plants; Galaxy Australia and UQ High-performance computing for data analysis.
The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Effect of elevated temperature on capsicum chlorosis virus (CaCV) infection in CaCV-resistant capsicum plants: (a) at 5 days post inoculation (dpi), a hypersensitive response (red arrows) was apparent in CaCV-inoculated leaves of capsicum plants grown at high temperature (HT) of 35 °C and ambient temperature (AT) of 25 °C; (b) at 10 dpi, necrotic spots (indicated by red arrows) were apparent on systemic leaves of some capsicum plants grown at HT but not on any plants grown at AT. Three representative plants are shown for each treatment.
Enlarge this image.Figure 2. Small RNA length distribution of 18-nucleotide to 26-nucleotide (million (M) reads) in 12 small RNA libraries. Each diagram represents three independent libraries constructed from samples collected from (a) capsicum chlorosis virus (CaCV)-inoculated plants grown at ambient temperature (AV), (b) mock-inoculated plants grown at ambient temperature (AM), (c) CaCV-inoculated plants grown at higher temperature (HV), and (d) mock-inoculated plants grown at higher temperature (HM).
Enlarge this image.Figure 3. Profiles of capsicum chlorosis virus (CaCV)-derived viral siRNAs. The 21-, 22-, and 24-nucleotide vsiRNA coverage across the CaCV genome L, M, and S segments in CaCV-infected capsicum plants grown at (a) higher temperature (HV), and (b) ambient temperature (AV). Standard error of 3 biological replicates is presented as the smoothed plots. The reads-per-million (RPM) scale is set at ±1000 for three datasets of the HV treatment and at ±20 for three datasets of the AV treatment.
Enlarge this image.Figure 4. Heat maps showing differentially expressed capsicum miRNAs in four pairwise comparisons: (a) capsicum chlorosis virus (CaCV)-infected capsicum grown at ambient temperature (AT) compared to mock-inoculated capsicum grown at AT (AV vs. AM); (b) CaCV-infected capsicum grown at higher temperature (HT) compared to mock-inoculated capsicum grown at HT (HV vs. HM); (c) mock-inoculated capsicum grown at HT compared to AT (HM vs. AM); (d) CaCV-infected capsicum grown at HT compared to AT (HV vs. AV); Columns represent independent biological replicates; rows represent different miRNAs. Clustering is based on Z-score hierarchical calculation. Red and blue indicate miRNAs with high and low expression, respectively, as shown on the color scale.
Enlarge this image.Figure 5. Acyclic graph showing hierarchical relationship of the enriched gene ontology (GO) terms associated with targets predicted from upregulated miRNAs in CaCV-infected capsicum grown at high temperature (HV) compared to ambient temperature (AV) in the biological process category. The color scale from yellow to red indicates an increasingly significant enrichment of GO terms.
Enlarge this image.Figure 6. Acyclic graph showing hierarchical relationship of the enriched gene ontology (GO) terms associated with targets predicted from upregulated miRNAs in CaCV-infected capsicum grown at high temperature (HV) compared to ambient temperature (AV) in the molecular function category. The color scale from yellow to red indicates an increasingly significant enrichment of GO terms.
Enlarge this image.Figure 7. Linear specific (S)-poly (A)-tailed quantitative real-time RT-PCR (qPCR) of miR408, miR397, and miR168 or Northern blot hybridization confirmation of miR164b and miR164c expression patterns compared to those obtained by sRNA high-throughput sequencing (HTS). The expression levels (log2-fold change) of five miRNAs in four pairwise comparisons: (a) CaCV-infected plants at high temperature (HV) vs. CaCV-infected plants at ambient temperature (AV); (b) HV vs. mock-inoculated plants at high temperature (HM); (c) AV vs. mock-inoculated plants at ambient temperature (AM); and (d) HM vs. AM are displayed for HTS and qPCR or Northern blot. U6 was used as internal reference for calibrating the expression of miRNAs. The error bars represent the mean (±standard error of the mean) of 4 biological replicates.
![View Image - Figure 8. MicroRNA-mediated regulation involved in capsicum resistance response to CaCV at high temperature (HV) or ambient temperature (AV). Expression patterns of miRNAs (in color) and targets (in black and white) in four pairwise comparisons (AV vs. mock-inoculated plants at ambient temperature (AM); HV vs. mock-inoculated plants at high temperature (HM); HM vs. AM; and HV vs. AV) were analyzed by real-time RT-PCR with the log2-fold change 2–ΔΔCt method or by Northern blot. Real-time RT-PCR of (a) Can-miR408a-3p and its targets: (b) CAN.G1061.9, (c) CAN.G671.1, and (d) CAN.G1305.35. Real-time RT-PCR analysis of (e) Can-miR397-5p and its targets: (f) CAN.G394.71, (g) CAN.G351.2, and (h) CAN.G355.8. Northern blot analysis of (i) Can- miR164b/c-5p, and real-time PCR analysis of its targets: (j) CAN.G394.28 and (k) CAN.G587.8. Real-time RT-PCR analysis of (l) Can-miR168-5p and its target (m) CAN.G637.6. Northern blot analysis was quantified through measuring signal strength using ibright. Actin and U6 were used as internal reference for calibrating the expression of targets and miRNAs, respectively. The error bars represent the mean (±standard error of the mean) of 4 biological replicates. Significant differences between treatments were assessed with Student’s t-test (* p [less than] 0.05; ** p [less than] 0.01).](https://media.proquest.com/media/hms/THUM/16/u6o4a?_a=ChgyMDI1MDcwODIwNDUzNzgzOTo5ODk0NTESBzEzODMzMjMaCk9ORV9TRUFSQ0giDTIxNi43My4yMTYuMzUqBzIwMzIzNTEyCjMxMTA2NDI4NjQ6C0lubGluZUltYWdlQgEwUgZPbmxpbmVaA0lMSWIEVEhVTWoKMjAyNC8wMS8wMXIKMjAyNC8xMi8zMXoAggEgUC0xMDAwMDg0LTIwMTM5NS1TVEFGRi1udWxsLW51bGySAQZPbmxpbmXKAWdNb3ppbGxhLzUuMCBBcHBsZVdlYktpdC81MzcuMzYgKEtIVE1MLCBsaWtlIEdlY2tvOyBjb21wYXRpYmxlOyBDbGF1ZGVCb3QvMS4wOyArY2xhdWRlYm90QGFudGhyb3BpYy5jb20p0gESU2Nob2xhcmx5IEpvdXJuYWxzmgIHUHJlUGFpZKoCG09TOkVNUy1GdWxsVGV4dC1nZXRGdWxsVGV4dMoCD0FydGljbGV8RmVhdHVyZeICAPICAPoCAVmCAwNXZWKKAxxDSUQ6MjAyNTA3MDgyMDQ1Mzc4Mzk6Mzc4MDI3&_s=x%2FHc1YdunXT0K4TzjF99dbVbohU%3D "View Image - Figure 8. MicroRNA-mediated regulation involved in capsicum resistance response to CaCV at high temperature (HV) or ambient temperature (AV). Expression patterns of miRNAs (in color) and targets (in black and white) in four pairwise comparisons (AV vs. mock-inoculated plants at ambient temperature (AM); HV vs. mock-inoculated plants at high temperature (HM); HM vs. AM; and HV vs. AV) were analyzed by real-time RT-PCR with the log2-fold change 2–ΔΔCt method or by Northern blot. Real-time RT-PCR of (a) Can-miR408a-3p and its targets: (b) CAN.G1061.9, (c) CAN.G671.1, and (d) CAN.G1305.35. Real-time RT-PCR analysis of (e) Can-miR397-5p and its targets: (f) CAN.G394.71, (g) CAN.G351.2, and (h) CAN.G355.8. Northern blot analysis of (i) Can- miR164b/c-5p, and real-time PCR analysis of its targets: (j) CAN.G394.28 and (k) CAN.G587.8. Real-time RT-PCR analysis of (l) Can-miR168-5p and its target (m) CAN.G637.6. Northern blot analysis was quantified through measuring signal strength using ibright. Actin and U6 were used as internal reference for calibrating the expression of targets and miRNAs, respectively. The error bars represent the mean (±standard error of the mean) of 4 biological replicates. Significant differences between treatments were assessed with Student’s t-test (* p [less than] 0.05; ** p [less than] 0.01).")
Enlarge this image.Figure 8. MicroRNA-mediated regulation involved in capsicum resistance response to CaCV at high temperature (HV) or ambient temperature (AV). Expression patterns of miRNAs (in color) and targets (in black and white) in four pairwise comparisons (AV vs. mock-inoculated plants at ambient temperature (AM); HV vs. mock-inoculated plants at high temperature (HM); HM vs. AM; and HV vs. AV) were analyzed by real-time RT-PCR with the log2-fold change 2–ΔΔCt method or by Northern blot. Real-time RT-PCR of (a) Can-miR408a-3p and its targets: (b) CAN.G1061.9, (c) CAN.G671.1, and (d) CAN.G1305.35. Real-time RT-PCR analysis of (e) Can-miR397-5p and its targets: (f) CAN.G394.71, (g) CAN.G351.2, and (h) CAN.G355.8. Northern blot analysis of (i) Can- miR164b/c-5p, and real-time PCR analysis of its targets: (j) CAN.G394.28 and (k) CAN.G587.8. Real-time RT-PCR analysis of (l) Can-miR168-5p and its target (m) CAN.G637.6. Northern blot analysis was quantified through measuring signal strength using ibright. Actin and U6 were used as internal reference for calibrating the expression of targets and miRNAs, respectively. The error bars represent the mean (±standard error of the mean) of 4 biological replicates. Significant differences between treatments were assessed with Student’s t-test (* p [less than] 0.05; ** p [less than] 0.01).
Enlarge this image.Figure 9. Validation of miRNA targets by 5′ RLM-RACE. The potential cleavage sites of (a) Can-miR408a-3p on CAN.G1061.9, (b) Can-miR397-5p on CAN.G355.8, (c) Can-miR168-5p on CAN.G637.6 in capsicums were mapped. The arrows indicate the cleavage sites and the numbers indicate clone frequencies. The dashes represent the standard base pairing rules and the dots indicate GU wobble base pairing.
Differentially expressed miRNAs in virus-infected capsicum versus mock-inoculated capsicum grown at high compared to ambient temperature.
| miRNA ID | HV/AV_log2 (FC) | AV/AM_log2 (FC) | HV/HM_log2 (FC) | HM/AM_log2 (FC) |
|---|---|---|---|---|
| Can-miR408a-5p | 4.94651183 | 2.15675264 | ||
| Can-miR408a-3p | 4.52786751 | −1.431929 | 1.87576962 | |
| Can-miR398b-3p | 4.3647467 | −1.1693631 | 2.7875221 | |
| Can-miR397-5p | 3.41052681 | −1.3643359 | 1.80648604 | |
| Can-miR397-3p | 3.09698569 | −1.1034217 | 1.41436391 | |
| Can-miRN492-5p | 2.4800554 | 2.07163552 | ||
| Can-miR_N7b-3p | 2.22044295 | 1.95841442 | ||
| Can-miR_N7c-3p | 2.22044295 | 1.95841442 | ||
| Can-miR_N7d-3p | 2.16627692 | 1.84961374 | ||
| Can-miR393-3p | 2.00053338 | |||
| Can-miRN37-Novel-5p | 1.97325814 | 1.2736888 | ||
| Can-miR_N2a-3p | 1.66675404 | |||
| Can-miR172a-m1-3p | 1.62135413 | |||
| Can-miR172d-3p | 1.61286777 | |||
| Can-miR_N71a-5p | 1.60761055 | |||
| Can-miR_N50a-5p | 1.51467255 | |||
| Can-miR403-3p | 1.29105902 | 1.04987498 | ||
| Can-miRN25-Novel-3p | 1.27226491 | 1.27136116 | ||
| Can-miR160a-m1-5p | 1.25709852 | |||
| Can-miRC13.1-m1-5p | 1.19626839 | |||
| Can-miRC13.2-m1-5p | 1.19626839 | |||
| Can-miR157c-5p | −1.0923326 | |||
| Can-miR157b-5p | −1.1162378 | |||
| Can-miR169h-3p | −1.1385036 | −1.8141374 | ||
| Can-miR_N80a-3p | −1.4595388 | |||
| Can-miR477a-5p | −1.6760649 | −3.4314472 | ||
| Can-miR164b-5p | −1.7200209 | |||
| Can-miR164c-5p | −1.7417164 | |||
| Can-miR477b-5p | −1.7912909 | −2.9432463 | ||
| Can-miRN482-5p | −1.8617146 | −2.2237655 | ||
| Can-miR399a-Known-5p | −1.9889216 | −5.5813609 | ||
| Can-miR399b-Known-5p | −1.9889216 | −5.5813609 | ||
| Can-miR399c-Known-5p | −1.9889216 | −5.5813609 | ||
| Can-miR399d-Known_1-5p | −1.9889216 | −5.5813609 | ||
| Can-miR399d-Known-5p | −1.9889216 | −5.5813609 |
HV/AV_log2 (FC)—log2-fold change when comparing miRNA expression in CaCV-infected capsicum grown at high temperature (HT) with infected plants grown at ambient temperature (AT); AV/AM_log2 (FC)—log2-fold change when comparing miRNA expression in CaCV-infected capsicum grown at AT to that in mock-inoculated capsicum grown at AT; HV/HM_log2 (FC)—log2-fold change when comparing miRNA expression in CaCV-infected capsicum grown at HT to that in mock-inoculated capsicum grown at HT; HM/AM_log2 (FC)—log2-fold change when comparing miRNA expression in mock-inoculated capsicum grown at HT and AT.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Tripodi, P.; Kumar, S. The Capsicum Crop: An Introduction. The Capsicum Genome; Ramchiary, N.; Kole, C. Springer International Publishing: Cham, Switzerland, 2019; pp. 1-8.
2. Zou, X.; Zhu, F. Origin, Evolution and Cultivation History of the Pepper. Acta Hortic. Sin.; 2022; 49, pp. 1371-1381. [DOI: https://dx.doi.org/10.16420/j.issn.0513-353x.2021-0853]
3. Khaitov, B.; Umurzokov, M.; Cho, K.-M.; Lee, Y.-J.; Park, K. Importance and production of chilli pepper; heat tolerance and efficient nutrient use under climate change conditions. Korean J. Agric. Sci.; 2019; 46, pp. 769-779. [DOI: https://dx.doi.org/10.7744/kjoas.20190059]
4. Penella, C.; Calatayud, A. Pepper Crop under Climate Change: Grafting as an Environmental Friendly Strategy. Climate Resilient Agriculture: Strategies and Perspectives; IntechOpen: London, UK, 2018; pp. 128-155.
5. Lu, J.; Guo, M.; Zhai, Y.; Gong, Z.; Lu, M. Differential Responses to the Combined Stress of Heat and Phytophthora capsici Infection Between Resistant and Susceptible Germplasms of Pepper (Capsicum annuum L.). J. Plant Growth Regul.; 2017; 36, pp. 161-173. [DOI: https://dx.doi.org/10.1007/s00344-016-9627-9]
6. Parisi, M.; Alioto, D.; Tripodi, P. Overview of Biotic Stresses in Pepper (Capsicum spp.): Sources of Genetic Resistance, Molecular Breeding and Genomics. Int. J. Mol. Sci.; 2020; 21, 2587. [DOI: https://dx.doi.org/10.3390/ijms21072587] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32276403]
7. Mandal, B.; Jain, R.K.; Krishnareddy, M.; Krishna Kumar, N.K.; Ravi, K.S.; Pappu, H.R. Emerging Problems of Tospoviruses (Bunyaviridae) and their Management in the Indian Subcontinent. Plant Dis.; 2012; 96, pp. 468-479. [DOI: https://dx.doi.org/10.1094/PDIS-06-11-0520]
8. Sharman, M.; Thomas, J.E.; Tree, D.; Persley, D.M. Natural host range and thrips transmission of capsicum chlorosis virus in Australia. Australas. Plant Pathol.; 2020; 49, pp. 45-51. [DOI: https://dx.doi.org/10.1007/s13313-019-00675-7]
9. Persley, D.M.; Thomas, J.E.; Sharman, M. Tospoviruses—An Australian perspective. Australas. Plant Pathol.; 2006; 35, pp. 161-180. [DOI: https://dx.doi.org/10.1071/AP06015]
10. Health, E.; Panel, o.P.; Bragard, C.; Baptista, P.; Chatzivassiliou, E.; Gonthier, P.; Jaques Miret, J.A.; Justesen, A.F.; MacLeod, A.; Magnusson, C.S. et al. Pest categorisation of Capsicum chlorosis virus. Eur. Food Saf. Auth. J.; 2022; 20, e07337. [DOI: https://dx.doi.org/10.2903/j.efsa.2022.7337]
11. Chen, K.; Xu, Z.; Yan, L.; Wang, G. Characterization of a New Strain of Capsicum chlorosis virus from Peanut (Arachis hypogaea L.) in China. J. Phytopathol.; 2007; 155, pp. 178-181. [DOI: https://dx.doi.org/10.1111/j.1439-0434.2007.01217.x]
12. Krishnareddy, M.; Rani, R.U.; Kumar, K.S.A.; Reddy, K.M.; Pappu, H.R. Capsicum chlorosis virus (Genus Tospovirus) Infecting Chili Pepper (Capsicum annuum) in India. Plant. Dis.; 2008; 92, 1469. [DOI: https://dx.doi.org/10.1094/PDIS-92-10-1469B]
13. Widana Gamage, S.M.K.; McGrath, D.J.; Persley, D.M.; Dietzgen, R.G. Transcriptome Analysis of Capsicum Chlorosis Virus-Induced Hypersensitive Resistance Response in Bell Capsicum. PLoS ONE; 2016; 11, e0159085. [DOI: https://dx.doi.org/10.1371/journal.pone.0159085] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27398596]
14. McGarth, D. Capsicum Breeding for Tospovirus Resistance; Horticultural Australia Ltd.: Sydney, NSW, Australia, 2006.
15. Abudurexiti, A.; Adkins, S.; Alioto, D.; Alkhovsky, S.V.; Avšič-Županc, T.; Ballinger, M.J.; Bente, D.A.; Beer, M.; Bergeron, É.; Blair, C.D. et al. Taxonomy of the order Bunyavirales: Update 2019. Arch. Virol.; 2019; 164, pp. 1949-1965. [DOI: https://dx.doi.org/10.1007/s00705-019-04253-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31065850]
16. Whitfield, A.E.; Ullman, D.E.; German, T.L. Tospovirus-thrips interactions. Annu. Rev. Phytopathol.; 2005; 43, pp. 459-489. [DOI: https://dx.doi.org/10.1146/annurev.phyto.43.040204.140017]
17. Widana Gamage, S.; Persley, D.M.; Higgins, C.M.; Dietzgen, R.G. First complete genome sequence of a capsicum chlorosis tospovirus isolate from Australia with an unusually large S RNA intergenic region. Arch. Virol.; 2015; 160, pp. 869-872. [DOI: https://dx.doi.org/10.1007/s00705-014-2324-8]
18. Honjo, M.N.; Emura, N.; Kawagoe, T.; Sugisaka, J.; Kamitani, M.; Nagano, A.J.; Kudoh, H. Seasonality of interactions between a plant virus and its host during persistent infection in a natural environment. Int. Soc. Microb. Ecol. J.; 2020; 14, pp. 506-518. [DOI: https://dx.doi.org/10.1038/s41396-019-0519-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31664159]
19. Makarova, S.; Makhotenko, A.; Spechenkova, N.; Love, A.J.; Kalinina, N.O.; Taliansky, M. Interactive responses of potato (Solanum tuberosum L.) plants to heat stress and infection with potato virus Y. Front. Microbiol.; 2018; 9, 2582. [DOI: https://dx.doi.org/10.3389/fmicb.2018.02582]
20. Wang, Y.; Bao, Z.; Zhu, Y.; Hua, J. Analysis of temperature modulation of plant defense against biotrophic microbes. Mol. Plant Microbe Interact.; 2009; 22, pp. 498-506. [DOI: https://dx.doi.org/10.1094/MPMI-22-5-0498]
21. Anfoka, G.; Moshe, A.; Fridman, L.; Amrani, L.; Rotem, O.; Kolot, M.; Zeidan, M.; Czosnek, H.; Gorovits, R. Tomato yellow leaf curl virus infection mitigates the heat stress response of plants grown at high temperatures. Sci. Rep.; 2016; 6, 19715. [DOI: https://dx.doi.org/10.1038/srep19715]
22. Canto, T.; Aranda, M.A.; Fereres, A. Climate change effects on physiology and population processes of hosts and vectors that influence the spread of hemipteran-borne plant viruses. Glob. Chang. Biol.; 2009; 15, pp. 1884-1894. [DOI: https://dx.doi.org/10.1111/j.1365-2486.2008.01820.x]
23. Velasquez, A.C.; Castroverde, C.D.M.; He, S.Y. Plant-Pathogen Warfare under Changing Climate Conditions. Curr. Biol.; 2018; 28, pp. R619-R634. [DOI: https://dx.doi.org/10.1016/j.cub.2018.03.054]
24. Tsai, W.-A.; Brosnan, C.A.; Mitter, N.; Dietzgen, R.G. Perspectives on plant virus diseases in a climate change scenario of elevated temperatures. Stress Biol.; 2022; 2, 37. [DOI: https://dx.doi.org/10.1007/s44154-022-00058-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37676437]
25. Gouveia, B.C.; Calil, I.P.; Machado, J.P.B.; Santos, A.A.; Fontes, E.P.B. Immune Receptors and Co-receptors in Antiviral Innate Immunity in Plants. Front. Microbiol; 2017; 7, 2139. [DOI: https://dx.doi.org/10.3389/fmicb.2016.02139]
26. de Ronde, D.; Lohuis, D.; Kormelink, R. Identification and characterization of a new class of Tomato spotted wilt virus isolates that break Tsw-based resistance in a temperature-dependent manner. Plant Pathol.; 2019; 68, pp. 60-71. [DOI: https://dx.doi.org/10.1111/ppa.12952]
27. Moury, B.; Selassie, K.G.; Marchoux, G.; Daubèze, A.-M.; Palloix, A. High temperature effects on hypersensitive resistance to Tomato Spotted Wilt Tospovirus (TSWV) in pepper (Capsicum chinense Jacq.). Eur. J. Plant Pathol.; 1998; 104, pp. 489-498. [DOI: https://dx.doi.org/10.1023/A:1008618022144]
28. Soler, S.; Diez, M.J.; Nuez, F. Effect of Temperature Regime and Growth Stage Interaction on Pattern of Virus Presence in TSWV-Resistant Accessions of Capsicum chinense. Plant Dis.; 1998; 82, pp. 1199-1204. [DOI: https://dx.doi.org/10.1094/PDIS.1998.82.11.1199]
29. Szajko, K.; Yin, Z.; Marczewski, W. Accumulation of miRNA and mRNA Targets in Potato Leaves Displaying Temperature-Dependent Responses to Potato Virus Y. Potato Res.; 2019; 62, pp. 379-392. [DOI: https://dx.doi.org/10.1007/s11540-019-9417-4]
30. Dexheimer, P.J.; Cochella, L. MicroRNAs: From Mechanism to Organism. Front. Cell Dev. Biol.; 2020; 8, 409. [DOI: https://dx.doi.org/10.3389/fcell.2020.00409]
31. Liu, S.-R.; Zhou, J.-J.; Hu, C.-G.; Wei, C.-L.; Zhang, J.-Z. MicroRNA-Mediated Gene Silencing in Plant Defense and Viral Counter-Defense. Front. Microbiol.; 2017; 8, 1801. [DOI: https://dx.doi.org/10.3389/fmicb.2017.01801]
32. D’Ario, M.; Griffiths-Jones, S.; Kim, M. Small RNAs: Big Impact on Plant Development. Trends. Plant Sci.; 2017; 22, pp. 1056-1068. [DOI: https://dx.doi.org/10.1016/j.tplants.2017.09.009]
33. Borges, F.; Martienssen, R.A. The expanding world of small RNAs in plants. Nat. Rev. Mol. Cell Biol.; 2015; 16, pp. 727-741. [DOI: https://dx.doi.org/10.1038/nrm4085]
34. Zhang, B.; Li, W.; Zhang, J.; Wang, L.; Wu, J. Roles of Small RNAs in Virus-Plant Interactions. Viruses; 2019; 11, 827. [DOI: https://dx.doi.org/10.3390/v11090827]
35. Križnik, M.; Baebler, Š.; Gruden, K. Roles of small RNAs in the establishment of tolerant interaction between plants and viruses. Curr. Opin. Virol.; 2020; 42, pp. 25-31. [DOI: https://dx.doi.org/10.1016/j.coviro.2020.04.006] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32480352]
36. Shivaprasad, P.V.; Chen, H.-M.; Patel, K.; Bond, D.M.; Santos, B.A.C.M.; Baulcombe, D.C. A MicroRNA Superfamily Regulates Nucleotide Binding Site–Leucine-Rich Repeats and Other mRNAs. Plant Cell; 2012; 24, pp. 859-874. [DOI: https://dx.doi.org/10.1105/tpc.111.095380]
37. Li, F.; Pignatta, D.; Bendix, C.; Brunkard, J.O.; Cohn, M.M.; Tung, J.; Sun, H.; Kumar, P.; Baker, B. MicroRNA regulation of plant innate immune receptors. Proc. Natl. Acad. Sci. USA; 2012; 109, pp. 1790-1795. [DOI: https://dx.doi.org/10.1073/pnas.1118282109] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22307647]
38. Sharma, N.; Sahu, P.P.; Prasad, A.; Muthamilarasan, M.; Waseem, M.; Khan, Y.; Thakur, J.K.; Chakraborty, S.; Prasad, M. The Sw5a gene confers resistance to ToLCNDV and triggers an HR response after direct AC4 effector recognition. Proc. Natl. Acad. Sci. USA; 2021; 118, e2101833118. [DOI: https://dx.doi.org/10.1073/pnas.2101833118] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34385303]
39. Zhang, C.; Ding, Z.; Wu, K.; Yang, L.; Li, Y.; Yang, Z.; Shi, S.; Liu, X.; Zhao, S.; Yang, Z. et al. Suppression of Jasmonic Acid-Mediated Defense by Viral-Inducible MicroRNA319 Facilitates Virus Infection in Rice. Mol. Plant; 2016; 9, pp. 1302-1314. [DOI: https://dx.doi.org/10.1016/j.molp.2016.06.014]
40. Yang, J.; Zhang, F.; Li, J.; Chen, J.P.; Zhang, H.M. Integrative Analysis of the microRNAome and Transcriptome Illuminates the Response of Susceptible Rice Plants to Rice Stripe Virus. PLoS ONE; 2016; 11, e0146946. [DOI: https://dx.doi.org/10.1371/journal.pone.0146946]
41. Tao, H.; Jia, Z.; Gao, X.; Gui, M.; Li, Y.; Liu, Y. Analysis of the miRNA expression profile involved in the tomato spotted wilt orthotospovirus–pepper interaction. Virus. Res.; 2022; 312, 198710. [DOI: https://dx.doi.org/10.1016/j.virusres.2022.198710]
42. Zhang, T.; Zhai, J.; Zhang, X.; Ling, L.; Li, M.; Xie, S.; Song, M.; Ma, C. Interactive Web-based Annotation of Plant MicroRNAs with iwa-miRNA. Genom. Proteom. Bioinform.; 2022; 20, pp. 557-567. [DOI: https://dx.doi.org/10.1016/j.gpb.2021.02.010]
43. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome. Biol.; 2009; 10, R25. [DOI: https://dx.doi.org/10.1186/gb-2009-10-3-r25]
44. Kuang, Z.; Wang, Y.; Li, L.; Yang, X. miRDeep-P2: Accurate and fast analysis of the microRNA transcriptome in plants. Bioinformatics; 2019; 35, pp. 2521-2522. [DOI: https://dx.doi.org/10.1093/bioinformatics/bty972]
45. Guo, Z.; Kuang, Z.; Wang, Y.; Zhao, Y.; Tao, Y.; Cheng, C.; Yang, J.; Lu, X.; Hao, C.; Wang, T. et al. PmiREN: A comprehensive encyclopedia of plant miRNAs. Nucleic. Acids. Res.; 2019; 48, pp. D1114-D1121. [DOI: https://dx.doi.org/10.1093/nar/gkz894] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31602478]
46. Chen, C.; Li, J.; Feng, J.; Liu, B.; Feng, L.; Yu, X.; Li, G.; Zhai, J.; Meyers, B.C.; Xia, R. sRNAanno—A database repository of uniformly annotated small RNAs in plants. Hortic. Res.; 2021; 8, 45. [DOI: https://dx.doi.org/10.1038/s41438-021-00480-8]
47. Lunardon, A.; Johnson, N.R.; Hagerott, E.; Phifer, T.; Polydore, S.; Coruh, C.; Axtell, M.J. Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants. Genome Res.; 2020; 30, pp. 497-513. [DOI: https://dx.doi.org/10.1101/gr.256750.119]
48. Meng, J.; Liu, D.; Sun, C.; Luan, Y. Prediction of plant pre-microRNAs and their microRNAs in genome-scale sequences using structure-sequence features and support vector machine. BMC Bioinform.; 2014; 15, 423. [DOI: https://dx.doi.org/10.1186/s12859-014-0423-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25547126]
49. Axtell, M.J.; Meyers, B.C. Revisiting Criteria for Plant MicroRNA Annotation in the Era of Big Data. Plant Cell; 2018; 30, pp. 272-284. [DOI: https://dx.doi.org/10.1105/tpc.17.00851] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29343505]
50. Heberle, H.; Meirelles, G.V.; da Silva, F.R.; Telles, G.P.; Minghim, R. InteractiVenn: A web-based tool for the analysis of sets through Venn diagrams. BMC Bioinform.; 2015; 16, 169. [DOI: https://dx.doi.org/10.1186/s12859-015-0611-3]
51. Friedländer, M.R.; Mackowiak, S.D.; Li, N.; Chen, W.; Rajewsky, N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic. Acids Res.; 2012; 40, pp. 37-52. [DOI: https://dx.doi.org/10.1093/nar/gkr688]
52. Fletcher, S.J.; Boden, M.; Mitter, N.; Carroll, B.J. SCRAM: A pipeline for fast index-free small RNA read alignment and visualization. Bioinformatics; 2018; 34, pp. 2670-2672. [DOI: https://dx.doi.org/10.1093/bioinformatics/bty161]
53. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.; 2014; 15, 550. [DOI: https://dx.doi.org/10.1186/s13059-014-0550-8]
54. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B; 1995; 57, pp. 289-300. [DOI: https://dx.doi.org/10.1111/j.2517-6161.1995.tb02031.x]
55. Bedre, R. Reneshbedre/Bioinfokit: Bioinformatics Data Analysis and Visualization Toolkit. Zenodo; 2020; [DOI: https://dx.doi.org/10.5281/zenodo.3747737]
56. Dai, X.; Zhuang, Z.; Zhao, P.X. psRNATarget: A plant small RNA target analysis server (2017 release). Nucleic. Acids Res.; 2018; 46, pp. W49-W54. [DOI: https://dx.doi.org/10.1093/nar/gky316]
57. Van Bel, M.; Diels, T.; Vancaester, E.; Kreft, L.; Botzki, A.; Van de Peer, Y.; Coppens, F.; Vandepoele, K. PLAZA 4.0: An integrative resource for functional, evolutionary and comparative plant genomics. Nucleic. Acids Res.; 2017; 46, pp. D1190-D1196. [DOI: https://dx.doi.org/10.1093/nar/gkx1002]
58. Tian, T.; Liu, Y.; Yan, H.; You, Q.; Yi, X.; Du, Z.; Xu, W.; Su, Z. agriGO v2.0: A GO analysis toolkit for the agricultural community, 2017 update. Nucleic. Acids Res.; 2017; 45, pp. W122-W129. [DOI: https://dx.doi.org/10.1093/nar/gkx382] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28472432]
59. Hong, G.; Zhang, W.; Li, H.; Shen, X.; Guo, Z. Separate enrichment analysis of pathways for up- and downregulated genes. J. R. Soc. Interface; 2014; 11, 20130950. [DOI: https://dx.doi.org/10.1098/rsif.2013.0950]
60. Benjamini, Y.; Yekutieli, D. The Control of the False Discovery Rate in Multiple Testing Under Dependency. Ann. Stat.; 2001; 29, pp. 1165-1188. [DOI: https://dx.doi.org/10.1214/aos/1013699998]
61. Kang, K.; Zhang, X.; Liu, H.; Wang, Z.; Zhong, J.; Huang, Z.; Peng, X.; Zeng, Y.; Wang, Y.; Yang, Y. et al. A novel real-time PCR assay of microRNAs using S-Poly(T), a specific oligo(dT) reverse transcription primer with excellent sensitivity and specificity. PLoS ONE; 2012; 7, e48536. [DOI: https://dx.doi.org/10.1371/journal.pone.0048536]
62. Niu, Y.; Zhang, L.; Qiu, H.; Wu, Y.; Wang, Z.; Zai, Y.; Liu, L.; Qu, J.; Kang, K.; Gou, D. An improved method for detecting circulating microRNAs with S-Poly(T) Plus real-time PCR. Sci. Rep.; 2015; 5, 15100. [DOI: https://dx.doi.org/10.1038/srep15100]
63. Xie, S.; Zhu, Q.; Qu, W.; Xu, Z.; Liu, X.; Li, X.; Li, S.; Ma, W.; Miao, Y.; Zhang, L. et al. sRNAPrimerDB: Comprehensive primer design and search web service for small non-coding RNAs. Bioinformatics; 2019; 35, pp. 1566-1572. [DOI: https://dx.doi.org/10.1093/bioinformatics/bty852]
64. Tang, Z.; Yu, J.; Xie, J.; Lyu, J.; Feng, Z.; Dawuda, M.M.; Liao, W.; Wu, Y.; Hu, L. Physiological and Growth Response of Pepper (Capsicum annum L.) Seedlings to Supplementary Red/Blue Light Revealed through Transcriptomic Analysis. Agronomy; 2019; 9, 139. [DOI: https://dx.doi.org/10.3390/agronomy9030139]
65. Tsai, W.A.; Shafiei-Peters, J.R.; Mitter, N.; Dietzgen, R.G. Effects of Elevated Temperature on the Susceptibility of Capsicum Plants to Capsicum Chlorosis Virus Infection. Pathogens; 2022; 11, 200. [DOI: https://dx.doi.org/10.3390/pathogens11020200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35215143]
66. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods; 2001; 25, pp. 402-408. [DOI: https://dx.doi.org/10.1006/meth.2001.1262]
67. Mitter, N.; Dietzgen, R.G. Use of hairpin RNA constructs for engineering plant virus resistance. Methods Mol. Biol.; 2012; 894, pp. 191-208. [DOI: https://dx.doi.org/10.1007/978-1-61779-882-5_13] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22678581]
68. Kasschau, K.D.; Xie, Z.; Allen, E.; Llave, C.; Chapman, E.J.; Krizan, K.A.; Carrington, J.C. P1/HC-Pro, a viral suppressor of RNA silencing, interferes with Arabidopsis development and miRNA function. Dev. Cell; 2003; 4, pp. 205-217. [DOI: https://dx.doi.org/10.1016/S1534-5807(03)00025-X]
69. Song, C.; Wang, C.; Zhang, C.; Korir, N.K.; Yu, H.; Ma, Z.; Fang, J. Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). BioMed Cent. Genom.; 2010; 11, 431. [DOI: https://dx.doi.org/10.1186/1471-2164-11-431]
70. Mallory, A.C.; Elmayan, T.; Vaucheret, H. MicroRNA maturation and action--the expanding roles of ARGONAUTEs. Curr. Opin. Plant Biol.; 2008; 11, pp. 560-566. [DOI: https://dx.doi.org/10.1016/j.pbi.2008.06.008]
71. Varallyay, E.; Valoczi, A.; Agyi, A.; Burgyan, J.; Havelda, Z. Plant virus-mediated induction of miR168 is associated with repression of ARGONAUTE1 accumulation. Eur. Mol. Biol. Organ. J.; 2010; 29, pp. 3507-3519. [DOI: https://dx.doi.org/10.1038/emboj.2010.215]
72. Varallyay, E.; Havelda, Z. Unrelated viral suppressors of RNA silencing mediate the control of ARGONAUTE1 level. Mol. Plant Pathol.; 2013; 14, pp. 567-575. [DOI: https://dx.doi.org/10.1111/mpp.12029]
73. Fang, X.; Qi, Y. RNAi in Plants: An Argonaute-Centered View. Plant Cell; 2016; 28, pp. 272-285. [DOI: https://dx.doi.org/10.1105/tpc.15.00920]
74. Fei, Y.; Mao, Y.; Shen, C.; Wang, R.; Zhang, H.; Huang, J. WPMIAS: Whole-degradome-based Plant MicroRNA-Target Interaction Analysis Server. Bioinformatics; 2019; 36, pp. 1937-1939. [DOI: https://dx.doi.org/10.1093/bioinformatics/btz820]
75. Zhang, L.; Qin, C.; Mei, J.; Chen, X.; Wu, Z.; Luo, X.; Cheng, J.; Tang, X.; Hu, K.; Li, S.C. Identification of MicroRNA Targets of Capsicum spp. Using MiRTrans-a Trans-Omics Approach. Front. Plant Sci.; 2017; 8, 495. [DOI: https://dx.doi.org/10.3389/fpls.2017.00495]
76. Marques de Carvalho, L.; Benda, N.D.; Vaughan, M.M.; Cabrera, A.R.; Hung, K.; Cox, T.; Abdo, Z.; Allen, L.H.; Teal, P.E. Mi-1-Mediated Nematode Resistance in Tomatoes is Broken by Short-Term Heat Stress but Recovers Over Time. J. Nematol.; 2015; 47, pp. 133-140. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26170475]
77. Romero, A.M.; Kousik, C.S.; Ritchie, D.F. Temperature Sensitivity of the Hypersensitive Response of Bell Pepper to Xanthomonas axonopodis pv. vesicatoria. Phytopathology; 2002; 92, pp. 197-203. [DOI: https://dx.doi.org/10.1094/PHYTO.2002.92.2.197] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18943094]
78. Whitham, S.; McCormick, S.; Baker, B. The N gene of tobacco confers resistance to tobacco mosaic virus in transgenic tomato. Proc. Natl. Acad. Sci. USA; 1996; 93, pp. 8776-8781. [DOI: https://dx.doi.org/10.1073/pnas.93.16.8776] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8710948]
79. Chung, B.N.; Lee, J.H.; Kang, B.C.; Koh, S.W.; Joa, J.H.; Choi, K.S.; Ahn, J.J. HR-Mediated Defense Response is Overcome at High Temperatures in Capsicum Species. Plant Pathol. J.; 2018; 34, pp. 71-77. [DOI: https://dx.doi.org/10.5423/PPJ.NT.06.2017.0120] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29422790]
80. Hwang, D.G.; Park, J.H.; Lim, J.Y.; Kim, D.; Choi, Y.; Kim, S.; Reeves, G.; Yeom, S.I.; Lee, J.S.; Park, M. et al. The hot pepper (Capsicum annuum) microRNA transcriptome reveals novel and conserved targets: A foundation for understanding MicroRNA functional roles in hot pepper. PLoS ONE; 2013; 8, e64238. [DOI: https://dx.doi.org/10.1371/journal.pone.0064238]
81. Seo, E.; Kim, T.; Park, J.H.; Yeom, S.-I.; Kim, S.; Seo, M.-K.; Shin, C.; Choi, D. Genome-wide comparative analysis in Solanaceous species reveals evolution of microRNAs targeting defense genes in Capsicum spp. DNA Res.; 2018; 25, pp. 561-575. [DOI: https://dx.doi.org/10.1093/dnares/dsy025]
82. Xu, X.-W.; Li, T.; Li, Y.; Li, Z.-X. Identification and Analysis of C. annuum microRNAs by High-Throughput Sequencing and Their Association with High Temperature and High Air Humidity Stress. Int. J. Bioautomat.; 2015; 19, pp. 459-472.
83. Yang, S.; Zhang, Z.; Chen, W.; Li, X.; Zhou, S.; Liang, C.; Li, X.; Yang, B.; Zou, X.; Liu, F. et al. Integration of mRNA and miRNA profiling reveals the heterosis of three hybrid combinations of Capsicum annuum varieties. GM Crops Food; 2021; 12, pp. 224-241. [DOI: https://dx.doi.org/10.1080/21645698.2020.1852064]
84. Apostolova, E.; Hadjieva, N.; Ivanova, D.P.; Yahubyan, G.; Baev, V.; Gozmanova, M. MicroRNA expression dynamics reshape the cultivar-specific response of pepper (Capsicum annuum L.) to potato spindle tuber viroid (PSTVd) infection. Sci. Hortic.; 2021; 278, 109845. [DOI: https://dx.doi.org/10.1016/j.scienta.2020.109845]
85. Lee, M.H.; Jeon, H.S.; Kim, S.H.; Chung, J.H.; Roppolo, D.; Lee, H.J.; Cho, H.J.; Tobimatsu, Y.; Ralph, J.; Park, O.K. Lignin-based barrier restricts pathogens to the infection site and confers resistance in plants. Eur. Mol. Biol. Organ. J; 2019; 38, e101948. [DOI: https://dx.doi.org/10.15252/embj.2019101948] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31559647]
86. Kimmins, W.C.; Wuddah, D. Hypersensitive Resistance: Determination of Lignin in Leaves with a Localized Virus Infection. Phytopathology; 1977; 67, pp. 1012-1016. [DOI: https://dx.doi.org/10.1094/Phyto-67-1012]
87. Li, C.; Liu, X.-X.; Abouelnasr, H.; Mohamed Hamed, A.; Kou, M.; Tang, W.; Yan, H.; Wang, X.; Wang, X.-X.; Zhang, Y.-G. et al. Inhibition of miR397 by STTM technology to increase sweetpotato resistance to SPVD. J. Integr. Agric.; 2022; 21, pp. 2865-2875. [DOI: https://dx.doi.org/10.1016/j.jia.2022.07.054]
88. Ma, C.; Burd, S.; Lers, A. miR408 is involved in abiotic stress responses in Arabidopsis. Plant J.; 2015; 84, pp. 169-187. [DOI: https://dx.doi.org/10.1111/tpj.12999]
89. Moxon, S.; Jing, R.; Szittya, G.; Schwach, F.; Rusholme Pilcher, R.L.; Moulton, V.; Dalmay, T. Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. Genome Res.; 2008; 18, pp. 1602-1609. [DOI: https://dx.doi.org/10.1101/gr.080127.108]
90. Mohorianu, I.; Schwach, F.; Jing, R.; Lopez-Gomollon, S.; Moxon, S.; Szittya, G.; Sorefan, K.; Moulton, V.; Dalmay, T. Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns. Plant J.; 2011; 67, pp. 232-246. [DOI: https://dx.doi.org/10.1111/j.1365-313X.2011.04586.x]
91. Linsen, S.E.; de Wit, E.; Janssens, G.; Heater, S.; Chapman, L.; Parkin, R.K.; Fritz, B.; Wyman, S.K.; de Bruijn, E.; Voest, E.E. et al. Limitations and possibilities of small RNA digital gene expression profiling. Nat. Methods; 2009; 6, pp. 474-476. [DOI: https://dx.doi.org/10.1038/nmeth0709-474]
92. Fuchs, R.T.; Sun, Z.; Zhuang, F.; Robb, G.B. Bias in ligation-based small RNA sequencing library construction is determined by adaptor and RNA structure. PLoS ONE; 2015; 10, e0126049. [DOI: https://dx.doi.org/10.1371/journal.pone.0126049] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25942392]
93. Dabney, J.; Meyer, M. Length and GC-biases during sequencing library amplification: A comparison of various polymerase-buffer systems with ancient and modern DNA sequencing libraries. BioTechniques; 2012; 52, pp. 87-94. [DOI: https://dx.doi.org/10.2144/000113809]
94. Abdel-Ghany, S.E.; Pilon, M. MicroRNA-mediated systemic down-regulation of copper protein expression in response to low copper availability in Arabidopsis. J. Biol. Chem.; 2008; 283, pp. 15932-15945. [DOI: https://dx.doi.org/10.1074/jbc.M801406200]
95. Trindade, I.; Capitão, C.; Dalmay, T.; Fevereiro, M.P.; Santos, D.M. miR398 and miR408 are up-regulated in response to water deficit in Medicago truncatula. Planta; 2010; 231, pp. 705-716. [DOI: https://dx.doi.org/10.1007/s00425-009-1078-0]
96. Li, C.; Li, D.; Zhou, H.; Li, J.; Lu, S. Analysis of the laccase gene family and miR397-/miR408-mediated posttranscriptional regulation in Salvia miltiorrhiza. PeerJ; 2019; 7, e7605. [DOI: https://dx.doi.org/10.7717/peerj.7605]
97. Lu, S.; Li, Q.; Wei, H.; Chang, M.J.; Tunlaya-Anukit, S.; Kim, H.; Liu, J.; Song, J.; Sun, Y.H.; Yuan, L. et al. Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa. Proc. Natl. Acad. Sci. USA; 2013; 110, pp. 10848-10853. [DOI: https://dx.doi.org/10.1073/pnas.1308936110] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23754401]
98. Wang, C.Y.; Zhang, S.; Yu, Y.; Luo, Y.C.; Liu, Q.; Ju, C.; Zhang, Y.C.; Qu, L.H.; Lucas, W.J.; Wang, X. et al. MiR397b regulates both lignin content and seed number in Arabidopsis via modulating a laccase involved in lignin biosynthesis. Plant Biotechnol. J.; 2014; 12, pp. 1132-1142. [DOI: https://dx.doi.org/10.1111/pbi.12222] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24975689]
99. Swetha, C.; Basu, D.; Pachamuthu, K.; Tirumalai, V.; Nair, A.; Prasad, M.; Shivaprasad, P.V. Major Domestication-Related Phenotypes in Indica Rice Are Due to Loss of miRNA-Mediated Laccase Silencing. Plant Cell; 2018; 30, pp. 2649-2662. [DOI: https://dx.doi.org/10.1105/tpc.18.00472] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30341147]
100. Huang, S.; Zhou, J.; Gao, L.; Tang, Y. Plant miR397 and its functions. Funct. Plant Biol.; 2021; 48, pp. 361-370. [DOI: https://dx.doi.org/10.1071/FP20342] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33333000]
101. Zhao, Q.; Tobimatsu, Y.; Zhou, R.; Pattathil, S.; Gallego-Giraldo, L.; Fu, C.; Jackson, L.A.; Hahn, M.G.; Kim, H.; Chen, F. et al. Loss of function of cinnamyl alcohol dehydrogenase 1 leads to unconventional lignin and a temperature-sensitive growth defect in Medicago truncatula. Proc. Natl. Acad. Sci. USA; 2013; 110, pp. 13660-13665. [DOI: https://dx.doi.org/10.1073/pnas.1312234110]
102. Liu, Q.; Luo, L.; Zheng, L. Lignins: Biosynthesis and Biological Functions in Plants. Int. J. Mol. Sci.; 2018; 19, 335. [DOI: https://dx.doi.org/10.3390/ijms19020335]
103. Rosas-Diaz, T.; Zhang, D.; Fan, P.; Wang, L.; Ding, X.; Jiang, Y.; Jimenez-Gongora, T.; Medina-Puche, L.; Zhao, X.; Feng, Z. et al. A virus-targeted plant receptor-like kinase promotes cell-to-cell spread of RNAi. Proc. Natl. Acad. Sci. USA; 2018; 115, pp. 1388-1393. [DOI: https://dx.doi.org/10.1073/pnas.1715556115]
104. Swida-Barteczka, A.; Szweykowska-Kulinska, Z. Micromanagement of Developmental and Stress-Induced Senescence: The Emerging Role of MicroRNAs. Genes; 2019; 10, 210. [DOI: https://dx.doi.org/10.3390/genes10030210] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30871088]
105. Lopez-Gomollon, S.; Mohorianu, I.; Szittya, G.; Moulton, V.; Dalmay, T. Diverse correlation patterns between microRNAs and their targets during tomato fruit development indicates different modes of microRNA actions. Planta; 2012; 236, pp. 1875-1887. [DOI: https://dx.doi.org/10.1007/s00425-012-1734-7]
106. Kidner, C.A.; Martienssen, R.A. Spatially restricted microRNA directs leaf polarity through ARGONAUTE1. Nature; 2004; 428, pp. 81-84. [DOI: https://dx.doi.org/10.1038/nature02366] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14999284]
107. Kawashima, C.G.; Yoshimoto, N.; Maruyama-Nakashita, A.; Tsuchiya, Y.N.; Saito, K.; Takahashi, H.; Dalmay, T. Sulphur starvation induces the expression of microRNA-395 and one of its target genes but in different cell types. Plant J.; 2009; 57, pp. 313-321. [DOI: https://dx.doi.org/10.1111/j.1365-313X.2008.03690.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18801012]
108. Singh, R.; Singh, S.; Parihar, P.; Mishra, R.K.; Tripathi, D.K.; Singh, V.P.; Chauhan, D.K.; Prasad, S.M. Reactive Oxygen Species (ROS): Beneficial Companions of Plants’ Developmental Processes. Front. Plant Sci.; 2016; 7, 1299. [DOI: https://dx.doi.org/10.3389/fpls.2016.01299] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27729914]
109. Pérez-Llorca, M.; Munné-Bosch, S. Aging, stress, and senescence in plants: What can biological diversity teach us?. Geroscience; 2021; 43, pp. 167-180. [DOI: https://dx.doi.org/10.1007/s11357-021-00336-y]
110. Singh, A.; Permar, V.; Basavaraj, A.; Bhoopal, S.T.; Praveen, S. Effect of Temperature on Symptoms Expression and Viral RNA Accumulation in Groundnut Bud Necrosis Virus Infected Vigna unguiculata. Iran. J. Biotechnol.; 2018; 16, pp. 227-234. [DOI: https://dx.doi.org/10.21859/ijb.1846]
111. Li, Z.; Peng, J.; Wen, X.; Guo, H. Ethylene-insensitive3 is a senescence-associated gene that accelerates age-dependent leaf senescence by directly repressing miR164 transcription in Arabidopsis. Plant Cell; 2013; 25, pp. 3311-3328. [DOI: https://dx.doi.org/10.1105/tpc.113.113340]
112. Li, T.; Gonzalez, N.; Inzé, D.; Dubois, M. Emerging Connections between Small RNAs and Phytohormones. Trends. Plant Sci.; 2020; 25, pp. 912-929. [DOI: https://dx.doi.org/10.1016/j.tplants.2020.04.004]
113. Podzimska-Sroka, D.; O’Shea, C.; Gregersen, P.L.; Skriver, K. NAC Transcription Factors in Senescence: From Molecular Structure to Function in Crops. Plants; 2015; 4, pp. 412-448. [DOI: https://dx.doi.org/10.3390/plants4030412]
114. Dalmadi, Á.; Miloro, F.; Bálint, J.; Várallyay, É.; Havelda, Z. Controlled RISC loading efficiency of miR168 defined by miRNA duplex structure adjusts ARGONAUTE1 homeostasis. Nucleic. Acids. Res.; 2021; 49, pp. 12912-12928. [DOI: https://dx.doi.org/10.1093/nar/gkab1138] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34850097]
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Brosnan, Christopher A; Mitter, Neena
; Dietzgen, Ralf G