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

NAC61 regulates late- and post-ripening osmotic, oxidative, and biotic stress responses in grapevine

During late- and post-ripening stages, grape berry undergoes profound biochemical and physiological changes whose molecular control is poorly understood. Here, we report the role of NAC61, a grapevine NAC transcription factor, in regulating different processes involved in berry ripening progression. NAC61 is highly expressed during post-harvest berry dehydration and its expression pattern is closely related to sugar concentration. The ectopic expression of NAC61 in Nicotiana benthamiana leaves resulted in low stomatal conductance, high leaf temperature, tissue collapse and a higher relative water content. Transcriptome analysis of grapevine leaves transiently overexpressing NAC61 and DNA affinity purification and sequencing analyses allowed us to narrow down a list of NAC61-regulated genes. Direct regulation of the stilbene synthase regulator MYB14, the osmotic stress-related gene DHN1b, the Botrytis cinerea susceptibility gene WRKY52, and NAC61 itself was validated. We also demonstrate that NAC61 interacts with NAC60, a proposed master regulator of grapevine organ maturation, in the activation of MYB14 and NAC61 expression. Overall, our findings establish NAC61 as a key player in a regulatory network that governs stilbenoid metabolism and osmotic, oxidative, and biotic stress responses that are the hallmark of late- and post-ripening grape stages.

Full-Length Transcriptome Characterization and Functional Analysis of Pathogenesis-Related Proteins in Lilium Oriental Hybrid 'Sorbonne' Infected with Botrytis elliptica

Gray mold (Botrytis elliptica) causes a deleterious fungal disease that decreases the ornamental value and yield of lilies. Lilium oriental hybrid 'Sorbonne' is a variety that is resistant to gray mold. Understanding the mechanism of resistance against B. elliptica infection in 'Sorbonne' can provide a basis for the genetic improvement in lily plants. In this study, a PacBio Sequel II system was used to sequence the full-length transcriptome of Lilium 'Sorbonne' after inoculation with B. elliptica. A total of 46.64 Gb subreads and 19,102 isoforms with an average length of 1598 bp were obtained. A prediction analysis revealed 263 lncRNAs, and 805 transcription factors, 4478 simple sequence repeats, and 17,752 coding sequences were identified. Pathogenesis-related proteins (PR), which may play important roles in resistance against B. elliptica infection, were identified based on the full-length transcriptome data and previously obtained second-generation transcriptome data. Nine non-redundant potential LhSorPR proteins were identified and assigned to two groups that were composed of two LhSorPR4 and seven LhSorPR10 proteins based on their genetic relatedness. The real-time quantitative reverse transcription PCR (qRT-PCR) results showed that the patterns of expression of nine differentially expressed PR genes under B. elliptica stress were basically consistent with the results of transcriptome sequencing. The pattern of expression of LhSorPR4s and LhSorPR10s genes in different tissues was analyzed, and the expression of each gene varied. Furthermore, we verified the function of LhSorPR4-2 gene in Lilium. The expression of LhSorPR4-2 was induced by phytohormones such as methyl jasmonate, salicylic acid, and ethephon. Moreover, the promoter region of LhSorPR4-2 was characterized by several functional domains associated with phytohormones and stress response. The overexpression of LhSorPR4-2 gene in 'Sorbonne' increased the resistance of the lily plant to B. elliptica and correlated with high chitinase activity. This study provides a full-length transcript database and functionally analyzed the resistance of PR gene to B. elliptica in Lilium, thereby introducing the candidate gene LhSorPR4-2 to breed resistance in Lilium.

Degradation of oxalic acid produced by Botrytis elliptica infection in two ploidy levels of Lilium rosthornii Diels

