Identification of novel and candidate miRNAs in rice by high throughput sequencing

Identification of miRNAs associated with Aspergillus flavus infection and their targets in groundnut (Arachis hypogaea L.)

BACKGROUND

The quality of groundnut produce is adversely impacted due to aflatoxin contamination by the fungus Aspergillus flavus. Although the transcriptomic control is not fully understood, the interaction between long non-coding RNAs and microRNAs in regulating A. flavus and aflatoxin contamination remains unclear. This study was carried out to identify microRNAs (miRNAs) to enhance the understanding of in vitro seed colonization (IVSC) resistance mechanism in groundnut.

RESULT

In this study, resistant (J 11) and susceptible (JL 24) varieties of groundnut were treated with toxigenic A. flavus (strain AF-11-4), and total RNA was extracted at 1 day after inoculation (1 DAI), 2 DAI, 3 DAI and 7 DAI. Seeds of JL 24 showed higher mycelial growth than J 11 at successive days after inoculation. A total of 208 known miRNAs belonging to 36 miRNA families, with length varying from 20-24 nucleotides, were identified, along with 27 novel miRNAs, with length varying from 20-22 nucleotides. Using psRNATarget server, 952 targets were identified for all the miRNAs. The targeted genes function as disease resistant proteins encoding, auxin responsive proteins, squamosa promoter binding like proteins, transcription factors, pentatricopeptide repeat-containing proteins and growth regulating factors. Through differential expression analysis, seven miRNAs (aly-miR156d-3p, csi-miR1515a, gma-miR396e, mtr-miR2118, novo-miR-n27, ptc-miR482d-3p and ppe-miR396a) were found common among 1 DAI, 2 DAI, 3 DAI and 7 DAI in J 11, whereas ten miRNAs (csi-miR159a-5p, csi-miR164a-3p, novo-miR-n17, novo-miR-n2, osa-miR162b, mtr-miR2118, ptc-miR482d-3p, ptc-miR167f-3p, stu-miR319-3p and zma-miR396b-3p) were found common among 1 DAI, 2 DAI, 3 DAI and 7 DAI in JL 24. Two miRNAs, ptc-miR482d-3p and mtr-miR2118, showed contrasting expression at different time intervals between J 11 and JL 24. These two miRNAs were found to target those genes with NBS-LRR function, making them potential candidates for marker development in groundnut breeding programs aimed at enhancing resistance against A. flavus infection.

CONCLUSION

This study enhances our understanding of the involvement of two miRNAs namely, ptc-miR482d-3p and mtr-miR2118, along with their NBS-LRR targets, in conferring resistance against A. flavus-induced aflatoxin contamination in groundnut under in vitro conditions.

Transcriptome and MicroRNA Analysis of Juglans regia in Response to Low-Temperature Stress

Walnuts are among the globally significant woody food and oil tree species. At high latitudes, they frequently experience late-frost damage, inducing low-temperature stress, which significantly affects walnut seedlings. The aim of this study was to investigate the physiological and biochemical alterations in walnut seedlings under low-temperature (LT) stress along with its underlying molecular mechanisms. Physiological indices were determined, and the transcriptome and miRNA were sequenced by sampling leaves (0 h, 6 h, 12 h, 24 h, and 48 h) of two-month-old live seedlings of walnuts treated with a low temperature of 4 °C. The results indicated that LT stress induced an increase in electrical conductivity and malondialdehyde content while simultaneously causing a reduction in Fv/Fm. From the transcriptome comparison between the control and treated groups, a total of 12,566 differentially expressed genes (DEGs) were identified, consisting of 6829 up-regulated and 5737 down-regulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the DEGs were primarily enriched in polysaccharide metabolic processes, responses to abscisic acid and phenylpropanoid biosynthesis pathways. Furthermore, the miRNA database identified 1052 miRNAs in response to low-temperature stress in walnuts; these miRNAs were found to target 7043 predicted genes. Through the integration and analysis of transcriptome and miRNA data, 244 differential DEGs were identified. Following GO and KEGG enrichment analyses of the differential target genes, we identified that these genes primarily regulate pathways involved in starch and sucrose metabolism, glyoxylate and dicarboxylate metabolism, and glycerophospholipid biosynthesis, as well as phenylalanine, tyrosine, and tryptophan biosynthesis, in walnut leaves under LT stress. Additionally, we conducted an in-depth analysis of the associations between differentially expressed genes (DEGs) and differentially expressed microRNAs (DEMs) within the starch and sucrose metabolism pathway. Real-time fluorescent quantitative PCR (qRT-PCR) validation of the expression patterns of a subset of differential genes confirmed the accuracy of the transcriptome data. This study unveils the potential molecular mechanisms underlying walnut's response to low-temperature stress, providing valuable genetic resources for future research on the cold tolerance mechanisms of walnut in response to late-frost damage.

MicroRNA164 Affects Plant Responses to UV Radiation in Perennial Ryegrass

Increasing the ultraviolet radiation (UV) level, particularly UV-B due to damage to the stratospheric ozone layer by human activities, has huge negative effects on plant and animal metabolism. As a widely grown cool-season forage grass and turfgrass in the world, perennial ryegrass (Lolium perenne) is UV-B-sensitive. To study the effects of miR164, a highly conserved microRNA in plants, on perennial ryegrass under UV stress, both OsmiR164a overexpression (OE164) and target mimicry (MIM164) transgenic perennial ryegrass plants were generated using agrobacterium-mediated transformation, and UV-B treatment (~600 μw cm-2) of 7 days was imposed. Morphological and physiological analysis showed that the miR164 gene affected perennial ryegrass UV tolerance negatively, demonstrated by the more scorching leaves, higher leaf electrolyte leakage, and lower relative water content in OE164 than the WT and MIM164 plants after UV stress. The increased UV sensitivity could be partially due to the reduction in antioxidative capacity and the accumulation of anthocyanins. This study indicated the potential of targeting miR164 and/or its targeted genes for the genetic manipulation of UV responses in forage grasses/turfgrasses; further research to reveal the molecular mechanism underlying how miR164 affects plant UV responses is needed.

Small RNA-Seq to Unveil the miRNA Expression Patterns and Identify the Target Genes in Panax ginseng

Panax ginseng, renowned for its medicinal properties, relies on adventitious roots and hairy roots as crucial sources for the production of ginsenosides. Despite the widespread utilization of ginseng, investigations into its miRNAs have remained scarce. To address this gap, two samples of ginseng adventitious roots and ginseng hairy roots were collected, and subsequent construction and sequencing of small RNA libraries of ginseng adventitious roots and hairy roots were performed using the Illumina HiSeq X Ten platform. The analysis of the sequencing data unveiled total miRNAs 2432. The miR166 and miR396 were the most highly expressed miRNA families in ginseng. The miRNA expression analysis results were used to validate the qRT-PCR. Target genes of miRNA were predicted and GO function annotation and KEGG pathway analysis were performed on target genes. It was found that miRNAs are mainly involved in synthetic pathways and biological processes in plants, which include metabolic and bioregulatory processes. The plant miRNAs enriched KEGG pathways are associated with some metabolism, especially amino acid metabolism and carbohydrate metabolism. These results provide valuable insights miRNAs and their roles in metabolic processes in ginseng.

Molecular mechanisms of miR172a and its target gene LbrTOE3 regulating maturation in Lilium

MAIN CONCLUSION

Lbr-miR172a could promote the growth phase transition and shorten maturation in Lilium, while LbrTOE3 inhibited this process and prolonged the growth period. Lilium is an ornamental flower with high economic value for both food and medicinal purposes. However, under natural conditions, Lilium bulbs take a long time and cost more to grow to commercial size. This research was conducted to shorten the maturation time by subjecting Lilium bulbs to alternating temperature treatment. To explore the molecular mechanism of the vegetative phase change (VPC) in Lilium after variable temperature treatment, the key module miR172a-TOE3 was selected based on a combined omics analysis. Gene cloning and transgene functional validation showed that overexpression of Lbr-mir172a promoted a phase change, while overexpression of LbrTOE3 inhibited this process. Subcellular localization and transcriptional activation assays indicated that LbrTOE3 was predominantly localized in the nucleus and showed transcriptional activity. In situ hybridization showed that LbrTOE3 expression was significantly downregulated after alternating temperature treatment. This study elucidates the molecular mechanisms of the phase transition of Lilium and provides a scientific basis for the phase transition in other plants.

