The role of small RNAs in resistant melon cultivar against Phelipanche aegyptiaca parasitization.

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'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's resistance to broomrape.

Identification of microRNAs involved in ectomycorrhizal formation in Populus tomentosa.

The majority of woody plants are able to form ectomycorrhizal (ECM) symbioses with fungi. During symbiotic development, plants undergo a complex re-programming process involving a series of physiological and morphological changes. MicroRNAs (miRNAs) are important components of the regulatory network underlying symbiotic development. To elucidate the mechanisms of miRNAs and miRNA-mediated mRNA cleavage during symbiotic development, we conducted high-throughput sequencing of small RNAs and degradome tags from roots of Populus tomentosa inoculated with Cenococcum geophilum. This process led to the annotation of 51 differentially expressed miRNAs between non-mycorrhizal and mycorrhizal roots of P. tomentosa, including 13 novel miRNAs. Increased or decreased accumulation of several novel and conserved miRNAs in ECM roots, including miR162, miR164, miR319, miR396, miR397, miR398, novel-miR44 and novel-miR47, suggests essential roles for these miRNAs in ECM formation. The degradome analysis identified root transcripts as miRNA-mediated mRNA cleavage targets, which was confirmed using real-time quantitative PCR. Several of the identified miRNAs and corresponding targets are involved in arbuscular mycorrhizal symbioses. In summary, increased or decreased accumulation of specific miRNAs and miRNA-mediated cleavage of symbiosis-related genes indicate that miRNAs play important roles in the regulatory network underlying symbiotic development.

Integrated analysis of smRNAome, transcriptome, and degradome data to decipher microRNAs regulating costunolide biosynthesis in Saussurea lappa.

Saussurea lappa (S. lappa) has been known to synthesize medicinally important, costunolide. Due to its immense therapeutic importance, understanding of regulatory mechanism associated with its biosynthesis is crucial. The identification of genes and transcription factors (TFs) in S. lappa, created a clear picture of costunolide biosynthesis pathways. Further to understand the regulation of costunolide biosynthesis by miRNAs, an integrated study of transcriptome, miRNAs, and degradome was performed. Identified candidate miRNAs and associated feed-forward loops (FFLs) illustrates their regulatory role in secondary metabolite biosynthesis. Small RNA and degradome sequencing were performed for leaf and root tissues to determine miRNAs-targets pairs. A total of 711 and 525 such targets were obtained for novel and known miRNAs respectively. This data was used to generate costunolide-specific miRNA-TF-gene interactome to perform systematic analyses through graph theoretical approach. Interestingly, miR171c.1 and sla-miR121 were identified as key regulators to connect and co-regulate both mevalonate and sesquiterpenoid pathways to bio-synthesize costunolide. Tissue-specific FFLs were identified to be involved in costunolide biosynthesis which further suggests the evolutionary co-relation of root-specific networks in synthesis of secondary metabolites in addition to leaf-specific networks. This integrative approach allowed us to determine candidate miRNAs and associated tissue-specific motifs involved in the diversification of secondary metabolites. MiRNAs identified in present study can provide alternatives for bioengineering tool to enhance the synthesis of costunolide and other secondary metabolites in S. lappa.

Novel insights on genes and pathways involved in Pinus elliottii response to resinosis.

Pinus elliottii, an important coniferous timber species, has recently become one of the most popular sources of resin in China. Resinosis is a common disease that may negatively affect pine tree growth and production. In this study, we used single-molecule real-time sequencing and Illumina RNA sequencing to generate an accurate transcriptome for P. elliottii. The transcriptome included 90,026 transcripts, 5160 long non-coding RNAs and 7710 transcription factors. We then analyzed RNA-sequencing, small RNA-sequencing and degradome data to identify genes, miRNAs and key miRNA-target pairs involved in response to resinosis in P. elliottii. We identified 1305 genes and 1151 miRNAs exhibiting significant differential expression in response to resinosis. According to the degradome sequencing analysis, 318 differentially expressed transcripts were targets of 14 differentially expressed miRNAs. Our study has provided resources for further functional characterization of genes and miRNAs involved in resinosis in P. elliottii, which should aid the future disease-resistance breeding of this species.

