The Function of miRNAs in Plants

Physiological and molecular response to drought in somatic plants from Pinus radiata embryonal masses induced at high temperatures

Drought and heat are among the major abiotic stresses in forest trees and are directly related with the consequences of climatic change. Many responses to abiotic stresses in plants have been associated with plant memory but mechanisms underlying this phenomenon remain unclear. Somatic embryogenesis, which is considered one of the most important methods for large-scale vegetative propagation of plants, is also used for stress induction and study the mechanisms involved in adaptation to abiotic stress. Specifically, heat stress during initiation stage of somatic embryogenesis has shown to have an impact in differential expression of stress related genes in pines. Modifications caused by a previous stress could eventually influence the stress tolerance of somatic plants years later. In this study we analysed the response to drought in 2-year-old radiata pine somatic plants, derived from embryonal masses initiated at 60 °C, at physiological, transcriptomic and amino acid accumulation level. Our results showed a more pronounce response to drought in plants coming from 60 °C treatment, which presented lower values in several physiological parameters as well as higher proline and tyrosine levels. Additionally, the transcriptomic response to drought was stronger in heat primed plants compared to control plants, suggesting a memory acquired two years before.

CsMPDB 1.0: An interactive web application for visualizing and exploring the microRNAs and phasiRNAs of tea plant (Camellia sinensis var. sinensis 'Shuchazao')

The microRNAs and phasiRNAs of plant are small non-coding RNAs with important functions through regulating gene expression at the post-transcriptional level. However, identifying miRNAs, phasiRNAs and their target genes from numerous sequencing raw data requires multiple software and command-line operations, which are time-consuming and labor-intensive for non-model plants. Therefore, we present CsMPDB (miRNAs and phasiRNAs database of Camellia sinensis), an interactive web application with multiple analysis modules developed to visualize and explore miRNA and phasiRNA in tea plants based on 259 sRNA-seq samples and 24 degradome-seq samples in NCBI. The source code for the CsMPDB was written in R/shiny. It is compatible, extendable, and portable to be easily set up on different operating systems, and can be accessed at http://myshiny.cpolar.io/CsMPDB. This application plays an important role in accelerating the functional study of sRNAs in the transcriptional regulation of tea plants, and has important reference value for the development of sRNA databases of other species.

MicroRNA-Mediated Post-Transcriptional Regulation of Enzymes Involved in Herbicide Resistance in Echinochloa oryzicola (Vasinger) Vasinger

Herbicide resistance is an emerging phytosanitary threat, causing serious yield and economic losses. Although this phenomenon has been widely studied, only recently has the role of epigenetic factors in its occurrence been considered. In the present study, we analyzed the microRNA-mediated regulation in Echinochloa oryzicola (Vasinger) Vasinger (late-watergrass) of the expression of cytochromes P450, glutathione S-transferase (GST), and eIF4B, all of which are enzymes involved in profoxydim (AURA®) detoxification. Before and after profoxydim application, the expression profiles of microRNAs (miRNAs) were selected for their ability to target the genes considered, and their targets were assessed by means of RT-qPCR. Susceptible and resistant biotypes showed different responses to this herbicide. After profoxydim application, in resistant biotypes, osa-miR2099-5p, ath-miR396b, osa-miR395f, osa-miR396a-5p, osa-miR166a-5p, osa-miR166d-5p, gra-miR8759, and gma-miR396f were not triggered, allowing the expression of CYP81A, GSTF1, and eIF4B genes and the herbicide's detoxification. Meanwhile, the transcription of ata-miR166c-5p, ath-miR847, osa-miR5538, and gra-miR7487c was triggered, down-regulating CYP71AK2, CYP72A254, CYP72A122, and EcGST expression. In susceptible biotypes, the herbicide stimulated ata-miR166c-5p, ath-miR847, osa-miR5538, gra-miR7487c, osa-miR166a-5p, and gra-miR8759, down-regulating their respective target genes (CYP72A122, CYP71AK2, EcGST, CYP72A254, CYP81A12, and eIF4B). A better understanding of the role of miRNA-mediated epigenetic regulation in herbicide resistance will be useful in planning more targeted and sustainable methods for controlling this phytosanitary threat.

Validation of suitable reference microRNAs for qRT-PCR in Osmanthus fragrans under abiotic stress, hormone and metal ion treatments

Introduction

Sweet osmanthus (Osmanthus fragrans) is a prominent woody ornamental plant extensively utilized in horticulture, the food industry, cosmetics, and traditional Chinese medicine. MicroRNAs (miRNAs) are crucial regulators of gene regulation, playing a vital role in enabling plants to adapt to environmental fluctuations. Despite their significance, research on miRNA expression in O. fragrans under adverse stress conditions remains limited. Therefore, the selection of appropriate reference miRNAs is essential to ensure accurate miRNA expression analysis.

Methods

In this study, qRT-PCR technology was combined with four algorithms (i.e., delta-Ct, geNorm, NormFinder, and BestKeeper) to systematically evaluate the expression stability of 14 candidate miRNAs across eleven environmental conditions, including under abiotic stress, under hormone and metal ion treatments, during flower opening and senescence, and across various tissues.

Results

The results revealed that under hormone treatments, ofr-miR159b-3p, novel8, and novel3 exhibited high expression stability; under abiotic stress, ofr-miR159b-3p, novel8, ofr-miR403-3p, and novel2 demonstrated considerable stability; during metal ion treatments, novel3, ofr-miR159b-3p, novel33, novel2, and ofr-miR395e were identified as stable miRNAs; in different tissues, novel2 and ofr-miR395e were relatively stable; and during flower opening and senescence, novel33 and ofr-miR395e maintained stable expression.

Discussion

This study represents the first comprehensive assessment of reference miRNA stability in O. fragrans, providing a reliable framework for miRNA expression analysis under diverse conditions, including flower development and senescence, abiotic stress, hormone treatments, and metal ion treatments. These findings carry significant implications for future research into the function of miRNAs.

In silico identification of chilli genome encoded MicroRNAs targeting the 16S rRNA and secA genes of "Candidatus phytoplasma trifolii"