High Daytime Temperature Responsive MicroRNA Profiles in Developing Grains of Rice Varieties with Contrasting Chalkiness

High temperature impairs starch biosynthesis in developing rice grains and thereby increases chalkiness, affecting the grain quality. Genome encoded microRNAs (miRNAs) fine-tune target transcript abundances in a spatio-temporal specific manner, and this mode of gene regulation is critical for a myriad of developmental processes as well as stress responses. However, the role of miRNAs in maintaining rice grain quality/chalkiness during high daytime temperature (HDT) stress is relatively unknown. To uncover the role of miRNAs in this process, we used five contrasting rice genotypes (low chalky lines Cyp, Ben, and KB and high chalky lines LaGrue and NB) and compared the miRNA profiles in the R6 stage caryopsis samples from plants subjected to prolonged HDT (from the onset of fertilization through R6 stage of caryopsis development). Our small RNA analysis has identified approximately 744 miRNAs that can be grouped into 291 families. Of these, 186 miRNAs belonging to 103 families are differentially regulated under HDT. Only two miRNAs, Osa-miR444f and Osa-miR1866-5p, were upregulated in all genotypes, implying that the regulations greatly varied between the genotypes. Furthermore, not even a single miRNA was commonly up/down regulated specifically in the three tolerant genotypes. However, three miRNAs (Osa-miR1866-3p, Osa-miR5150-3p and canH-miR9774a,b-3p) were commonly upregulated and onemiRNA (Osa-miR393b-5p) was commonly downregulated specifically in the sensitive genotypes (LaGrue and NB). These observations suggest that few similarities exist within the low chalky or high chalky genotypes, possibly due to high genetic variation. Among the five genotypes used, Cypress and LaGrue are genetically closely related, but exhibit contrasting chalkiness under HDT, and thus, a comparison between them is most relevant. This comparison revealed a general tendency for Cypress to display miRNA regulations that could decrease chalkiness under HDT compared with LaGrue. This study suggests that miRNAs could play an important role in maintaining grain quality in HDT-stressed rice.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

The Integrated mRNA and miRNA Approach Reveals Potential Regulators of Flowering Time in Arundina graminifolia

Orchids are among the most precious flowers in the world. Regulation of flowering time is one of the most important targets to enhance their ornamental value. The beauty of Arundina graminifolia is its year-round flowering, although the molecular mechanism of this flowering ability remains masked. Therefore, we performed a comprehensive assessment to integrate transcriptome and miRNA sequencing to disentangle the genetic regulation of flowering in this valuable species. Clustering analyses provided a set of molecular regulators of floral transition and floral morphogenesis. We mined candidate floral homeotic genes, including FCA, FPA, GI, FT, FLC, AP2, SOC1, SVP, GI, TCP, and CO, which were targeted by a variety of miRNAs. MiR11091 targeted the highest number of genes, including candidate regulators of phase transition and hormonal control. The conserved miR156-miR172 pathway of floral time regulation was evident in our data, and we found important targets of these miRNAs in the transcriptome. Moreover, endogenous hormone levels were determined to decipher the hormonal control of floral buds in A. graminifolia. The qRT-PCR analysis of floral and hormonal integrators validated the transcriptome expression. Therefore, miRNA-mediated mining of candidate genes with hormonal regulation forms the basis for comprehending the complex regulatory network of perpetual flowering in precious orchids. The findings of this study can do a great deal to broaden the breeding programs for flowering time manipulation of orchids.

Identification of Conserved and Novel MicroRNAs with their Targets in Garden Pea (Pisum Sativum L.) Leaves by High-Throughput Sequencing

MicroRNAs (miRNAs) are single-stranded, endogenous, non-coding RNAs of 20–24 nucleotides that play a significant role in post-transcriptional gene regulation. Various conserved and novel miRNAs have been characterized, especially from the plant species whose genomes were well-characterized; however, information on miRNA in economically important plants such as pea ( Pisum sativum L.) is limited. In this study, I have identified conserved and novel miRNA in garden pea plant leaves samples along with their targets by analyzing the next generation sequencing (NGS) data. The raw data obtained from NGS were processed and 1.38 million high-quality non-redundant reads were retained for analysis, this tremendous quantity of reads indicates a large and diverse small RNA population in pea leaves. After analyzing the deep sequencing data, 255 conserved and 11 novel miRNAs were identified in the garden pea leaves sample. Utilizing psRNATarget tool, the miRNA targets of conserved and novel miRNA were predicted. Further, the functional annotation of the miRNA targets were performed using blast2Go software and the target gene products were predicted. The miRNA target gene products along with GO_ID (Gene Ontology Identifier) were categorized into biological processes, cellular components, and molecular functions. The information obtained from this study will provide genomic resources that will help in understanding miRNA-mediated post-transcriptional gene regulation in garden peas.

Identification of differentially expressed miRNAs and their target genes in response to brassinolide treatment on flowering of tree peony (Paeonia ostii)

Tree peony is a famous flower plant in China, but the short and concentrated flowering period limits its ornamental value and economic value. Brassinolide (BR) plays an important role in plant growth and development including flowering. There have been a large number of reports on the molecular aspects of the flowering process, but the genetic mechanism that was responsible for miRNA-guided regulation of tree peony is almost unclear. In this study, the leaves of tree peony cultivar, 'Feng Dan', were sprayed with different concentrations of BR, and the obvious bloom delay was found at the treatment with BR 50 μg/L. The small RNA sequencing and transcriptome sequencing were performed on the petals of tree peony under an untreated control (CK) and the treatment with BR 50 μg/L during four consecutive flowering development stages. A total of 22 known miRNAs belonging to 12 families were identified and 84 novel miRNAs were predicted. Combined with transcriptome data, a total of 376 target genes were predicted for the 18 differentially expressed known miRNAs and 177 target genes were predicted for the 23 differentially expressed novel miRNAs. Additionally, the potential miRNAs and their target genes were identified, including miR156b targeting SPL, miR172a_4 targeting AP2 and four novel miRNAs targeting SPA1, and revealed that they might affect the flowering time in tree peony. Collectively, these results would provide a theoretical basis for further analysis of miRNA-guided regulation on flowering period in tree peony.

High-Throughput Sequencing Reveals Novel microRNAs Involved in the Continuous Flowering Trait of Longan (Dimocarpus longan Lour.)

A major determinant of fruit production in longan (Dimocarpus longan Lour.) is the difficulty of blossoming. In this study, high-throughput microRNA sequencing (miRNA-Seq) was carried out to compare differentially expressed miRNAs (DEmiRNAs) and their target genes between a continuous flowering cultivar 'Sijimi' (SJ), and a unique cultivar 'Lidongben' (LD), which blossoms only once in the season. Over the course of our study, 1662 known miRNAs and 235 novel miRNAs were identified and 13,334 genes were predicted to be the target of 1868 miRNAs. One conserved miRNA and 29 new novel miRNAs were identified as differently expressed; among them, 16 were upregulated and 14 were downregulated. Through the KEGG pathway and cluster analysis of DEmiRNA target genes, three critical regulatory pathways, plant-pathogen interaction, plant hormone signal transduction, and photosynthesis-antenna protein, were discovered to be strongly associated with the continuous flowering trait of the SJ. The integrated correlation analysis of DEmiRNAs and their target mRNAs revealed fourteen important flowering-related genes, including COP1-like, Casein kinase II, and TCP20. These fourteen flowering-related genes were targeted by five miRNAs, which were novel-miR137, novel-miR76, novel-miR101, novel-miR37, and csi-miR3954, suggesting these miRNAs might play vital regulatory roles in flower regulation in longan. Furthermore, novel-miR137 was cloned based on small RNA sequencing data analysis. The pSAK277-miR137 transgenic Arabidopsis plants showed delayed flowering phenotypes. This study provides new insight into molecular regulation mechanisms of longan flowering.

Genome-Wide Identification and In Silico Analysis of ZF-HD Transcription Factor Genes in Zea mays L

Zinc finger-homeodomain proteins are amongst the most prominent transcription factors (TFs) involved in biological processes, such as growth, development, and morphogenesis, and assist plants in alleviating the adverse effects of abiotic and biotic stresses. In the present study, genome-wide identification and expression analyses of the maize ZHD gene family were conducted. A total of 21 ZHD genes with different physicochemical properties were found distributed on nine chromosomes in maize. Through sequence alignment and phylogenetic analysis, we divided ZHD proteins into eight groups that have variations in gene structure, motif distribution, and a conserved ZF domain. Synteny analysis indicated duplication in four pairs of genes and the presence of orthologues of maize in monocots. Ka/Ks ratios suggested that strong pure selection occurred during evolution. Expression profiling revealed that the genes are evenly expressed in different tissues. Most of the genes were found to make a contribution to abiotic stress response, plant growth, and development. Overall, the evolutionary research on exons and introns, motif distributions, and cis-acting regions suggests that these genes play distinct roles in biological processes which may provide a basis for further study of these genes' functions in other crops.