Identification of miRNAs Mediating Seed Storability of Maize during Germination Stage by High-Throughput Sequencing, Transcriptome and Degradome Sequencing.

Seed storability is an important trait for improving grain quality and germplasm conservation, but little is known about the regulatory mechanisms and gene networks involved. MicroRNAs (miRNAs) are small non-coding RNAs regulating the translation and accumulation of their target mRNAs by means of sequence complementarity and have recently emerged as critical regulators of seed germination. Here, we used the germinating embryos of two maize inbred lines with significant differences in seed storability to identify the miRNAs and target genes involved. We identified a total of 218 previously known and 448 novel miRNAs by miRNA sequencing and degradome analysis, of which 27 known and 11 newly predicted miRNAs are differentially expressed in two maize inbred lines, as measured by Gene Ontology (GO) enrichment analysis. We then combined transcriptome sequencing and real-time quantitative polymerase chain reaction (RT-PCR) to screen and confirm six pairs of differentially expressed miRNAs associated with seed storability, along with their negative regulatory target genes. The enrichment analysis suggested that the miRNAs/target gene mediation of seed storability occurs via the ethylene activation signaling pathway, hormone synthesis and signal transduction, as well as plant organ morphogenesis. Our results should help elucidate the mechanisms through which miRNAs are involved in seed storability in maize.

Identification of miRNAs and their target genes associated with improved maize seed vigor induced by gibberellin.

High seed vigor is crucial for agricultural production owing to its potential in high quality and yield of crops and a better understanding of the molecular mechanism associated with maize seed vigor is highly necessary. To better understand the involvement and regulatory mechanism of miRNAs correlated with maize seed vigor, small RNAs and degradome sequencing of two inbred lines Yu537A and Yu82 were performed. A total of 791 mature miRNAs were obtained with different expressions, among of which 505 miRNAs were newly identified and the rest miRNAs have been reported before by comparing the miRNAs with the sequences in miRbase database. Analysis of miRNA families showed maize seeds contain fewer miRNA families and larger miRNA families compared with animals, indicating that functions of miRNAs in maize seeds were more synergistic than animals. Degradome sequencing was used to identify the targets of miRNAs and the results showed a total of 6,196 targets were obtained. Function analysis of differentially expressed miRNAs and targets showed Glycan degradation and galactose metabolism were closely correlated with improved maize seed vigor. These findings provide valuable information to understand the involvement of miRNAs with maize seed vigor and these putative genes will be valuable resources for improving the seed vigor in future maize breeding.

Comparative Analysis of Salt Responsive MicroRNAs in Two Sweetpotato [Ipomoea batatas (L.) Lam.] Cultivars With Different Salt Stress Resistance.

Sweetpotato [Ipomoea batatas (L.) Lam.] is an important food, vegetable and economic crop, but its productivity is remarkably affected by soil salinity. MiRNAs are a class of endogenous non-coding small RNAs that play an important role in plant resistance to salt stress. However, the function of miRNAs still remains largely unknown in sweetpotato under salt stress. Previously, we identified salt-responsive miRNAs in one salt-sensitive sweetpotato cultivar "Xushu 32." In this study, we identified miRNAs in another salt-tolerant cultivar "Xushu 22" by high-throughput deep sequencing and compared the salt-responsive miRNAs between these two cultivars with different salt sensitivity. We identified 687 miRNAs in "Xushu 22," including 514 known miRNAs and 173 novel miRNAs. Among the 759 miRNAs from the two cultivars, 72 and 109 miRNAs were specifically expressed in "Xushu 32" and "Xushu 22," respectively, and 578 miRNAs were co-expressed. The comparison of "Xushu 32" and "Xushu 22" genotypes showed a total of 235 miRNAs with obvious differential expression and 177 salt-responsive miRNAs that were obviously differently expressed between "Xushu 32" and "Xushu 22" under salt stress. The target genes of the miRNAs were predicted and identified using the Target Finder tool and degradome sequencing. The results showed that most of the targets were transcription factors and proteins related to metabolism and stress response. Gene Ontology analysis revealed that these target genes are involved in key pathways related to salt stress response and secondary redox metabolism. The comparative analysis of salt-responsive miRNAs in sweetpotato cultivars with different salt sensitivity is helpful for understanding the regulatory pattern of miRNA in different sweetpotato genotypes and improving the agronomic traits of sweetpotato by miRNA manipulation in the future.