Phytoplasma, a potentially hazardous pathogen associated with witches' broom, is an economically harmful disease-producing bacteria that damages chilli cultivation. Phytoplasma-infected plants display various symptoms that indicate significant disruptions in normal plant physiology and behaviour. Diseases caused by phytoplasma are widespread and have a major economic impact on crop quality and yield. This work focuses on identifying and examining chilli microRNAs (miRNAs) as potential targets against the 16S rRNA and secA gene of "Candidatus Phytoplasma trifolii" ("Ca. P. trifolii") through plant miRNA prediction algorithms. Mature chilli miRNAs (CA-miRNAs) were collected and used to hybridise the 16S rRNA and secA genes. A total of four common CA-miRNAs were picked according to genetic consensus. Three algorithms applied in the present study suggested that the physiologically relevant, top-ranked miR169b_2 has a possibly specific site at nucleotide position 1,006 for targeting the 'Ca. P. trifolii' 16S rRNA gene. The circos algorithm was then utilised to create the miRNA-mRNA regulatory network. The free energy between the miRNA:mRNA duplex was also computed, and the best value of -17.46 kcal/mol was obtained for CA-miR166c_2. Currently, there are no suitable commercial 'Ca. P. trifolii'-resistant chilli crops. As a result, the expected biological data provide useful evidence for developing 'Ca. P. trifolii'-resistant chilli plants.

DNA self-assembly-boosted transcription amplification coupled with CRISPR/Cas13a system for plant microRNA analysis

MicroRNAs (miRNAs) play important roles in the growth process of plants, and some food-originated plant miRNAs have potential impacts on human health, which makes the detection of plant miRNAs of great significance. However, plant miRNAs are naturally modified with 2'-O-methyl at the 3'-terminal, which is difficult to be directly quantified by enzyme-catalyzed terminal polymerization protocols. Herein, we have proposed a simple strategy by coupling DNA self-assembly-boosted transcription amplification with CRISPR/Cas13a platform (termed as Cas13a-SATA) for the specific and sensitive detection of plant miRNA. In the Cas13a-SATA, the plant miRNA will mediate DNA self-assembly on the surface of microbeads and then trigger efficient transcription amplification to yield numerous single-stranded RNA (ssRNA) molecules, which can effectively activate the Cas13a trans-cleavage activity to generate intense fluorescence signal in a plant miRNA dosage-responsive manner. Using the Cas13a-SATA, we have realized the sensitive detection of plant miR156a with the limit of detection (LOD) down to 3.8 fM. Furthermore, Cas13a-SATA has been successfully applied to the accurate quantification of miR156a in Arabidopsis and maize, demonstrating its feasibility in analyzing plant miRNAs in real biological samples.

Machine learning-aided microRNA discovery for olive oil quality

MicroRNAs (miRNAs) are key regulators of gene expression in plants, influencing various biological processes such as oil quality and seed development. Although, our knowledge about miRNAs in olive (Olea europaea L.) is progressing, with several miRNAs being identified in previous studies, but most of these reported miRNAs have been predicted without the aid of a reference genome, primarily due to limited genome accessibility at the time. However, significant knowledge gaps still need to be improved in this area. This study addresses the complexities of miRNA detection in olive, using a high quality reference genome and a combination of genomics and machine learning-based methods. By leveraging random forest and support vector machine algorithms, we successfully identified 56 novel miRNAs in olive, surpassing the limitations of conventional homology-based methods. Our subsequent analysis revealed that some of these miRNAs are implicated in the regulation of key genes involved in oil quality. Within the context of oil biosynthesis pathways, the novel miRNA Oeu124369 regulates fatty acid biosynthesis by targeting acetyl-CoA acyltransferase 1 and palmitoyl-protein thioesterase, thereby influencing the production of acetyl-CoA and palmitic acid, respectively. These findings underscore the power of machine learning in unraveling the complex miRNA regulatory network in olive and provide a high quality miRNA resource for future research aimed at improving olive oil production by exploring the target genes of the identified miRNAs to understand their role and their biological processes.

Identification of Black Cumin (Nigella sativa) MicroRNAs by Next-Generation Sequencing and Their Implications in Secondary Metabolite Biosynthesis

Secondary metabolites are bioactive compounds believed to contribute to the pharmacological properties of plants. MicroRNAs (miRNAs) are small non-coding RNA molecules involved in post-transcriptional regulation and are thought to play an important role in regulating secondary metabolism biosynthesis. Nevertheless, the extent of miRNA involvement in secondary metabolism remains minimal. Nigella sativa (black cumin/black seed) is a popular medicinal and culinary plant known for its pharmaceutical properties; however, its genomic information is scarce. In this study, next-generation sequencing (NGS) technology was employed to obtain the miRNA profile of N. sativa, and their involvement in secondary metabolite biosynthesis was explored. A total of 25,139,003 unique reads ranging from 16 to 40 nucleotides were attained, out of which 240 conserved and 34 novel miRNAs were identified. Moreover, 6083 potential target genes were recognized in this study. Several conserved and novel black cumin miRNAs were found to target enzymes involved in the terpenoid, diterpenoid, phenylpropanoid, carotenoid, flavonoid, steroid, and ubiquinone biosynthetic pathways, among others, for example, beta-carotene 3-hydroxylase, gibberellin 3 beta-dioxygenase, trimethyltridecatetraene synthase, carboxylic ester hydrolases, acetyl-CoA C-acetyltransferase, isoprene synthase, peroxidase, shikimate O-hydroxycinnamoyltransferase, etc. Furthermore, sequencing data were validated through qPCR by checking the relative expression of eleven randomly selected conserved and novel miRNAs (nsa-miR164d, nsa-miR166a, nsa-miR167b, nsa-miR171a, nsa-miR390b, nsa-miR396, nsa-miR159a, nsa-miRN1, nsa-miRN29, nsa-miRN32, and nsa-miRN34) and their expression patterns were found to be corroborated with the sequencing data. We anticipate that this work will assist in clarifying the implications of miRNAs in plant secondary metabolism and aid in the generation of artificial miRNA-based strategies to overproduce highly valuable secondary metabolites from N. sativa.

Silencing Osa-miR827 via CRISPR/Cas9 protects rice against the blast fungus Magnaporthe oryzae