Oxalic acid (OA) is a crucial pathogenic factor for Sclerotinia spp. fungi, which is closely related to Botrytis spp. fungi. Whether OA is a pathogenic factor for the causal agent of grey mould in lily, Botrytis elliptica, and the response of lily to OA are poorly understood. To address these questions, lesion tissues and deposition of calcium oxalate (CaOX) and callose were observed in diploid and tetraploid leaves of L. rosthornii after inoculation with B. elliptica. Oxalate oxidase (OXO) activity and the transcript levels of some genes related to OA degradation (LrGLP1, LrGLP2 and LrWRKY4), reactive oxygen species (ROS) production/scavenging systems (LrRBOHD, LrGST, LrPOD and LrAPX1) and pathogen-related protein (PR) synthesis (LrCHI, LrBGL and LrPR10) were compared. After diploid and tetraploid leaves inoculation, lesion tissue and callose and CaOX were separately observed around in guard cells and stomata rather than the epidermis in the infected area. OXO activity was triggered at 2 h post-inoculation (hpi) in both ploidy leaves, and it was higher in the latter from 12-48 hpi. Expression of LrGLP1, LrGLP2, LrRBOHD, LrGST, LrPOD, LrCHI, LrBGL and LrPR10 was higher in tetraploids than in diploids from 24(12)-36(48) hpi. In conclusion, for B. elliptica, OA mainly chelates Ca²⁺ from the stomata cell wall. The strong capability to degrade OA and higher expression levels of some genes related to ROS accumulation/scavenging and PR synthesis may partially explain the relatively higher grey mould resistance of tetraploid L. rosthornii.

Soybean miR159-GmMYB33 Regulatory Network Involved in Gibberellin-Modulated Resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is an obligate sedentary biotroph that poses major threats to soybean production globally. Recently, multiple miRNAome studies revealed that miRNAs participate in complicated soybean-SCN interactions by regulating their target genes. However, the functional roles of miRNA and target genes regulatory network are still poorly understood. In present study, we firstly investigated the expression patterns of miR159 and targeted GmMYB33 genes. The results showed miR159-3p downregulation during SCN infection; conversely, GmMYB33 genes upregulated. Furthermore, miR159 overexpressing and silencing soybean hairy roots exhibited strong resistance and susceptibility to H. glycines, respectively. In particular, miR159-GAMYB genes are reported to be involve in GA signaling and metabolism. Therefore, we then investigated the effects of GA application on the expression of miR159-GAMYB module and the development of H. glycines. We found that GA directly controls the miR159-GAMYB module, and exogenous GA application enhanced endogenous biologically active GA1 and GA3, the abundance of miR159, lowered the expression of GmMYB33 genes and delayed the development of H. glycines. Moreover, SCN infection also results in endogenous GA content decreased in soybean roots. In summary, the soybean miR159-GmMYB33 module was directly involved in the GA-modulated soybean resistance to H. glycines.

Molecular mechanism of modulating miR482b level in tomato with botrytis cinerea infection

BACKGROUND

Plant miRNAs are involved in the response to biotic and abiotic stresses by altering their expression levels, and they play an important role in the regulation of plant resistance to stress. However, the molecular mechanism that regulates the expression levels of miRNAs in plants with biotic and abiotic stress still needs to be explored. Previously, we found that the expression of the miR482 family was changed in tomato infected by Botrytis cinerea. In this study, we investigated and uncovered the mechanism underlying the response of miR482 to B. cinerea infection in tomato.

RESULTS

First, RT-qPCR was employed to detect the expression patterns of miR482b in tomato infected by B. cinerea, and results showed that miR482b primary transcripts (pri-miR482b) were up-regulated in B. cinerea-infected leaves, but the mature miR482b was down-regulated. Subsequently, we used rapid amplification cDNA end method to amplify the full-length of pri-miR482b. Result showed that the pri-miR482b had two isoforms, with the longer one (consisting 300 bp) having an extra fragment of 53 bp in the 3'-end compared with the shorter one. In vitro Dicer assay indicated that the longer isoform pri-miR482b-x1 had higher efficiency in the post-transcriptional splicing of miRNA than the shorter isoform pri-miR482b-x2. In addition, the transcription level of mature miR482b was much higher in transgenic Arabidopsis overexpressing pri-miR482b-x1 than that in OE pri-miR482b-x2 Arabidopsis. These results confirmed that this extra 53 bp in pri-miR482b-x1 might play a key role in the miR482b biogenesis of post-transcription processing.

CONCLUSIONS

Extra 53 bp in pri-miR482b-x1 enhanced miR482b biogenesis, which elevated the transcription level of miR482b. This study clarified the response of miR482 to B. cinerea infection in tomato, thereby helping us further understand the molecular mechanisms that regulate the expression levels of other miRNAs.

Integrated mRNA and small RNA sequencing reveals a regulatory network associated with flower color in oriental hybrid lily