Methods for detecting RNA degradation intermediates in plants.

RNA degradation is an important process for controlling gene expression and is mediated by decapping / deadenylation-dependent or endonucleolytic cleavage-dependent RNA degradation mechanisms. High-throughput sequencing of RNA degradation intermediates was initially developed in Arabidopsis thaliana and similar RNA degradome sequencing methods were conducted in other eukaryotes. However, interpreting results obtained by these sequencing methods is fragmented, and an overview is needed. Here we review the findings and limitations of these sequencing methods and discuss the missing experiments needed to understand RNA degradation intermediates accurately. This review provides direction for future research on RNA degradation and is a reference for RNA degradome studies in other species.

Laccase Directed Lignification Is One of the Major Processes Associated With the Defense Response Against Pythium ultimum Infection in Apple Roots.

Apple replant disease (ARD), incited by a pathogen complex including Pythium ultimum, causes stunted growth or death of newly planted trees at replant sites. Development and deployment of resistant or tolerant rootstocks offers a cost-effective, ecologically friendly, and durable approach for ARD management. Maximized exploitation of natural resistance requires integrated efforts to identify key regulatory mechanisms underlying resistance traits in apple. In this study, miRNA profiling and degradome sequencing identified major miRNA pathways and candidate genes using six apple rootstock genotypes with contrasting phenotypes to P. ultimum infection. The comprehensive RNA-seq dataset offered an expansive view of post-transcriptional regulation of apple root defense activation in response to infection from P. ultimum. Several pairs of miRNA families and their corresponding targets were identified for their roles in defense response in apple roots, including miR397-laccase, miR398-superoxide dismutase, miR10986-polyphenol oxidase, miR482-resistance genes, and miR160-auxin response factor. Of these families, the genotype-specific expression patterns of miR397 indicated its fundamental role in developing defense response patterns to P. ultimum infection. Combined with other identified copper proteins, the importance of cellular fortification, such as lignification of root tissues by the action of laccase, may critically contribute to genotype-specific resistance traits. Our findings suggest that quick and enhanced lignification of apple roots may significantly impede pathogen penetration and minimize the disruption of effective defense activation in roots of resistant genotypes. The identified target miRNA species and target genes consist of a valuable resource for subsequent functional analysis of their roles during interaction between apple roots and P. ultimum.

Feature selection for RNA cleavage efficiency at specific sites using the LASSO regression model in Arabidopsis thaliana.

BACKGROUND: RNA degradation is important for the regulation of gene expression. Despite the identification of proteins and sequences related to deadenylation-dependent RNA degradation in plants, endonucleolytic cleavage-dependent RNA degradation has not been studied in detail. Here, we developed truncated RNA end sequencing in Arabidopsis thaliana to identify cleavage sites and evaluate the efficiency of cleavage at each site. Although several features are related to RNA cleavage efficiency, the effect of each feature on cleavage efficiency has not been evaluated by considering multiple putative determinants in A. thaliana. RESULTS: Cleavage site information was acquired from a previous study, and cleavage efficiency at the site level (CSsite value), which indicates the number of reads at each cleavage site normalized to RNA abundance, was calculated. To identify features related to cleavage efficiency at the site level, multiple putative determinants (features) were used to perform feature selection using the Least Absolute Shrinkage and Selection Operator (LASSO) regression model. The results indicated that whole RNA features were important for the CSsite value, in addition to features around cleavage sites. Whole RNA features related to the translation process and nucleotide frequency around cleavage sites were major determinants of cleavage efficiency. The results were verified in a model constructed using only sequence features, which showed that the prediction accuracy was similar to that determined using all features including the translation process, suggesting that cleavage efficiency can be predicted using only sequence information. The LASSO regression model was validated in exogenous genes, which showed that the model constructed using only sequence information can predict cleavage efficiency in both endogenous and exogenous genes. CONCLUSIONS: Feature selection using the LASSO regression model in A. thaliana identified 155 features. Correlation coefficients revealed that whole RNA features are important for determining cleavage efficiency in addition to features around the cleavage sites. The LASSO regression model can predict cleavage efficiency in endogenous and exogenous genes using only sequence information. The model revealed the significance of the effect of multiple determinants on cleavage efficiency, suggesting that sequence features are important for RNA degradation mechanisms in A. thaliana.