MicroRNAs (miRNAs) are short, non-coding RNAs that regulate gene expression at the post-transcriptional level. In plants, miRNAs participate in diverse developmental processes and adaptive responses to biotic and abiotic stress. MiR827 has long been recognized to be involved in plant responses to phosphate starvation. In rice, the miR827 regulates the expression of OsSPX-MFS1 and OsSPX-MFS2, these genes encoding vacuolar phosphate transporters. In this study, we demonstrated that miR827 plays a role in resistance to infection by the fungus Magnaporthe oryzae in rice. We show that MIR827 overexpression enhances susceptibility to infection by M. oryzae which is associated to a weaker induction of defense gene expression during pathogen infection. Conversely, CRISPR/Cas9-induced mutations in the MIR827 gene completely abolish miR827 production and confer resistance to M. oryzae infection. This resistance is accompanied by a reduction of leaf Pi content compared to wild-type plants, whereas Pi levels increase in leaves of the blast-susceptible miR827 overexpressor plants. In wild-type plants, miR827 accumulation in leaves decreases during the biotrophic phase of the infection process. Taken together, our data indicates that silencing MIR827 confers resistance to M. oryzae infection in rice while further supporting interconnections between Pi signaling and immune signaling in plants. Unravelling the role of miR827 during M. oryzae infection provides knowledge to improve blast resistance in rice by CRISPR/Cas9-editing of MIR827.

Insights into the Early Steps of the Symbiotic Interaction between Soybean (Glycine max) and Sinorhizobium fredii Symbiosis Using Transcriptome, Small RNA, and Degradome Sequencing

Symbiotic nitrogen fixation carried out by the soybean-rhizobia symbiosis increases soybean yield and reduces the amount of nitrogen fertilizer that has been applied. MicroRNAs (miRNAs) are crucial in plant growth and development, prompting an investigation into their role in the symbiotic interaction of soybean with partner rhizobia. Through integrated small RNA, transcriptome, and degradome sequencing analysis, 1215 known miRNAs, 314 of them conserved, and 187 novel miRNAs were identified, with 44 differentially expressed miRNAs in soybean roots inoculated with Sinorhizobium fredii HH103 and a ttsI mutant. The study unveiled that the known miRNA gma-MIR398a-p5 was downregulated in the presence of the ttsI mutation, while the target gene of gma-MIR398a-p5, Glyma.06G007500, associated with nitrogen metabolism, was upregulated. The results of this study offer insights for breeding high-efficiency nitrogen-fixing soybean varieties, enhancing crop yield and quality.

Deciphering miRNA-lncRNA-mRNA interaction through experimental validation of miRNAs, lncRNAs, and miRNA targets on mRNAs in Cajanus cajan

Pigeon pea (Cajanus cajan) is widely cultivated for its nutritional and medicinal value yet remains an orphan crop as productivity has not been improved because of a lack of genome and non-coding genome information. Non-coding RNAs, like miRNAs and long non-coding RNAs (lncRNAs), are involved in regulation of growth, metabolism, development, and stress response, and have a critical role in post-transcriptional gene regulation (PTGR). We attempted to elucidate the roles of miRNAs and lncRNAs in pigeon pea through experimental validation of computationally predicted miRNAs and lncRNAs and targets of miRNAs on mRNAs. We experimentally validated 20 miRNAs and 11 lncRNAs. We predicted cleavage sites of three miRNA targets: serine/threonine-protein kinase, polygalacturonase, beta-galactosidase. We identified 469 targets of 265 miRNAs and their functional annotations using computational methods. We built a miRNA-mRNA-lncRNA network model, with the miRNAs targeting both mRNAs and lncRNAs, to obtain information on the interplay of these three molecules. A confirmed interaction through experimental validation was established between miRNA, namely cca-miR1535a targeting the mRNA for beta-galactosidase, as well as the lncRNA cca-lnc-020033. Our findings increase knowledge of the non-coding genome of pigeon pea and their roles in PTGR and in improving agronomic traits of this pulse crop.

Species-Specific miRNAs Contribute to the Divergence between Deciduous and Evergreen Species in Ilex

MicroRNAs (miRNAs) are pivotal regulators of gene expression, playing crucial roles in plant developmental processes and environmental responses. However, the function of miRNAs in influencing deciduous traits has been little explored. Here, we utilized sRNA-seq on two deciduous species, Ilex polyneura (Hand.-Mazz.) S. Y. Hu and Ilex asprella Champ. ex Benth., along with an evergreen species, Ilex latifolia Thunb., to identify and annotate miRNAs within these species. Our analysis revealed 162 species-specific miRNAs (termed SS-miRNAs) from 120 families, underscoring the fundamental roles and potential influence of SS-miRNAs on plant phenotypic diversity and adaptation. Notably, three SS-miRNAs in I. latifolia were found to target crucial genes within the abscission signaling pathway. Analysis of cis-regulatory elements suggested a novel regulatory relationship that may contribute to the evergreen phenotype of I. latifolia by modulating the abscission process in a light-independent manner. These findings propose a potential mechanism by which SS-miRNAs can influence the conserved abscission pathway, contributing to the phenotypic divergence between deciduous and evergreen species within the genus Ilex.

Maize miRNAs and their putative target genes involved in chilling stress response in 5-day old seedlings

BACKGROUND

In the context of early sowing of maize as a promising adaptation strategy that could significantly reduce the negative effects of climate change, an in-depth understanding of mechanisms underlying plant response to low-temperature stress is demanded. Although microRNAs (miRNAs) have been recognized as key regulators of plant stress response, research on their role in chilling tolerance of maize during early seedling stages is scarce. Therefore, it is of great significance to explore chilling-responsive miRNAs, reveal their expression patterns and associated target genes, as well as to examine the possible functions of the conserved and novel miRNAs. In this study, the role of miRNAs was examined in 5d-old maize seedlings of one tolerant and one sensitive inbred line exposed to chilling (10/8 °C) stress for 6 h and 24 h, by applying high throughput sequencing.

RESULTS

A total of 145 annotated known miRNAs belonging to 30 families and 876 potentially novel miRNAs were identified. Differential expression (DE) analysis between control and stress conditions identified 98 common miRNAs for both genotypes at one time point and eight miRNAs at both time points. Target prediction and enrichment analysis showed that the DE zma-miR396, zma-miR156, zma-miR319, and zma-miR159 miRNAs modulate growth and development. Furthermore, it was found that several other DE miRNAs were involved in abiotic stress response: antioxidative mechanisms (zma-miR398), signal transduction (zma-miR156, zma-miR167, zma-miR169) and regulation of water content (zma-miR164, zma-miR394, zma-miR396). The results underline the zma-miRNAs involvement in the modulation of their target genes expression as an important aspect of the plant's survival strategy and acclimation to chilling stress conditions.

