Genome-Wide Analysis of MicroRNA Responses to the Phytohormone Abscisic Acid in Populus euphratica

Comparative Analysis of miRNA Expression Profiles under Salt Stress in Wheat

Salt stress is one of the important environmental factors that inhibit the normal growth and development of plants. Plants have evolved various mechanisms, including signal transduction regulation, physiological regulation, and gene transcription regulation, to adapt to environmental stress. MicroRNAs (miRNAs) play a role in regulating mRNA expression. Nevertheless, miRNAs related to salt stress are rarely reported in bread wheat (Triticum aestivum L.). In this study, using high-throughput sequencing, we analyzed the miRNA expression profile of wheat under salt stress. We identified 360 conserved and 859 novel miRNAs, of which 49 showed considerable changes in transcription levels after salt treatment. Among them, 25 were dramatically upregulated and 24 were downregulated. Using real-time quantitative PCR, we detected significant changes in the relative expression of miRNAs, and the results showed the same trend as the sequencing data. In the salt-treated group, miR109 had a higher expression level, while miR60 and miR202 had lower expression levels. Furthermore, 21 miRNAs with significant changes were selected from the differentially expressed miRNAs, and 1023 candidate target genes were obtained through the prediction of the website psRNATarget. Gene ontology (GO) analysis of the candidate target genes showed that the expressed miRNA may be involved in the response to biological processes, molecular functions, and cellular components. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis confirmed their important functions in RNA degradation, metabolic pathways, synthesis pathways, peroxisome, environmental adaptation, global and overview maps, and stress adaptation and the MAPK signal pathway. These findings provide a basis for further exploring the function of miRNA in wheat salt tolerance.

Comparative bioinformatics analysis and abiotic stress responses of expansin proteins in Cucurbitaceae members: watermelon and melon

Watermelon and melon are members of the Cucurbitaceae family including economically significant crops in the world. The expansin protein family, which is one of the members of the cell wall, breaks down the non-covalent bonds between cell wall polysaccharides, causing pressure-dependent cell expansion. Comparative bioinformatics and molecular characterization analysis of the expansin protein family were carried out in the watermelon (Citrullus lanatus) and melon (Cucumis melo) plants in the study. Gene expression levels of expansin family members were analyzed in leaf and root tissues of watermelon and melon under ABA, drought, heat, cold, and salt stress conditions by quantitative real-time PCR analysis. After comprehensive searches, 40 expansin proteins (22 ClaEXPA, 14 ClaEXPLA, and 4 ClaEXPB) in watermelon and 43 expansin proteins (19 CmEXPA, 15 CmEXPLA, 3 CmEXPB, and 6 CmEXPLB) in melon were identified. The greatest orthologous genes were identified with soybean expansin genes for watermelon and melon. However, the latest divergence time between orthologous genes was determined with poplar expansin genes for watermelon and melon expansin genes. ClaEXPA-04, ClaEXPA-09, ClaEXPB-01, ClaEXPB-03, and ClaEXPLA-13 genes in watermelon and CmEXPA-12, CmEXPA-10, and CmEXPLA-01 genes in melon can be involved in tissue development and abiotic stress response of the plant. The current study combining bioinformatics and experimental analysis can provide a detailed characterization of the expansin superfamily which has roles in growth and reaction to the stress of the plant. The study ensures detailed data for future studies examining gene functions including the roles in plant growth and stress conditions.

Plant responses to drought stress: microRNAs in action

Drought is common in most regions of the world, and it has a significant impact on plant growth and development. Plants, on the other hand, have evolved their own defense systems to deal with the extreme weather. The reprogramming of gene expression by microRNAs (miRNAs) is one of these defense mechanisms. miRNAs are short noncoding RNAs that have emerged as key post-transcriptional gene regulators in a variety of species. Drought stress modulates the expression of certain miRNAs that are functionally conserved across plant species. These characteristics imply that miRNA-based genetic changes might improve drought resistance in plants. This study highlights current knowledge of plant miRNA biogenesis, regulatory mechanisms and their role in drought stress responses. miRNAs functions and their adaptations by plants during drought stress has also been explained that can be exploited to promote drought-resistance among economically important crops.

Distinguishing friends from foes: Can smRNAs modulate plant interactions with beneficial and pathogenic organisms?