An overview of salinity stress, mechanism of salinity tolerance and strategies for its management in cotton

Salinity stress is one of the primary threats to agricultural crops resulting in impaired crop growth and development. Although cotton is considered as reasonably salt tolerant, it is sensitive to salt stress at some critical stages like germination, flowering, boll formation, resulting in reduced biomass and fiber production. The mechanism of partial ion exclusion (exclusion of Na+ and/or Cl-) in cotton appears to be responsible for the pattern of uptake and accumulation of harmful ions (Na+ and Cl) in tissues of plants exposed to saline conditions. Maintaining high tissue K+/Na+ and Ca2+/Na+ ratios has been proposed as a key selection factor for salt tolerance in cotton. The key adaptation mechanism in cotton under salt stress is excessive sodium exclusion or compartmentation. Among the cultivated species of cotton, Egyptian cotton (Gossypium barbadense L.) exhibit better salt tolerance with good fiber quality traits as compared to most cultivated cotton and it can be used to improve five quality traits and transfer salt tolerance into Upland or American cotton (Gossypium hirsutum L.) by interspecific introgression. Cotton genetic studies on salt tolerance revealed that the majority of growth, yield, and fiber traits are genetically determined, and controlled by quantitative trait loci (QTLs). Molecular markers linked to genes or QTLs affecting key traits have been identified, and they could be utilized as an indirect selection criterion to enhance breeding efficiency through marker-assisted selection (MAS). Transfer of genes for compatible solute, which are an important aspect of ion compartmentation, into salt-sensitive species is, theoretically, a simple strategy to improve tolerance. The expression of particular stress-related genes is involved in plant adaptation to environmental stressors. As a result, enhancing tolerance to salt stress can be achieved by marker assisted selection added with modern gene editing tools can boost the breeding strategies that defend and uphold the structure and function of cellular components. The intent of this review was to recapitulate the advancements in salt screening methods, tolerant germplasm sources and their inheritance, biochemical, morpho-physiological, and molecular characteristics, transgenic approaches, and QTLs for salt tolerance in cotton.

Versatile Roles of Microbes and Small RNAs in Rice and Planthopper Interactions

Planthopper infestation in rice causes direct and indirect damage through feeding and viral transmission. Host microbes and small RNAs (sRNAs) play essential roles in regulating biological processes, such as metabolism, development, immunity, and stress responses in eukaryotic organisms, including plants and insects. Recently, advanced metagenomic approaches have facilitated investigations on microbial diversity and its function in insects and plants, highlighting the significance of microbiota in sustaining host life and regulating their interactions with the environment. Recent research has also suggested significant roles for sRNA-regulated genes during rice-planthopper interactions. The response and behavior of the rice plant to planthopper feeding are determined by changes in the host transcriptome, which might be regulated by sRNAs. In addition, the roles of microbial symbionts and sRNAs in the host response to viral infection are complex and involve defense-related changes in the host transcriptomic profile. This review reviews the structure and potential functions of microbes and sRNAs in rice and the associated planthopper species. In addition, the involvement of the microbiota and sRNAs in the rice-planthopper-virus interactions during planthopper infestation and viral infection are discussed.

Excess nitrogen responsive HvMADS27 transcription factor controls barley root architecture by regulating abscisic acid level

Nitrogen (N) is an important element for plant growth and development. Although several studies have examined plants' response to N deficiency, studies on plants' response to excess N, which is common in fertilizer-based agrosystems, are limited. Therefore, the aim of this study was to examine the response of barley to excess N conditions, specifically the root response. Additionally, genomic mechanism of excess N response in barley was elucidated using transcriptomic technologies. The results of the study showed that barley MADS27 transcription factor was mainly expressed in the roots and its gene contained N-responsive cis-regulatory elements in the promoter region. Additionally, there was a significant decrease in HvMADS27 expression under excess N condition; however, its expression was not significantly affected under low N condition. Phenotypic analysis of the root system of HvMADS27 knockdown and overexpressing barley plants revealed that HvMADS27 regulates barley root architecture under excess N stress. Further analysis of wild-type (WT) and transgenic barley plants (hvmads27 kd and hvmads27 c-Myc OE) revealed that HvMADS27 regulates the expression of HvBG1 β-glucosidase, which in turn regulates abscisic acid (ABA) level in roots. Overall, the findings of this study showed that HvMADS27 expression is downregulated in barley roots under excess N stress, which induces HvBG1 expression, leading to the release of ABA from ABA-glucose conjugate, and consequent shortening of the roots.

Rapid sequence and functional diversification of a miRNA superfamily targeting calcium signaling components in seed plants

MicroRNA (miRNA)-directed posttranscriptional gene silencing (miR-PTGS) is an integral component of gene regulatory networks governing plant development and responses to the environment. The sequence homology between Sly-miR4376, a miRNA common to Solanaceae and reported to target autoinhibited Ca²⁺ -ATPase 10 (ACA10) messenger RNA (mRNA) in tomato, and Arabidopsis miR391 (Ath-miR391), previously annotated as a nonconserved member of the deeply conserved miR390 family, has prompted us to revisit the function of Ath-miR391, as well as its regulatory conservation. A combination of genetic, molecular, and bioinformatic analyses revealed a hidden conservation for miR-PTGS of ACA10 homologs in spermatophytes. We found that the Arabidopsis ACA10 mRNA undergoes miR391-directed cleavage in vivo. Furthermore, transgenic overexpression of miR391 recapitulated the compact inflorescence (cif) phenotypes characteristic of ACA10 loss-of-function mutants, due to miR391-directed PTGS of ACA10. Significantly, comprehensive data mining revealed robust evidence for widespread PTGS of ACA10 homologs directed by a superfamily of related miRNAs sharing a conserved sequence core. Intriguingly, the ACA-targeting miRNAs in Poaceae also direct PTGS for calmodulin-like proteins which are putative Ca²⁺ sensors. The PTGS of ACA10 homologs is therefore directed by a miRNA superfamily that is of ancient origin and has undergone rapid sequence diversification associated with functional innovation.

Pan-genomic, transcriptomic, and miRNA analyses to decipher genetic diversity and anthocyanin pathway genes among the traditional rice landraces

Black rice is famous for containing high anthocyanin while Joha rice is aromatic with low anthocyanin containing rice from the North Eastern Region (NER) of India. However, there are limited reports on the anthocyanin biosynthesis in Manipur Black rice. Therefore, the present study was aimed to understand the origin, domestication and anthocyanin biosynthesis pathways in Black rice using the next generation sequencing of approaches. With the sequencing data, various analyses were carried out for differential expression and construction of a pan-genome. Protein coding RNA and small RNA sequencing analysis aided in determining 7415 and 131 differentially expressed transcripts and miRNAs, respectively in NER rice. This is the first extensive study on identification and expression analysis of miRNAs and their target genes in regulating anthocyanin biosynthesis in NER rice. This study will aid in better understanding for decoding the theory of high or low anthocyanin content in different rice genotypes.

An Insight Into Pentatricopeptide-Mediated Chloroplast Necrosis via microRNA395a During Rhizoctonia solani Infection

Sheath blight (ShB) disease, caused by Rhizoctonia solani, is one of the major biotic stress-oriented diseases that adversely affect the rice productivity worldwide. However, the regulatory mechanisms are not understood yet comprehensively. In the current study, we had investigated the potential roles of miRNAs in economically important indica rice variety Pusa Basmati-1 upon R. solani infection by carrying out in-depth, high-throughput small RNA sequencing with a total data size of 435 million paired-end raw reads from rice leaf RNA samples collected at different time points. Detailed data analysis revealed a total of 468 known mature miRNAs and 747 putative novel miRNAs across all the libraries. Target prediction and Gene Ontology functional analysis of these miRNAs were found to be unraveling various cellular, molecular, and biological functions by targeting various plant defense-related genes. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed to validate the miRNAs and their putative target genes. Out of the selected miRNA-specific putative target genes, miR395a binding and its cleavage site on pentatricopeptide were determined by 5’ RACE-PCR. It might be possible that R. solani instigated chloroplast degradation by modulating the pentatricopeptide which led to increased susceptibility to fungal infection.