CONCLUSIONS

To our understanding, this is the first study on miRNAs in 5-d old seedlings' response to chilling stress, providing data on the role of known and novel miRNAs post-transcriptional regulation of expressed genes and contributing a possible platform for further network and functional analysis.

Insights into cucumber (Cucumis sativus) genetics: Genome-wide discovery and computational analysis of the Calreticulin Domain-Encoding gene (CDEG) family

Cucumber is an essential vegetable crop throughout the world. Cucumber development is vital for accomplishing both quality and productivity requirements. Meanwhile, numerous factors have resulted in substantial cucumber losses. However, the calreticulin domain-encoding genes (CDEGs) in cucumber were not well-characterized and had little function. In the genome-wide association study (GWAS), we recognized and characterized the CDEGs in Cucumis sativus (cucumber). Through a comprehensive study of C. sativus, our research has unveiled the presence of three unique genes, denoted as CsCRTb, CsCRT3, and CsCNX1, unevenly distributed on three chromosomes in the genome of C. sativus. In accordance to the phylogenetic investigation, these genes may be categorized into three subfamilies. Based on the resemblance with AtCDE genes, we reorganized the all CsCDE genes in accordance with international nomenclature. The expression analysis and cis-acting components revealed that each of CsCDE gene promoter region enclosed number of cis-elements connected with hormone and stress response. According to subcellular localization studies demonstrated that, they were found in deferent locations of the cell such as endoplasmic reticulum, plasma membrane, golgi apparatus, and vacuole, according to subcellular localization studies. Chromosomal distribution analysis and synteny analysis demonstrated the probability of segmental or tandem duplications within the cucumber CDEG gene family. Additionally, miRNAs displayed diverse modes of action, including mRNA cleavage and translational inhibition. We used the RNA seq data to analyze the expression of CDEG genes in response to cold stress and also improved cold tolerance, which was brought on by treating cucumber plants to an exogenous chitosan oligosaccharide spray. Our investigation revealed that these genes responded to this stress in a variety of ways, demonstrating that they may adapt quickly to environmental changes in cucumber plants. This study provides a base for further understanding in reference to CDE gene family and reveals that genes play significant functions in cucumber stress responses.

Evidence of microRNAs origination from chloroplast genome and their role in regulating Photosystem II protein N (psbN) mRNA

The microRNAs are endogenous, regulating gene expression either at the DNA or RNA level. Despite the availability of extensive studies on microRNA generation in plants, reports on their abundance, biogenesis, and consequent gene regulation in plant organelles remain naVve. Building on previous studies involving pre-miRNA sequencing in Abelmoschus esculentus, we demonstrated that three putative microRNAs were raised from the chloroplast genome. In the current study, we have characterized the genesis of these three microRNAs through a combination of bioinformatics and experimental approaches. The gene sequence for a miRNA, designated as AecpmiRNA1 (A. esculentus chloroplast miRNA), is potentially located in both the genomic DNA, i.e., nuclear and chloroplast genome. In contrast, the gene sequences for the other two miRNAs (AecpmiRNA2 and AecpmiRNA3) are exclusively present in the chloroplast genome. Target prediction revealed many potential mRNAs as targets for AecpmiRNAs. Further analysis using 5' RACE-PCR determined the AecpmiRNA3 binding and cleavage site at the photosystem II protein N (psbN). These results indicate that AecpmiRNAs are generated from the chloroplast genome, possessing the potential to regulate mRNAs arising from chloroplast gene(s). On the other side, the possibility of nuclear genome-derived mRNA regulation by AecpmiRNAs cannot be ruled out.

Non-coding RNAs and leaf senescence: Small molecules with important roles

Non-coding RNAs (ncRNAs) are a special class of functional RNA molecules that are not translated into proteins. ncRNAs have emerged as pivotal regulators of diverse developmental processes in plants. Recent investigations have revealed the association of ncRNAs with the regulation of leaf senescence, a complex and tightly regulated developmental process. However, a comprehensive review of the involvement of ncRNAs in the regulation of leaf senescence is still lacking. This manuscript aims to summarize the molecular mechanisms underlying ncRNAs-mediated leaf senescence and the potential applications of ncRNAs to manipulate the onset and progression of leaf senescence. Various classes of ncRNAs, including microRNAs (miRNAs), small interfering RNAs (siRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs), are discussed in terms of their regulatory mechanisms in leaf senescence. Furthermore, we explore the interactions between ncRNA and the key regulators of senescence, including transcription factors as well as core components in phytohormone signaling pathways. We also discuss the possible challenges and approaches related to ncRNA-mediated leaf senescence. This review contributes to a further understanding of the intricate regulatory network involving ncRNAs in leaf senescence.

Comprehensive analysis of OsJAZ gene family deciphers rhizobacteria-mediated nutrient stress modulation in rice

The JASMONATE-ZIM DOMAIN (JAZ) repressors are crucial proteins in jasmonic acid signaling pathway that are critical for plant growth. Therefore, the present study aimed to identify and characterize OsJAZs in the rice genome, revealing their structural attributes, regulatory elements, miRNA interactions, and subcellular localization. 23 JAZ transcripts across the 6 chromosomes of rice genome were identified having conserved domains and different physiochemical characteristics. Phylogenetically classified into five clades, they showed highest syntenic relationship with P. virgatum. The non-synonymous/synonymous values ranged from 0.44 to 1.21 suggesting purifying/stabilizing selection in OsJAZs. The study examined the 1.5 kb promoter region for cis-regulatory elements, and also identified 92 miRNAs targets. Furthermore, homology modeling provided insights into the 3D-structures of JAZ proteins while in-silico gene expression analysis revealed their functional diversity in various tissues and developmental stages. Additionally, qRT-PCR analysis highlighted their involvement in stress adaptation to sub-optimum nutrient conditions induced by plant-beneficial rhizobacteria Bacillus amyloliquefaciens (SN13) in two rice varieties. Distinct OsJAZ expression patterns in the two varieties correlated with altered root architecture, xylem structure, and lignification. These findings affirmed that specific up-or down-regulation of OsJAZs might play critical role in SN13 induced changes in the two varieties that enabled them to survive under stress.

