<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2024.1408926</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>The role of small RNAs in resistant melon cultivar against <italic>Phelipanche aegyptiaca</italic> parasitization</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes">
<name><surname>Mao</surname> <given-names>Jian-Cai</given-names></name>
<xref ref-type="author-notes" rid="fn0001"><sup>&#x2020;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2327013/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name><surname>Yan</surname> <given-names>Miao</given-names></name>
<xref ref-type="author-notes" rid="fn0001"><sup>&#x2020;</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Li</surname> <given-names>Jun-Hua</given-names></name>
<uri xlink:href="https://loop.frontiersin.org/people/2330760/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Yang</surname> <given-names>Jun-Yan</given-names></name>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Wang</surname> <given-names>Hao-Jie</given-names></name>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2330776/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff><institution>Hami Melon Research Center, Xinjiang Academy of Agricultural Sciences</institution>, <addr-line>Urumqi</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by" id="fn0002">
<p>Edited by: Hongju Zhu, Chinese Academy of Agricultural Sciences, China</p>
</fn>
<fn fn-type="edited-by" id="fn0003">
<p>Reviewed by: Tao Zhang, Chinese Academy of Sciences (CAS), China</p>
<p>Xingwang Liu, China Agricultural University, China</p>
<p>Shakil Ahmad, Hainan University, China</p>
</fn>
<corresp id="c001">&#x002A;Correspondence: Hao-Jie Wang, <email>wanghaojie@xaas.ac.cn</email></corresp>
<fn fn-type="equal" id="fn0001">
<p><sup>&#x2020;</sup>These authors have contributed equally to this work</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>07</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>15</volume>
<elocation-id>1408926</elocation-id>
<history>
<date date-type="received">
<day>29</day>
<month>03</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>22</day>
<month>04</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2024 Mao, Yan, Li, Yang and Wang.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Mao, Yan, Li, Yang and Wang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<p>Bidirectional trans-kingdom RNA silencing, a pivotal factor in plant-pathogen interactions, remains less explored in plant host-parasite dynamics. Here, using small RNA sequencing in melon root systems, we investigated microRNA (miRNA) expression variation in resistant and susceptible cultivars pre-and post-infection by the parasitic plant, broomrape. This approach revealed 979 known miRNAs and 110 novel miRNAs across 110 families. When comparing susceptible (F0) and resistant (R0) melon lines with broomrape infection (F25 and R25), 39 significantly differentially expressed miRNAs were observed in F25 vs. F0, 35 in R25 vs. R0, and 5 in R25 vs. F25. Notably, two miRNAs consistently exhibited differential expression across all comparisons, targeting genes linked to plant disease resistance. This suggests their pivotal role in melon&#x2019;s defense against broomrape. The target genes of these miRNAs were confirmed via degradome sequencing and validated by qRT-PCR, ensuring reliable sequencing outcomes. GO and KEGG analyses shed light on the molecular functions and pathways of these differential miRNAs. Furthermore, our study unveiled four trans-kingdom miRNAs, forming a foundation for exploring melon&#x2019;s resistance to broomrape.</p>
</abstract>
<kwd-group>
<kwd>trans-kingdom RNAi</kwd>
<kwd>melon</kwd>
<kwd>broomrape</kwd>
<kwd>small RNA sequencing</kwd>
<kwd>degradome sequencing</kwd>
</kwd-group>
<contract-num rid="cn1">KY2022002</contract-num>
<contract-num rid="cn2">32260758</contract-num>
<contract-sponsor id="cn1">Public Welfare Basic Scientific Research Funds of Xinjiang Academy of Agricultural Sciences</contract-sponsor>
<contract-sponsor id="cn2">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content></contract-sponsor>
<counts>
<fig-count count="6"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="38"/>
<page-count count="10"/>
<word-count count="5735"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Microbe and Virus Interactions with Plants</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="sec1">
<label>1</label>
<title>Introduction</title>
<p>Broomrape (<italic>Phelipanche aegyptiaca</italic> Pers) is a parasitic plant known for inflicting significant economic losses in melon production. This holoparasite attaches itself to the root of the melon plants and establishes a connection with the vascular tissue to absorb nutrients. This parasitic interaction profoundly impairs host development and substantially diminishes melon yields (<xref ref-type="bibr" rid="ref2">Bao et al., 2017</xref>). The attachment between broomrape and its host plant is facilitated through specialized structures called haustorium, which enable bidirectional movement of photosynthates and other essential macromolecules between the host and parasite, including microRNAs (miRNAs) (<xref ref-type="bibr" rid="ref2">Bao et al., 2017</xref>). miRNAs, a class of endogenous single-stranded non-coding RNAs spanning 20&#x2013;40&#x2009;nt, are widespread in plants (<xref ref-type="bibr" rid="ref22">Nizampatnam et al., 2015</xref>; <xref ref-type="bibr" rid="ref27">Song et al., 2019</xref>; <xref ref-type="bibr" rid="ref8">Cui et al., 2020</xref>). They play an important regulatory role in plant growth, development, hormone signal transduction, and stress response. Remarkably, miRNAs can even exert regulatory influence over gene expression in different species, highlighting their capacity for cross-species regulation (<xref ref-type="bibr" rid="ref12">Guo et al., 2017</xref>, <xref ref-type="bibr" rid="ref11">2021</xref>). In addition, some miRNAs exhibit the ability to traverse between host and parasite plants, resulting in the reverse silencing of target genes, leading to the phenomenon of &#x201C;trans-kingdom RNA silencing.&#x201D; This intriguing process involves the transfer of miRNAs and subsequent reverse silencing of target genes, thus influencing gene expression across the plant-parasite interface (<xref ref-type="bibr" rid="ref34">Zhang et al., 2016a</xref>,<xref ref-type="bibr" rid="ref36">b</xref>).</p>
<p>The phenomenon of trans-kingdom RNA silencing was first discovered in the interaction between plants and fungi. Notably, <italic>Botrytis cinerea</italic> demonstrated the ability to transfer a repertoire of small RNAs into plant cells, thereby triggering the silencing of immune response genes across diverse hosts, including <italic>Arabidopsis</italic> and tomato (<xref ref-type="bibr" rid="ref33">Weiberg et al., 2013</xref>). Subsequent investigations unveiled the mechanism by which pathogenic small RNAs, interfacing with <italic>Arabidopsis</italic> AGO1 and harnessing the host RNAi machinery, silence pivotal host genes such as MAPK kinases related to cell wall integrity and genes involved in reactive oxygen species accumulation (<xref ref-type="bibr" rid="ref33">Weiberg et al., 2013</xref>). This strategic employment of small RNAs enabled the pathogen to seize control of host immune responses, thereby facilitating successful invasion (<xref ref-type="bibr" rid="ref34">Zhang et al., 2016a</xref>,<xref ref-type="bibr" rid="ref36">b</xref>, <xref ref-type="bibr" rid="ref35">2022</xref>). The trans-kingdom miRNAs represent an innovative class of effector molecules adept at quelling host immune response. Notably, the miRNA (Bc-siR37) from <italic>Bostalospora cinerea</italic> was significantly enriched in the host <italic>Arabidopsis</italic> during the early stage of invasion. Bc-siR37 targeted a minimum of 15 genes, including WRKY transcription factors and cell wall modifying enzymes in <italic>Arabidopsis</italic> (<xref ref-type="bibr" rid="ref30">Wang B. et al., 2017</xref>). Trans-kingdom miRNAs regulate host target genes to help pathogens successfully invade and colonize. Similarly, in the context of wheat stripe rust, Pst-mi1R1 was identified as a key pathogenic factor of the pathogen, which perturbs the host by targeting the gene <italic>SM638</italic> in the genome of the host wheat. Partial silencing of <italic>SM638</italic> in wheat significantly reduced the resistance to wheat stripe rust in resistant varieties (<xref ref-type="bibr" rid="ref31">Wang M. et al., 2017</xref>). Recent results show have extended the scope of trans-kingdom RNA silencing to parasitic plants, exemplified by <italic>Cuscuta campestris</italic>. Endogenous miRNAs from this parasitic plant can be transported to the host <italic>Arabidopsis thaliana</italic>, guiding and modulating the parasitization process by regulating the expression of different target genes in the host (<xref ref-type="bibr" rid="ref13">Hudzik et al., 2023</xref>). This study underscore the intricate involvement of miRNAs in the host immune pathway, shedding light on a previously uncharted aspect of parasitic interactions.</p>
<p>The trans-kingdom RNAi constitutes a pivotal mechanism in bolstering plant defenses against pathogens, exhibiting bidirectional transmission capabilities between host and pathogen. Trans-kingdom RNAi vectors were constructed using Host-Induced Gene Silencing (HIGS) to silence pathogen virulence-associated target genes, thereby allowing the host to acquire resistance against the pathogen. <xref ref-type="bibr" rid="ref12">Guo et al. (2017</xref>, <xref ref-type="bibr" rid="ref11">2021)</xref> used HIGS technology to produce genetically stable transgenic cotton. Herein, the <italic>VdH1</italic> gene of <italic>Verticillium dahlia</italic>, pivotal for nuclear integrity, was effectively silenced, causing a reduction in virulence due to the loss of nucleus formation in the pathogen. This study is the first to demonstrate the potential of HIGS technology in plant disease resistance research (<xref ref-type="bibr" rid="ref34">Zhang et al., 2016a</xref>,<xref ref-type="bibr" rid="ref36">b</xref>). In a similar vein, <xref ref-type="bibr" rid="ref10">Govindarajulu et al. (2015)</xref> harnessed HIGS technology to produce transgenic lettuce resistant to <italic>Bremia lactucae</italic>. By attenuating a key lettuce gene through RNAi, resistance against the pathogen was realized (<xref ref-type="bibr" rid="ref34">Zhang et al., 2016a</xref>,<xref ref-type="bibr" rid="ref36">b</xref>). Likewise, <xref ref-type="bibr" rid="ref6">Chen et al. (2016)</xref> yielded transgenic wheat with enhanced resistance against <italic>Fusarium graminearum</italic> through HIGS-mediated silencing of the <italic>FcGls1</italic> gene, impairing mycelial cell wall integrity. These remarkable instances offer valuable insights for crop breeding via inter-species RNAi applications.</p>
<p>Evidently, trans-kingdom miRNAs exhibit a propensity for horizontal transfer across several species, such as from humans to malignant malaria-causing insects (<italic>Plasmodium falciparum</italic>) (<xref ref-type="bibr" rid="ref17">LaMonte et al., 2012</xref>); from bacteria to insects (<xref ref-type="bibr" rid="ref17">LaMonte et al., 2012</xref>); from parasitic plants to host plants (<xref ref-type="bibr" rid="ref25">Shahid et al., 2018</xref>; <xref ref-type="bibr" rid="ref13">Hudzik et al., 2023</xref>); from plant to insect (<xref ref-type="bibr" rid="ref14">Ibrahim et al., 2011</xref>); from fungal pathogens to plants (<xref ref-type="bibr" rid="ref33">Weiberg et al., 2013</xref>), from plant-to-fungal pathogens (<xref ref-type="bibr" rid="ref5">Cai et al., 2018</xref>); and even reciprocal transfers between plants and insects (<xref ref-type="bibr" rid="ref32">Wang et al., 2016</xref>). However, a notable gap exists in the exploration of transferable miRNAs in melon plants resistante to the parasitic broomrape. This research aims to unravel the intriguing implications of this novel defense strategy within the interaction between melon and broomrape. By delving into the dynamics of trans-kingdom RNAi, it aspires to shed light on the broader landscape of cross-border RNAi application in plant disease management. Furthermore, this research holds the potential to guide the evolution of pathogenic bacteria as miRNA targets, thereby advancing innovative strategies for combatting broomrape in melon cultivation.</p>
</sec>
<sec sec-type="materials|methods" id="sec2">
<label>2</label>
<title>Materials and methods</title>
<sec id="sec3">
<label>2.1</label>
<title>Plant materials and manipulations</title>
<p>The resistant line K1-6 and sensitive line K2-7 melon cultivars are preserved in Xinjiang Academy of Agricultural Sciences (<xref ref-type="fig" rid="fig1">Figure 1</xref>). Broomrape seeds, identified as <italic>Phelipanche aegyptiaca</italic> Pers, were collected in the Kashgar region of Xinjiang. The pots (Diameter: 16&#x2009;cm, Height: 12.5&#x2009;cm) were filled with a vermiculite-sand mixture (1:1) and infused with 0.02&#x2009;g of broomrape seeds. Each pot was planted with 1&#x2013;2 melon seeds, and this setup was repeated 10 times for each treatment. The pots were placed randomly in a greenhouse (Temperature 28&#x00B0;C, light/darkness&#x2009;=&#x2009;16: 8), receiving regular watering. Root samples from the parasitized melon were collected at 0&#x2009;days (non-infectious treatment) and 25&#x2009;days after infection for the resistant line K1-6 (R0 and R25) and sensitive line K2-7 (F0 and F25). Parasitic broomrape samples were taken for degradome sequencing to find target genes. Each sample was repeated 3 times, rapidly frozen in liquid nitrogen and stored at-80&#x00B0;C. The blank control was induced by strigolactone in Petri dishes and sampled after broomrape germination.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>Resistant and susceptible lines were inoculated for 25&#x2009;days. K1-6 is resistant line and K2-7 is susceptible line, the red box indicates the parasitic broomrape.</p>
</caption>
<graphic xlink:href="fmicb-15-1408926-g001.tif"/>
</fig>
</sec>
<sec id="sec4">
<label>2.2</label>
<title>Construction of small RNA library</title>
<p>The broomrape samples were sent to Genepioneer Biotechnologies for small RNA sequencing. Total RNA was extracted using the RNA Easy Fast (DP452), and small RNA libraries were prepared using the TruSeq Small RNA Sample Prep Kits (Illumina, San Diego, United States) according to standard procedures. After library preparation, the constructed libraries were sequenced using Illumina Hiseq 2000/2500 with a sequence reading length of 1&#x2009;&#x00D7;&#x2009;50&#x2009;bp.</p>
</sec>
<sec id="sec5">
<label>2.3</label>
<title>miRNA sequencing and data processing</title>
<p>The initial sequencing data was aligned with the reference genome, with subsequent filtration to compare with the broomrape genomes. The remaining data were analyzed using the miRNA data analysis software ACGT101-miR (LC Sciences, Houston, Texas, United States). The analysis process of the software is as follows: elimination of 3&#x2032; junction and irrelevant sequences to get clean data; Length screening, filtering out data with baseline length at 18&#x2013;25&#x2009;nt in the clean data; assessment against various RNA databases (mRNA, RFam, Repbase) to exclude non-miRNA data, retaining valid sequences for further analysis.</p>
</sec>
<sec id="sec6">
<label>2.4</label>
<title>miRNA identification</title>
<p>miRNA identification includes both known and novel miRNAs. We compared effective data with plant miRNA precursors and maturing sequences in the miRbase database. This comparison involved BLAST analysis of unique sequences spanning 18&#x2013;25 nucleotides against a specific species precursor in miRBase (miRBase 21.0; <ext-link xlink:href="http://www.mirbase.org/" ext-link-type="uri">http://www.mirbase.org/</ext-link>), in order to identify known miRNA as well as novel miRNA.</p>
</sec>
<sec id="sec7">
<label>2.5</label>
<title>Degradome sequencing identifies cross-border miRNA target genes</title>
<sec id="sec8">
<label>2.5.1</label>
<title>Destruction group building</title>
<p>The samples were sent to Genepioneer Biotechnologies for degradome sequencing. The proposed library steps are as follows: miRNA capture by magnetic beads, adaptor ligation; mixed reverse transcription of Biotinylated Random Primers and mRNA; PCR amplification, after completing the whole library preparation, the constructed library was sequenced with Illumina Hiseq 2000/2500, and the sequencing read length was 1&#x2009;&#x00D7;&#x2009;50&#x2009;bp at the single end.</p>
</sec>
</sec>
<sec id="sec9">
<label>2.6</label>
<title>Degradome data analysis and target gene prediction</title>
<p>The raw data obtained from sequencing were processed through a series of data processing to obtain comparable sequenced pairs for subsequent analysis. The comparable pair sequences were compared with the melon cDNA database sequences to generate the degradome density files. The target gene mRNA sequences paired with the small RNA sequences of sequenced species were predicted by the Target finder target gene prediction software. The Allen Score was used to evaluate the degree of base complementarity between miRNA and target gene binding sites, indicating the binding site pairing penalty (where lower values indicate better pairing).</p>
</sec>
<sec id="sec10">
<label>2.7</label>
<title>Target gene prediction</title>
<p>The target genes corresponding to melon miRNAs predicted by Target finder software based on the principle of base complementary pairing. The mRNAs in the generated degradome density files were combined and analyzed to identify the target genes of the miRNAs. Functional annotation of the target genes was performed by GO enrichment analysis and KEGG enrichment analysis.</p>
</sec>
<sec id="sec11">
<label>2.8</label>
<title>qRT-PCR</title>
<p>Total RNA was digested by DNase I and cDNA were synthesized using Rrans Script Green miRNA Two-Step qRT-PCR SuperMix. Differential miRNAs were randomly chosen for validation, and primers were designed accordingly. Each sample underwent triplicate repetitions to determine relative expression of the target miRNA, H<sub>2</sub>O was used as a negative control and CmADP was used as a reference gene. The primers used in RT-PCR are listed in <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S1</xref>.</p>
</sec>
</sec>
<sec sec-type="results" id="sec12">
<label>3</label>
<title>Results</title>
<sec id="sec13">
<label>3.1</label>
<title>sRNA sequencing in melon varieties resistant or sensitive to broomrape</title>
<p>Here, using small RNA sequencing in melon root systems, we investigated miRNA expression variation in resistant and susceptible cultivars pre-and post-infection by the parasitic plant, broomrape (<xref ref-type="fig" rid="fig1">Figure 1</xref>). The results showed that a total of 12 melon sRNA libraries were prepared for sequencing on the Illumina SE50 platform. The integrity of the sequencing data was notably robust, with each of the four samples yielding a range of clean data reads (ReadSum) spanning from 11,023,148 to 11,977,704, and a clean date base (BaseSum) ranging from 562,180,582 to 610,862,904. These data were characterized by outstanding quality, as evidenced by Q20&#x2009;&#x003E;&#x2009;98% and Q30&#x2009;&#x003E;&#x2009;96% for each sample and GC content with a consistent range from 52.99 to 54.27% (<xref ref-type="table" rid="tab1">Table 1</xref>). The above results indicate that the quality of the raw data obtained from sequencing meets the requirements for subsequent analysis.</p>
<table-wrap position="float" id="tab1">
<label>Table 1</label>
<caption>
<p>sRNA classification annotation statistics.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Type</th>
<th align="center" valign="top">F0</th>
<th align="center" valign="top">F25</th>
<th align="center" valign="top">R0</th>
<th align="center" valign="top">R25</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">rRNA</td>
<td align="char" valign="middle" char="(">1,847,784 (40.7%)</td>
<td align="char" valign="middle" char="(">2,563,363 (45.9%)</td>
<td align="char" valign="middle" char="(">2,545,300 (49.6%)</td>
<td align="char" valign="middle" char="(">2,357,108 (46.7%)</td>
</tr>
<tr>
<td align="left" valign="middle">tRNA</td>
<td align="char" valign="middle" char="(">252,608 (5.5%)</td>
<td align="char" valign="middle" char="(">293,954 (5.2%)</td>
<td align="char" valign="middle" char="(">288,017 (5.6%)</td>
<td align="char" valign="middle" char="(">242,178 (4.9%)</td>
</tr>
<tr>
<td align="left" valign="middle">snRNA</td>
<td align="char" valign="middle" char="(">75 (0.0%)</td>
<td align="char" valign="middle" char="(">115 (0.0%)</td>
<td align="char" valign="middle" char="(">113 (0.0%)</td>
<td align="char" valign="middle" char="(">83 (0.0%)</td>
</tr>
<tr>
<td align="left" valign="middle">snoRNA</td>
<td align="char" valign="middle" char="(">45,619 (1.1%)</td>
<td align="char" valign="middle" char="(">50,316 (0.9%)</td>
<td align="char" valign="middle" char="(">65,543 (1.3%)</td>
<td align="char" valign="middle" char="(">39,611 (0.8%)</td>
</tr>
<tr>
<td align="left" valign="middle">Repeat</td>
<td align="char" valign="middle" char="(">1,586,690 (34.57%)</td>
<td align="char" valign="middle" char="(">1,602,387 (28.6%)</td>
<td align="char" valign="middle" char="(">1,380,110 (26.9%)</td>
<td align="char" valign="middle" char="(">1,436,123 (28.6%)</td>
</tr>
<tr>
<td align="left" valign="middle">miRNA</td>
<td align="char" valign="middle" char="(">19,303 (0.4%)</td>
<td align="char" valign="middle" char="(">7,219 (0.1%)</td>
<td align="char" valign="middle" char="(">7,936 (0.2%)</td>
<td align="char" valign="middle" char="(">2,285 (0.4%)</td>
</tr>
<tr>
<td align="left" valign="middle">Other</td>
<td align="char" valign="middle" char="(">811,715 (1.8%)</td>
<td align="char" valign="middle" char="(">1,099,407 (19.2%)</td>
<td align="char" valign="middle" char="(">842,163 (16.4%)</td>
<td align="char" valign="middle" char="(">964,461 (10.1%)</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The length distribution statistics were performed on the clean reads. The distribution of fragment lengths serves as an illuminating marker to identify small RNA species. miRNAs were usually concentrated at 21&#x2013;22&#x2009;nt, as shown in <xref ref-type="fig" rid="fig2">Figure 2A</xref>. The small RNA lengths of the samples were significantly concentrated at 21&#x2009;nt. In addition, an insightful examination of the first point base preference in clean reads was conducted. A robust preference for the base pair U at the initial position was pronounced in the context of 21&#x2009;nt fragments (<xref ref-type="fig" rid="fig2">Figure 2B</xref>). The results with the annotation of small RNAs are shown in <xref ref-type="table" rid="tab1">Table 1</xref>.</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>miRNA length distribution and nucleotide bias at each position. <bold>(A)</bold> The length distribution chart of identified miRNA. <bold>(B)</bold> Distribution of bases at each miRNA site.</p>
</caption>
<graphic xlink:href="fmicb-15-1408926-g002.tif"/>
</fig>
</sec>
<sec id="sec14">
<label>3.2</label>
<title>miRNA identification</title>
<p>The clean reads after removing the sequences of non-miRNAs were compared with plant miRNA precursors and mature sequences in the miRBase database. A total of 979 known miRNAs and 110 novel miRNAs were identified, which belong to 111 miRNA families. These miRNA families contain miRNA members ranging from 1&#x2013;10. Notably, the MIR156 family emerged as the most populous (<xref ref-type="fig" rid="fig3">Figure 3</xref>), harboring an array of members such as cme-MIR156a, cme-MIR156b and cme-MIR156c. These miRNAs, closely linked with defense mechanisms and the accumulation of secondary metabolites in host plants, underscore their significance within this context.</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption>
<p>Top 20 miRNA family distribution.</p>
</caption>
<graphic xlink:href="fmicb-15-1408926-g003.tif"/>
</fig>
</sec>
<sec id="sec15">
<label>3.3</label>
<title>miRNA differential expression analysis</title>
<p>Employing high-throughput sequencing, the systematic exploration of differential gene expression during biotic stress affords a rapid avenue for the annotation, analysis, and prediction of target genes, thereby facilitating a deeper comprehension of the intricate molecular mechanisms underpinning plant resistance to biotic stress. Our differential expression miRNA analysis for F0 and F25 revealed a total of 39 differentially expressed miRNAs, comprising 32 miRNAs significantly up-regulated and 7 significantly down-regulated (<xref ref-type="fig" rid="fig4">Figures 4A</xref>,<xref ref-type="fig" rid="fig4">D</xref>). For R0 and R25, we observed 35 differentially expressed miRNAs, with 29 displaying upregulated expression and 6 evincing downregulated expression (<xref ref-type="fig" rid="fig4">Figures 4B</xref>,<xref ref-type="fig" rid="fig4">D</xref>). A comparative analysis between R25 and F25 unveiled 5 differentially expressed miRNAs, comprising 4 up-regulated and 1 down-regulated (<xref ref-type="fig" rid="fig4">Figures 4C</xref>,<xref ref-type="fig" rid="fig4">D</xref>).</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption>
<p>Differential expression of miRNA volcanic plot and Venn diagram. <bold>(A)</bold> miRNA volcano plot of sensitive line at 0 and 25&#x2009;days after inoculation. <bold>(B)</bold> miRNA volcano plot of resistant line at 0 and 25&#x2009;days after inoculation. <bold>(C)</bold> miRNA volcano plot of differences between sensitive and resistant lines after 25&#x2009;days of inoculation. <bold>(D)</bold> Histogram of differential miRNA statistics. <bold>(E)</bold> Venn diagram of differentially expressed miRNAs. <bold>(F)</bold> Correlation analysis between qRT-PCR and miRNA sequencing results. The horizontal coordinate of the <bold>A&#x2013;C</bold> plots indicates the change of miRNA expression ploidy in melon after being parasitized by broomrape; the vertical coordinate indicates the significance level of expression difference. The expression of up-regulated miRNA is represented by red dots, and the down-regulated miRNA by green dots indicates that the blue dots have not changed significantly.</p>
</caption>
<graphic xlink:href="fmicb-15-1408926-g004.tif"/>
</fig>
<p>Furthermore, in order to elucidate the mechanism of resistance to Broomrape, this study focused on the differences in miRNA expression before and after broomrape infection in K1-6 (R25 VS R0) and in K2-7 (F25 VS F0). The results showed 22 miRNAs exhibiting consistent significant differential expression trends in R25 VS R0 and F25 VS F0, comprising 4 miRNAs down-regulated and 18 up-regulated in both comparisons. Interestingly, two miRNAs, namely cme-miR408 and cme-miR398a, displayed up-regulation across R25 VS R0, F25 VS F0 and R25 VS F25 comparisons (<xref ref-type="fig" rid="fig4">Figure 4E</xref>).</p>
<p>To verify the accuracy of miRNA sequencing data, a random selection of 6 miRNAs, including cme-miR408, cme-miR398a, cme-miR394a, cme-miR397, cme-miR477a and cme-miR160a, was subject to quantitative real-time PCR (qRT-PCR) analysis. The results showed a high degree of correlation between the qRT-PCR and sequencing data of the selected miRNAs, indicating that the miRNA sequencing data in this experiment were authentic and reliable (<xref ref-type="fig" rid="fig4">Figure 4F</xref>).</p>
</sec>
<sec id="sec16">
<label>3.4</label>
<title>Degradome sequencing for identification of miRNA target genes</title>
<p>In order to gain deeper insights into the role of melon miRNAs in response to broomrape infestation, this study used degradome sequencing to identify the target genes of miRNAs. The degradome library was constructed from a mixed sample of F0, F25, R0 and R25. The data was compared to a quality control to predict the target genes of miRNAs. A total of 172 miRNA-targeted pairs were detected. In a bid to validate the predictive accuracy of these miRNA-target interactions, a judicious selection of 5 pairs of miRNA-targets for subsequent validation, including cme-miR398a-MELO3C026431, cme-miR408-MELO3C008424, cme-miR477b-MELO3C026237, cme-miR160a-MELO3C011372, and cme-miR162-MELO3C031048. The results showed a compelling negative correlation between the expression trends of miRNAs and their corresponding target genes, indicating that the miRNA target genes found by degradome were reliable (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p>
<fig position="float" id="fig5">
<label>Figure 5</label>
<caption>
<p>The verify of miRNA-target gene expression by qRT-PCR.</p>
</caption>
<graphic xlink:href="fmicb-15-1408926-g005.tif"/>
</fig>
</sec>
<sec id="sec17">
<label>3.5</label>
<title>Functional enrichment of differential miRNA target genes</title>
<p>The miRNA target genes identified by degradome sequencing were subject to GO enrichment analysis. The corresponding GO annotations were classified into three categories according to biological process, cellular component and molecular functions, enriching 432, 341, and 156 target genes, respectively. The target genes were significantly enriched in cellular and metabolic processes in the biological process category. In addition, the target genes were significantly enriched in cells, cell parts and binding within the cellular component and molecular functions domain (<xref ref-type="fig" rid="fig6">Figure 6A</xref>).</p>
<fig position="float" id="fig6">
<label>Figure 6</label>
<caption>
<p>Functional annotation of the target genes of the differential expressed miRNAs. <bold>(A)</bold> GO enrichment analysis of differentially expressed miRNA target genes. The horizontal coordinate is the GO function entry and the vertical coordinate is the number of differential genes. <bold>(B)</bold> Analysis of KEEG enrichment pathway for differentially expressed miRNA target genes. Each row in the graph represents a KEGG pathway. The horizontal coordinate is the enrichment factor.</p>
</caption>
<graphic xlink:href="fmicb-15-1408926-g006.tif"/>
</fig>
<p>To glean further insights into the intricate pathways associated with melon miRNAs&#x2019; response to broomrape infestation, a comparative analysis was pursued by aligning the miRNA-target genes with the KEGG database. Target genes are enriched in phosphatidylinositol signaling and purine metabolism, which are important pathways involved in the defense response of melon. The results underscore the importance of miRNA regulation in orchestrating melon&#x2019;s defense response to broomrape infestation (<xref ref-type="fig" rid="fig6">Figure 6B</xref>).</p>
</sec>
<sec id="sec18">
<label>3.6</label>
<title>Transferable miRNA screening</title>
<p>We identified a total of four known miRNAs transferred from melon to broomrape through screening of the transferable miRNA, which were cme-miR156a, cme-m0671-5p, cme-m0738-3p, and cme-miR398a. These miRNAs are mainly involved in structural constituent of cell wall, transposon protein, glutamate receptor and cysteine-rich repeat secretory protein (<xref ref-type="table" rid="tab2">Table 2</xref>). From the functions of these transferred miRNAs, it is speculated that the miRNAs transferred from melon to broomrape may inhibit lepidote growth by affecting the biosynthesis of lepidote miRNAs.</p>
<table-wrap position="float" id="tab2">
<label>Table 2</label>
<caption>
<p>Transferable miRNA and target gene annotation information.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Micro RNA</th>
<th align="left" valign="top" rowspan="2">Target gene ID</th>
<th align="center" valign="top" colspan="2">Gene position</th>
<th align="center" valign="top" rowspan="2">CDs (bp)</th>
<th align="left" valign="top" rowspan="2">Target gene annotation</th>
<th align="left" valign="top" rowspan="2">Arabidopsis homology</th>
</tr>
<tr>
<th align="center" valign="top">Start</th>
<th align="center" valign="top">End (+/&#x2212;)</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top">cme-miR156a</td>
<td align="left" valign="top">MELO3C019024</td>
<td align="center" valign="middle">15,870,317</td>
<td align="center" valign="middle">15,872,012 (&#x2212;)</td>
<td align="center" valign="middle">1,146</td>
<td align="left" valign="top">At5g64090</td>
<td align="left" valign="middle">AT5G64090.1</td>
</tr>
<tr>
<td align="left" valign="top">cme-m0671-5p</td>
<td align="left" valign="top">MELO3C035153</td>
<td align="center" valign="middle">26,029,432</td>
<td align="center" valign="middle">26,037,682 (&#x2212;)</td>
<td align="center" valign="middle">1,156</td>
<td align="left" valign="top">CACTA En-Spm transposon protein</td>
<td align="left" valign="middle">AT3G56600.3</td>
</tr>
<tr>
<td align="left" valign="top">cme-m0738-3p</td>
<td align="left" valign="top">MELO3C010853</td>
<td align="center" valign="middle">28,110,891</td>
<td align="center" valign="middle">28,112,794 (&#x2212;)</td>
<td align="center" valign="middle">1,668</td>
<td align="left" valign="top">Glutamate receptor 2.5-like</td>
<td align="left" valign="middle">AT5G11210.9</td>
</tr>
<tr>
<td align="left" valign="top">cme-miR398a</td>
<td align="left" valign="top">MELO3C026431</td>
<td align="center" valign="middle">21,484,312</td>
<td align="center" valign="middle">21,487,759 (&#x2212;)</td>
<td align="center" valign="middle">906</td>
<td align="left" valign="top">Cysteine-rich repeat secretory 12-like protein</td>
<td align="left" valign="middle">AT2G01660.3</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec sec-type="discussion" id="sec19">
<label>4</label>
<title>Discussion</title>
<p>miRNAs represent a class of non-coding RNAs pivotal for regulatory functions that exert critical influences in almost all aspects of plant growth and development. Additionally, they play important regulatory roles in response to environmental variables such as light and nutrients, and both biotic and abiotic stresses (<xref ref-type="bibr" rid="ref4">Budak et al., 2015</xref>; <xref ref-type="bibr" rid="ref26">Shriram et al., 2016</xref>; <xref ref-type="bibr" rid="ref1">Aslam et al., 2022</xref>; <xref ref-type="bibr" rid="ref15">Kajla et al., 2023</xref>). The advent of small RNA sequencing, based on high-throughput sequencing, has effectively dismantled barriers in miRNA research, facilitating a more incisive exploration of their biological functions. While some reports have used small RNA sequencing to study melon miRNAs, the context of broomrape resistance in melon remains elusive. In this study, a total of 979 known miRNAs and 110 novel miRNAs were identified, belonging to 111 miRNA families. Intriguingly, the MIR156 family emerged as the most prolific, associated with plant resistance and accumulation of secondary metabolites (<xref ref-type="bibr" rid="ref7">Cui et al., 2014</xref>; <xref ref-type="bibr" rid="ref29">Wang, 2014</xref>; <xref ref-type="bibr" rid="ref23">Rahimi et al., 2022</xref>). The results imply a potentially pivotal role for this miRNA family in the resistance of melon to broomrape infestation.</p>
<p>In this study, small RNA sequencing was performed to identify and quantify the expression of miRNAs in the root systems of susceptible and resistant melon lines before and after being infected by broomrape. The differential expression profiles were constructed by comparing the miRNA expression of resistant and susceptible strains. The Venn diagrams unveiled noteworthy distinctions, with cme-miR398a and cme-miR408 standing out as significantly divergent across all examined combinations. It is worth stating that while we identified the differentially expressed miRNAs in resistant and susceptible lines by small RNA sequencing, the target genes of these miRNAs and their regulatory functions remain unknown, which can be well addressed by degradome sequencing. In comparison, predictive target genes based on sequence complimentarity often entail high rates of false positives. Alternately, experimental approaches, though comprehensive, suffer from limitations of expense and time. In this context, degradome sequencing presents an elegant solution. The degradome sequencing technology is based on the principle that plant miRNAs can shear target mRNAs. The intersection of the sheared mRNA fragments are collected and sequenced by high-throughput sequencing technology, thereby allowing the prediction of miRNA target genes by Target finder software. The degradome sequencing effectively overcome the drawbacks and defects of the previous methods for identifying miRNA targets. The target genes of these two miRNAs encode cysteine-rich repeat secretory protein and copper-transporting ATPase. A recent study showed that both proteins are associated with the defense response of the organisms (<xref ref-type="bibr" rid="ref28">Tao et al., 2021</xref>). Therefore, we speculate that these two miRNAs and their target genes play an important role in the resistance of melon to broomrape infestation.</p>
<p>It is well documented that there is a complex transfer of substances between the host and the parasitic organism, including miRNA, and that these substances have essential effects on both parties after the transfer. miRNAs can be transferred through the host to the parasitic plant and improves the host resistance by RNA interference. It has been shown that <italic>Cuscuta campestris</italic> transfers a large number of miRNAs to the host Arabidopsis and Benthamiana, which may function as virulence factors during parasitization as cross-species regulators of host gene expression. However, the extent of long-distance miRNA movement during broomrape parasitization has remained a question. The results of this study show that certain miRNAs are transferred from resistant melon varieties to broomrape. One of the screened transfer miRNAs, cme-miR156a, is a very conserved small non-coding RNAs in plants, and its target gene family is the SPL (SQUAOSA promoter binding like) family (<xref ref-type="bibr" rid="ref20">Matthew et al., 2002</xref>; <xref ref-type="bibr" rid="ref3">Bonnet et al., 2010</xref>). The miR156-SPL pathway plays a regulatory role in several plant growth and developmental processes, including phytohormone signaling and stress response (<xref ref-type="bibr" rid="ref29">Wang, 2014</xref>; <xref ref-type="bibr" rid="ref18">Lei et al., 2016</xref>; <xref ref-type="bibr" rid="ref21">Miao et al., 2019</xref>). Therefore, it is hypothesized that the melon can initiate defense response through the overexpression and transfer of miR156 to broomrape, thus affecting the infestation process of broomrape. In addition, miR398a is mainly involved in the negative regulation of dismutase genes (<xref ref-type="bibr" rid="ref24">Raichur et al., 2007</xref>; <xref ref-type="bibr" rid="ref38">Zidenga et al., 2012</xref>). When melon was infected by broomrape, the expression of cme-miR398a in the roots was elevated, indicating that cme-miR398a may increase ROS content by inhibiting the activity of reactive oxygen species (ROS) enzymes in the cells at the infected site, which eventually resulted in an allergic necrotic response and prevented broomrape infestation.</p>
<p>In this study, a combination of small RNA sequencing and degradome sequencing was used to identify the differentially expressed miRNAs and their target genes between susceptible and resistant lines when infested by broomrape and to explore the defense response of melon against broomrape. In addition, 4 transferable miRNAs were identified in the analysis of melon and broomrape interactions. This comprehensive investigation into the novel defense strategy that underpins melon&#x2019;s interaction with broomrape furnishes valuable reference insights, with implications extending to the broader application of cross-border RNA interference in the realm of plant disease control.</p>
</sec>
<sec sec-type="data-availability" id="sec20">
<title>Data availability statement</title>
<p>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/<xref ref-type="supplementary-material" rid="SM1">Supplementary material</xref>.</p>
</sec>
<sec sec-type="author-contributions" id="sec21">
<title>Author contributions</title>
<p>J-CM: Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing. MY: Writing &#x2013; review &#x0026; editing. J-HL: Funding acquisition, Writing &#x2013; review &#x0026; editing. J-YY: Writing &#x2013; original draft. H-JW: Funding acquisition, Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing.</p>
</sec>
</body>
<back>
<sec sec-type="funding-information" id="sec22">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the Public Welfare Basic Scientific Research Funds of Xinjiang Academy of Agricultural Sciences (Grant no. KY2022002), and National Natural Science Foundation of China (Grant no. 32260758).</p>
</sec>
<ack>
<p>The authors are grateful to BL (Xinjiang Academy of Agricultural Sciences) for his suggestions to this article.</p>
</ack>
<sec sec-type="COI-statement" id="sec23">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>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.</p>
</sec>
<sec sec-type="supplementary-material" id="sec24">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2024.1408926/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmicb.2024.1408926/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Table_1.DOCX" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aslam</surname> <given-names>M.</given-names></name> <name><surname>Fakher</surname> <given-names>B.</given-names></name> <name><surname>Qin</surname> <given-names>Y.</given-names></name></person-group> (<year>2022</year>). <article-title>Big role of small RNAs in female gametophyte development</article-title>. <source>Int. J. Mol. Sci.</source> <volume>23</volume>:<fpage>1979</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ijms23041979</pub-id>, PMID: <pub-id pub-id-type="pmid">35216096</pub-id></citation></ref>
<ref id="ref2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bao</surname> <given-names>Y. Z.</given-names></name> <name><surname>Yao</surname> <given-names>Z. Q.</given-names></name> <name><surname>Cao</surname> <given-names>X. L.</given-names></name> <name><surname>Peng</surname> <given-names>J. F.</given-names></name> <name><surname>Xu</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>M. X.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Transcriptome analysis of <italic>Phelipanche aegyptiaca</italic> seed germination mechanisms stimulated by fluridone, TIS108, and GR24</article-title>. <source>PLoS One</source> <volume>12</volume>:<fpage>e187539</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0187539</pub-id>, PMID: <pub-id pub-id-type="pmid">29099877</pub-id></citation></ref>
<ref id="ref3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bonnet</surname> <given-names>E.</given-names></name> <name><surname>He</surname> <given-names>Y.</given-names></name> <name><surname>Billiau</surname> <given-names>K.</given-names></name> <name><surname>Van de Peer</surname> <given-names>Y.</given-names></name></person-group> (<year>2010</year>). <article-title>TAPIR, a web server for the prediction of plant microRNA targets, including target mimics</article-title>. <source>Bioinformatics</source> <volume>26</volume>, <fpage>1566</fpage>&#x2013;<lpage>1568</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btq233</pub-id></citation></ref>
<ref id="ref4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Budak</surname> <given-names>H.</given-names></name> <name><surname>Kantar</surname> <given-names>M.</given-names></name> <name><surname>Bulut</surname> <given-names>R.</given-names></name> <name><surname>Akpinar</surname> <given-names>B. A.</given-names></name></person-group> (<year>2015</year>). <article-title>Stress responsive miRNAs and isomiRs in cereals</article-title>. <source>Plant Sci.</source> <volume>235</volume>, <fpage>1</fpage>&#x2013;<lpage>13</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.plantsci.2015.02.008</pub-id>, PMID: <pub-id pub-id-type="pmid">25900561</pub-id></citation></ref>
<ref id="ref5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cai</surname> <given-names>Q.</given-names></name> <name><surname>Qiao</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>M.</given-names></name> <name><surname>He</surname> <given-names>B.</given-names></name> <name><surname>Lin</surname> <given-names>F. M.</given-names></name> <name><surname>Palmquist</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Plants send small RNAs in extracellular vesicles to fungal pathogen to silence virulence genes</article-title>. <source>Science</source> <volume>360</volume>, <fpage>1126</fpage>&#x2013;<lpage>1129</lpage>. doi: <pub-id pub-id-type="doi">10.1126/science.aar4142</pub-id>, PMID: <pub-id pub-id-type="pmid">29773668</pub-id></citation></ref>
<ref id="ref6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>W.</given-names></name> <name><surname>Kastner</surname> <given-names>C.</given-names></name> <name><surname>Nowara</surname> <given-names>D.</given-names></name> <name><surname>Oliveira-Garcia</surname> <given-names>E.</given-names></name> <name><surname>Rutten</surname> <given-names>T.</given-names></name> <name><surname>Zhao</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Host-induced silencing of <italic>Fusarium culmorum</italic> genes protects wheat from infection</article-title>. <source>J. Exp. Bot.</source> <volume>67</volume>, <fpage>4979</fpage>&#x2013;<lpage>4991</lpage>. doi: <pub-id pub-id-type="doi">10.1093/jxb/erw263</pub-id>, PMID: <pub-id pub-id-type="pmid">27540093</pub-id></citation></ref>
<ref id="ref7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cui</surname> <given-names>L. G.</given-names></name> <name><surname>Shan</surname> <given-names>J. X.</given-names></name> <name><surname>Shi</surname> <given-names>M.</given-names></name> <name><surname>Gao</surname> <given-names>J. P.</given-names></name> <name><surname>Lin</surname> <given-names>H. X.</given-names></name></person-group> (<year>2014</year>). <article-title>The miR156-SPL9-DFR pathway coordinates the relationship between development and abiotic stress tolerance in plants</article-title>. <source>Plant J.</source> <volume>80</volume>, <fpage>1108</fpage>&#x2013;<lpage>1117</lpage>. doi: <pub-id pub-id-type="doi">10.1111/tpj.12712</pub-id></citation></ref>
<ref id="ref8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cui</surname> <given-names>L.</given-names></name> <name><surname>Zheng</surname> <given-names>F.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Zhang</surname> <given-names>C.</given-names></name> <name><surname>Xiao</surname> <given-names>F.</given-names></name> <name><surname>Ye</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>MiR156a-targeted SBP-box transcription factor SlSPL13 regulates inflorescence morphogenesis by directly activating SFT in tomato</article-title>. <source>Plant Biotechnol. J.</source> <volume>18</volume>, <fpage>1670</fpage>&#x2013;<lpage>1682</lpage>. doi: <pub-id pub-id-type="doi">10.1111/pbi.13331</pub-id>, PMID: <pub-id pub-id-type="pmid">31916387</pub-id></citation></ref>
<ref id="ref9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>David-Schwartz</surname> <given-names>R.</given-names></name> <name><surname>Runo</surname> <given-names>S.</given-names></name> <name><surname>Townsley</surname> <given-names>B.</given-names></name> <name><surname>Machuka</surname> <given-names>J.</given-names></name> <name><surname>Sinha</surname> <given-names>N.</given-names></name></person-group> (<year>2008</year>). <article-title>Long-distance transport of mRNA via parenchyma cells and phloem across the host-parasite junction in <italic>Cuscuta</italic></article-title>. <source>New Phytol.</source> <volume>179</volume>, <fpage>1133</fpage>&#x2013;<lpage>1141</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1469-8137.2008.02540.x</pub-id>, PMID: <pub-id pub-id-type="pmid">18631294</pub-id></citation></ref>
<ref id="ref10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Govindarajulu</surname> <given-names>M.</given-names></name> <name><surname>Epstein</surname> <given-names>L.</given-names></name> <name><surname>Wroblewski</surname> <given-names>T.</given-names></name> <name><surname>Michelmore</surname> <given-names>R. W.</given-names></name></person-group> (<year>2015</year>). <article-title>Host-induced gene silencing inhibits the biotrophic pathogen causing downy mildew of lettuce</article-title>. <source>Plant Biotechnol. J.</source> <volume>13</volume>, <fpage>875</fpage>&#x2013;<lpage>883</lpage>. doi: <pub-id pub-id-type="doi">10.1111/pbi.12307</pub-id>, PMID: <pub-id pub-id-type="pmid">25487781</pub-id></citation></ref>
<ref id="ref11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guo</surname> <given-names>C.</given-names></name> <name><surname>Jiang</surname> <given-names>Y.</given-names></name> <name><surname>Shi</surname> <given-names>M.</given-names></name> <name><surname>Wu</surname> <given-names>X.</given-names></name> <name><surname>Wu</surname> <given-names>G.</given-names></name></person-group> (<year>2021</year>). <article-title>ABI5 acts downstream of miR159 to delay vegetative phase change in <italic>Arabidopsis</italic></article-title>. <source>New Phytol.</source> <volume>231</volume>, <fpage>339</fpage>&#x2013;<lpage>350</lpage>. doi: <pub-id pub-id-type="doi">10.1111/nph.17371</pub-id>, PMID: <pub-id pub-id-type="pmid">33774835</pub-id></citation></ref>
<ref id="ref12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guo</surname> <given-names>C.</given-names></name> <name><surname>Xu</surname> <given-names>Y.</given-names></name> <name><surname>Shi</surname> <given-names>M.</given-names></name> <name><surname>Lai</surname> <given-names>Y.</given-names></name> <name><surname>Wu</surname> <given-names>X.</given-names></name> <name><surname>Wang</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Repression of miR156 by miR159 regulates the timing of the juvenile-to-adult transition in <italic>Arabidopsis</italic></article-title>. <source>Plant Cell</source> <volume>29</volume>, <fpage>1293</fpage>&#x2013;<lpage>1304</lpage>. doi: <pub-id pub-id-type="doi">10.1105/tpc.16.00975</pub-id>, PMID: <pub-id pub-id-type="pmid">28536099</pub-id></citation></ref>
<ref id="ref13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hudzik</surname> <given-names>C.</given-names></name> <name><surname>Maguire</surname> <given-names>S.</given-names></name> <name><surname>Guan</surname> <given-names>S.</given-names></name> <name><surname>Held</surname> <given-names>J.</given-names></name> <name><surname>Axtell</surname> <given-names>M. J.</given-names></name></person-group> (<year>2023</year>). <article-title>Trans-species microRNA loci in the parasitic plant <italic>Cuscuta campestris</italic> have a U6-like snRNA promoter</article-title>. <source>Plant Cell</source> <volume>35</volume>, <fpage>1834</fpage>&#x2013;<lpage>1847</lpage>. doi: <pub-id pub-id-type="doi">10.1093/plcell/koad076</pub-id>, PMID: <pub-id pub-id-type="pmid">36896651</pub-id></citation></ref>
<ref id="ref14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ibrahim</surname> <given-names>H. M.</given-names></name> <name><surname>Alkharouf</surname> <given-names>N. W.</given-names></name> <name><surname>Meyer</surname> <given-names>S. L.</given-names></name> <name><surname>Aly</surname> <given-names>M. A.</given-names></name> <name><surname>Gamal El-Din Ael</surname> <given-names>K.</given-names></name> <name><surname>Hussein</surname> <given-names>E. H.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Post-transcriptional gene silencing of root-knot nematode in transformed soybean roots</article-title>. <source>Exp. Parasitol.</source> <volume>127</volume>, <fpage>90</fpage>&#x2013;<lpage>99</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.exppara.2010.06.037</pub-id>, PMID: <pub-id pub-id-type="pmid">20599433</pub-id></citation></ref>
<ref id="ref15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kajla</surname> <given-names>M.</given-names></name> <name><surname>Roy</surname> <given-names>A.</given-names></name> <name><surname>Singh</surname> <given-names>I. K.</given-names></name> <name><surname>Singh</surname> <given-names>A.</given-names></name></person-group> (<year>2023</year>). <article-title>Regulation of the regulators: transcription factors controlling biosynthesis of plant secondary metabolites during biotic stresses and their regulation by miRNAs</article-title>. <source>Front. Plant Sci.</source> <volume>14</volume>:<fpage>1126567</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fpls.2023.1126567</pub-id>, PMID: <pub-id pub-id-type="pmid">36938003</pub-id></citation></ref>
<ref id="ref16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>G.</given-names></name> <name><surname>Westwood</surname> <given-names>J. H.</given-names></name></person-group> (<year>2015</year>). <article-title>Macromolecule exchange in <italic>Cuscuta</italic>-host plant interactions</article-title>. <source>Curr. Opin. Plant Biol.</source> <volume>26</volume>, <fpage>20</fpage>&#x2013;<lpage>25</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.pbi.2015.05.012</pub-id>, PMID: <pub-id pub-id-type="pmid">26051214</pub-id></citation></ref>
<ref id="ref17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>LaMonte</surname> <given-names>G.</given-names></name> <name><surname>Philip</surname> <given-names>N.</given-names></name> <name><surname>Reardon</surname> <given-names>J.</given-names></name> <name><surname>Lacsina</surname> <given-names>J. R.</given-names></name> <name><surname>Majoros</surname> <given-names>W.</given-names></name> <name><surname>Chapman</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Translocation of sickle cell erythrocyte microRNAs into <italic>Plasmodium falciparum</italic> inhibits parasite translation and contributes to malaria resistance</article-title>. <source>Cell Host Microbe</source> <volume>12</volume>, <fpage>187</fpage>&#x2013;<lpage>199</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.chom.2012.06.007</pub-id>, PMID: <pub-id pub-id-type="pmid">22901539</pub-id></citation></ref>
<ref id="ref18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lei</surname> <given-names>K. J.</given-names></name> <name><surname>Lin</surname> <given-names>Y. M.</given-names></name> <name><surname>An</surname> <given-names>G. Y.</given-names></name></person-group> (<year>2016</year>). <article-title>MiR156 modulates rhizosphere acidification in response to phosphate limitation in <italic>Arabidopsis</italic></article-title>. <source>J. Plant Res.</source> <volume>129</volume>, <fpage>275</fpage>&#x2013;<lpage>284</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s10265-015-0778-8</pub-id>, PMID: <pub-id pub-id-type="pmid">26659856</pub-id></citation></ref>
<ref id="ref19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>H.</given-names></name> <name><surname>Wang</surname> <given-names>X.</given-names></name> <name><surname>Wang</surname> <given-names>H. D.</given-names></name> <name><surname>Wu</surname> <given-names>J.</given-names></name> <name><surname>Ren</surname> <given-names>J.</given-names></name> <name><surname>Meng</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title><italic>Escherichia coli</italic> noncoding RNAs can affect gene expression and physiology of <italic>Caenorhabditis elegans</italic></article-title>. <source>Nat. Commun.</source> <volume>3</volume>:<fpage>1073</fpage>. doi: <pub-id pub-id-type="doi">10.1038/ncomms2071</pub-id>, PMID: <pub-id pub-id-type="pmid">23011127</pub-id></citation></ref>
<ref id="ref20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Matthew</surname> <given-names>W. R.</given-names></name> <name><surname>Brenda</surname> <given-names>J. R.</given-names></name> <name><surname>Lee</surname> <given-names>P. L.</given-names></name> <name><surname>Christopher</surname> <given-names>B. B.</given-names></name> <name><surname>Bonnie</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2002</year>). <article-title>Prediction of plant microRNA targets</article-title>. <source>Cell</source> <volume>110</volume>, <fpage>513</fpage>&#x2013;<lpage>520</lpage>. doi: <pub-id pub-id-type="doi">10.1016/s0092-8674(02)00863-2</pub-id></citation></ref>
<ref id="ref21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Miao</surname> <given-names>C.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Yao</surname> <given-names>J.</given-names></name> <name><surname>Hua</surname> <given-names>K.</given-names></name> <name><surname>Liu</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>The grain yield modulator miR156 regulates seed dormancy through the gibberellin pathway in rice</article-title>. <source>Nat. Commun.</source> <volume>10</volume>:<fpage>3822</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41467-019-11830-5</pub-id>, PMID: <pub-id pub-id-type="pmid">31444356</pub-id></citation></ref>
<ref id="ref22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nizampatnam</surname> <given-names>N. R.</given-names></name> <name><surname>Schreier</surname> <given-names>S. J.</given-names></name> <name><surname>Damodaran</surname> <given-names>S.</given-names></name> <name><surname>Adhikari</surname> <given-names>S.</given-names></name> <name><surname>Subramanian</surname> <given-names>S.</given-names></name></person-group> (<year>2015</year>). <article-title>microRNA160 dictates stage-specific auxin and cytokinin sensitivities and directs soybean nodule development</article-title>. <source>Plant J.</source> <volume>84</volume>, <fpage>140</fpage>&#x2013;<lpage>153</lpage>. doi: <pub-id pub-id-type="doi">10.1111/tpj.12965</pub-id>, PMID: <pub-id pub-id-type="pmid">26287653</pub-id></citation></ref>
<ref id="ref23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rahimi</surname> <given-names>A.</given-names></name> <name><surname>Karami</surname> <given-names>O.</given-names></name> <name><surname>Balazadeh</surname> <given-names>S.</given-names></name> <name><surname>Offringa</surname> <given-names>R.</given-names></name></person-group> (<year>2022</year>). <article-title>MiR156-independent repression of the ageing pathway by longevity-promoting AHL proteins in <italic>Arabidopsis</italic></article-title>. <source>New Phytol.</source> <volume>235</volume>, <fpage>2424</fpage>&#x2013;<lpage>2438</lpage>. doi: <pub-id pub-id-type="doi">10.1111/nph.18292</pub-id>, PMID: <pub-id pub-id-type="pmid">35642455</pub-id></citation></ref>
<ref id="ref24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Raichur</surname> <given-names>S.</given-names></name> <name><surname>Lau</surname> <given-names>P.</given-names></name> <name><surname>Staels</surname> <given-names>B.</given-names></name> <name><surname>Muscat</surname> <given-names>G. E.</given-names></name></person-group> (<year>2007</year>). <article-title>Retinoid-related orphan receptor gamma regulates several genes that control metabolism in skeletal muscle cells: links to modulation of reactive oxygen species production</article-title>. <source>J. Mol. Endocrinol.</source> <volume>39</volume>, <fpage>29</fpage>&#x2013;<lpage>44</lpage>. doi: <pub-id pub-id-type="doi">10.1677/jme.1.00010</pub-id>, PMID: <pub-id pub-id-type="pmid">17601883</pub-id></citation></ref>
<ref id="ref25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shahid</surname> <given-names>S.</given-names></name> <name><surname>Kim</surname> <given-names>G.</given-names></name> <name><surname>Johnson</surname> <given-names>N. R.</given-names></name> <name><surname>Wafula</surname> <given-names>E.</given-names></name> <name><surname>Wang</surname> <given-names>F.</given-names></name> <name><surname>Coruh</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>MicroRNAs from the parasitic plant <italic>Cuscuta campestris</italic> target host messenger RNAs</article-title>. <source>Nature</source> <volume>553</volume>, <fpage>82</fpage>&#x2013;<lpage>85</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nature25027</pub-id>, PMID: <pub-id pub-id-type="pmid">29300014</pub-id></citation></ref>
<ref id="ref26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shriram</surname> <given-names>V.</given-names></name> <name><surname>Kumar</surname> <given-names>V.</given-names></name> <name><surname>Devarumath</surname> <given-names>R. M.</given-names></name> <name><surname>Khare</surname> <given-names>T. S.</given-names></name> <name><surname>Wani</surname> <given-names>S. H.</given-names></name></person-group> (<year>2016</year>). <article-title>MicroRNAs as potential targets for abiotic stress tolerance in plants</article-title>. <source>Front. Plant Sci.</source> <volume>7</volume>:<fpage>817</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fpls.2016.00817</pub-id>, PMID: <pub-id pub-id-type="pmid">27379117</pub-id></citation></ref>
<ref id="ref27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>X.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Cao</surname> <given-names>X.</given-names></name> <name><surname>Qi</surname> <given-names>Y.</given-names></name></person-group> (<year>2019</year>). <article-title>microRNAs and their regulatory roles in plant-environment interactions</article-title>. <source>Annu. Rev. Plant Biol.</source> <volume>70</volume>, <fpage>489</fpage>&#x2013;<lpage>525</lpage>. doi: <pub-id pub-id-type="doi">10.1146/annurev-arplant-050718-100334</pub-id>, PMID: <pub-id pub-id-type="pmid">30848930</pub-id></citation></ref>
<ref id="ref28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tao</surname> <given-names>Z.</given-names></name> <name><surname>Yalong</surname> <given-names>Z.</given-names></name> <name><surname>Dewei</surname> <given-names>T.</given-names></name> <name><surname>Xiaozhong</surname> <given-names>Z.</given-names></name> <name><surname>Xueshi</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Rice blast resistance-associated genes based on different RNA-seq resources</article-title>. <source>Agric. Biotechnol.</source> <volume>10</volume>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>,</citation></ref>
<ref id="ref29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>J. W.</given-names></name></person-group> (<year>2014</year>). <article-title>Regulation of flowering time by the miR156-mediated age pathway</article-title>. <source>J. Exp. Bot.</source> <volume>65</volume>, <fpage>4723</fpage>&#x2013;<lpage>4730</lpage>. doi: <pub-id pub-id-type="doi">10.1093/jxb/eru246</pub-id>, PMID: <pub-id pub-id-type="pmid">24958896</pub-id></citation></ref>
<ref id="ref30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>B.</given-names></name> <name><surname>Sun</surname> <given-names>Y.</given-names></name> <name><surname>Song</surname> <given-names>N.</given-names></name> <name><surname>Zhao</surname> <given-names>M.</given-names></name> <name><surname>Liu</surname> <given-names>R.</given-names></name> <name><surname>Feng</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title><italic>Puccinia striiformis</italic> f. sp. <italic>tritici</italic> microRNA-like RNA 1 (<italic>Pst</italic>-milR1), an important pathogenicity factor of <italic>Pst</italic>, impairs wheat resistance to <italic>Pst</italic> by suppressing the wheat pathogenesis-related 2 gene</article-title>. <source>New Phytol.</source> <volume>215</volume>, <fpage>338</fpage>&#x2013;<lpage>350</lpage>. doi: <pub-id pub-id-type="doi">10.1111/nph.14577</pub-id>, PMID: <pub-id pub-id-type="pmid">28464281</pub-id></citation></ref>
<ref id="ref31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>M.</given-names></name> <name><surname>Weiberg</surname> <given-names>A.</given-names></name> <name><surname>Dellota</surname> <given-names>E. J.</given-names></name> <name><surname>Yamane</surname> <given-names>D.</given-names></name> <name><surname>Jin</surname> <given-names>H.</given-names></name></person-group> (<year>2017</year>). <article-title>Botrytis small RNA Bc-siR37 suppresses plant defense genes by cross-kingdom RNAi</article-title>. <source>RNA Biol.</source> <volume>14</volume>, <fpage>421</fpage>&#x2013;<lpage>428</lpage>. doi: <pub-id pub-id-type="doi">10.1080/15476286.2017.1291112</pub-id>, PMID: <pub-id pub-id-type="pmid">28267415</pub-id></citation></ref>
<ref id="ref32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>M.</given-names></name> <name><surname>Weiberg</surname> <given-names>A.</given-names></name> <name><surname>Lin</surname> <given-names>F. M.</given-names></name> <name><surname>Thomma</surname> <given-names>B. P.</given-names></name> <name><surname>Huang</surname> <given-names>H. D.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Bidirectional cross-kingdom RNAi and fungal uptake of external RNAs confer plant protection</article-title>. <source>Nat. Plants</source> <volume>2</volume>:<fpage>16151</fpage>. doi: <pub-id pub-id-type="doi">10.1038/nplants.2016.151</pub-id>, PMID: <pub-id pub-id-type="pmid">27643635</pub-id></citation></ref>
<ref id="ref33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Weiberg</surname> <given-names>A.</given-names></name> <name><surname>Wang</surname> <given-names>M.</given-names></name> <name><surname>Lin</surname> <given-names>F. M.</given-names></name> <name><surname>Zhao</surname> <given-names>H.</given-names></name> <name><surname>Zhang</surname> <given-names>Z.</given-names></name> <name><surname>Kaloshian</surname> <given-names>I.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>Fungal small RNAs suppress plant immunity by hijacking host RNA interference pathways</article-title>. <source>Science</source> <volume>342</volume>, <fpage>118</fpage>&#x2013;<lpage>123</lpage>. doi: <pub-id pub-id-type="doi">10.1126/science.1239705</pub-id>, PMID: <pub-id pub-id-type="pmid">24092744</pub-id></citation></ref>
<ref id="ref34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>T.</given-names></name> <name><surname>Jin</surname> <given-names>Y.</given-names></name> <name><surname>Zhao</surname> <given-names>J. H.</given-names></name> <name><surname>Gao</surname> <given-names>F.</given-names></name> <name><surname>Zhou</surname> <given-names>B. J.</given-names></name> <name><surname>Fang</surname> <given-names>Y. Y.</given-names></name> <etal/></person-group>. (<year>2016a</year>). <article-title>Host-induced gene silencing of the target gene in fungal cells confers effective resistance to the cotton wilt disease pathogen <italic>Verticillium dahliae</italic></article-title>. <source>Mol. Plant</source> <volume>9</volume>, <fpage>939</fpage>&#x2013;<lpage>942</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.molp.2016.02.008</pub-id>, PMID: <pub-id pub-id-type="pmid">26925819</pub-id></citation></ref>
<ref id="ref35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>B. S.</given-names></name> <name><surname>Li</surname> <given-names>Y. C.</given-names></name> <name><surname>Guo</surname> <given-names>H. S.</given-names></name> <name><surname>Zhao</surname> <given-names>J. H.</given-names></name></person-group> (<year>2022</year>). <article-title><italic>Verticillium dahliae</italic> secretes small RNA to target host MIR157d and retard plant floral transition during infection</article-title>. <source>Front. Plant Sci.</source> <volume>13</volume>:<fpage>847086</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fpls.2022.847086</pub-id>, PMID: <pub-id pub-id-type="pmid">35519822</pub-id></citation></ref>
<ref id="ref36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>T.</given-names></name> <name><surname>Zhao</surname> <given-names>Y. L.</given-names></name> <name><surname>Zhao</surname> <given-names>J. H.</given-names></name> <name><surname>Wang</surname> <given-names>S.</given-names></name> <name><surname>Jin</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>Z. Q.</given-names></name> <etal/></person-group>. (<year>2016b</year>). <article-title>Cotton plants export microRNAs to inhibit virulence gene expression in a fungal pathogen</article-title>. <source>Nat. Plants</source> <volume>2</volume>:<fpage>16153</fpage>. doi: <pub-id pub-id-type="doi">10.1038/nplants.2016.153</pub-id>, PMID: <pub-id pub-id-type="pmid">27668926</pub-id></citation></ref>
<ref id="ref37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>C.</given-names></name> <name><surname>Liu</surname> <given-names>T.</given-names></name> <name><surname>Chang</surname> <given-names>Y. N.</given-names></name> <name><surname>Duan</surname> <given-names>C. G.</given-names></name></person-group> (<year>2019</year>). <article-title>Small RNA functions as a trafficking effector in plant immunity</article-title>. <source>Int. J. Mol. Sci.</source> <volume>20</volume>:<fpage>2816</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ijms20112816</pub-id>, PMID: <pub-id pub-id-type="pmid">31181829</pub-id></citation></ref>
<ref id="ref38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zidenga</surname> <given-names>T.</given-names></name> <name><surname>Leyva-Guerrero</surname> <given-names>E.</given-names></name> <name><surname>Moon</surname> <given-names>H.</given-names></name> <name><surname>Siritunga</surname> <given-names>D.</given-names></name> <name><surname>Sayre</surname> <given-names>R.</given-names></name></person-group> (<year>2012</year>). <article-title>Extending cassava root shelf life via reduction of reactive oxygen species production</article-title>. <source>Plant Physiol.</source> <volume>159</volume>, <fpage>1396</fpage>&#x2013;<lpage>1407</lpage>. doi: <pub-id pub-id-type="doi">10.1104/pp.112.200345</pub-id>, PMID: <pub-id pub-id-type="pmid">22711743</pub-id></citation></ref>
</ref-list>
</back>
</article>