Degradome sequencing-based identification of phasiRNAs biogenesis pathways in Oryza sativa.

BACKGROUND: The microRNAs(miRNA)-derived secondary phased small interfering RNAs (phasiRNAs) participate in post-transcriptional gene silencing and play important roles in various bio-processes in plants. In rice, two miRNAs, miR2118 and miR2275, were mainly responsible for triggering of 21-nt and 24-nt phasiRNAs biogenesis, respectively. However, relative fewer phasiRNA biogenesis pathways have been discovered in rice compared to other plant species, which limits the comprehensive understanding of phasiRNA biogenesis and the miRNA-derived regulatory network. RESULTS: In this study, we performed a systematical searching for phasiRNA biogenesis pathways in rice. As a result, five novel 21-nt phasiRNA biogenesis pathways and five novel 24-nt phasiRNA biogenesis pathways were identified. Further investigation of their regulatory function revealed that eleven novel phasiRNAs in 21-nt length recognized forty-one target genes. Most of these genes were involved in the growth and development of rice. In addition, five novel 24-nt phasiRNAs targeted to the promoter of an OsCKI1 gene and thereafter resulted in higher level of methylation in panicle, which implied their regulatory function in transcription of OsCKI1,which acted as a regulator of rice development. CONCLUSIONS: These results substantially extended the information of phasiRNA biogenesis pathways and their regulatory function in rice.

Changes in mRNA Degradation Efficiencies under Varying Conditions Are Regulated by Multiple Determinants in Arabidopsis thaliana.

Multiple mechanisms are involved in gene expression, with mRNA degradation being critical for the control of mRNA accumulation. In plants, although some trans-acting factors and motif sequences have been identified in deadenylation-dependent mRNA degradation, endonucleolytic cleavage-dependent mRNA degradation has not been studied in detail. Previously, we developed truncated RNA-end sequencing (TREseq) in Arabidopsis thaliana and detected G-rich sequence motifs around 5' degradation intermediates. However, it remained to be elucidated whether degradation efficiencies of 5' degradation intermediates in A. thaliana vary among growth conditions and developmental stages. To address this issue, we conducted TREseq of cultured cells under heat stress and at three developmental stages (seedlings, expanding leaves and expanded leaves) and compared 5' degradation intermediates data among the samples. Although some 5' degradation intermediates had almost identical degradation efficiencies, others differed among conditions. We focused on the genes and sites whose degradation efficiencies differed. Changes in degradation efficiencies at the gene and site levels revealed an effect on mRNA accumulation in all comparisons. These changes in degradation efficiencies involved multiple determinants, including mRNA length and translation efficiency. These results suggest that several determinants govern the efficiency of mRNA degradation in plants, helping the organism to adapt to varying conditions by controlling mRNA accumulation.