In their agro-ecological habitats, plants are constantly challenged by fungal interactions that might be pathogenic or beneficial in nature, and thus, plants need to exhibit appropriate responses to discriminate between them. Such interactions involve sophisticated molecular mechanism of signal exchange, signal transduction and regulation of gene expression. Small RNAs (smRNAs), including the microRNAs (miRNAs), form an essential layer of regulation in plant developmental processes as well as in plant adaptation to environmental stresses, being key for the outcome during plant-microbial interactions. Further, smRNAs are mobile signals that can go across kingdoms from one interacting partner to the other and hence can be used as communication as well as regulatory tools not only by the host plant but also by the colonising fungus. Here, largely with a focus on plant-fungal interactions and miRNAs, we will discuss the role of smRNAs, and how they might help plants to discriminate between friends and foes.

Identification of Novel miRNAs and Their Target Genes in the Response to Abscisic Acid in Arabidopsis

miRNAs are involved in various biological processes, including adaptive responses to abiotic stress. To understand the role of miRNAs in the response to ABA, ABA-responsive miRNAs were identified by small RNA sequencing in wild-type Arabidopsis, as well as in abi1td, mkkk17, and mkkk18 mutants. We identified 10 novel miRNAs in WT after ABA treatment, while in abi1td, mkkk17, and mkkk18 mutants, three, seven, and nine known miRNAs, respectively, were differentially expressed after ABA treatment. One novel miRNA (miRn-8) was differentially expressed in the mkkk17 mutant. Potential target genes of the miRNA panel were identified using psRNATarget. Sequencing results were validated by quantitative RT-PCR of several known and novel miRNAs in all genotypes. Of the predicted targets of novel miRNAs, seven target genes of six novel miRNAs were further validated by 5' RLM-RACE. Gene ontology analyses showed the potential target genes of ABA-responsive known and novel miRNAs to be involved in diverse cellular processes in plants, including development and stomatal movement. These outcomes suggest that a number of the identified miRNAs have crucial roles in plant responses to environmental stress, as well as in plant development, and might have common regulatory roles in the core ABA signaling pathway.

Identification and molecular characterization of miRNAs and their target genes associated with seed development through small RNA sequencing in chickpea

Multiple studies have attempted to dissect the molecular mechanism underlying seed development in chickpea (Cicer arietinum L.). These studies highlight the need to focus on the role of miRNAs in regulating storage protein accumulation in seeds. Therefore, a total of 8,856,691 short-read sequences were generated from a small RNA library of developing chickpea seeds and were analyzed using miRDeep-P to identify 74 known and 26 novel miRNA sequences. Known miRNAs were classified into 22 miRNA families with miRNA156 family being most abundant. Of the 26 putative novel miRNAs identified, only 22 could be experimentally validated using stem loop end point PCR. Differential expression analyses led to the identification of known as well as novel miRNAs that could regulate various stages of chickpea seed development. In silico target prediction revealed several important target genes and transcription factors like SPL, mediator of RNA Polymerase II transcription subunit 12, aspartic proteinase and NACs, which were further validated by real-time PCR analysis. A comparative expression analysis in chickpea genotypes with contrasting seed protein content revealed one known (Car-miR156h) and two novel miRNA (Car-novmiR7 and Car-novmiR23) candidates to be highly expressed in the LPC (low protein content) chickpea genotypes, targets of which are known to regulate seed storage protein accumulation. Therefore, this study provides a useful resource in the form of miRNA and their targets which can be further utilized to understand and manipulate various regulatory mechanisms involved in seed development with the overall aim of improving yield and nutrition attributes in chickpea.

Emerging Connections between Small RNAs and Phytohormones

Small RNAs (sRNAs), mainly including miRNAs and siRNAs, are ubiquitous in eukaryotes. sRNAs mostly negatively regulate gene expression via (post-)transcriptional gene silencing through DNA methylation, mRNA cleavage, or translation inhibition. The mechanisms of sRNA biogenesis and function in diverse biological processes, as well as the interactions between sRNAs and environmental factors, like (a)biotic stress, have been deeply explored. Phytohormones are central in the plant's response to stress, and multiple recent studies highlight an emerging role for sRNAs in the direct response to, or the regulation of, plant hormonal pathways. In this review, we discuss recent progress on the unraveling of crossregulation between sRNAs and nine plant hormones.

MicroRNA6443-mediated regulation of FERULATE 5-HYDROXYLASE gene alters lignin composition and enhances saccharification in Populus tomentosa