Machine learning in photosynthesis: Prospects on sustainable crop development

Improving photosynthesis is a promising avenue to increase food security. Studying photosynthetic traits with the aim to improve efficiency has been one of many strategies to increase crop yield but analyzing large data sets presents an ongoing challenge. Machine learning (ML) represents a ubiquitous tool that can provide a more elaborate data analysis. Here we review the application of ML in various domains of photosynthetic research, as well as in photosynthetic pigment studies. We highlight how correlating hyperspectral data with photosynthetic parameters to improve crop yield could be achieved through various ML algorithms. We also propose strategies to employ ML in promoting photosynthetic pigment research for furthering crop yield.

Comparative expression analysis of sucrose phosphate synthase gene family in a low and high sucrose Pakistani sugarcane cultivars

Sugarcane is the world's largest cultivated crop by biomass and is the main source of sugar and biofuel. Sucrose phosphate synthase (SPS) enzymes are directly involved in the synthesis of sucrose. Here, we analyzed and compared one of the important gene families involved in sucrose metabolism in a high and low sucrose sugarcane cultivar. A comprehensive in silico analysis of the SoSPS family displayed their phylogenetic relationship, gene and protein structure, miRNA targets, protein interaction network (PPI), gene ontology and collinearity. This was followed by a spatial expression analysis in two different sugarcane varieties. The phylogenetic reconstruction distributed AtSPS, ZmSPS, OsSPS, SoSPS and SbSPS into three main groups (A, B, C). The regulatory region of SoSPS genes carries ABRE, ARE, G-box, and MYC as the most dominant cis-regulatory elements. The PPI analysis predicted a total of 14 unique proteins interacting with SPS. The predominant expression of SPS in chloroplast clearly indicates that they are the most active in the organelle which is the hub of photosynthesis. Similarly, gene ontology attributed SPS to sucrose phosphate synthase and glucosyl transferase molecular functions, as well as sucrose biosynthetic and disaccharide biological processes. Overall, the expression of SPS in CPF252 (high sucrose variety) was higher in leaf and culm as compared to that of CPF 251 (low sucrose variety). In brief, this study adds to the present literature about sugarcane, sucrose metabolism and role of SPS in sucrose metabolism thereby opening up further avenues of research in crop improvement.

microRNA408 and its encoded peptide regulate sulfur assimilation and arsenic stress response in Arabidopsis

MicroRNAs (miRNAs) are small non-coding RNAs that play a central role in regulating various developmental and biological processes. The expression of miRNAs is differentially modulated in response to various biotic and abiotic stresses. Recent findings have shown that some pri-miRNAs encode small regulatory peptides known as microRNA-encoded peptides (miPEPs). miPEPs regulate the growth and development of plants by modulating corresponding miRNA expression; however, the role of these peptides under different stress conditions remains unexplored. Here, we report that pri-miR408 encodes a small peptide, miPEP408, that regulates the expression of miR408, its targets, and associated phenotype in Arabidopsis. We also report that miR408, apart from Plantacyanin (ARPN) and Laccase3 (LAC3), targets a glutathione S-transferase (GSTU25) that plays a role in sulfur assimilation and exhibits a range of detoxification activities with the environmental pollutant. Plants overexpressing miR408 showed severe sensitivity under low sulfur (LS), arsenite As(III), and LS + As(III) stress, while miR408 mutants developed using the CRISPR/Cas9 approach showed tolerance. Transgenic lines showed phenotypic alteration and modulation in the expression of genes involved in the sulfur reduction pathway and affect sulfate and glutathione accumulation. Similar to miR408 overexpressing lines, the exogenous application of synthetic miPEP408 and miPEP408OX lines led to sensitivity in plants under LS, As(III), and combined LS + As(III) stress compared to the control. This study suggests the involvement of miR408 and miPEP408 in heavy metal and nutrient deficiency responses through modulation of the sulfur assimilation pathway.

Potential role of the regulatory miR1119-MYC2 module in wheat (Triticum aestivum L.) drought tolerance

MicroRNA (miRNA)-target gene modules are essential components of plants' abiotic stress signalling pathways Little is known about the drought-responsive miRNA-target modules in wheat, but systems biology approaches have enabled the prediction of these regulatory modules and systematic study of their roles in responses to abiotic stresses. Using such an approach, we sought miRNA-target module(s) that may be differentially expressed under drought and non-stressed conditions by mining Expressed Sequence Tag (EST) libraries of wheat roots and identified a strong candidate (miR1119-MYC2). We then assessed molecular and physiochemical differences between two wheat genotypes with contrasting drought tolerance in a controlled drought experiment and assessed possible relationships between their tolerance and evaluated traits. We found that the miR1119-MYC2 module significantly responds to drought stress in wheat roots. It is differentially expressed between the contrasting wheat genotypes and under drought versus non-stressed conditions. We also found significant associations between the module's expression profiles and ABA hormone content, water relations, photosynthetic activities, H₂O₂ levels, plasma membrane damage, and antioxidant enzyme activities in wheat. Collectively, our results suggest that a regulatory module consisting of miR1119 and MYC2 may play an important role in wheat's drought tolerance.

What Do We Know about Barley miRNAs?

Plant miRNAs are powerful regulators of gene expression at the post-transcriptional level, which was repeatedly proved in several model plant species. miRNAs are considered to be key regulators of many developmental, homeostatic, and immune processes in plants. However, our understanding of plant miRNAs is still limited, despite the fact that an increasing number of studies have appeared. This systematic review aims to summarize our current knowledge about miRNAs in spring barley (Hordeum vulgare), which is an important agronomical crop worldwide and serves as a common monocot model for studying abiotic stress responses as well. This can help us to understand the connection between plant miRNAs and (not only) abiotic stresses in general. In the end, some future perspectives and open questions are summarized.

Genome-wide and molecular characterization of the DNA replication helicase 2 (DNA2) gene family in rice under drought and salt stress