Identification of microRNAs and their Endonucleolytic Cleavaged target mRNAs in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) ranks the third among the most common malignancies globally. It is well known that microRNAs (miRNAs) play vital roles in destabilizing mRNAs and repressing their translations in this disease. However, the mechanism of miRNA-induced mRNA cleavage remains to be investigated. METHOD: In this study, high-throughput small RNA (sRNA) sequencing was utilized to identify and profile miRNAs from six pairs of colorectal cancer tissues (CTs) and adjacent tissues (CNs). Degradome sequencing (DS) was employed to detect the cleaved target genes. The Database for Annotation, Visualization and Integrated Discovery (DAVID) software was used for GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis. RESULTS: In total, 1278 known miRNAs (clustered into 337 families) and 131 novel miRNAs were characterized in the CT and CN libraries, respectively. Of those, 420 known and eight novel miRNAs were defined as differentially expressed miRNAs (DEmiRNAs) by comparing the expression levels observed in the CT and CN libraries. Furthermore, through DS, 9685 and 202 potential target transcripts were characterized as target genes for 268 known and 33 novel miRNAs, respectively. It was further predicted that a total of 264 targeted genes for the 85 DEmiRNAs are involved in proteoglycans in cancer and the AMP-activated protein kinase signaling pathway. After systemic analysis of prognosis-related miRNA targets in those cancer-related signal pathways, we found that two targets ezrin (EZR) and hematopoietic cell-specific Lyn substrate 1 (HCLS1) had the potential prognostic characteristics with CRC regarding over survival (OS) or recurrence. CONCLUSION: In total, we found that endonucleolytic miRNA-directed mRNA cleavage occurs in CRC. A number of potential genes targeted by CRC-related miRNAs were identified and some may have the potential as prognosis markers of CRC. The present findings may lead to an improved better appreciation of the novel interaction mode between miRNAs and target genes in CRC.

Identification and integrated analysis of glyphosate stress-responsive microRNAs, lncRNAs, and mRNAs in rice using genome-wide high-throughput sequencing.

BACKGROUND: Glyphosate has become the most widely used herbicide in the world. Therefore, the development of new varieties of glyphosate-tolerant crops is a research focus of seed companies and researchers. The glyphosate stress-responsive genes were used for the development of genetically modified crops, while only the EPSPS gene has been used currently in the study on glyphosate-tolerance in rice. Therefore, it is essential and crucial to intensify the exploration of glyphosate stress-responsive genes, to not only acquire other glyphosate stress-responsive genes with clean intellectual property rights but also obtain non-transgenic glyphosate-tolerant rice varieties. This study is expected to elucidate the responses of miRNAs, lncRNAs, and mRNAs to glyphosate applications and the potential regulatory mechanisms in response to glyphosate stress in rice. RESULTS: Leaves of the non-transgenic glyphosate-tolerant germplasm CA21 sprayed with 2 mg·ml- 1 glyphosate (GLY) and CA21 plants with no spray (CK) were collected for high-throughput sequencing analysis. A total of 1197 DEGs, 131 DELs, and 52 DEMs were identified in the GLY samples in relation to CK samples. Genes were significantly enriched for various biological processes involved in detoxification of plant response to stress. A total of 385 known miRNAs from 59 miRNA families and 94 novel miRNAs were identified. Degradome analysis led to the identification of 32 target genes, of which, the squamosa promoter-binding-like protein 12 (SPL12) was identified as a target of osa-miR156a_L + 1. The lncRNA-miRNA-mRNA regulatory network consisted of osa-miR156a_L + 1, two transcripts of SPL12 (LOC_Os06g49010.3 and LOC_Os06g49010.5), and 13 lncRNAs (e.g., MSTRG.244.1 and MSTRG.16577.1). CONCLUSION: Large-scale expression changes in coding and noncoding RNA were observed in rice mainly due to its response to glyphosate. SPL12, osa-miR156, and lncRNAs (e.g., MSTRG.244.1 and MSTRG.16577.1) could be a novel ceRNA mechanism in response to glyphosate in rice by regulating transcription and metal ions binding. These findings provide a theoretical basis for breeding glyphosate-tolerant rice varieties and for further research on the biogenesis of glyphosate- tolerance in rice.

Adaptive regulation of virulence genes by microRNA-like RNAs in Valsa mali.

MicroRNAs play important roles in the regulation of gene expression in plants and animals. However, little information is known about the action mechanism and function of fungal microRNA-like RNAs (milRNAs). In this study, combining deep sequencing, molecular and histological assays, milRNAs and their targets in the phytopathogenic fungus Valsa mali were isolated and identified. A critical milRNA, Vm-milR16, was identified to adaptively regulate the expression of virulence genes. Fourteen isolated milRNAs showed high expression abundance. Based on the assessment of a pathogenicity function of these milRNAs, Vm-milR16 was found to be a critical milRNA in V. mali by regulating sucrose non-fermenting 1 (VmSNF1), 4,5-DOPA dioxygenase extradiol (VmDODA), and a hypothetical protein (VmHy1). During V. mali infection, Vm-milR16 is downregulated, while its targets are upregulated. Overexpression of Vm-milR16, but not mutated Vm-milR16, significantly reduces the expression of targets and virulence of V. mali. Furthermore, deletion of VmSNF1, VmDODA and VmHy1 significantly reduce virulence of V. mali. All three targets seem to be essential for oxidative stress response and VmSNF1 is required for expression of pectinase genes during V. mali-host interaction. Our results demonstrate Vm-milRNAs contributing to the infection of V. mali on apple trees by adaptively regulating virulence genes.