Ferulate 5-hydroxylase (F5H) is a limiting enzyme involved in biosynthesizing sinapyl (S) monolignol in angiosperms. Genetic regulation of F5H can influence S monolignol synthesis and therefore improve saccharification efficiency and biofuel production. To date, little is known about whether F5H is post-transcriptionally regulated by endogenous microRNAs (miRNAs) in woody plants. Here, we report that a microRNA, miR6443, specifically regulates S lignin biosynthesis during stem development in Populus tomentosa. In situ hybridization showed that miR6443 is preferentially expressed in vascular tissues. We further identified that F5H2 is the direct target of miR6443. Overexpression of miR6443 decreased the transcript level of F5H2 in transgenic plants, resulting in a significant reduction in S lignin content. Conversely, reduced miR6443 expression by short tandem target mimics (STTM) elevated F5H2 transcripts, therefore increasing S lignin composition. Introduction of a miR6443-resistant form of F5H2 into miR6443-overexpression plants restored lignin ectopic composition, supporting that miR6443 specifically regulated S lignin biosynthesis by repressing F5H2 in P. tomentosa. Furthermore, saccharification assays revealed decreased hexose yields by 7.5-24.5% in miR6443-overexpression plants compared with the wild-type control, and increased hexoses yields by 13.2-14.6% in STTM6443-overexpression plants. Collectively, we demonstrate that miR6443 modulates S lignin biosynthesis by specially regulating F5H2 in P. tomentosa.

Improved genome assembly provides new insights into genome evolution in a desert poplar (Populus euphratica)

Populus euphratica is well adapted to extreme desert environments and is an important model species for elucidating the mechanisms of abiotic stress resistance in trees. The current assembly of P. euphratica genome is highly fragmented with many gaps and errors, thereby impeding downstream applications. Here, we report an improved chromosome-level reference genome of P. euphratica (v2.0) using single-molecule sequencing and chromosome conformation capture (Hi-C) technologies. Relative to the previous reference genome, our assembly represents a nearly 60-fold improvement in contiguity, with a scaffold N50 size of 28.59 Mb. Using this genome, we have found that extensive expansion of Gypsy elements in P. euphratica led to its rapid increase in genome size compared to any other Salicaceae species studied to date, and potentially contributed to adaptive divergence driven by insertions near genes involved in stress tolerance. We also detected a wide range of unique structural rearrangements in P. euphratica, including 2,549 translocations, 454 inversions, 121 tandem and 14 segmental duplications. Several key genes likely to be involved in tolerance to abiotic stress were identified within these regions. This high-quality genome represents a valuable resource for poplar breeding and genetic improvement in the future, as well as comparative genomic analysis with other Salicaceae species.

MiR319 mediated salt tolerance by ethylene

Salinity-induced accumulation of certain microRNAs accompanied by gaseous phytohormone ethylene production has been recognized as a mechanism of plant salt tolerance. MicroRNA319 (miR319) has been characterized as an important player in abiotic stress resistance in some C3 plants, such as Arabidopsis thaliana and rice. However, its role in the dedicated biomass plant switchgrass (Panicum virgatum L.), a C4 plant, has not been reported. Here, we show crosstalk between miR319 and ethylene (ET) for increasing salt tolerance. By overexpressing Osa-MIR319b and a target mimicry form of miR319 (MIM319), we showed that miR319 positively regulated ET synthesis and salt tolerance in switchgrass. By experimental treatments, we demonstrated that ET-mediated salt tolerance in switchgrass was dose-dependent, and miR319 regulated the switchgrass salt response by fine-tuning ET synthesis. Further experiments showed that the repression of a miR319 target, PvPCF5, in switchgrass also led to enhanced ethylene accumulation and salt tolerance in transgenic plants. Genome-wide transcriptome analysis demonstrated that overexpression of miR319 (OE-miR319) down-regulated the expression of key genes in the methionine (Met) cycle but promoted the expression of genes in ethylene synthesis. The results enrich our understanding of the synergistic effects of the miR319-PvPCF5 module and ethylene synthesis in the salt tolerance of switchgrass, a C4 bioenergy plant.

Exploration of miRNA-mediated fertility regulation network of cytoplasmic male sterility during flower bud development in soybean