Rice plants experience various biotic (such as insect and pest attack) and abiotic (such as drought, salt, heat, and cold etc.) stresses during the growing season, resulting in DNA damage and the subsequent losses in rice production. DNA Replication Helicase/Nuclease2 (DNA2) is known to be involved in DNA replication and repair. In animals and yeast DNA2 are well characterized because it has the abilities of both helicase and nuclease, it plays a crucial role in DNA replication in the nucleus and mitochondrial genomes. However; they are not fully examined in plants due to less focused on plants damage repair. To fill this research gap, the current study focused on the genome-wide identification and characterization of OsDNA2 genes, along with analyses of their transcriptional expression, duplication, and phylogeny in rice. Overall, 17 OsDNA2 members were reported to be found on eight different chromosomes (2, 3, 4, 6, 7, 9, 10, and 11). Among these chromosomes (Chr), Chr4 contained a maximum of six OsDNA2 genes. Based on phylogenetic analysis, the OsDNA2 gene members were clustered into three different groups. Furthermore, the conserved domains, gene structures, and cis-regulatory elements were systematically investigated. Gene duplication analysis revealed that OsDNA2_2 had an evolutionary relationship with OsDNA2_14, OsDNA2_5 with OsDNA2_6, and OsDNA2_1 with OsDNA2_8. Moreover, results showed that the conserved domain (AAA_11 superfamily) were present in the OsDNA2 genes, which belongs to the DEAD-like helicase superfamily. In addition, to understand the post-transcriptional modification of OsDNA2 genes, miRNAs were predicted, where 653 miRNAs were reported to target 17 OsDNA2 genes. The results indicated that at the maximum, OsDNA2_1 and OsDNA2_4 were targeted by 74 miRNAs each, and OsDNA2_9 was less targeted (20 miRNAs). The three-dimensional (3D) structures of 17 OsDNA2 proteins were also predicted. Expression of OsDNA2 members was also carried out under drought and salt stresses, and conclusively their induction indicated the possible involvement of OsDNA2 in DNA repair under stress when compared with the control. Further studies are recommended to confirm where this study will offer valuable basic data on the functioning of DNA2 genes in rice and other crop plants.

Integrated analyses reveal the response of peanut to phosphorus deficiency on phenotype, transcriptome and metabolome

BACKGROUND

Phosphorus (P) is one of the most essential macronutrients for crops. The growth and yield of peanut (Arachis hypogaea L.) are always limited by P deficiency. However, the transcriptional and metabolic regulatory mechanisms were less studied. In this study, valuable phenotype, transcriptome and metabolome data were analyzed to illustrate the regulatory mechanisms of peanut under P deficiency stress.

RESULT

In present study, two treatments of P level in deficiency with no P application (-P) and in sufficiency with 0.6 mM P application (+ P) were used to investigate the response of peanut on morphology, physiology, transcriptome, microRNAs (miRNAs), and metabolome characterizations. The growth and development of plants were significantly inhibited under -P treatment. A total of 6088 differentially expressed genes (DEGs) were identified including several transcription factor family genes, phosphate transporter genes, hormone metabolism related genes and antioxidant enzyme related genes that highly related to P deficiency stress. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that 117 genes were annotated in the phenylpropanoid biosynthesis pathway under P deficiency stress. A total of 6 miRNAs have been identified significantly differential expression between + P and -P group by high-throughput sequencing of miRNAs, including two up-regulated miRNAs (ahy-miR160-5p and ahy-miR3518) and four down-regulated miRNAs (ahy-miR408-5p, ahy-miR408-3p, ahy-miR398, and ahy-miR3515). Further, the predicted 22 target genes for 6 miRNAs and cis-elements in 2000 bp promoter region of miRNA genes were analyzed. A total of 439 differentially accumulated metabolites (DAMs) showed obviously differences in two experimental conditions.

CONCLUSIONS

According to the result of transcripome and metabolome analyses, we can draw a conclusion that by increasing the content of lignin, amino acids, and levan combining with decreasing the content of LPC, cell reduced permeability, maintained stability, raised the antioxidant capacity, and increased the P uptake in struggling for survival under P deficiency stress.

Identification and differential analysis of noncoding RNAs in response to drought in Phyllostachys aureosulcata f. spectabilis

The role of noncoding RNAs (ncRNAs) in plant resistance to abiotic stresses is increasingly being discovered. Drought stress is one of the most common stresses that affecting plant growth, and high intensity drought has a significant impact on the normal growth of plants. In this study, a high-throughput sequencing was performed on plant tissue samples of Phyllostachys aureosulcata f. spectabilis C. D. Chu et C. S. Chao by drought treatment for 0, 2, 4 and 6 days. The sequencing results were analysed bioinformatically. We detected 336,946 RNAs among all 12 samples, including 192,098 message RNAs (mRNAs), 142,761 long noncoding RNAs (lncRNAs), 1,670 circular RNAs (circRNAs), and 417 microRNAs (miRNAs). We detected 2,419 differentially expressed (DE) ncRNAs, including 213 DE circRNAs, 2,088 DE lncRNAs and 118 DE miRNAs. Then, we used Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to functionally predict DE ncRNAs. The results showed that most DE ncRNAs are involved in the response to drought stress, mainly in biochemical reactions involved in some metabolites, as well as in organelle activities. In addition, we validated two random circRNAs and demonstrated their circularity. We also found a stable internal reference gene available for Phyllostachys aureosulcata f. spectabilis and validated the accuracy of this experiment by quantitative real-time polymerase chain reaction (qRT-PCR).

De novo assembly and characterization of the draft genome of the cashew (Anacardium occidentale L.)

Cashew is the second most important tree nut crop in the global market. Cashew is a diploid and heterozygous species closely related to the mango and pistachio. Its improvement by conventional breeding is slow due to the long juvenile phase. Despite the economic importance, very little genomics/transcriptomics information is available for cashew. In this study, the Oxford nanopore reads and Illumina reads were used for de novo assembly of the cashew genome. The hybrid assembly yielded a 356.6 Mb genome corresponding to 85% of the estimated genome size (419 Mb). The BUSCO analysis showed 91.8% of genome completeness. Transcriptome mapping showed 92.75% transcripts aligned with the assembled genome. Gene predictions resulted in the identification of 31,263 genes coding for a total of 35,000 gene isoforms. About 46% (165 Mb) of the cashew genome comprised of repetitive sequences. Phylogenetic analyses of the cashew with nine species showed that it was closely related to Mangifera indica. Analysis of cashew genome revealed 3104 putative R-genes. The first draft assembly of the genome, transcriptome and R gene information generated in this study would be the foundation for understanding the molecular basis of economic traits and genomics-assisted breeding in cashew.