Anthocyanins are one of the main components of pigments, that are responsible for a wide range of colors in plants. To clarify the regulatory mechanism of anthocyanin biosynthesis in oriental hybrid lily, UPLC/MS analysis was performed to identify the pigments in two cultivars (white and pink). Four major anthocyanins were identified in pink cultivar, and no anthocyanins were detected in white cultivar. Transcriptome and small RNA sequencing (sRNAseq) analyses were performed using tepal tissues at two floral developmental stages from the two cultivars. In total, 55,698 transcripts were assembled, among which 233 were annotated as putative anthocyanin-related transcripts. Differential expression analysis and qRT-PCR results confirmed that most of the anthocyanin-related structural genes had higher expression levels in pink cultivar than in white cultivar. Conversely, LhANR showed a significantly high expression level in white cultivar. Annotated transcription factors (TFs), including MYB activators, MYB repressors and bHLHs, that putatively inhibit or enhance the expression of anthocyanin-related genes were identified. LhMYBA1, an anthocyanin activator, was isolated, and its heterologous expression resulted in a remarkably high level of anthocyanin accumulation. Additionally, 73 differentially expressed microRNAs (miRNAs), including 23 known miRNAs, were detected through sRNAseq. The miRNA target prediction showed that several anthocyanin-related genes might be targeted by miRNAs. Expression profile analysis revealed that these miRNAs showed higher expression levels at later floral developmental stages in white cultivar than in pink cultivar. The results indicated that anthocyanin deficiency in white cultivar might be influenced by multiple levels of suppressive mechanisms, including mRNAs and sRNAs.

Metabolic and Transcriptomic Profiling of Lilium Leaves Infected With Botrytis elliptica Reveals Different Stages of Plant Defense Mechanisms

Botrytis elliptica, the causal agent of gray mold disease, poses a major threat to commercial Lilium production, limiting its ornamental value and yield. The molecular and metabolic regulation mechanisms of Lilium's defense response to B. elliptica infection have not been completely elucidated. Here, we performed transcriptomic and metabolomic analyses of B. elliptica resistant Lilium oriental hybrid "Sorbonne" to understand the molecular basis of gray mold disease resistance in gray mold disease. A total of 115 differentially accumulated metabolites (DAMs) were detected by comparing the different temporal stages of pathogen infection. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed the differentially expressed genes (DEGs) and DAMs were enriched in the phenylpropanoid and flavonoid pathways at all stages of infection, demonstrating the prominence of these pathways in the defense response of "Sorbonne" to B. elliptica. Network analysis revealed high interconnectivity of the induced defense response. Furthermore, time-course analysis of the transcriptome and a weighted gene coexpression network analysis (WGCNA) led to the identification of a number of hub genes at different stages, revealing that jasmonic acid (JA), salicylic acid (SA), brassinolide (BR), and calcium ions (Ca²⁺) play a crucial role in the response of "Sorbonne" to fungal infection. Our work provides a comprehensive perspective on the defense response of Lilium to B. elliptica infection, along with a potential transcriptional regulatory network underlying the defense response, thereby offering gene candidates for resistance breeding and metabolic engineering of Lilium.

The lre-miR159a-LrGAMYB pathway mediates resistance to grey mould infection in Lilium regale

Grey mould is one of the most determinative factors of lily growth and plays a major role in limiting lily productivity. MicroRNA159 (miR159) is a highly conserved microRNA in plants, and participates in the regulation of plant development and stress responses. Our previous studies revealed that lre-miR159a participates in the response of Lilium regale to Botrytis elliptica according to deep sequencing analyses; however, the response mechanism remains unknown. Here, lre-miR159a and its target LrGAMYB gene were isolated from L. regale. Transgenic Arabidopsis overexpressing lre-MIR159a exhibited larger leaves and smaller necrotic spots on inoculation with Botrytis than those of wild-type and overexpressing LrGAMYB plants. The lre-MIR159a overexpression also led to repressed expression of two targets of miR159, AtMYB33 and AtMYB65, and enhanced accumulation of hormone-related genes, including AtPR1, AtPR2, AtNPR1, AtPDF1.2, and AtLOX for both the jasmonic acid and salicylic acid pathways. Moreover, lower levels of H₂ O₂ and O2- were observed in lre-MIR159a transgenic Arabidopsis, which reduced the damage from reactive oxygen species accumulation. Taken together, these results indicate that lre-miR159a positively regulates resistance to grey mould by repressing the expression of its target LrGAMYB gene and activating a defence response.

Comparative transcriptome profiling uncovers a Lilium regale NAC transcription factor, LrNAC35, contributing to defence response against cucumber mosaic virus and tobacco mosaic virus