Cytoplasmic male sterility (CMS) plays an important role in the production of soybean hybrid seeds. MicroRNAs (miRNAs) are a class of non-coding endogenous ~ 21 nt small RNAs that play crucial roles in flower and pollen development by targeting genes in plants. To dissect the function of miRNAs in soybean CMS, a total of 558 known miRNAs, 10 novel miRNAs, and 466 target genes were identified in flower buds of the soybean CMS line NJCMS1A and its restorer line NJCMS1C through small RNA sequencing and degradome analysis. In addition, miRNA-mediated editing events were also observed, and the two most frequently observed editing types (A → G and C → U) were validated by cloning and sequencing. And as the base editing occurred, some targets were filtered, such as gma-miR2118b-P6GT with Glyma.08G122000.2. Further integrated analysis of transcriptome and small RNA found some miRNAs and their targets' expression patterns showing a negative correlation, such as gma-miR156b/GmSPL9a and gma-miR4413b/GmPPR. Furthermore, opposite expression pattern was observed between gma-miR156b and GmSPL9 during early stage of flower bud development. Taken together, the regulatory network of gma-miR156b/GmSPL9 and gma-miR4413b/GmPPR with flower bud development in soybean CMS was developed. Most importantly, previous reports showed that these targets might be related to pollen development and male sterility, suggesting that both conserved and species-specific miRNAs might act as regulators of flower bud development in soybean CMS. These findings may provide a better understanding of the miRNA-mediated regulatory networks of CMS mechanisms in soybean.

Exploration of ABA Responsive miRNAs Reveals a New Hormone Signaling Crosstalk Pathway Regulating Root Growth of Populus euphratica

Abscisic acid (ABA) plays an important role in the regulation of plant adaptation, seed germination, and root development in plants. However, the mechanism of ABA regulation of root development is still poorly understood, especially through the miRNA-mediated pathway. Here, small RNA (sRNA)-seq and degradome-seq were used to analyze the miRNAs’ responsive to ABA in the stems and roots of P. euphratica, a model tree species for abiotic stress-resistance research. In total, 255 unique mature sequences, containing 154 known miRNAs and 101 novel miRNAs were identified, among which 33 miRNAs and 54 miRNAs were responsive to ABA in the roots and stems, respectively. Furthermore, the analysis of these miRNAs and their targets revealed a new hormone signaling crosstalk model of ABA regulation of root growth through miRNA-mediated pathways, such as peu-miR-n68 mediation of the crosstalk between ABA and the brassinosteroid (BR) signaling pathway and peu-miR477b mediation of the crosstalk between ABA and Gibberellic acid (GA) signaling. Taken together, our genome-wide analysis of the miRNAs provides a new insight into the mechanism of ABA regulation of root growth in Populus.

Plant microRNAs: Front line players against invading pathogens

Plants are attacked by a large number of pathogens. To defend against these pathogens, plants activate or repress a vast array of genes. For genetic expression and reprogramming, host endogenous small RNAs (sRNAs) are the key factors. Among these sRNAs, microRNAs (miRNAs) mediate gene regulation through RNA silencing at the post-transcriptional level and play an essential role in the defense responses to biotic and abiotic stress. In the recent years, high-throughput sequencing has enabled the researchers to uncover the role of plant miRNAs during pathogen invasion. So here we have reviewed the recent research findings illustrating the plant miRNAs active involvement in various defense processes during fungal, bacterial, viral and nematode infections. However, rapid validation of direct targets of miRNAs is the dire need of time, which can be very helpful in improving the plant resistance against various pathogenic diseases.

Characterization of Conserved and Novel microRNAs in Lilium lancifolium Thunb. by High-Throughput Sequencing

MicroRNAs (miRNAs) are among the class of noncoding small RNA molecules and play a crucial role in post-transcriptional regulation in plants. Although Lilium is one of the most popular ornamental flowers worldwide, however, there is no report on miRNAs identification. In the present study, therefore, miRNAs and their targets were identified from flower, leaf, bulblet and bulb of Lilium lancifolium Thunb. by high-throughput sequencing and bioinformatics analysis. In this study, a total of 38 conserved miRNAs belonging to 17 miRNA families and 44 novel miRNAs were identified. In total, 366 target genes for conserved miRNAs and 415 target genes for novel miRNAs were predicted. The majority of the target genes for conserved miRNAs were transcriptional factors and novel miRNAs targeted mainly protein coding genes. A total of 53 cleavage sites belonging to 6 conserved miRNAs families and 14 novel miRNAs were identified using degradome sequencing. Twenty-three miRNAs were randomly selected, then, their credibility was confirmed using northern blot or stem-loop qRT-PCR. The results from qRT-PCR analysis showed the expression pattern of 4 LL-miRNAs was opposite to their targets. Therefore, our finding provides an important basis to understand the biological functions of miRNAs in Lilium.