Integrated Multi-Omics Perspective to Strengthen the Understanding of Salt Tolerance in Rice

Salt stress is one of the major constraints to rice cultivation worldwide. Thus, the development of salt-tolerant rice cultivars becomes a hotspot of current rice breeding. Achieving this goal depends in part on understanding how rice responds to salt stress and uncovering the molecular mechanism underlying this trait. Over the past decade, great efforts have been made to understand the mechanism of salt tolerance in rice through genomics, transcriptomics, proteomics, metabolomics, and epigenetics. However, there are few reviews on this aspect. Therefore, we review the research progress of omics related to salt tolerance in rice and discuss how these advances will promote the innovations of salt-tolerant rice breeding. In the future, we expect that the integration of multi-omics salt tolerance data can accelerate the solution of the response mechanism of rice to salt stress, and lay a molecular foundation for precise breeding of salt tolerance.

An Integration of MicroRNA and Transcriptome Sequencing Analysis Reveal Regulatory Roles of miRNAs in Response to Chilling Stress in Wild Rice

A chromosome single segment substitution line (CSSL) DC90, which was generated by introgressing CTS-12, a locus derived from common wild rice (Oryza rufipogon Griff.), into the 9311 (Oryza sativa L. ssp. indica) background, exhibits a chilling tolerance phenotype under chilling stress. Here, an integration of microRNA (miRNA) deep sequencing and transcriptomic sequencing analysis was performed to explore the expression profiles of miRNAs and their target genes mediated by CTS-12 under chilling stress, and to reveal the possible regulatory mechanisms of miRNAs that are involved in chilling tolerance. Integration analysis revealed that a number of differentially expressed miRNAs (DEMs) and putative target genes with different expression patterns and levels were identified in 9311 and DC90 under chilling stress. KEGG enrichment analysis revealed that the target genes that are regulated by chilling-induced miRNAs are involved in the regulation of various biological processes/pathways, including protein biosynthesis, redox process, photosynthetic process, and chloroplast development in two genotypes. CRISPR/Cas9 editing of the target genes of the key DEMs in a chilling tolerant rice variety Zhonghua 11 (ZH11) found that LOC_Os11g48020 (OsGL1-11), one of the putative target genes of osa-miR1846a/b-5p and encoding a wax synthesis protein, is correlated with a chilling stress tolerance phenotype, implying osa-miR1846a/b-5p/OsGL1-11 plays an important role in CTS-12-mediated chilling stress tolerance regulatory pathway(s). Therefore, we speculate that the CTS-12 may regulate the key miRNA target genes in response to chilling stress by differential regulation of miRNAs in wild rice, thereby resulting in the variation of chilling tolerance phenotype between 9311 and DC90.

The Multiverse of Plant Small RNAs: How Can We Explore It?

Plant small RNAs (sRNAs) are a heterogeneous group of noncoding RNAs with a length of 20–24 nucleotides that are widely studied due to their importance as major regulators in various biological processes. sRNAs are divided into two main classes—microRNAs (miRNAs) and small interfering RNAs (siRNAs)—which differ in their biogenesis and functional pathways. Their identification and enrichment with new structural variants would not be possible without the use of various high-throughput sequencing (NGS) techniques, allowing for the detection of the total population of sRNAs in plants. Classifying sRNAs and predicting their functional role based on such high-performance datasets is a nontrivial bioinformatics task, as plants can generate millions of sRNAs from a variety of biosynthetic pathways. Over the years, many computing tools have been developed to meet this challenge. Here, we review more than 35 tools developed specifically for plant sRNAs over the past few years and explore some of their basic algorithms for performing tasks related to predicting, identifying, categorizing, and quantifying individual sRNAs in plant samples, as well as visualizing the results of these analyzes. We believe that this review will be practical for biologists who want to analyze their plant sRNA datasets but are overwhelmed by the number of tools available, thus answering the basic question of how to choose the right one for a particular study.

The Genetic and Hormonal Inducers of Continuous Flowering in Orchids: An Emerging View

Orchids are the flowers of magnetic beauty. Vivid and attractive flowers with magnificent shapes make them the king of the floriculture industry. However, the long-awaited flowering is a drawback to their market success, and therefore, flowering time regulation is the key to studies about orchid flower development. Although there are some rare orchids with a continuous flowering pattern, the molecular regulatory mechanisms are yet to be elucidated to find applicable solutions to other orchid species. Multiple regulatory pathways, such as photoperiod, vernalization, circadian clock, temperature and hormonal pathways are thought to signalize flower timing using a group of floral integrators. This mini review, thus, organizes the current knowledge of floral time regulators to suggest future perspectives on the continuous flowering mechanism that may help to plan functional studies to induce flowering revolution in precious orchid species.

The rice EP3 and OsFBK1 E3 ligases alter plant architecture and flower development, and affect transcript accumulation of microRNA pathway genes and their targets

ERECTA PANICLE 3 (EP3) and ORYZA SATIVA F-BOX KELCH 1 (OsFBK1) proteins share 57% and 54% sequence identity with the Arabidopsis F-box protein HAWAIIAN SKIRT (HWS). Previously we showed that EP3 is a functional orthologue of HWS. Here we demonstrate that OsFBK1 is another functional orthologue of HWS and show the complexity of interaction between EP3 and OsFBK1 genes at different developmental stages of the plant. qRT-PCR expression analyses and studies of EP3-GFP and OsFBK1-RFP promoter reporter lines demonstrate that although EP3 and OsFBK1 expression can be detected in the same tissues some cells exclusively express EP3 or OsFBK1 whilst others co-express both genes. Loss, reduction or gain-of-function lines for EP3 and OsFBK1, show that EP3 and OsFBK1 affect plant architecture, organ size, floral organ number and size, floral morphology, pollen viability, grain size and weight. We have identified the putative orthologue genes of the rice microRNA pathway for ORYZA SATIVA DAWDLE (OsDDL) and ORYZA SATIVA SERRATE (OsSE), and demonstrated that EP3 and OsFBK1 affect their transcript levels as well as those of CROWN ROOT DEFECT 1/ORYZA SATIVA Exportin-5 HASTY (CRD1/OsHST), ORYZA SATIVA DICER-LIKE 1 (OsDCL) and ORYZA SATIVA WEAVY LEAF1 (OsWAF1). We show that EP3 affects OsPri-MIR164, OsNAM1 and OsNAC1 transcript levels. OsNAC1 transcripts are modified by OsFBK1, suggesting two independent regulatory pathways, one via EP3 and OsMIR164 and the other via OsFBK1. Our data propose that EP3 and OsFBK1 conjointly play similar roles in rice to how HWS does in Arabidopsis.

MicroRNAs Are Involved in Regulating Plant Development and Stress Response through Fine-Tuning of TIR1/AFB-Dependent Auxin Signaling

Auxin, primarily indole-3-acetic acid (IAA), is a versatile signal molecule that regulates many aspects of plant growth, development, and stress response. Recently, microRNAs (miRNAs), a type of short non-coding RNA, have emerged as master regulators of the auxin response pathways by affecting auxin homeostasis and perception in plants. The combination of these miRNAs and the autoregulation of the auxin signaling pathways, as well as the interaction with other hormones, creates a regulatory network that controls the level of auxin perception and signal transduction to maintain signaling homeostasis. In this review, we will detail the miRNAs involved in auxin signaling to illustrate its in planta complex regulation.

Computational screening of miRNAs and their targets in saffron (Crocus sativus L.) by transcriptome mining