The Involvement of microRNAs in Plant Lignan Biosynthesis—Current View

Lignans, as secondary metabolites synthesized within a phenylpropanoid pathway, play various roles in plants, including their involvement in growth and plant defense processes. The health and nutritional benefits of lignans are unquestionable, and many studies have been devoted to these attributes. Although the regulatory role of miRNAs in the biosynthesis of secondary metabolites has been widely reported, there is no systematic review available on the miRNA-based regulatory mechanism of lignans biosynthesis. However, the genetic background of lignan biosynthesis in plants is well characterized. We attempted to put together a regulatory mosaic based on current knowledge describing miRNA-mediated regulation of genes, enzymes, or transcription factors involved in this biosynthesis process. At the same time, we would like to underline the fact that further research is necessary to improve our understanding of the miRNAs regulating plant lignan biosynthesis by exploitation of current approaches for functional identification of miRNAs.

MicroRNA: a new signal in plant-to-plant communication

Plants can communicate inter- and intraspecifically using signals transmitted via root exudate and volatiles released into the atmosphere. A recent study by Betti et al. discovered that miRNA is one of the signals used during plant communication. MiRNAs are secreted by plants and change the gene expression in neighbouring plants.

Novel miRNA-SSRs for Improving Seed Hardness Trait of Pomegranate (Punica granatum L.)

Present research discovered novel miRNA-SSRs for seed type trait from 761 potential precursor miRNA sequences of pomegranate. SSR mining and BLASTx of the unique sequences identified 69 non-coding pre-miRNA sequences, which were then searched for BLASTn homology against Dabenzi genome. Sixty three true pri-miRNA contigs encoding 213 pre-miRNAs were predicted. Analysis of the resulting sequences enabled discovery of SSRs within pri-miRNA (227) and pre-miRNA sequences (79). A total of 132 miRNA-SSRs were developed for seed type trait from 63 true pri-miRNAs, of which 46 were specific to pre-miRNAs. Through ePCR, 123 primers were validated and mapped on eight Tunisia chromosomes. Further, 80 SSRs producing specific amplicons were ePCR-confirmed on multiple genomes i.e. Dabenzi, Taishanhong, AG2017 and Tunisia, yielding a set of 63 polymorphic SSRs (polymorphism information content ≥0.5). Of these, 32 miRNA-SSRs revealed higher polymorphism level (89.29%) when assayed on six pomegranate genotypes. Furthermore, target prediction and network analysis suggested a possible association of miRNA-SSRs i.e. miRNA_SH_SSR69, miRNA_SH_SSR36, miRNA_SH_SSR103, miRNA_SH_SSR35 and miRNA_SH_SSR53 with seed type trait. These miRNA-SSRs would serve as important genomic resource for rapid and targeted improvement of seed type trait of pomegranate.

Transcriptome and Small RNA Profiling of Potato Virus Y Infected Potato Cultivars, Including Systemically Infected Russet Burbank

Potatoes are the world’s most produced non-grain crops and an important food source for billions of people. Potatoes are susceptible to numerous pathogens that reduce yield, including Potato virus Y (PVY). Genetic resistance to PVY is a sustainable way to limit yield and quality losses due to PVY infection. Potato cultivars vary in their susceptibility to PVY and include susceptible varieties such as Russet Burbank, and resistant varieties such as Payette Russet. Although the loci and genes associated with PVY-resistance have been identified, the genes and mechanisms involved in limiting PVY during the development of systemic infections have yet to be fully elucidated. To increase our understanding of PVY infection, potato antiviral responses, and resistance, we utilized RNA sequencing to characterize the transcriptomes of two potato cultivars. Since transcriptional responses associated with the extreme resistance response occur soon after PVY contact, we analyzed the transcriptome and small RNA profile of both the PVY-resistant Payette Russet cultivar and PVY-susceptible Russet Burbank cultivar 24 h post-inoculation. While hundreds of genes, including terpene synthase and protein kinase encoding genes, exhibited increased expression, the majority, including numerous genes involved in plant pathogen interactions, were downregulated. To gain greater understanding of the transcriptional changes that occur during the development of systemic PVY-infection, we analyzed Russet Burbank leaf samples one week and four weeks post-inoculation and identified similarities and differences, including higher expression of genes involved in chloroplast function, photosynthesis, and secondary metabolite production, and lower expression of defense response genes at those time points. Small RNA sequencing identified different populations of 21- and 24-nucleotide RNAs and revealed that the miRNA profiles in PVY-infected Russet Burbank plants were similar to those observed in other PVY-tolerant cultivars and that during systemic infection ~32% of the NLR-type disease resistance genes were targeted by 21-nt small RNAs. Analysis of alternative splicing in PVY-infected potato plants identified splice variants of several hundred genes, including isoforms that were more dominant in PVY-infected plants. The description of the PVYN-Wi-associated transcriptome and small RNA profiles of two potato cultivars described herein adds to the body of knowledge regarding differential outcomes of infection for specific PVY strain and host cultivar pairs, which will help further understanding of the mechanisms governing genetic resistance and/or virus-limiting responses in potato plants.

A prescient evolutionary model for genesis, duplication and differentiation of MIR160 homologs in Brassicaceae