High-throughput deep sequencing reveals the important role that microRNAs play in the salt response in sweet potato (Ipomoea batatas L.).

BACKGROUND: MicroRNAs (miRNAs), a class of small regulatory RNAs, have been proven to play important roles in plant growth, development and stress responses. Sweet potato (Ipomoea batatas L.) is an important food and industrial crop that ranks seventh in staple food production. However, the regulatory mechanism of miRNA-mediated abiotic stress response in sweet potato remains unclear. RESULTS: In this study, we employed deep sequencing to identify both conserved and novel miRNAs from salinity-exposed sweet potato cultivars and its untreated control. Twelve small non-coding RNA libraries from NaCl-free (CK) and NaCl-treated (Na150) sweet potato leaves and roots were constructed for salt-responsive miRNA identification in sweet potatoes. A total of 475 known miRNAs (belonging to 66 miRNA families) and 175 novel miRNAs were identified. Among them, 51 (22 known miRNAs and 29 novel miRNAs) were significantly up-regulated and 76 (61 known miRNAs and 15 novel miRNAs) were significantly down-regulated by salinity stress in sweet potato leaves; 13 (12 known miRNAs and 1 novel miRNAs) were significantly up-regulated and 9 (7 known miRNAs and 2 novel miRNAs) were significantly down-regulated in sweet potato roots. Furthermore, 636 target genes of 314 miRNAs were validated by degradome sequencing. Deep sequencing results confirmed by qRT-PCR experiments indicated that the expression of most miRNAs exhibit a negative correlation with the expression of their targets under salt stress. CONCLUSIONS: This study provides insights into the regulatory mechanism of miRNA-mediated salt response and molecular breeding of sweet potatoes though miRNA manipulation.

Identification of milRNAs and their target genes in Ganoderma lucidum by high-throughput sequencing and degradome analysis.

MicroRNAs (miRNAs in animals and plants or milRNAs in fungi) are endogenous noncoding RNAs that can regulate gene expression. However, little information is known about milRNAs and their target genes in Ganoderma lucidum. Here, we systematically predicted and characterised the milRNAs and their target genes across the three developmental stages of G. lucidum. A total of 168 unique milRNAs were predicted using a small RNA sequencing method. For them, 1612 target sequences corresponding to 1311 unique genes were predicted by degradome sequencing. We selected 42 predicted milRNAs and performed RT-PCR amplification and Sanger sequencing of the products. Five products were found to have sequences similar to those predicted, confirming the presence of milRNAs in G. lucidum, and demonstrating the difficulty in their validation. Among the 168 milRNAs, 111 were found to be significantly differentially expressed across the three developmental stages (q ≤ 0.05). The expression levels of 12 milRNAs were measured by stem-loop quantitative real-time polymerase chain reaction. Eight of them were in line with the sequencing results (r ≥ 0.9, p ≤ 0.05). These 12 milRNAs and their target genes form 16 milRNA-target gene pairs. The expression profiles of 8 of these 16 miRNA-target pairs were negatively correlated, according to real-time quantitative analysis, whereas the other eight pairs were positively correlated. Furthermore, the results of functional enrichment analysis showed that the target genes of milRNAs mapped to the Gene Ontology terms 'GTP binding' and 'FAD binding' were enriched in specific developmental stages. These target genes were related to the biosynthesis of triterpenes and polysaccharides and lignin degradation pathway in G. lucidum. In summary, this study indicates that milRNAs may play crucial regulatory roles in various biological processes of G. lucidum and open up new avenues for research on milRNAs' biosyntheses and functions in basidiomycetes.