A robust workflow for the identification of miRNAs and their targets in saffron was developed. MicroRNA-mediated gene regulation in saffron is potentially involved in several biological processes, including the biosynthesis of highly valuable apocarotenoids. Saffron (Crocus sativus L.) is the most expensive spice in the world and a major source of apocarotenoids. Even though miRNAs (20-24 nt non-coding small RNAs) are important regulators of gene expression at transcriptional and post-transcriptional levels, their role in saffron has not been thoroughly investigated. As a result, a workflow for computational identification of miRNAs and their targets can be useful to uncover the regulatory networks underlying biological processes in this valuable plant. The efficiency of several assembly tools such as Trans-ABySS, Trinity, Bridger, rnaSPAdes, and EvidentialGene was evaluated based on both reference-based and reference-free metrics using transcriptome data. A reliable workflow for computational identification of miRNAs and their targets in saffron was described. The EvidentialGene was found to be the most efficient de novo transcriptome assembler for saffron as a complex triploid model, followed by the Trinity. In total, 66 miRNAs from 19 different families that target 2880 genes, including several transcription factors involved in the flowering transition, were identified. Three of the identified targets were involved in the terpenoids backbone biosynthesis. CsCCD and CsUGT genes involved in the apocarotenoids biosynthetic pathway were targeted by csa-miR156g and csa-miR156b-3p, revealing a unique post-transcriptional regulation dynamic in saffron. The identified miRNAs and their targets add to our understanding of the many biological roles of miRNAs in saffron and shed new light on the control of the apocarotenoid biosynthetic pathway in this valuable plant.

miRPreM and tiRPreM: Improved methodologies for the prediction of miRNAs and tRNA-induced small non-coding RNAs for model and non-model organisms

In recent years, microRNAs (miRNAs) and tRNA-derived RNA fragments (tRFs) have been reported extensively following different approaches of identification and analysis. Comprehensively analyzing the present approaches to overcome the existing variations, we developed a benchmarking methodology each for the identification of miRNAs and tRFs, termed as miRNA Prediction Methodology (miRPreM) and tRNA-induced small non-coding RNA Prediction Methodology (tiRPreM), respectively. We emphasized the use of respective genome of organism under study for mapping reads, sample data with at least two biological replicates, normalized read count support and novel miRNA prediction by two standard tools with multiple runs. The performance of these methodologies was evaluated by using Oryza coarctata, a wild rice species as a case study for model and non-model organisms. With organism-specific reference genome approach, 98 miRNAs and 60 tRFs were exclusively found. We observed high accuracy (13 out of 15) when tested these genome-specific miRNAs in support of analyzing the data with respective organism. Such a strong impact of miRPreM, we have predicted more than double number of miRNAs (186) as compared with the traditional approaches (79) and with tiRPreM, we have predicted all known classes of tRFs within the same small RNA data. Moreover, the methodologies presented here are in standard form in order to extend its applicability to different organisms rather than restricting to plants. Hence, miRPreM and tiRPreM can fulfill the need of a comprehensive methodology for miRNA prediction and tRF identification, respectively, for model and non-model organisms.

Investigations on Regulation of MicroRNAs in Rice Reveal [Ca2+]cyt Signal Transduction Regulated MicroRNAs

MicroRNAs (miRNAs) are critical components of the multidimensional regulatory networks in eukaryotic systems. Given their diverse spectrum of function, it is apparent that the transcription, processing, and activity of the miRNAs themselves, is very dynamically regulated. One of the most important and universally implicated signaling molecules is [Ca2+]cyt. It is known to regulate a plethora of developmental and metabolic processes in both plants and animals; however, its impact on the regulation of miRNA expression is relatively less explored. The current study employed a combination of internal and external calcium channel inhibitors to establishing that [Ca2+]cyt signatures actively regulate miRNA expression in rice. Involvement of [Ca2+]cyt in the regulation of miRNA expression was further confirmed by treatment with calcimycin, the calcium ionophore. Modulation of the cytosolic calcium levels was also found to regulate the drought-responsive expression as well as ABA-mediated response of miRNA genes in rice seedlings. The study further establishes the role of calmodulins and Calmodulin-binding Transcription Activators (CAMTAs) as important components of the signal transduction schema that regulates miRNA expression. Yeast one-hybrid assay established that OsCAMTA4 & 6 are involved in the transcriptional regulation of miR156a and miR167h. Thus, the study was able to establish that [Ca2+]cyt is actively involved in regulating the expression of miRNA genes both under control and stress conditions.

Genome and transcriptome analysis to understand the role diversification of cytochrome P450 gene under excess nitrogen treatment

BACKGROUND

Panax notoginseng (Burk.) F. H. Chen (P. notoginseng) is a medicinal plant. Cytochrome P450 (CYP450) monooxygenase superfamily is involved in the synthesis of a variety of plant hormones. Studies have shown that CYP450 is involved in the synthesis of saponins, which are the main medicinal component of P. notoginseng. To date, the P. notoginseng CYP450 family has not been systematically studied, and its gene functions remain unclear.

RESULTS

In this study, a total of 188 PnCYP genes were identified, these genes were divided into 41 subfamilies and clustered into 9 clans. Moreover, we identified 40 paralogous pairs, of which only two had Ka/Ks ratio greater than 1, demonstrating that most PnCYPs underwent purification selection during evolution. In chromosome mapping and gene replication analysis, 8 tandem duplication and 11 segmental duplication events demonstrated that PnCYP genes were continuously replicating during their evolution. Gene ontology (GO) analysis annotated the functions of 188 PnCYPs into 21 functional subclasses, suggesting the functional diversity of these gene families. Functional divergence analyzed the members of the three primitive branches of CYP51, CYP74 and CYP97 at the amino acid level, and found some critical amino acid sites. The expression pattern of PnCYP450 related to nitrogen treatment was studied using transcriptome sequencing data, 10 genes were significantly up-regulated and 37 genes were significantly down-regulated. Combined with transcriptome sequencing analysis, five potential functional genes were screened. Quantitative real-time PCR (qRT-PCR) indicated that these five genes were responded to methyl jasmonate (MEJA) and abscisic acid (ABA) treatment.

CONCLUSIONS

These results provide a valuable basis for comprehending the classification and biological functions of PnCYPs, and offer clues to study their biological functions in response to nitrogen treatment.

Transcriptome and Small RNA Combined Sequencing Analysis of Cold Tolerance in Non-heading Chinese Cabbage

Background

Non-heading Chinese cabbage (Brassica rapa ssp. chinensis) is an important leaf vegetable grown worldwide. However, there has currently been not enough transcriptome and small RNA combined sequencing analysis of cold tolerance, which hinders further functional genomics research.

Results

In this study, 63.43 Gb of clean data was obtained from the transcriptome analysis. The clean data of each sample reached 6.99 Gb, and the basic percentage of Q30 was 93.68% and above. The clean reads of each sample were sequence aligned with the designated reference genome (Brassica rapa, IVFCAASv1), and the efficiency of the alignment varied from 81.54 to 87.24%. According to the comparison results, 1,860 new genes were discovered in Pak-choi, of which 1,613 were functionally annotated. Among them, 13 common differentially expressed genes were detected in all materials, including seven upregulated and six downregulated. At the same time, we used quantitative real-time PCR to confirm the changes of these gene expression levels. In addition, we sequenced miRNA of the same material. Our findings revealed a total of 34,182,333 small RNA reads, 88,604,604 kinds of small RNAs, among which the most common size was 24 nt. In all materials, the number of common differential miRNAs is eight. According to the corresponding relationship between miRNA and its target genes, we carried out Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis on the set of target genes on each group of differentially expressed miRNAs. Through the analysis, it is found that the distributions of candidate target genes in different materials are different. We not only used transcriptome sequencing and small RNA sequencing but also used experiments to prove the expression levels of differentially expressed genes that were obtained by sequencing. Sequencing combined with experiments proved the mechanism of some differential gene expression levels after low-temperature treatment.

Conclusion

In all, this study provides a resource for genetic and genomic research under abiotic stress in Pak-choi.

Genome-Wide Identification and Evolutionary Analysis of Argonaute Genes in Hexaploid Bread Wheat

Argonaute (AGO) proteins play a pivotal role in plant growth and development as the core components of RNA-induced silencing complex (RISC). However, no systematic characterization of AGO genes in wheat has been reported to date. In this study, a total number of 69 TaAGO genes in the hexaploid bread wheat (Triticum aestivum cv. Chinese Spring) genome, divided into 10 subfamilies, were identified. Compared to all wheat genes, TaAGOs showed a significantly lower evolutionary rate, which is consistent with their high conservation in eukaryotes. However, the homoeolog retention was remarkably higher than the average, implying the nonredundant biological importance of TaAGO genes in bread wheat. Further homoeologous gene expression bias analyses revealed that TaAGOs may have undergone neofunctionalization after polyploidization and duplication through the divergent expression of homoeologous gene copies, to provide new opportunities for the generation of adaptive traits. Moreover, quantitative real-time polymerase chain reaction (qRT-PCR) analyses indicated that TaAGO gene expression was involved in response to heat, drought, and salt stresses. Our results would provide a theoretical basis for future studies on the biological functions of TaAGO genes in wheat and other gramineous species.