Cucumber mosaic virus (CMV) is a highly prevalent viral pathogen causing substantial damage to the bulb and cut-flower production of Lilium spp. Here, we performed an Illumina RNA sequencing (RNA-Seq) study on the leaf tissues of a virus-resistant species Lilium regale inoculated with mock control and CMV. A total of 1346 differentially expressed genes (DEGs) were identified in the leaves of L. regale upon CMV inoculation, which contained 34 up-regulated and 40 down-regulated DEGs that encode putative transcription factors (TFs). One up-regulated TF, LrNAC35, belonging to the NAM/ATAF/CUC (NAC) superfamily, was selected for further functional characterization. Aside from CMV, lily mottle virus and lily symptomless virus infections provoked a striking increase in LrNAC35 transcripts in both resistant and susceptible Lilium species. The treatments with low temperature and several stress-related hormones activated LrNAC35 expression, contrary to its reduced expression under salt stress. Ectopic overexpression of LrNAC35 in petunia (Petunia hybrida) resulted in reduced susceptibility to CMV and Tobacco mosaic virus infections, and enhanced accumulation of lignin in the cell walls. Four lignin biosynthetic genes, including PhC4H, Ph4CL, PhHCT and PhCCR, were found to be up-regulated in CMV-infected petunia lines overexpressing LrNAC35. In vivo promoter-binding tests showed that LrNAC35 specifically regulated the expression of Ph4CL. Taken together, our results suggest a positive role of transcriptome-derived LrNAC35 in transcriptional modulation of host defence against viral attack.

Identification and characterization of known and novel microRNAs in strawberry fruits induced by Botrytis cinerea

MicroRNAs are endogenous small non-coding RNAs that negatively regulate mRNAs, mainly at the post-transcriptional level, and play an important role in resistance response of plants. To date, there are few reports on resistance response of strawberry miRNAs to pathogens. In this study, using high-throughput sequencing, 134 conserved and 35 novel miRNAs were identified in six libraries within the treatment of Botrytis cinerea. A total 497 potential target genes were predicted using Fragaria vesca genome. Most of the differential expressed miRNAs in strawberry fruits were up-regulated in early libraries and down-regulated in late libraries. PIRL, the target gene of miR5290a, showed the opposite expressed trend compared with miR5290 from T1 to T3 libraries, and functional analysis of the PIRL gene shows that it has obvious resistance to B. cinerea in the strawberry fruits with overexpressed PIRL gene. We speculate that miR5290a negatively regulates its target gene PIRL to increase resistance to pathogen infection, and further analysis of PIRL function is meaningful for studying the plant-pathogen relationship and improving strawberry fruit quality and yield.

Analysis of WRKY transcription factors and characterization of two Botrytis cinerea-responsive LrWRKY genes from Lilium regale

A major constraint in producing lilies is gray mold caused by Botrytis elliptica and B. cinerea. WRKY transcription factors play important roles in plant immune responses. However, limited information is available about the WRKY gene family in lily plants. In this study, 23 LrWRKY genes with complete WRKY domains were identified from the Botrytis-resistant species Lilium regale. The putative WRKY genes were divided into seven subgroups (Group I, IIa-e, and III) according to their structural features. Sequence alignment revealed that LrWRKY proteins have a highly conserved WRKYGQK domain and a variant, the WRKYGKK domain, and these proteins generally contained similar motif compositions throughout the same subgroup. Functional annotation predicted they might be involved in biological processes related to abiotic and biotic stresses. A qRT-PCR analysis confirmed that expression of six LrWRKY genes in L. regale or the susceptible Asian hybrid 'Yale' was induced by B. cinerea infection. Among these genes, LrWRKY4, LrWRKY8 and LrWRKY10 were expressed at a higher level in L. regale than 'Yale', while the expression of LrWRKY6 and LrWRKY12 was lower in L. regale. Furthermore, LrWRKY4 and LrWRKY12 genes, which also respond to salicylic acid (SA) and methyl jasmonate (MeJA) treatments, were isolated from L. regale. Subcellular localization analysis determined that they were targeted to the nucleus. Constitutive expression of LrWRKY4 and LrWRKY12 in Arabidopsis resulted in plants that were more resistant to B. cinerea than wild-type plants. This resistance was coupled with the transcriptional changes of SA and JA-responsive genes. Overall, our study provides valuable information about the structural and functional characterization of LrWRKY genes that will not only deepen our understanding of the molecular mechanisms underlying the defense of lily against B. cinerea but also offer potential targets for cultivar improvement via biotechnology.

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.