Different Plant Species Have Common Sequence Features Related to mRNA Degradation Intermediates.

mRNA degradation is an important cellular mechanism involved in the control of gene expression. Several genome-wide profiling methods have been developed for detecting mRNA degradation in plants and animals. However, because many of these techniques use poly (A) mRNA for library preparation, degradation intermediates are often only detected near the 3'-ends of transcripts. Previously, we developed the Truncated RNA End Sequencing (TREseq) method using Arabidopsis thaliana, and demonstrated that this method ameliorates 3'-end bias. In analyses using TREseq, we observed G-rich sequences near the 5'-ends of degradation intermediates. However, this finding remained to be confirmed in other plant species. Hence, in this study, we conducted TREseq analyses in Lactuca sativa (lettuce), Oryza sativa (rice) and Rosa hybrida (rose). These species including A. thaliana were selected to encompass a diverse range in the angiosperm phylogeny. The results revealed similar sequence features near the 5'-ends of degradation intermediates, and involvement of translation process in all four species. In addition, homologous genes have similar efficiencies of mRNA degradation in different plants, suggesting that similar mechanisms of mRNA degradation are conserved across plant species. These strong sequence features were not observed in previous degradome analyses among different species in plants.

Identification of miRNAs Involved in Bacillus velezensis FZB42-Activated Induced Systemic Resistance in Maize.

Bacillus velezensis FZB42 is able to activate induced systemic resistance (ISR) to enhance plant defense response against pathogen infections. Though the roles of microRNAs (miRNAs) in Bacillus-triggered ISR have been reported in Arabidopsis, the maize miRNAs responsible for the Bacillus-activated ISR process have not been discovered. To explore the maize miRNAs involved in ISR, maize miRNAs in response to FZB42 (ISR activating), FZB42△sfp△alss (deficient in triggering ISR), and a control for 12 h were sequenced. A total of 146 known miRNAs belonging to 30 miRNA families and 217 novel miRNAs were identified. Four miRNAs specifically repressed in FZB42-treatment were selected as candidate ISR-associated miRNAs. All of them contained at least one defense response-related cis-element, suggesting their potential roles in activating the ISR process. Interestingly, three of the four candidate ISR-associated miRNAs belong to the conserved miR169 family, which has previously been confirmed to play roles in abiotic stress response. Moreover, 52 mRNAs were predicted as potential targets of these candidate ISR-associated miRNAs through TargetFinder software and degradome sequencing. Gene Ontology (GO) and network analyses of target genes showed that these differentially expressed miRNA might participate in the ISR process by regulating nuclear factor Y transcription factor. This study is helpful in better understanding the regulatory roles of maize miRNAs in the Bacillus-activated ISR process.

Identification of miRNAs and Their Target Genes Involved in Cucumber Fruit Expansion Using Small RNA and Degradome Sequencing.

Fruit expansion is an essential and very complex biological process. Regulatory roles of microRNAs (miRNAs) and miRNA-mRNA modules in the cucumber fruit expansion are not yet to be investigated. In this work, 1253 known and 1269 novel miRNAs were identified from nine cucumber fruit small RNA (sRNA) libraries through high-throughput sequencing. A total of 105 highly differentially expressed miRNAs were recognized in the fruit on five days post anthesis with pollination (EXP_5d) sRNA library. Further, expression patterns of 11 differentially expressed miRNAs were validated by quantitative real-time PCR (qRT-PCR). The expression patterns were similar to sRNAs sequencing data. Transcripts of 1155 sequences were predicted as target genes of differentially expressed miRNAs by degradome sequencing. Gene Ontology (GO) enrichment showed that these target genes were involved in 24 biological processes, 15 cell components and nine molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that these target genes were significantly enriched in 19 pathways and the enriched KEGG pathways were associated with environmental adaptation, signal transduction and translation. Based on the functional prediction of miRNAs and target genes, our findings suggest that miRNAs have a potential regulatory role in cucumber fruit expansion by targeting their target genes, which provide important data for understanding the miRNA-mediated regulatory networks controlling fruit expansion in cucumber. Specific miRNAs could be selected for further functional research and molecular breeding in cucumber.