Genome-wide analysis of microRNA156 and its targets, the genes encoding SQUAMOSA promoter-binding protein-like (SPL) transcription factors, in the grass family Poaceae

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play an important role in post-transcriptional gene regulation in plants and animals by targeting messenger RNAs (mRNAs) for cleavage or repressing translation of specific mRNAs. The first miRNA identified in plants, miRNA156 (miR156), targets the SQUAMOSA promoter-binding protein-like (SPL) transcription factors, which play critical roles in plant phase transition, flower and plant architecture, and fruit development. We identified multiple copies of MIR156 and SPL in the rice, Brachypodium, sorghum, maize, and foxtail millet genomes. Sequence and chromosomal synteny analysis showed that both MIR156s and SPLs are conserved across species in the grass family. Analysis of expression data of the SPLs in eleven juvenile and adult rice tissues revealed that four non-miR156-targeted genes were highly expressed and three miR156-targeted genes were only slightly expressed in all tissues/developmental stages. The remaining SPLs were highly expressed in the juvenile stage, but their expression was lower in the adult stage. It has been proposed that under strong selective pressure, non-miR156-targeted mRNA may be able to re-structure to form a miRNA-responsive element. In our analysis, some non-miR156-targeted SPLs (SPL5/8/10) had gene structure and gene expression patterns similar to those of miR156-targeted genes, suggesting that they could diversify into miR156-targeted genes. DNA methylation profiles of SPLs and MIR156s in different rice tissues showed diverse methylation patterns, and hypomethylation of non-CG sites was observed in rice endosperm. Our findings suggested that MIR156s and SPLs had different origination and evolutionary mechanisms: the SPLs appear to have resulted from vertical evolution, whereas MIR156s appear to have resulted from strong evolutionary selection on mature sequences.

Priming by High Temperature Stress Induces MicroRNA Regulated Heat Shock Modules Indicating Their Involvement in Thermopriming Response in Rice

Rice plants often encounter high temperature stress, but the associated coping strategies are poorly understood. It is known that a prior shorter exposure to high temperature, called thermo-priming, generally results in better adaptation of the plants to subsequent exposure to high temperature stress. High throughput sequencing of transcript and small RNA libraries of rice seedlings primed with short exposure to high temperature followed by high temperature stress and from plants exposed to high temperature without priming was performed. This identified a number of transcripts and microRNAs (miRs) that are induced or down regulated. Among them osa-miR531b, osa-miR5149, osa-miR168a-5p, osa-miR1846d-5p, osa-miR5077, osa-miR156b-3p, osa-miR167e-3p and their respective targets, coding for heat shock activators and repressors, showed differential expression between primed and non-primed plants. These findings were further validated by qRT-PCR. The results indicate that the miR-regulated heat shock proteins (HSPs)/heat shock transcription factors (HSFs) may serve as important regulatory nodes which are induced during thermo-priming for plant survival and development under high temperatures.

Comprehensive analysis of the isomiRome in the vegetative organs of the conifer Pinus pinaster under contrasting water availability

An increasing number of microRNAs (miRNAs) and miRNA-related sequences produced during miRNA biogenesis, comprising the isomiRome, have been recently highlighted in different species as critical mediators of environmental stress responses. Conifers have some of the largest known genomes but an extensive characterization of the isomiRome from any conifer species has been lacking. We provide here a comprehensive overview of the Pinus pinaster isomiRome expressed in roots, stem and needles under well-watered and drought conditions. From the 13,441 unique small RNA sequences identified, 2,980 were annotated as canonical miRNAs or miRNA* and the remaining were classified as isomiRNA or miRNA-like sequences. A survey of their expression patterns highlighted roots as the most responsive organ under drought, where specific sequences of which a 24-nt novel miRNA stood out, were strongly down-regulated. Given the putative roles of the miRNA-targeted transcripts validated specifically in root tissues, some of the miRNAs, conserved and novel, are shortlisted as potential regulators of drought response. These results provide a valuable resource for comparative studies between gymnosperms and angiosperms. Furthermore, it evidences high transferability of the isomiRome between pine species being a useful basis for further molecular regulation and physiological studies, and especially those focused on adaptation to drought conditions.

Deciphering the role of microRNAs during Pi54 gene mediated Magnaporthe oryzae resistance response in rice

The broad-spectrum resistance gene Pi54 confers resistance to multiple isolates of Magnaporthe oryzae in rice. In order to decipher the molecular mechanism underlying the Pi54 mediated resistance in rice line Taipei309 ^Pi54 (carrying Pi54), miRNAome study was performed at 24 h post-inoculation (hpi) with M. oryzae. A total of 222 known miRNAs representing 101 miRNA families were found in this study. Of these, 29 and 24 miRNAs were respectively up- and down-regulated in the resistant Taipei309 ^Pi54 . Defence response (DR) genes, like, NBSGO35, and OsWAK129b, and genes related to transcription factors were up-regulated in Taipei309 ^Pi54 line. The vast array of miRNA candidates identified here are miR159c, miR167c, miR2100, miR2118o, miR2118l, miR319a, miR393, miR395l, miR397a, miR397b, miR398, miR439g, miR531b, miR812f, and miR815c, and they manifest their role in balancing the interplay between various DR genes during Pi54 mediated resistance. We also validated miRNA/target gene pairs involved in hormone signalling, and cross-talk among hormone pathways regulating the rice immunity. This study suggests that the Pi54 gene mediated blast resistance is influenced by several microRNAs through PTI and ETI components in the rice line Taipei309 ^Pi54 , leading to incompatible host-pathogen interaction.

Genome-wide miRNA analysis and integrated network for flavonoid biosynthesis in Osmanthus fragrans

Background

Osmanthus fragrans is an important economical plant containing multiple secondary metabolites including flavonoids and anthocyanins. During the past years, the roles of miRNAs in regulating the biosynthesis of secondary metabolites in plants have been widely investigated. However, few studies on miRNA expression profiles and the potential roles in regulating flavonoid biosynthesis have been reported in O. fragrans.

Results

In this study, we used high-throughput sequencing technology to analyze the expression profiles of miRNAs in leaf and flower tissues of O. fragrans. As a result, 106 conserved miRNAs distributed in 47 families and 88 novel miRNAs were identified. Further analysis showed there were 133 miRNAs differentially expressed in leaves and flowers. Additionally, the potential target genes of miRNAs as well as the related metabolic pathways were predicted. In the end, flavonoid content was measured in flower and leaf tissues and potential role of miR858 in regulating flavonoid synthesis was illustrated in O. fragrans.

Conclusions

This study not only provided the genome-wide miRNA profiles in the flower and leaf tissue of O. fragrans, but also investigated the potential regulatory role of miR858a in flavonoid synthesis in O. fragrans. The results specifically indicated the connection of miRNAs to the regulation of secondary metabolite biosynthesis in non-model economical plant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07439-y.

Physiological, Transcriptomic, and Metabolic Responses of Ginkgo biloba L. to Drought, Salt, and Heat Stresses

Ginkgo biloba L. is highly adaptable and resistant to a range of abiotic stressors, allowing its growth in various environments. However, it is unclear how G. biloba responds to common environmental stresses. We explored the physiological, transcriptomic, and metabolic responses of G. biloba to short-term drought, salt, and heat stresses. Proline, H2O2, and ABA contents, along with CAT activity, increased under all three types of stress. SOD activity increased under salt and heat stresses, while soluble protein and IAA contents decreased under drought and salt stresses. With respect to metabolites, D-glyceric acid increased in response to drought and salt stresses, whereas isomaltose 1, oxalamide, and threonine 2 increased under drought. Piceatannol 2,4-hydroxybutyrate and 1,3-diaminopropane increased under salt stress, whereas 4-aminobutyric acid 1 and galactonic acid increased in response to heat stress. Genes regulating nitrogen assimilation were upregulated only under drought, while the GRAS gene was upregulated under all three types of stressors. ARF genes were downregulated under heat stress, whereas genes encoding HSF and SPL were upregulated. Additionally, we predicted that miR156, miR160, miR172, and their target genes participate in stress responses. Our study provides valuable data for studying the multilevel response to drought, salinity, and heat in G. biloba.

Identification of genuine and novel miRNAs in Amaranthus hypochondriacus from high-throughput sequencing data