MicroRNA expression changes following synthesis of three full-sib Populus triploid populations with different heterozygosities

Through high-throughput sequencing, we compared the relative expression levels of miRNA in three full-sib Populus triploid populations with that in their parents and one diploid hybrid population. We found similar numbers of miRNAs differentially expressed between the parents and the four progeny hybrid populations. In addition, unbalanced parental expression level dominance of miRNAs were found in the three allotriploid and interspecific hybrid populations, which may reprogram gene expression networks and contribute to the growth of Populus hybrids. These results indicated that hybridization has a great impact on the miRNA expression variation in the newly synthesized Populus triploid and diploid hybrid populations. However, we also found no significant differences in miRNA expression among one diploid and three triploid hybrid populations, hinting that miRNA abundances do not increase with the genome content. No dosage effect of miRNA expression could lead to dosage-dependent negative effects on target genes and their downstream pathway in polyploids. We speculate that polyploids may gain advantages from the slight decrease in miRNA regulation, suggesting an important molecular mechanism of polyploid advantage. Hybridization with three types of induced 2n gametes transmitted different parental heterozygosities has been proven as an efficient method for Populus triploid production. Several researches have shown that miRNA could be non-additively expressed in allopolyploids. However, it is still unclear whether the non-additively expressed miRNAs result from the effect of hybridization or polyploidization, and whether a dose response to the additional genomic content exists for the expression of miRNA. Toward this end, through high-throughput sequencing, we compared the expression levels of miRNA in three full-sib Populus triploid populations with that in their parents and one interspecific hybrid population. We found similar numbers of miRNAs differentially expressed between the parents and the four progeny hybrid populations. Unbalanced parental expression level dominance of miRNAs were found in the three triploid and diploid hybrid populations, which may reprogram gene expression networks and affect the growth of Populus hybrids. These results indicated that hybridization has a great impact on the miRNA expression variation in the newly synthesized Populus triploid and diploid hybrid populations. However, we also found no significant differences in miRNA expression among the three triploid populations and the diploid hybrid population. No dosage effect of miRNA expression could lead to dosage-dependent negative effects on target genes and their downstream pathway in polyploids. We speculate that polyploids may gain advantages from the decrease in miRNA negative regulation, suggesting an important molecular mechanism of polyploid advantage.

Physiological and transcriptome analysis of heteromorphic leaves and hydrophilic roots in response to soil drying in desert Populus euphratica

Populus euphratica Olivier, which has been considered as a tree model for the study of higher plant response to abiotic stresses, survive in the desert ecosystem characterized by extreme drought stress. To survive in the harsh environmental condition the plant species have developed some plasticity such as the development of heteromorphic leaves and well-developed roots system. We investigated the physiological and molecular mechanisms enabling this species to cope with severe stress caused by drought. The heterophylly, evolved from linear to toothed-ovate shape, showed the significant difference in cuticle thickness, stomata densities, and sizes. Physiological parameters, SOD, POD, PPO, CAT activity, free proline, soluble protein and MDA contents fluctuated in response to soil drying. Gene expression profile of roots monitored at control and 4 moisture gradients regimes showed the up-regulation of 124, 130, 126 and 162 and down-regulation of 138, 251, 314, 168 DEGs, respectively. Xyloglucan endotransglucosylase/ hydrolase gene (XET) up-regulated at different moisture gradients, was cloned and expressed in tobacco. The XET promoter sequence harbors the drought signaling responsive cis-elements. The promoter expression activity varies in different organs. Over-expression and knocked down transgenic tobacco plant analysis confirmed the role of XET gene in roots growth and drought resistance.

The role of peu-miR164 and its target PeNAC genes in response to abiotic stress in Populus euphratica