MicroRNA160 is a class of nitrogen-starvation responsive genes which governs establishment of root system architecture by down-regulating AUXIN RESPONSE FACTOR genes (ARF10, ARF16 and ARF17) in plants. The high copy number of MIR160 variants discovered by us from land plants, especially polyploid crop Brassicas, posed questions regarding genesis, duplication, evolution and function. Absence of studies on impact of whole genome and segmental duplication on retention and evolution of MIR160 homologs in descendent plant lineages prompted us to undertake the current study. Herein, we describe ancestry and fate of MIR160 homologs in Brassicaceae in context of polyploidy driven genome re-organization, copy number and differentiation. Paralogy amongst Brassicaceae MIR160a, MIR160b and MIR160c was inferred using phylogenetic analysis of 468 MIR160 homologs from land plants. The evolutionarily distinct MIR160a was found to represent ancestral form and progenitor of MIR160b and MIR160c. Chronology of evolutionary events resulting in origin and diversification of genomic loci containing MIR160 homologs was delineated using derivatives of comparative synteny. A prescient model for causality of segmental duplications in establishment of paralogy in Brassicaceae MIR160, with whole genome duplication accentuating the copy number increase, is being posited in which post-segmental duplication events viz. differential gene fractionation, gene duplications and inversions are shown to drive divergence of chromosome segments. While mutations caused the diversification of MIR160a, MIR160b and MIR160c, duplicated segments containing these diversified genes suffered gene rearrangements via gene loss, duplications and inversions. Yet the topology of phylogenetic and phenetic trees were found congruent suggesting similar evolutionary trajectory. Over 80% of Brassicaceae genomes and subgenomes showed a preferential retention of single copy each of MIR160a, MIR160b and MIR160c suggesting functional relevance. Thus, our study provides a blue-print for reconstructing ancestry and phylogeny of MIRNA gene families at genomics level and analyzing the impact of polyploidy on organismal complexity. Such studies are critical for understanding the molecular basis of agronomic traits and deploying appropriate candidates for crop improvement.

Sequence and functional analysis of MIR319 promoter homologs from Brassica juncea reveals regulatory diversification and altered expression under stress

KEY MESSAGE

Extensive regulatory divergence during development, abiotic stress and ABA regime observed amongst promoter homologs and homeologs of MIR319 from Brassica juncea. Gene duplication followed by sub-functionalization, neo-functionalization, and pseudogenization are routes to functional and adaptive diversification. The influence of polyploidy on protein-coding genes is well investigated but little is known about their impact on transcriptional regulation of MIRNA gene family. The present study was therefore performed with an aim to uncover regulatory diversification of MIR319 homologs and homeologs in Brassica juncea. We employed comparative genomics to identify and isolate six promoter homologs of MIR319 from B. juncea. Regulatory diversification was studied using analysis of reporter activity driven by BjMIR319 promoters in a heterologous system employing promoter-reporter fusion constructs. MIR319 is known to play important roles in leaf and flower development, and multiple stress responses. Reporter activity was therefore monitored during development, hormonal and stress regimes. In-silico analyses revealed differential distribution of cis-regulatory motifs and functional analysis revealed distinct spatiotemporal expression patterns. The significance of presence of selected cis-regulatory motifs corresponding to heat, cold, salt and ABA stress were further functionally validated. It was observed that promoter of Bj -MIR319a-A01 was upregulated in response to cold and salt stress, while promoter of Bj -MIR319c-A04 (D1) and Bj -MIR319c-A05 (FL) were downregulated in response to high temperature. In summary, comparative analysis of homologous promoters from Brassica juncea, an allopolyploid revealed extensive sequence and functional diversity. Spatiotemporal activity of reporter gene driven by BjMIR319 promoter was distinct, and partially overlapping with from those reported previously for A. thaliana. The present study clearly demonstrates regulatory divergence amongst promoter homologs of MIR319 in Brassica juncea during development and stress response, and underlines the urgent need for dissection of promoter function and detailed characterization including identification of interacting trans-factors. Genbank accession numbers: MT379853-MT379858.

Metal and Metalloid Toxicity in Plants: An Overview on Molecular Aspects

Worldwide, the effects of metal and metalloid toxicity are increasing, mainly due to anthropogenic causes. Soil contamination ranks among the most important factors, since it affects crop yield, and the metals/metalloids can enter the food chain and undergo biomagnification, having concomitant effects on human health and alterations to the environment. Plants have developed complex mechanisms to overcome these biotic and abiotic stresses during evolution. Metals and metalloids exert several effects on plants generated by elements such as Zn, Cu, Al, Pb, Cd, and As, among others. The main strategies involve hyperaccumulation, tolerance, exclusion, and chelation with organic molecules. Recent studies in the omics era have increased knowledge on the plant genome and transcriptome plasticity to defend against these stimuli. The aim of the present review is to summarize relevant findings on the mechanisms by which plants take up, accumulate, transport, tolerate, and respond to this metal/metalloid stress. We also address some of the potential applications of biotechnology to improve plant tolerance or increase accumulation.

Identification of Small RNAs During High Light Acclimation in Arabidopsis thaliana

The biological significance of non-coding RNAs (ncRNAs) has been firmly established to be important for the regulation of genes involved in stress acclimation. Light plays an important role for the growth of plants providing the energy for photosynthesis; however, excessive light conditions can also cause substantial defects. Small RNAs (sRNAs) are a class of non-coding RNAs that regulate transcript levels of protein-coding genes and mediate epigenetic silencing. Next generation sequencing facilitates the identification of small non-coding RNA classes such as miRNAs (microRNAs) and small-interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), but changes in the ncRNA transcriptome in response to high light are poorly understood. We subjected Arabidopsis plants to high light conditions and performed a temporal in-depth study of the transcriptome data after 3 h, 6 h, and 2 days of high light treatment. We identified a large number of high light responsive miRNAs and sRNAs derived from NAT gene pairs, lncRNAs and TAS transcripts. We performed target predictions for differentially expressed miRNAs and correlated their expression levels through mRNA sequencing data. GO analysis of the targets revealed an overrepresentation of genes involved in transcriptional regulation. In A. thaliana, sRNA-mediated regulation of gene expression in response to high light treatment is mainly carried out by miRNAs and sRNAs derived from NAT gene pairs, and from lncRNAs. This study provides a deeper understanding of sRNA-dependent regulatory networks in high light acclimation.

Characterizing the Role of the miR156-SPL Network in Plant Development and Stress Response

MicroRNA (miRNA) is a short, single-stranded, non-coding RNA found in eukaryotic cells that can regulate the expression of many genes at the post-transcriptional level. Among various plant miRNAs with diverse functions, miR156 plays a key role in biological processes, including developmental regulation, immune response, metabolic regulation, and abiotic stress. MiRNAs have become the regulatory center for plant growth and development. MicroRNA156 (miR156) is a highly conserved and emerging tool for the improvement of plant traits, including crop productivity and stress tolerance. Fine-tuning of squamosa promoter biding-like (SPL) gene expression might be a useful strategy for crop improvement. Here, we studied the regulation of the miR156 module and its interaction with SPL factors to understand the developmental transition of various plant species. Furthermore, this review provides a strong background for plant biotechnology and is an important source of information for further molecular breeding to optimize farming productivity.