Amaranth has been proposed as an exceptional alternative for food security and climate change mitigation. Information about the distribution, abundance, or specificity of miRNAs in amaranth species is scare. Here, small RNAs from seedlings under control, drought, heat, and cold stress conditions of the Amaranthus hypocondriacus variety "Gabriela" were sequenced and miRNA loci identified in the amaranth genome using the ShortStack software. Fifty-three genuine miRNA clustersthirty-nine belonging to conserved families, and fourteen novel, were identified. Identification of their target genes suggests that conserved amaranth miRNAs are involved in growth and developmental processes, as well as stress responses. MiR0005, an amaranth-specific miRNA, exhibited an unusual high level of expression, akin to that of conserved miRNAs. Overall, our results broaden our knowledge regarding the distribution, abundance and expression of miRNAs in amaranth, providing the basis for future research on miRNAs and their functions in this important species.

Identification and expression analysis of miRNAs and elucidation of their role in salt tolerance in rice varieties susceptible and tolerant to salinity

Soil salinization is a serious problem for cultivation of rice, as among cereals rice is the most salt sensitive crop, and more than 40% of the total agricultural land amounting to approximately 80 million ha the world over is salt affected. Salinity affects a plant in a varieties of ways, including ion toxicity, osmotic stress and oxidative damage. Since miRNAs occupy the top place in biochemical events determining a trait, understanding their role in salt tolerance is highly desirable, which may allow introduction of the trait in the rice cultivars of choice through biotechnological interventions. High throughput sequencing of sRNAs in the root and shoot tissues of the seedlings of the control and NaCl treated Pokkali, a salt-tolerant rice variety, identified 75 conserved miRNAs and mapped 200 sRNAs to the rice genome as novel miRNAs. Expression of nine novel miRNAs and two conserved miRNAs were confirmed by Northern blotting. Several of both conserved and novel miRNAs that expressed differentially in root and/or shoot tissues targeted transcription factors like AP2/EREBP domain protein, ARF, NAC, MYB, NF-YA, HD-Zip III, TCP and SBP reported to be involved in salt tolerance or in abiotic stress tolerance in general. Most of the novel miRNAs expressed in the salt tolerant wild rice Oryza coarctata, suggesting conservation of miRNAs in taxonomically related species. One of the novel miRNAs, osa-miR12477, also targeted L-ascorbate oxidase (LAO), indicating build-up of oxidative stress in the plant upon salt treatment, which was confirmed by DAB staining. Thus, salt tolerance might involve miRNA-mediated regulation of 1) cellular abundance of the hormone signaling components like EREBP and ARF, 2) synthesis of abiotic stress related transcription factors, and 3) antioxidative component like LAO for mitigation of oxidative damage. The study clearly indicated importance of osa-miR12477 regulated expression of LAO in salt tolerance in the plant.

Identification of drought-responsive miRNAs in Hippophae tibetana using high-throughput sequencing

MicroRNAs (miRNAs) play an important role in abiotic stress response in plants. However, the total miRNA profiles (miRNome) and drought-responsive miRNAs in H. tibetana have not been identified. In this study, we present the first report on the miRNome profiles of H. tibetana by high-throughput sequencing technology. 116 known and 4 predicted novel miRNAs were all identified in six H. tibetana samples. Moreover, to reveal the drought-responsive miRNAs in H. tibetana, we compared the miRNA profiles of H. tibetana grown under water sufficiency and drought stress. The results showed that 39 known miRNAs were up-regulated, while 34 miRNAs were downregulated under drought stress. Moreover, the expression of two novel miRNAs (novel_mir_24 and novel_mir_87) showed notable changes in response to drought stress. The target genes of these differentially expressed miRNAs were mainly enriched in cellular process, metabolic process, cell part, and response to stimulus. The identified drought-responsive miRNAs might be used for improving drought tolerance in H. tibetana and other plateau plants.

Independent origin of MIRNA genes controlling homologous target genes by partial inverted duplication of antisense-transcribed sequences

Some microRNAs (miRNAs) are key regulators of developmental processes, mainly by controlling the accumulation of transcripts encoding transcription factors that are important for morphogenesis. MADS-box genes encode a family of transcription factors which control diverse developmental processes in flowering plants. Here we study the convergent evolution of two MIRNA (MIR) gene families, named MIR444 and MIR824, targeting members of the same clade of MIKC^C -group MADS-box genes. We show that these two MIR genes most likely originated independently in monocots (MIR444) and in Brassicales (eudicots, MIR824). We provide evidence that, in both cases, the future target gene was transcribed in antisense prior to the evolution of the MIR genes. Both MIR genes then likely originated by a partial inverted duplication of their target genes, resulting in natural antisense organization of the newly evolved MIR gene and its target gene at birth. We thus propose a model for the origin of MIR genes, MEPIDAS (MicroRNA Evolution by Partial Inverted Duplication of Antisense-transcribed Sequences). MEPIDAS is a refinement of the inverted duplication hypothesis. According to MEPIDAS, a MIR gene evolves at a genomic locus at which the future target gene is also transcribed in the antisense direction. A partial inverted duplication at this locus causes the antisense transcript to fold into a stem-loop structure that is recognized by the miRNA biogenesis machinery to produce a miRNA that regulates the gene at this locus. Our analyses exemplify how to elucidate the origin of conserved miRNAs by comparative genomics and will guide future studies. OPEN RESEARCH BADGE: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/genbank/.

Deep-sequencing of Solanum commersonii small RNA libraries reveals riboregulators involved in cold stress response

Among wild species used in potato breeding, Solanum commersonii displays the highest tolerance to low temperatures under both acclimated (ACC) and non-acclimated (NACC) conditions. It is also the first wild potato relative with a known whole genome sequence. Recent studies have shown that abiotic stresses induce changes in the expression of many small non-coding RNA (sncRNA). We determined the small non-coding RNA (sncRNAome) of two clones of S. commersonii contrasting in their cold response phenotypes via smRNAseq. Differential analysis provided evidence that expression of several miRNAs changed in response to cold stress conditions. Conserved miR408a and miR408b changed their expression under NACC conditions, whereas miR156 and miR169 were differentially expressed only under ACC conditions. We also report changes in tasiRNA and secondary siRNA expression under both stress conditions. Our results reveal possible roles of sncRNA in the regulatory networks associated with tolerance to low temperatures and provide useful information for a more strategic use of genomic resources in potato breeding.

Applications and Trends of Machine Learning in Genomics and Phenomics for Next-Generation Breeding

Crops are the major source of food supply and raw materials for the processing industry. A balance between crop production and food consumption is continually threatened by plant diseases and adverse environmental conditions. This leads to serious losses every year and results in food shortages, particularly in developing countries. Presently, cutting-edge technologies for genome sequencing and phenotyping of crops combined with progress in computational sciences are leading a revolution in plant breeding, boosting the identification of the genetic basis of traits at a precision never reached before. In this frame, machine learning (ML) plays a pivotal role in data-mining and analysis, providing relevant information for decision-making towards achieving breeding targets. To this end, we summarize the recent progress in next-generation sequencing and the role of phenotyping technologies in genomics-assisted breeding toward the exploitation of the natural variation and the identification of target genes. We also explore the application of ML in managing big data and predictive models, reporting a case study using microRNAs (miRNAs) to identify genes related to stress conditions.

Expressivity of the key genes associated with seed and pod development is highly regulated via lncRNAs and miRNAs in Pigeonpea

Non-coding RNA’s like miRNA, lncRNA, have gained immense importance as a significant regulatory factor in different physiological and developmental processes in plants. In an effort to understand the molecular role of these regulatory agents, in the present study, 3019 lncRNAs and 227 miRNAs were identified from different seed and pod developmental stages in Pigeonpea, a major grain legume of Southeast Asia and Africa. Target analysis revealed that 3768 mRNAs, including 83 TFs were targeted by lncRNAs; whereas 3060 mRNA, including 154 TFs, were targeted by miRNAs. The targeted transcription factors majorly belong to WRKY, MYB, bHLH, etc. families; whereas the targeted genes were associated with the embryo, seed, and flower development. Total 302 lncRNAs interact with miRNAs and formed endogenous target mimics (eTMs) which leads to sequestering of the miRNAs present in the cell. Expression analysis showed that notably, Cc_lncRNA-2830 expression is up-regulated and sequestrates miR160h in pod leading to higher expression of the miR160h target gene, Auxin responsive factor-18. A similar pattern was observed for SPIKE, Auxin signaling F-box-2, Bidirectional sugar transporter, and Starch synthetase-2 eTMs. All the identified target mRNAs code for transcription factor and genes are involved in the processes like cell division, plant growth and development, starch synthesis, sugar transportation and accumulation of storage proteins which are essential for seed and pod development. On a combinatorial basis, our study provides a lncRNA and miRNA based regulatory insight into the genes governing seed and pod development in Pigeonpea.