Plant miR164 family is highly conserved and miR164 members regulate conserved targets belonging to NAC transcription factors. Our previous studies have revealed that peu-miR164a-e and its target gene POPTR_0007s08420 participate in abiotic stress response in Populus euphratica according to deep sequencing and degradome sequencing. In this study, miR164 family comprises six members that generate two mature products (miR164a-e and miR164f) and target seven NAC genes in P. euphratica. Co-expression in Nicotiana benthamiana and 5' RACE confirmed that peu-miR164 directs PeNAC070, PeNAC012 and PeNAC028 mRNAs cleavage. Expression profiles of primary peu-miR164 a/b/c/d/e bear similarity to those of peu-miR164a-e, whereas PeNAC070 and PeNAC081 showed inverse expression patterns with peu-miR164a-e under abiotic stresses. Existence of cis-acting elements in PeNAC070 promoter (ABRE,MBs, Box-W1, GC-motif, and W-box) and in peu-MIR164b promoter (HSE) further confirmed different responses of peu-miR164 and PeNAC070 to abiotic stresses. Histochemical β-glucuronidase (GUS) staining revealed that GUS activities increased when ProPeNAC070::GUS transgenic Arabidopsis plants were exposed to NaCl, mannitol and abscisic acid (ABA), whereas GUS activity of Propeu-MIR164b::GUS plants decreased under ABA treatment. Subcellular localization and transactivation assays showed that PeNAC070 protein was localized to the nucleus and exhibited transactivation activity at the C-terminal. Overexpression of PeNAC070 in Arabidopsis promoted lateral root development, delayed stem elongation, and increased sensitivity of transgenic plants to drought and salt stresses. This study aids in understanding the adaptability of P. euphratica to extreme drought and salt environment by analysing tissue-specific expression patterns of miR164-regulated and specific promoter-regulated PeNAC genes.

MicroRNA156 improves drought stress tolerance in alfalfa (Medicago sativa) by silencing SPL13

Alfalfa (Medicago sativa) is an important forage crop that is often grown in areas that frequently experience drought and water shortage. MicroRNA156 (miR156) is an emerging tool for improving various traits in plants. We tested the role of miR156d in drought response of alfalfa, and observed a significant improvement in drought tolerance of miR156 overexpression (miR156OE) alfalfa genotypes compared to the wild type control (WT). In addition to higher survival and reduced water loss, miR156OE genotypes also maintained higher stomatal conductance compared to WT during drought stress. Furthermore, we observed an enhanced accumulation of compatible solute (proline) and increased levels of abscisic acid (ABA) and antioxidants in miR156OE genotypes. Similarly, alfalfa plants with reduced expression of miR156-targeted SPL13 showed reduced water loss and enhanced stomatal conductance, chlorophyll content and photosynthetic assimilation. Several genes known to be involved in drought tolerance were differentially expressed in leaf and root of miR156 overexpression plants. Taken together, our findings reveal that miR156 improves drought tolerance in alfalfa at least partially by silencing SPL13.

Transcriptome-Wide Analysis of Botrytis elliptica Responsive microRNAs and Their Targets in Lilium Regale Wilson by High-Throughput Sequencing and Degradome Analysis

MicroRNAs, as master regulators of gene expression, have been widely identified and play crucial roles in plant-pathogen interactions. A fatal pathogen, Botrytis elliptica, causes the serious folia disease of lily, which reduces production because of the high susceptibility of most cultivated species. However, the miRNAs related to Botrytis infection of lily, and the miRNA-mediated gene regulatory networks providing resistance to B. elliptica in lily remain largely unexplored. To systematically dissect B. elliptica-responsive miRNAs and their target genes, three small RNA libraries were constructed from the leaves of Lilium regale, a promising Chinese wild Lilium species, which had been subjected to mock B. elliptica treatment or B. elliptica infection for 6 and 24 h. By high-throughput sequencing, 71 known miRNAs belonging to 47 conserved families and 24 novel miRNA were identified, of which 18 miRNAs were downreguleted and 13 were upregulated in response to B. elliptica. Moreover, based on the lily mRNA transcriptome, 22 targets for 9 known and 1 novel miRNAs were identified by the degradome sequencing approach. Most target genes for elliptica-responsive miRNAs were involved in metabolic processes, few encoding different transcription factors, including ELONGATION FACTOR 1 ALPHA (EF1a) and TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR 2 (TCP2). Furthermore, the expression patterns of a set of elliptica-responsive miRNAs and their targets were validated by quantitative real-time PCR. This study represents the first transcriptome-based analysis of miRNAs responsive to B. elliptica and their targets in lily. The results reveal the possible regulatory roles of miRNAs and their targets in B. elliptica interaction, which will extend our understanding of the mechanisms of this disease in lily.