Computational methods for the ab initio identification of novel microRNA in plants: a systematic review

BACKGROUND

MicroRNAs (miRNAs) play a vital role as post-transcriptional regulators in gene expression. Experimental determination of miRNA sequence and structure is both expensive and time consuming. The next-generation sequencing revolution, which facilitated the rapid accumulation of biological data has brought biology into the "big data" domain. As such, developing computational methods to predict miRNAs has become an active area of inter-disciplinary research.

OBJECTIVE

The objective of this systematic review is to focus on the developments of ab initio plant miRNA identification methods over the last decade.

DATA SOURCES

Five databases were searched for relevant articles, according to a well-defined review protocol.

STUDY SELECTION

The search results were further filtered using the selection criteria that only included studies on novel plant miRNA identification using machine learning.

DATA EXTRACTION

Relevant data from each study were extracted in order to carry out an analysis on their methodologies and findings.

RESULTS

Results depict that in the last decade, there were 20 articles published on novel miRNA identification methods in plants of which only 11 of them were primarily focused on plant microRNA identification. Our findings suggest a need for more stringent plant-focused miRNA identification studies.

CONCLUSION

Overall, the study accuracies are of a satisfactory level, although they may generate a considerable number of false negatives. In future, attention must be paid to the biological plausibility of computationally identified miRNAs to prevent further propagation of biologically questionable miRNA sequences.

Identification of microRNAs regulated by tobacco curly shoot virus co-infection with its betasatellite in Nicotiana benthamiana

BACKGROUND

MicroRNAs (miRNAs) are a class of 21-24 nucleotide endogenous non-coding small RNAs that play important roles in plant development and defense responses to biotic and abiotic stresses. Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus, cause leaf curling and plant stunting symptoms in many Solanaceae plants. The betasatellite of TbCSV (TbCSB) induces more severe symptoms and enhances virus accumulation when co-infect the plants with TbCSV.

METHODS

In this study, miRNAs regulated by TbCSV and TbCSB co-infection in Nicotiana benthamiana were characterized using high-throughput sequencing technology.

RESULTS

Small RNA sequencing analysis revealed that a total of 13 known miRNAs and 42 novel miRNAs were differentially expressed in TbCSV and TbCSB co-infected N. benthamiana plants. Several potential miRNA-targeted genes were identified through data mining and were involved in both catalytic and metabolic processes, in addition to plant defense mechanisms against virus infections according to Gene Ontology (GO) analyses. In addition, the expressions of several differentially expressed miRNAs and their miRNA-targeted gene were validated through quantitative real time polymerase chain reaction (qRT-PCR) approach.

CONCLUSIONS

A large number of miRNAs are identified, and their target genes, functional annotations also have been explored. Our results provide the information on N. benthamiana miRNAs and would be useful to further understand miRNA regulatory mechanisms after TbCSV and TbCSB co-infection.

Identification and Temporal Expression Analysis of Conserved and Novel MicroRNAs in the Leaves of Winter Wheat Grown in the Field

Cold acclimation and vegetative/reproductive transition are two important evolutionary adaptive mechanisms for winter wheat surviving the freezing temperature in winter and successful seeds setting in the next year. MicroRNA (miRNA) is a class of regulatory small RNAs (sRNAs), which plays critical roles in the growth and development of plants. However, the regulation mechanism of miRNAs during cold acclimation and vegetative/reproductive transition of winter wheat is not much understood. In this study, four sRNA libraries from leaves of winter wheat grown in the field at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage were analyzed to identify known and novel miRNAs and to understand their potential roles in the growth and development of winter wheat. We examined miRNA expression using a high-throughput sequencing technique. A total of 373 known, 55 novel, and 27 putative novel miRNAs were identified. Ninety-one miRNAs were found to be differentially expressed at the four stages. Among them, the expression of six known and eight novel miRNAs was significantly suppressed at the winter dormancy stage, whereas the expression levels of seven known and eight novel miRNAs were induced at this stage; three known miRNAs and three novel miRNAs were significantly induced at the spring green-up stage; six known miRNAs were induced at the spring green-up stage and reached the highest expression level at the jointing stage; and 20 known miRNAs and 10 novel miRNAs were significantly induced at the jointing stage. Expression of a number of representative differentially expressed miRNAs was verified using quantitative real-time polymerase chain reaction (qRT-PCR). Potential target genes for known and novel miRNAs were predicted. Moreover, six novel target genes for four Pooideae species-specific miRNAs and two novel miRNAs were verified using the RNA ligase-mediated 5'-rapid amplification of cDNA ends (RLM-5'RACE) technique. These results indicate that miRNAs are key non-coding regulatory factors modulating the growth and development of wheat. Our study provides valuable information for in-depth understanding of the regulatory mechanism of miRNAs in cold acclimation and vegetative/reproductive transition of winter wheat grown in the field.

Comparative profiling of roots small RNA expression and corresponding gene ontology and pathway analyses for low- and high-cadmium-accumulating genotypes of wheat in response to cadmium stress

MicroRNAs (miRNAs) participate in multiple biological processes in plant. Cd accumulation ability differs among varieties in wheat, but little is known about miRNAs and their function in Cd accumulation of wheat under Cd stress. Therefore, the present study detected small RNAs responsible for differential Cd accumulation between two contrasting wheat genotypes (low-Cd accumulation one L17 and high-Cd accumulation one H17) to identify novel targets to further study Cd stress in wheat. Under normal conditions, 139 miRNAs were differentially expressed between L17 and H17, while this value reached 142 after Cd exposure. For Cd-induced DEMs, total 25 miRNAs were differentially expressed in L17 after Cd treatment, while, 70 Cd-induced DEMs were found in H17. Moreover, GO analysis revealed that target genes of DEMs related to lipid biosynthetic process and chlorophyll binding are uniquely enriched in L17, target genes of DEMs related to ribosome biogenesis and sucrose alpha-glucosidase activity are uniquely enriched in H17. By pathway analysis, target genes of DEMs related to PI3K-Akt signaling pathway was specifically enriched in L17, target genes of DEMs related to carbohydrate digestion and absorption pathway was uniquely enriched in H17. In addition, miRNA-gene co-expression showed that tae-miR9774 was uniquely expressed between L17Cd and L17CK, while tae-miR398 was specially expressed between H17Cd and H17CK. Our results suggested that Cd-accumulating ability of L17 and H17 varied from the expression of induced unique miRNA, such as expression of tae-miR-9774 and tae-miR-398. Our study not only provide the foundation for further exploring the miRNAs-induced molecular mechanisms of Cd accumulation in wheat but also supply novel strategies for improving phytoremediation ability of food plants through genetic engineering.

Post-transcriptional adaptation of the aquatic plant Spirodela polyrhiza under stress and hormonal stimuli

The Lemnaceae family comprises aquatic plants of angiosperms gaining attention due to their utility in wastewater treatment, and rapid production of biomass that can be used as feed, fuel, or food. Moreover, it can serve as a model species for neotenous growth and environmental adaptation. The latter properties are subject to post-transcriptional regulation of gene expression, meriting investigation of how miRNAs in Spirodela polyrhiza, the most basal and most thoroughly sequenced member of the family, are expressed under different growth conditions. To further scientific understanding of its capacity to adapt to environmental cues, we measured miRNA expression and processing of their target sequences under different temperatures, and in the presence of abscisic acid, copper, kinetin, nitrate, and sucrose. Using two small RNA sequencing experiments and one degradome sequencing experiment, we provide evidence for 108 miRNAs. Sequencing cleaved mRNAs validated 42 conserved miRNAs with 83 targets and 24 novel miRNAs regulating 66 targets and created a list of 575 predicted and verified targets. These analyses revealed condition-induced changes in miRNA expression and cleavage activity, and resulted in the addition of stringently reviewed miRNAs to miRBase. This combination of small RNA and degradome sequencing provided not only high confidence predictions of conserved and novel miRNAs and targets, but also a view of the post-transcriptional regulation of adaptations. A unique aspect is the role of miR156 and miR172 expression and activity in its clonal propagation and neoteny. Additionally, low levels of 24 nt sRNAs were observed, despite the lack of recent retrotransposition.