Genome-wide identification of soybean microRNAs and their targets reveals their organ-specificity and responses to phosphate starvation

microRNAs: Key Regulators in Plant Responses to Abiotic and Biotic Stresses via Endogenous and Cross-Kingdom Mechanisms

Dramatic shifts in global climate have intensified abiotic and biotic stress faced by plants. Plant microRNAs (miRNAs)-20-24 nucleotide non-coding RNA molecules-form a key regulatory system of plant gene expression; playing crucial roles in plant growth; development; and defense against abiotic and biotic stress. Moreover, they participate in cross-kingdom communication. This communication encompasses interactions with other plants, microorganisms, and insect species, collectively exerting a profound influence on the agronomic traits of crops. This article comprehensively reviews the biosynthesis of plant miRNAs and explores their impact on plant growth, development, and stress resistance through endogenous, non-transboundary mechanisms. Furthermore, this review delves into the cross-kingdom regulatory effects of plant miRNAs on plants, microorganisms, and pests. It proceeds to specifically discuss the design and modification strategies for artificial miRNAs (amiRNAs), as well as the protection and transport of miRNAs by exosome-like nanovesicles (ELNVs), expanding the potential applications of plant miRNAs in crop breeding. Finally, the current limitations associated with harnessing plant miRNAs are addressed, and the utilization of synthetic biology is proposed to facilitate the heterologous expression and large-scale production of miRNAs. This novel approach suggests a plant-based solution to address future biosafety concerns in agriculture.

miR2118 Negatively Regulates Bacterial Blight Resistance through Targeting Several Disease Resistance Genes in Rice

Plant miRNAs are a class of noncoding RNA with a length of 21-24 nt that play an important role in plant responses to biotic and abiotic stresses. Bacterial blight (BB) caused by Xanthomonas oryzae pv. oryzae (Xoo) is one of the most serious bacterial diseases in rice. Our previous work showed that osa-miR2118b/n was induced by Xoo infection. However, the biological function of miR2118 has not yet been characterized in experiments. Herein, we constructed MIR2118b OE, as well as single and double mutants of MIR2118b/n using CRISPR/Cas9. Further results showed that osa-MIR2118b OE plants exhibited longer lesion lengths than the wild type after Xoo inoculation, while MIR2118 CRISPR plants exhibited shorter lesion lengths than the wild type after Xoo inoculation. Co-transformation experiments in rice protoplasts indicated that osa-miR2118 negatively regulated the transcripts of three nucleotide-binding sites and leucine-rich repeat (NLR) genes (LOC_Os08g42700.1, LOC_Os01g05600.1, and LOC_Os12g37290.1) which are predicted target genes of miR2118, but not the mutated NLR genes with a 3 bp insertion at the center of the binding sites. The transcriptional level of the three NLR genes was reversed relative to osa-miR2118 in the MIR2118b OE and MIR2118b CRISPR plants. The above results demonstrate that osa-miR2118b/n negatively regulates the resistance to bacterial blight through negatively regulating several NLR genes.

Mechanistic basis for mitigating drought tolerance by selenium application in tobacco (Nicotiana tabacum L.): a multi-omics approach

The lack of irrigation water in agricultural soils poses a significant constraint on global crop production. In-depth investigation into microRNAs (miRNAs) has been widely used to achieve a comprehensive understanding of plant defense mechanisms. However, there is limited knowledge on the association of miRNAs with drought tolerance in cigar tobacco. In this study, a hydroponic experiment was carried out to identify changes in plant physiological characteristics, miRNA expression and metabolite profile under drought stress, and examine the mitigating effects of selenium (Se) application. The shoot dry weight of drought-stressed plants was approximately half (50.3%) of that in non-stressed (control) conditions. However, plants supplied with Se attained 38.8% greater shoot dry weight as compared to plants with no Se supply under drought stress. Thirteen miRNAs were identified to be associated with drought tolerance. These included 7 known (such as nta-miR156b and nta-miR166a) and 6 novel miRNAs (such as novel-nta-miR156-5p and novel-nta-miR209-5p) with the target genes of squamosa promoter-binding-like protein 4 (SPL4), serine/threonine protein phosphatase 2A (PPP2A), cation/calcium exchanger 4-like (CCX4), extensin-1-like (EXT1) and reduced wall acetylation 2 (RWA2). Further investigation revealed that the expression levels of Ext1 and RWA2 were significantly decreased under drought stress but increased with Se addition. Moreover, key metabolites such as catechin and N-acetylneuraminic acid were identified, which may play a role in the regulation of drought tolerance. The integrated analysis of miRNA sequencing and metabolome highlighted the significance of the novel-nta-miR97-5p- LRR-RLK- catechin pathway in regulating drought tolerance. Our findings provide valuable insights into the molecular mechanisms underlying drought tolerance and Se-induced stress alleviation in cigar tobacco.

Medicago truncatula PHO2 genes have distinct roles in phosphorus homeostasis and symbiotic nitrogen fixation

Three PHO2-like genes encoding putative ubiquitin-conjugating E2 enzymes of Medicago truncatula were characterized for potential roles in phosphorous (P) homeostasis and symbiotic nitrogen fixation (SNF). All three genes, MtPHO2A, B and C, contain miR399-binding sites characteristic of PHO2 genes in other plant species. Distinct spatiotemporal expression patterns and responsiveness of gene expression to P- and N-deprivation in roots and shoots indicated potential roles, especially for MtPHO2B, in P and N homeostasis. Phenotypic analysis of pho2 mutants revealed that MtPHO2B is integral to Pi homeostasis, affecting Pi allocation during plant growth under nutrient-replete conditions, while MtPHO2C had a limited role in controlling Pi homeostasis. Genetic analysis also revealed a connection between Pi allocation, plant growth and SNF performance. Under N-limited, SNF conditions, Pi allocation to different organs was dependent on MtPHO2B and, to a lesser extent, MtPHO2C and MtPHO2A. MtPHO2A also affected Pi homeostasis associated with nodule formation. Thus, MtPHO2 genes play roles in systemic and localized, i.e., nodule, P homeostasis affecting SNF.

INDETERMINATE1 autonomously regulates phosphate homeostasis upstream of the miR399-ZmPHO2 signaling module in maize

The macronutrient phosphorus is essential for plant growth and development. Plants have evolved multiple strategies to increase the efficiency of phosphate (Pi) acquisition to protect themselves from Pi starvation. However, the crosstalk between Pi homeostasis and plant development remains to be explored. Here, we report that overexpressing microRNA399 (miR399) in maize (Zea mays) is associated with premature senescence after pollination. Knockout of ZmPHO2 (Phosphate 2), a miR399 target, resulted in a similar premature senescence phenotype. Strikingly, we discovered that INDETERMINATE1 (ID1), a floral transition regulator, inhibits the transcription of ZmMIR399 genes by directly binding to their promoters, alleviating the repression of ZmPHO2 by miR399 and ultimately contributing to the maintenance of Pi homeostasis in maize. Unlike ZmMIR399 genes, whose expression is induced by Pi deficiency, ID1 expression was independent of the external inorganic orthophosphate status, indicating that ID1 is an autonomous regulator of Pi homeostasis. Furthermore, we show that ZmPHO2 was under selection during maize domestication and cultivation, resulting in a more sensitive response to Pi starvation in temperate maize than in tropical maize. Our study reveals a direct functional link between Pi-deprivation sensing by the miR399-ZmPHO2 regulatory module and plant developmental regulation by ID1.

Identification of microRNAs involved in the Phosphate starvation response in Oil Palm (Elaeis guineensis Jacq.)

BACKGROUND

Plant microRNA, often known as miRNA, is a novel form of gene expression regulator that is known to play a significant role in phosphate starvation. The identification of microRNAs involved in the response to phosphate starvation in oil palms is beneficial for breeding programs.

METHOD

The main nursery stage seedlings of two oil palm progenies were treated with three different fertiliser namely: complete fertiliser with urea, P2O5, K2O, and MgO based on the standard procedure as a control (C); fertiliser with urea, K2O, MgO without P2O5 (P0); and no fertiliser (F0) for 24 weeks. A total of six oil palm roots were subjected to RNA isolation, followed by miRNA sequencing using the Illumina HiSeq 4000 platform, and all reads were computationally analysed.

RESULTS

In total, 119 potential miRNAs related to 5,891 genes were identified. The P-specific miRNAs were assumed based on the miRNAs that identified without P fertilizer treatment, resulted of twenty miRNA sequences in the treatment comparison of (C vs P0) vs (C vs F0). Those 20 miRNA sequences were grouped into 9 families, namely EgmiR319; EgmiR399; EgmiR396; EgmiR172; EgmiR156; EgmiR157; miR5648; miR5645; and EgmiRNA_unidentified. Two miRNAs were selected for RT-qPCR validation, namely EgMir399 and EgMir172. Their expression pattern was similar with the RNA sequencing results and shown opposite expression pattern with their target genes, UBC E2 24 and APETALA2, respectively.

CONCLUSIONS

The nine micro RNA families was identified in oil palm root tissue at phosphate starvation.

Plant phosphate status influences root biotic interactions

Phosphorus (P) deficiency stress in combination with biotic stress(es) severely impacts crop yield. Plant responses to P deficiency overlapping with that of other stresses exhibit a high degree of complexity involving different signaling pathways. On the one hand, plants engage with rhizosphere microbiome/arbuscular mycorrhizal fungi for improved phosphate (Pi) acquisition and plant stress response upon Pi deficiency; on the other hand, this association is gets disturbed under Pi sufficiency. This nutrient-dependent response is highly regulated by the phosphate starvation response (PSR) mediated by the master regulator, PHR1, and its homolog, PHL. It is interesting to note that Pi status (deficiency/sufficiency) has a varying response (positive/negative) to different biotic encounters (beneficial microbes/opportunistic pathogens/insect herbivory) through a coupled PSR-PHR1 immune system. This also involves crosstalk among multiple players including transcription factors, defense hormones, miRNAs, and Pi transporters, among others influencing the plant-biotic-phosphate interactions. We provide a comprehensive view of these key players involved in maintaining a delicate balance between Pi homeostasis and plant immunity. Finally, we propose strategies to utilize this information to improve crop resilience to Pi deficiency in combination with biotic stresses.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

Cooperative interactions between nitrogen fixation and phosphorus nutrition in legumes

Legumes such as soybean are considered important crops as they provide proteins and oils for humans and livestock around the world. Different from other crops, leguminous crops accumulate nitrogen (N) for plant growth through symbiotic nitrogen fixation (SNF) in coordination with rhizobia. A number of studies have shown that efficient SNF requires the cooperation of other nutrients, especially phosphorus (P), a nutrient deficient in most soils. During the last decades, great progress has been made in understanding the molecular mechanisms underlying the interactions between SNF and P nutrition, specifically through the identification of transporters involved in P transport to nodules and bacteroids, signal transduction, and regulation of P homeostasis in nodules. These studies revealed a distinct N-P interaction in leguminous crops, which is characterized by specific signaling cross talk between P and SNF. This review aimed to present an updated picture of the cross talk between N fixation and P nutrition in legumes, focusing on soybean as a model crop, and Medicago truncatula and Lotus japonicus as model plants. We also discuss the possibilities for enhancing SNF through improving P nutrition, which are important for high and sustainable production of leguminous crops.

Functional Analysis of Bna-miR399c-PHO2 Regulatory Module Involved in Phosphorus Stress in Brassica napus

Phosphorus stress is one of the important factors restricting plant growth and development, and the microRNA (miRNA) family is involved in the regulation of the response to plant nutrient stress by repressing the expression of target genes at the post-transcriptional or translational level. miR399 is involved in the transportation of phosphate in multiple plants by improving tolerance to low Pi conditions. However, the effect of miR399 on the response of low Pi stress in rapeseed (Brassica napus L.) is unclear. The present study showed a significant increase in taproot length and lateral root number of plants overexpressing Bna-miR399c, while the biomass and Pi accumulation in shoots and roots increased, and the anthocyanin content decreased and chlorophyll content improved under low Pi stress. The results illustrate that Bna-miR399c could enhance the uptake and transportation of Pi in soil, thus making B. napus more tolerant to low Pi stress. Furthermore, we confirmed that BnPHO2 is one of the targets of Bna-miR399c, and the rejection of Pi in rapeseed seedlings increased due to the overexpression of BnPHO2. Hence, we suggest that miR399c-PHO2 module can effectively regulate the homeostasis of Pi in B. napus. Our study can also provide the theoretical basis for germplasm innovation and the design of intelligent crops with low nutrient input and high yield to achieve the dual objectives of income and yield increase and environmental protection in B. napus.

Ectopic expression of GmNF-YA8 in Arabidopsis delays flowering via modulating the expression of gibberellic acid biosynthesis- and flowering-related genes and promotes lateral root emergence in low phosphorus conditions

NUCLEAR FACTOR Y subunit alpha (NF-YA), together with NF-YB and NF-YC, regulates plant growth and development, as well as plant responses to biotic and abiotic stresses. Although extensive studies have examined the functions of NF-YAs in Arabidopsis thaliana, the roles of NF- YAs in Glycinme max are poorly understood. In this study, we identified a phosphorus (P) starvation-responsive NF-YA8 in soybean. The expression of GmNF-YA8 is induced by low P or low nitrogen in leaves, but not by potassium or iron starvation, respectively. GmNF-YA8 is localized in the nucleus and plasma membrane. Ectopic expression of GmNF-YA8 inhibits plant growth and delayed flowering in Arabidopsis. Exogenous application of gibberellic acid (GA) rescues the delayed flowering phenotype in Arabidopsis overexpressing GmNF-YA8 lines GmNF-YA8OE-05 and GmNF-YA8OE-20. Moreover, quantitative real time PCR (qRT-PCR) verified that overexpression of GmNF-YA8 downregulates GA20ox2 and GA3ox2 expression, but upregulates GA2ox2 and GA2ox3 that encode enzymes, which inactive bioactive GAs. Consistent with the late flowering phenotype of Arabidopsis trangenic lines that overexpress GmNF-YA8, the transcript levels of flowering-promoting genes AP1, CO, LFY, and SOC1 are reduced. In addition, overexpression of GmNF-YA8 promotes the emergence of lateral root (LR) primordium from epidermis rather than the initiation of LR in low P, and increases the LR density in low nitrogen. Our results provide insights into the roles of GmNF-YA8.

Roles of non-coding RNAs in the hormonal and nutritional regulation in nodulation and nitrogen fixation

Symbiotic nitrogen fixation is an important component in the nitrogen cycle and is a potential solution for sustainable agriculture. It is the result of the interactions between the plant host, mostly restricted to legume species, and the rhizobial symbiont. From the first encounter between the host and the symbiont to eventual successful nitrogen fixation, there are delicate processes involved, such as nodule organogenesis, rhizobial infection thread progression, differentiation of the bacteroid, deregulation of the host defense systems, and reallocation of resources. All these processes are tightly regulated at different levels. Recent evidence revealed that non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), participate in these processes by controlling the transcription and translation of effector genes. In general, ncRNAs are functional transcripts without translation potential and are important gene regulators. MiRNAs, negative gene regulators, bind to the target mRNAs and repress protein production by causing the cleavage of mRNA and translational silencing. LncRNAs affect the formation of chromosomal loops, DNA methylation, histone modification, and alternative splicing to modulate gene expression. Both lncRNAs and circRNAs could serve as target mimics of miRNA to inhibit miRNA functions. In this review, we summarized and discussed the current understanding of the roles of ncRNAs in legume nodulation and nitrogen fixation in the root nodule, mainly focusing on their regulation of hormone signal transduction, the autoregulation of nodulation (AON) pathway and nutrient homeostasis in nodules. Unraveling the mediation of legume nodulation by ncRNAs will give us new insights into designing higher-performance leguminous crops for sustainable agriculture.

Integrated mRNA and microRNA expression analysis of root response to phosphate deficiency in Medicago sativa

The deficiency of available phosphate significantly limits plant growth and development. This study sought to investigate how alfalfa (Medicago sativa), a high-yielding and high-quality forage widely cultivated worldwide, responds to phosphate deficiency stress by integrating transcriptional and post-transcriptional data. In this study, 6,041 differentially expressed genes (DEGs) were identified in alfalfa roots under phosphate deficiency conditions. Furthermore, psRNATarget, RNAhybrid, and TargetFinder were used to predict the target genes of 137 differentially expressed miRNAs (DEMs) in the root. In total, 3,912 DEGs were predicted as target genes. Pearson correlation analysis revealed 423 pairs of miRNA-mRNA regulatory relationships. MiRNA negatively regulates mRNA involved in regulatory pathways of phosphate deficiency responses in alfalfa. miR156e targeted squamosa promoter-binding-like protein 13A (SPL13), miR160c targeted auxin response factor 18 (ARF18), and miR2587a controlled glycolysis and citrate cycle via Phosphoenolpyruvate carboxykinase (ATP) (PCKA). Novel-miR27 regulated SPX domain-containing protein that controls phosphate transport in alfalfa root, novel-miR3-targeted sulfoquinovosyl transferase SQD2 controlled sulfolipid synthesis and glutathione S-transferase (GST; mediated by miR169j/k and novel-miR159) regulated glutathione metabolism. miR399l regulated auxin-responsive protein SAUR72 involved in IAA signal transduction, while abscisic acid receptor PYL4 (regulated by novel-miR205 and novel-miR83) participated in ABA signal transduction. Combined miRNA-mRNA enrichment analysis showed that most miRNAs regulate the phosphate starvation response of alfalfa by modulating target genes involved in carbohydrate metabolism, sulfolipid metabolism, glutathione metabolism, and hormone signal transduction. Therefore, this study provides new insights into the post-transcriptional regulation mechanism of phosphate deficiency responses and new perspectives on phosphate assimilation pathways in alfalfa and other legumes.

Epigenomics as Potential Tools for Enhancing Magnitude of Breeding Approaches for Developing Climate Resilient Chickpea

Epigenomics has become a significant research interest at a time when rapid environmental changes are occurring. Epigenetic mechanisms mainly result from systems like DNA methylation, histone modification, and RNA interference. Epigenetic mechanisms are gaining importance in classical genetics, developmental biology, molecular biology, cancer biology, epidemiology, and evolution. Epigenetic mechanisms play important role in the action and interaction of plant genes during development, and also have an impact on classical plant breeding programs, inclusive of novel variation, single plant heritability, hybrid vigor, plant-environment interactions, stress tolerance, and performance stability. The epigenetics and epigenomics may be significant for crop adaptability and pliability to ambient alterations, directing to the creation of stout climate-resilient elegant crop cultivars. In this review, we have summarized recent progress made in understanding the epigenetic mechanisms in plant responses to biotic and abiotic stresses and have also tried to provide the ways for the efficient utilization of epigenomic mechanisms in developing climate-resilient crop cultivars, especially in chickpea, and other legume crops.

Innovation and appropriation in mycorrhizal and rhizobial Symbioses

Most land plants benefit from endosymbiotic interactions with mycorrhizal fungi, including legumes and some nonlegumes that also interact with endosymbiotic nitrogen (N)-fixing bacteria to form nodules. In addition to these helpful interactions, plants are continuously exposed to would-be pathogenic microbes: discriminating between friends and foes is a major determinant of plant survival. Recent breakthroughs have revealed how some key signals from pathogens and symbionts are distinguished. Once this checkpoint has been passed and a compatible symbiont is recognized, the plant coordinates the sequential development of two types of specialized structures in the host. The first serves to mediate infection, and the second, which appears later, serves as sophisticated intracellular nutrient exchange interfaces. The overlap in both the signaling pathways and downstream infection components of these symbioses reflects their evolutionary relatedness and the common requirements of these two interactions. However, the different outputs of the symbioses, phosphate uptake versus N fixation, require fundamentally different components and physical environments and necessitated the recruitment of different master regulators, NODULE INCEPTION-LIKE PROTEINS, and PHOSPHATE STARVATION RESPONSES, for nodulation and mycorrhization, respectively.

A Multi-Level Iterative Bi-Clustering Method for Discovering miRNA Co-regulation Network of Abiotic Stress Tolerance in Soybeans

Although growing evidence shows that microRNA (miRNA) regulates plant growth and development, miRNA regulatory networks in plants are not well understood. Current experimental studies cannot characterize miRNA regulatory networks on a large scale. This information gap provides an excellent opportunity to employ computational methods for global analysis and generate valuable models and hypotheses. To address this opportunity, we collected miRNA-target interactions (MTIs) and used MTIs from Arabidopsis thaliana and Medicago truncatula to predict homologous MTIs in soybeans, resulting in 80,235 soybean MTIs in total. A multi-level iterative bi-clustering method was developed to identify 483 soybean miRNA-target regulatory modules (MTRMs). Furthermore, we collected soybean miRNA expression data and corresponding gene expression data in response to abiotic stresses. By clustering these data, 37 MTRMs related to abiotic stresses were identified, including stress-specific MTRMs and shared MTRMs. These MTRMs have gene ontology (GO) enrichment in resistance response, iron transport, positive growth regulation, etc. Our study predicts soybean MTRMs and miRNA-GO networks under different stresses, and provides miRNA targeting hypotheses for experimental analyses. The method can be applied to other biological processes and other plants to elucidate miRNA co-regulation mechanisms.

Differential microRNA expression, microRNA arm switching, and microRNA:long noncoding RNA interaction in response to salinity stress in soybean

BACKGROUND

Soybean is a major legume crop with high nutritional and environmental values suitable for sustainable agriculture. Noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), are important regulators of gene functions in eukaryotes. However, the interactions between these two types of ncRNAs in the context of plant physiology, especially in response to salinity stress, are poorly understood.

RESULTS

Here, we challenged a cultivated soybean accession (C08) and a wild one (W05) with salt treatment and obtained their small RNA transcriptomes at six time points from both root and leaf tissues. In addition to thoroughly analyzing the differentially expressed miRNAs, we also documented the first case of miRNA arm-switching (miR166m), the swapping of dominant miRNA arm expression, in soybean in different tissues. Two arms of miR166m target different genes related to salinity stress (chloroplastic beta-amylase 1 targeted by miR166m-5p and calcium-dependent protein kinase 1 targeted by miR166m-3p), suggesting arm-switching of miR166m play roles in soybean in response to salinity stress. Furthermore, two pairs of miRNA:lncRNA interacting partners (miR166i-5p and lncRNA Gmax_MSTRG.35921.1; and miR394a-3p and lncRNA Gmax_MSTRG.18616.1) were also discovered in reaction to salinity stress.

CONCLUSIONS

This study demonstrates how ncRNA involves in salinity stress responses in soybean by miRNA arm switching and miRNA:lncRNA interactions. The behaviors of ncRNAs revealed in this study will shed new light on molecular regulatory mechanisms of stress responses in plants, and hence provide potential new strategies for crop improvement.

The Evolution and Functional Roles of miR408 and Its Targets in Plants

MicroRNA408 (miR408) is an ancient and highly conserved miRNA, which is involved in the regulation of plant growth, development and stress response. However, previous research results on the evolution and functional roles of miR408 and its targets are relatively scattered, and there is a lack of a systematic comparison and comprehensive summary of the detailed evolutionary pathways and regulatory mechanisms of miR408 and its targets in plants. Here, we analyzed the evolutionary pathway of miR408 in plants, and summarized the functions of miR408 and its targets in regulating plant growth and development and plant responses to various abiotic and biotic stresses. The evolutionary analysis shows that miR408 is an ancient and highly conserved microRNA, which is widely distributed in different plants. miR408 regulates the growth and development of different plants by down-regulating its targets, encoding blue copper (Cu) proteins, and by transporting Cu to plastocyanin (PC), which affects photosynthesis and ultimately promotes grain yield. In addition, miR408 improves tolerance to stress by down-regulating target genes and enhancing cellular antioxidants, thereby increasing the antioxidant capacity of plants. This review expands and promotes an in-depth understanding of the evolutionary and regulatory roles of miR408 and its targets in plants.

Integration of small RNA, degradome, and transcriptome sequencing data illustrates the mechanism of low phosphorus adaptation in Camellia oleifera

Phosphorus (P) is an indispensable macronutrient for plant growth and development, and it is involved in various cellular biological activities in plants. Camellia oleifera is a unique high-quality woody oil plant that grows in the hills and mountains of southern China. However, the available P content is deficient in southern woodland soil. Until now, few studies focused on the regulatory functions of microRNAs (miRNAs) and their target genes under low inorganic phosphate (Pi) stress. In this study, we integrated small RNA, degradome, and transcriptome sequencing data to investigate the mechanism of low Pi adaptation in C. oleifera. We identified 40,689 unigenes and 386 miRNAs by the deep sequencing technology and divided the miRNAs into four different groups. We found 32 miRNAs which were differentially expressed under low Pi treatment. A total of 414 target genes of 108 miRNAs were verified by degradome sequencing. Gene ontology (GO) functional analysis of target genes found that they were related to the signal response to the stimulus and transporter activity, indicating that they may respond to low Pi stress. The integrated analysis revealed that 31 miRNA-target pairs had negatively correlated expression patterns. A co-expression regulatory network was established based on the profiles of differentially expressed genes. In total, three hub genes (ARF22, WRKY53, and SCL6), which were the targets of differentially expressed miRNAs, were discovered. Our results showed that integrated analyses of the small RNA, degradome, and transcriptome sequencing data provided a valuable basis for investigating low Pi in C. oleifera and offer new perspectives on the mechanism of low Pi tolerance in woody oil plants.

Noncoding-RNA-Mediated Regulation in Response to Macronutrient Stress in Plants

Macronutrient elements including nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), and sulfur (S) are required in relatively large and steady amounts for plant growth and development. Deficient or excessive supply of macronutrients from external environments may trigger a series of plant responses at phenotypic and molecular levels during the entire life cycle. Among the intertwined molecular networks underlying plant responses to macronutrient stress, noncoding RNAs (ncRNAs), mainly microRNAs (miRNAs) and long ncRNAs (lncRNAs), may serve as pivotal regulators for the coordination between nutrient supply and plant demand, while the responsive ncRNA-target module and the interactive mechanism vary among elements and species. Towards a comprehensive identification and functional characterization of nutrient-responsive ncRNAs and their downstream molecules, high-throughput sequencing has produced massive omics data for comparative expression profiling as a first step. In this review, we highlight the recent findings of ncRNA-mediated regulation in response to macronutrient stress, with special emphasis on the large-scale sequencing efforts for screening out candidate nutrient-responsive ncRNAs in plants, and discuss potential improvements in theoretical study to provide better guidance for crop breeding practices.

Soybean responds to phosphate starvation through reversible protein phosphorylation

Phosphorus (P) deficiency is considered as a major constraint on crop production. Although a set of adaptative strategies are extensively suggested in soybean (Glycine max) to phosphate (Pi) deprivation, molecular mechanisms underlying reversible protein phosphorylation in soybean responses to P deficiency remains largely unclear. In this study, isobaric tags for relative and absolute quantitation, combined with liquid chromatography and tandem mass spectrometry analysis was performed to identify differential phosphoproteins in soybean roots under Pi sufficient and deficient conditions. A total of 427 phosphoproteins were found to exhibit differential accumulations, with 213 up-regulated and 214 down-regulated. Among them, a nitrate reductase, GmNR4 exhibiting increased phosphorylation levels under low Pi conditions, was further selected to evaluate the effects of phosphorylation on its nitrate reductase activity and subcellular localization. Mutations of GmNR4 phosphorylation levels significantly influenced its activity in vitro, but not for its subcellular localization. Taken together, identification of differential phosphoproteins reveled the complex regulatory pathways for soybean adaptation to Pi starvation through reversible protein phosphorylation.

MicroRNAs: Tiny, powerful players of metal stress responses in plants

Metal contamination of the environment is a widespread problem threatening sustainable and safe crop production. Physio-biochemical and molecular mechanisms of plant responses to metal exposure have been studied to establish the best possible agronomical or biotechnological methods to tackle metal contamination. Metal stress tolerance is regulated by several molecular effectors among which microRNAs are one of the key master regulators of plant growth and stress responses in plants. MicroRNAs are known to coordinate multitude of plant responses to metal stress through antioxidant functions, root growth, hormonal signalling, transcription factors and metal transporters. The present review discusses integrative functions of microRNAs in the regulation of metal stress in plants, which will be useful for engineering stress tolerance traits for improved plant growth and productivity in metal stressed situations.

Differentially expressed microRNAs that target functional genes in mature soybean nodules

MicroRNAs (miRNAs) are important regulators of biological functions in plants. To find out what roles miRNAs play in regulating symbiotic nitrogen fixation (SNF) in soybean [Glycine max (L.) Merr.], we identified high-confidence differentially expressed (DE) miRNAs from uninoculated roots (UR), rhizobium-inoculated roots (IR), and nodules (NODs) of soybean by robust small RNA sequencing (sRNA-seq). Based on their predicted target messenger RNAs (mRNAs), the expression profiles of some of these DE miRNAs could be linked to nodule functions. In particular, several miRNAs associated with nutrient transportation genes were differentially expressed in IRs and mature NODs. MiR399b, specifically, was highly induced in IRs and NODs, as well as by inorganic phosphate (Pi) starvation. In composite soybean plants overexpressing miR399b, PHOSPHATE2 (PHO2), a known target of miR399b that inhibits the activities of high-affinity Pi transporters, was strongly repressed. In addition, the overexpression of miR399b in the roots of transgenic composite plants significantly improved whole-plant Pi and ureide concentrations and the overall growth in terms of leaf node numbers and whole-plant dry weight. Our findings suggest that the induction of miR399b in NODs could enhance nitrogen fixation and soybean growth, possibly via improving Pi uptake to achieve a better Pi-nitrogen balance to promote SNF in nodules.

Root-omics for drought tolerance in cool-season grain legumes

Root traits can be exploited to increase the physiological efficiency of crop water use under drought. Root length, root hairs, root branching, root diameter, and root proliferation rate are genetically defined traits that can help to improve the water productivity potential of crops. Recently, high-throughput phenotyping techniques/platforms have been used to screen the germplasm of major cool-season grain legumes for root traits and their impact on different physiological processes, including nutrient uptake and yield potential. Advances in omics approaches have led to the dissection of genomic, proteomic, and metabolomic structures of these traits. This knowledge facilitates breeders to improve the water productivity and nutrient uptake of cultivars under limited soil moisture conditions in major cool-season grain legumes that usually face terminal drought. This review discusses the advances in root traits and their potential for developing drought-tolerant cultivars in cool-season grain legumes.

Shoot-derived miR2111 controls legume root and nodule development

Legumes control their nodule numbers through the autoregulation of nodulation (AON). Rhizobia infection stimulates the production of root-derived CLE peptide hormones that are translocated to the shoot where they regulate a new signal. We used soybean to demonstrate that this shoot-derived signal is miR2111, which is transported via phloem to the root where it targets transcripts of Too Much Love (TML), a negative regulator of nodulation. Shoot perception of rhizobia-induced CLE peptides suppresses miR2111 expression, resulting in TML accumulation in roots and subsequent inhibition of nodule organogenesis. Feeding synthetic mature miR2111 via the petiole increased nodule numbers per plant. Likewise, elevating miR2111 availability by over-expression promoted nodulation, while target mimicry of TML induced the opposite effect on nodule development in wild-type plants and alleviated the supernodulating and stunted root growth phenotypes of AON-defective mutants. Additionally, in non-nodulating wild-type plants, ectopic expression of miR2111 significantly enhanced lateral root emergence with a decrease in lateral root length and average root diameter. In contrast, hairy roots constitutively expressing the target mimic construct exhibited reduced lateral root density. Overall, these findings demonstrate that miR2111 is both the critical shoot-to-root factor that positively regulates root nodule development and also acts to shape root system architecture.

PHO1 family members transport phosphate from infected nodule cells to bacteroids in Medicago truncatula

Legumes play an important role in the soil nitrogen availability via symbiotic nitrogen fixation (SNF). Phosphate (Pi) deficiency severely impacts SNF because of the high Pi requirement of symbiosis. Whereas PHT1 transporters are involved in Pi uptake into nodules, it is unknown how Pi is transferred from the plant infected cells to nitrogen-fixing bacteroids. We hypothesized that Medicago truncatula genes homologous to Arabidopsis PHO1, encoding a vascular apoplastic Pi exporter, are involved in Pi transfer to bacteroids. Among the seven MtPHO1 genes present in M. truncatula, we found that two genes, namely MtPHO1.1 and MtPHO1.2, were broadly expressed across the various nodule zones in addition to the root vascular system. Expressions of MtPHO1.1 and MtPHO1.2 in Nicotiana benthamiana mediated specific Pi export. Plants with nodule-specific downregulation of both MtPHO1.1 and MtPHO1.2 were generated by RNA interference (RNAi) to examine their roles in nodule Pi homeostasis. Nodules of RNAi plants had lower Pi content and a three-fold reduction in SNF, resulting in reduced shoot growth. Whereas the rate of 33Pi uptake into nodules of RNAi plants was similar to control, transfer of 33Pi from nodule cells into bacteroids was reduced and bacteroids activated their Pi-deficiency response. Our results implicate plant MtPHO1 genes in bacteroid Pi homeostasis and SNF via the transfer of Pi from nodule infected cells to bacteroids.

Profiling of MicroRNAs and Their Targets in Roots and Shoots Reveals a Potential MiRNA-Mediated Interaction Network in Response to Phosphate Deficiency in the Forestry Tree Betula luminifera

Inorganic phosphate (Pi) is often lacking in natural and agro-climatic environments, which impedes the growth of economically important woody species. Plants have developed strategies to cope with low Pi (LP) availability. MicroRNAs (miRNAs) play important roles in responses to abiotic stresses, including nutrition stress, by regulating target gene expression. However, the miRNA-mediated regulation of these adaptive responses and their underlying coordinating signals are still poorly understood in forestry trees such as Betula luminifera. Transcriptomic libraries, small RNA (sRNA) libraries, and a mixed degradome cDNA library of B. luminifera roots and shoots treated under LP and normal conditions (CK) were constructed and sequenced using next-generation deep sequencing. A comprehensive B. luminifera transcriptome derived from its roots and shoots was constructed, and a total of 76,899 unigenes were generated. Analysis of the transcriptome identified 8,095 and 5,584 differentially expressed genes in roots and shoots, respectively, under LP conditions. sRNA sequencing analyses indicated that 66 and 60 miRNAs were differentially expressed in roots and shoots, respectively, under LP conditions. A total of 109 and 112 miRNA-target pairs were further validated in the roots and shoots, respectively, using degradome sequencing. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of differential miRNA targets indicated that the "ascorbate and aldarate metabolism" pathway responded to LP, suggesting miRNA-target pairs might participating in the removing of reactive oxidative species under LP stress. Moreover, a putative network of miRNA-target interactions involved in responses to LP stress in B. luminifera is proposed. Taken together, these findings provide useful information to decipher miRNA functions and establish a framework for exploring P signaling networks regulated by miRNAs in B. luminifera and other woody plants. It may provide new insights into the genetic engineering of high use efficiency of Pi in forestry trees.

Quantitative Proteomics Reveals that GmENO2 Proteins Are Involved in Response to Phosphate Starvation in the Leaves of Glycine max L

Soybean (Glycine max L.) is a major crop providing important source for protein and oil for human life. Low phosphate (LP) availability is a critical limiting factor affecting soybean production. Soybean plants develop a series of strategies to adapt to phosphate (Pi) limitation condition. However, the underlying molecular mechanisms responsible for LP stress response remain largely unknown. Here, we performed a label-free quantification (LFQ) analysis of soybean leaves grown under low and high phosphate conditions. We identified 267 induced and 440 reduced differential proteins from phosphate-starved leaves. Almost a quarter of the LP decreased proteins are involved in translation processes, while the LP increased proteins are accumulated in chlorophyll biosynthetic and carbon metabolic processes. Among these induced proteins, an enolase protein, GmENO2a was found to be mostly induced protein. On the transcriptional level, GmENO2a and GmENO2b, but not GmENO2c or GmENO2d, were dramatically induced by phosphate starvation. Among 14 enolase genes, only GmENO2a and GmENO2b genes contain the P1BS motif in their promoter regions. Furthermore, GmENO2b was specifically induced in the GmPHR31 overexpressing soybean plants. Our findings provide molecular insights into how soybean plants tune basic carbon metabolic pathway to adapt to Pi deprivation through the ENO2 enzymes.

Transcriptional, metabolic, physiological and developmental responses of switchgrass to phosphorus limitation

Knowing how switchgrass (Panicum virgatum L.) responds and adapts to phosphorus (P)-limitation will aid efforts to optimize P acquisition and use in this species for sustainable biomass production. This integrative study investigated the impacts of mild, moderate, and severe P-stress on genome transcription and whole-plant metabolism, physiology and development in switchgrass. P-limitation reduced overall plant growth, increased root/shoot ratio, increased root branching at moderate P-stress, and decreased root diameter with increased density and length of root hairs at severe P-stress. RNA-seq analysis revealed thousands of genes that were differentially expressed under moderate and severe P-stress in roots and/or shoots compared to P-replete plants, with many stress-induced genes involved in transcriptional and other forms of regulation, primary and secondary metabolism, transport, and other processes involved in P-acquisition and homeostasis. Amongst the latter were multiple miRNA399 genes and putative targets of these. Metabolite profiling showed that levels of most sugars and sugar alcohols decreased with increasing P stress, while organic and amino acids increased under mild and moderate P-stress in shoots and roots, although this trend reversed under severe P-stress, especially in shoots.

Strontium stress disrupts miRNA biogenesis by reducing HYL1 protein levels in Arabidopsis

Strontium (Sr) is an emerging environmental pollutant that has become a major global concern after the nuclear accident at the Fukushima Daiichi Nuclear Power Plant in 2011. Although many studies have demonstrated the harmful effects of Sr on plant growth and development at the physiological level, knowledge regarding how plants sense and respond to Sr stress at the molecular level is limited. Recent studies have suggested that microRNAs (miRNAs) function as key regulators of plant growth and development as well as in the responses of plants to environmental stresses, including salinity, drought, cold, nutrient starvation, and heavy metals. In this study, we examined the global expression profile of miRNAs under Sr stress using small RNA sequencing analysis in Arabidopsis to better understand the molecular basis of plant responses to Sr stress. To identify specific Sr-responsive miRNAs, we performed comparative miRNA expression profiling analysis using control, CaCl₂-, and SrCl₂-treated seedlings. Compared to the control treatment, the expressions of most miRNAs were considerably decreased in the Sr-treated seedlings. However, under Sr stress, the expressions of primary miRNAs (pri-miRNAs) and their target genes were significantly increased; the protein levels of HYPONASTIC LEAVES 1 (HYL1), one of the core components of the microprocessor complex, were strongly reduced despite the increased HYL1 mRNA expression. In addition, hyl1-2 mutant plants were shown to be more sensitive to Sr stress than wild-type plants. Collectively, our results strongly suggested that Sr stress may be associated with the disruption of miRNA biogenesis by reducing the protein level of HYL1, which is required to maintain proper growth and development for plants. Our findings further indicated that some miRNAs may play important roles in plant responses to Sr stress.

Genome-wide identification of low phosphorus responsive microRNAs in two soybean genotypes by high-throughput sequencing

MicroRNAs (miRNAs) have been reported to be correlated with various stress responses in soybean, but only a few miRNAs have been demonstrated to respond to low phosphorus (LP) stress. To unravel the response mechanisms of miRNAs to low-P stress, the roots of two representative soybean genotypes with different P efficiency, Nannong94-156 (a LP-tolerant genotype) and Bogao (a LP-sensitive genotype), were used for the construction of RNA sequencing (RNA-seq) libraries under low/normal-P treatment by high-throughput sequencing. In total, 603 existing miRNAs and 1699 novel miRNAs belonging to 248 and 1582 families in all samples were identified, respectively. Among these miRNAs, 777 miRNAs were differentially expressed (DE) across different P levels and genotypes. Furthermore, putative targets of DE miRNAs were predicted, and these miRNAs mainly targeted ERF (ethylene responsive factor), auxin response factors (ARF), zinc finger protein, MYB, and NAC domain transcription factors. Gene ontology (GO) analysis showed that targets of DE miRNAs were significantly enriched in binding, metabolic processes, biological regulation, response to stress, and phosphorus metabolic processes. In addition, the expression profiles of chosen P-responsive miRNAs and target genes were validated by quantitative real-time PCR (qRT-PCR). Our study focused on genome-wide miRNA identification in two representative soybean genotypes under low-P stress. Overall, the DE miRNAs across different P levels and genotypes and their putative target genes will provide useful information for further study of miRNAs mediating low-P response and facilitate improvements in soybean breeding.

Global identification and characterization of miRNA family members responsive to potassium deprivation in wheat (Triticum aestivum L.)

Potassium (K) is essential for plant growth and stress responses. MicroRNAs (miRNAs) are involved in adaptation to nutrient deprivation through modulating gene expression. Here, we identified the miRNAs responsive to K deficiency in Triticum aestivum based on high-throughput small RNA sequencing analyses. Eighty-nine miRNAs, including 68 previously reported ones and 21 novel ones, displayed differential expression under K deficiency. In Gene Ontology and Kyoto Encyclopedia and Genome analyses, the putative target genes of the differentially expressed miRNAs were categorized into functional groups associated with ADP-binding activity, secondary metabolic pathways, and biosynthesis and metabolism. Functional characterization of tae-miR408, an miRNA significantly down-regulated under K deficiency, revealed its important role in mediating low-K tolerance. Compared with wild type, transgenic tobacco lines overexpressing tae-miR408 showed significantly improved K uptake, biomass, photosynthesis, and reactive oxygen species scavenging under K deficiency. These results show that distinct miRNAs function in the plant response to K deficiency through regulating target genes involved in energy metabolism and various secondary metabolic pathways. Our findings shed light on the plant response to K deficiency mediated by miRNAs in T. aestivum. Distinct miRNAs, such as tae-miR408, are valuable targets for generating crop varieties with improved K-use efficiency.

Genome wide identification, characterization and validation of novel miRNA-based SSR markers in pomegranate (Punica granatum L.)

A total of 17,439 mature miRNAs (~ 21 nt) earlier generated through RNA seq in the pomegranate were used for in silico analysis. After complexity reduction, a total of 1922 representative mature miRNAs were selected and used as query sequences against pomegranate genome to retrieve 2540 homologous contigs with flanking regions (~ 800). By using pre-miRNA prediction web server, a total of 1028 true contigs harbouring pri-miRNAs encoding 1162 pre-miRNAs were identified. Survey of these sequences for SSRs yielded a total of 1358 and 238 SSRs specific to pri-miRNA and pre-miRNAs, respectively. Of these, primer pairs were designed for 897 pri-miRNA and 168 pre-miRNA SSRs. In pri-miRNA sequences, hexa-nucleotides repeats were found to be most abundant (44.18%) followed by mono- (18.41%) and di-nucleotide (17.01%), which is also observed in pre-miRNA sequences. Further, a set of 51 randomly selected pre-miRNA-SSRs was examined for marker polymorphism. The experimental validation of these markers on eight pomegranate genotypes demonstrated 92.15% polymorphism. Utility of these functional markers was confirmed via examination of genetic diversity of 18 pomegranate genotypes using 15 miRNA-SSRs. Further, potential application of miRNA-SSRs for discovery of trait specific candidate genes was showed by validating 51 mature miRNA against publically available 2047 EST sequences of pomegranate by target and network analysis. In summary, the current study offers novel functional molecular markers for pomegranate genetic improvement.

Novel insights into expansion and functional diversification of MIR169 family in tomato

MAIN CONCLUSION

Expansion of MIR169 members by duplication and new mature forms, acquisition of new promoters, differential precursor-miRNA processivity and engaging novel targets increase the functional diversification of MIR169 in tomato. MIR169 family is an evolutionarily conserved miRNA family in plants. A systematic in-depth analysis of MIR169 family in tomato is lacking. We report 18 miR169 precursors, annotating new loci for MIR169a, b and d, as well as 3 novel mature isoforms (MIR169f/g/h). The family has expanded by both tandem- and segmental-duplication events during evolution. A tandem-pair MIR169b/b-1 and MIR169b-2/h is polycistronic in nature coding for three MIR169b isoforms and a new variant miR169h, that is evidently absent in the wild relatives S. pennellii and S. pimpinellifolium. Seven novel miR169 targets including RNA-binding protein, protein-phosphatase, aminotransferase, chaperone, tetratricopeptide-repeat-protein, and transcription factors ARF-9B and SEPELLATA-3 were established by efficient target cleavage in the presence of specific precursors as well as increased target abundance upon miR169 chelation by short-tandem-target-mimic construct in transient assays. Comparative antagonistic expression profiles of MIR169:target pairs suggest MIR169 family as ubiquitous regulator of various abiotic stresses (heat, cold, dehydration and salt) and developmental pathways. This regulation is partly brought about by acquisition of new promoters as demonstrated by promoter MIR169:GUS reporter assays as well as differential processivity of different precursors and miRNA cleavage efficiencies. Thus, the current study augments the functional horizon of MIR169 family with applications for stress tolerance in crops.

Transcriptome-wide identification and profiling of miRNAs in a stress-tolerant conifer Sabina chinensis

miRNAs are important regulatory components involving in many biological processes, including plant development, vegetative and reproductive growth, and stress response. However, identification and characterization of miRNAs still remain limited for conifer species. In this study, with deep sequencing, we obtained 1,314,450 unique reads with 18-30 nt length from a stress-tolerant conifer, Sabina chinensis. We identified 37 conserved and 103 novel miRNAs, their unique characteristics were further analyzed, and 10 randomly selected were validated by qRT-PCR. Through miRNA target predictions and annotations, we found miRNA may have several targets as well a target could be regulated by several miRNAs, and a total of 2,397 mRNAs were predicted to be targets of the 140 miRNAs. These targets included not only important transcription factors such as auxin response factors, but also indispensable non-transcriptional factor proteins. Pathway-based analysis showed that S. chinensis miRNAs are involved in 172 metabolic pathways, of which 3 were discovered in adaptation-related pathways, indicating their possible relevance to the species' stress-tolerance characteristics. This study is expected to lay the foundation for exploring the regulative roles of miRNAs in development, growth, and response to environmental stresses of S. chinensis.

Determination of the MiRNAs Related to Bean Pyralid Larvae Resistance in Soybean Using Small RNA and Transcriptome Sequencing

Soybean is one of the most important oil crops in the world. Bean pyralid is a major leaf-feeding insect of soybean. In order to screen out the functional genes and regulatory pathways related to the resistance for bean pyralid larvae, the small RNA and transcriptome sequencing were performed based on the highly resistant material (Gantai-2-2) and highly susceptible material (Wan 82-178) of soybean. The results showed that, when comparing 48 h feeding with 0 h feeding, 55 differentially expressed miRNAs were identified in Gantai-2-2 and 58 differentially expressed miRNAs were identified in Wan82-178. When comparing Gantai-2-2 with Wan82-178, 77 differentially expressed miRNAs were identified at 0 h feeding, and 70 differentially expressed miRNAs were identified at 48 h feeding. The pathway analysis of the predicted target genes revealed that the plant hormone signal transduction, RNA transport, protein processing in the endoplasmic reticulum, zeatin biosynthesis, ubiquinone and other terpenoid-quinone biosynthesis, and isoquinoline alkaloid biosynthesis may play important roles in soybean's defense against the stress caused by bean pyralid larvae. According to conjoint analysis of the miRNA/mRNA, a total of 20 differentially expressed miRNAs were negatively correlated with 26 differentially expressed target genes. The qRT-PCR analysis verified that the small RNA sequencing results were credible. According to the analyses of the differentially expressed miRNAs, we speculated that miRNAs are more likely to play key roles in the resistance to insects. Gma-miR156q, Gma-miR166u, Gma-miR166b, Gma-miR166j-3p, Gma-miR319d, Gma-miR394a-3p, Gma-miR396e, and so on-as well as their negatively regulated differentially expressed target genes-may be involved in the regulation of soybean resistance to bean pyralid larvae. These results laid a foundation for further in-depth research regarding the action mechanisms of insect resistance.

Transcriptome-wide identification and characterization of microRNAs responsive to phosphate starvation in Populus tomentosa

miRNAs (microRNAs) are ~ 21-nt non-coding small RNAs (sRNAs) that play crucial regulatory roles in plant biotic and abiotic stress responses. Phosphorus (Pi) deficiency constrains plant growth and reduces yields worldwide. To identify tree miRNAs and evaluate their functions in the response to low Pi, we identified 261 known and 31 candidate novel miRNA families from three sRNA libraries constructed from Populus tomentosa subjected to sufficient or Pi deficiency condition or to restoration of a sufficient Pi level after Pi deficiency. Pi deficiency resulted in significant changes in the abundance of TPM (transcript per million) of 65 known and 3 novel miRNAs. Interestingly, four miRNAs responsive to low N-miR167, miR394, miR171, and miR857-were found to be involved in the response to low Pi. Thirty-five known and one novel miRNAs responded dynamically to Pi fluctuations, suggesting their involvement in the response to Pi deficiency. miRNA clusters comprising 36 miRNAs were identified in 10 chromosomes. Intriguingly, nine pairs of sense and antisense miRNAs transcribed from the same loci were detected in P. tomentosa, which is the first such report in woody plants. Moreover, target genes of the known miRNAs and novel miRNA candidates with significantly changed abundance were predicted, and their functions were annotated. Degradome sequencing supported the identified targets of miRNAs in P. tomentosa. These findings will enhance our understanding of universal and specific molecular regulatory mechanisms of trees under nutrition stress and may facilitate improvement of the Pi utilization efficiency of woody plants.

Integration of metabolome and transcriptome analyses highlights soybean roots responding to phosphorus deficiency by modulating phosphorylated metabolite processes

Phosphorus (P) is a major constituent of biomolecules in plant cells, and is an essential plant macronutrient. Low phosphate (Pi) availability in soils is a major constraint on plant growth. Although a complex variety of plant responses to Pi starvation has been well documented, few studies have integrated both global transcriptome and metabolome analyses to shed light on molecular mechanisms underlying metabolic responses to P deficiency. This study is the first time to investigate global profiles of metabolites and transcripts in soybean (Glycine max) roots subjected to Pi starvation through targeted liquid chromatography electrospray ionization mass spectrometry (LC-ESI-MS/MS) and RNA-sequencing analyses. This integrated analysis allows for assessing coordinated transcriptomic and metabolic responses in terms of both pathway enzyme expression and regulatory levels. Between two Pi availability treatments, a total of 155 metabolites differentially accumulated in soybean roots, of which were phosphorylated metabolites, flavonoids and amino acids. Meanwhile, a total of 1644 differentially expressed genes (DEGs) were identified in soybean roots, including 1199 up-regulated and 445 down-regulated genes. Integration of metabolome and transcriptome analyses revealed Pi starvation responsive connection between specific metabolic processes in soybean roots, especially metabolic processes involving phosphorylated metabolites (e.g., phosphorylated lipids and nucleic acids). Taken together, this study suggests that complex molecular responses scavenging internal Pi from phosphorylated metabolites are typical adaptive strategies soybean roots employ as responses to Pi starvation. Identified DEGs will provide potential target region for future efforts to develop P-efficient soybean cultivars.

Differential Expression of MicroRNAs During Root Formation in Taxus Chinensis Var. mairei Cultivars

MicroRNAs (miRNAs) have been shown to play key roles in the regulation of plant growth and development by modifying the expression of their target genes. However, the influence of miRNAs on root formation and development in woody plants, such as Taxus chinensis, remains largely unknown. In the current study, we explored the phytohormone-response and nutrition-response miRNA expression profiles during T. chinensis rooting by quantitative real-time PCR (qPCR). We identified six phytohormone-response miRNAs, namely, miR164a, miR165, miR167a, miR171b, miR319, and miR391, and eight nutrition-response miRNAs, namely, miR169b, miR395a, miR399c, miR408, miR826, miR827, miR857, and miR2111a, that were differentially expressed at different rooting phases of T. chinensis. Using northern blot analysis of the putative target genes of these miRNAs, we detected the relative gene expression changes of the target genes. Taken together, our results suggest that miRNAs are involved in root formation of T. chinensis and that miRNAs may play important regulatory roles in primary root, crown root, and root hair formation by targeting phytohormone and/or nutrition response genes in T. chinensis. For the first time, these results expand our understanding of the molecular mechanisms of plant root formation and development in a conifer species.

Homeostasis in the soybean miRNA396-GRF network is essential for productive soybean cyst nematode infections

Heterodera glycines, the soybean cyst nematode, penetrates soybean roots and migrates to the vascular cylinder where it forms a feeding site called the syncytium. MiRNA396 (miR396) targets growth-regulating factor (GRF) genes, and the miR396-GRF1/3 module is a master regulator of syncytium development in model cyst nematode H. schachtii infection of Arabidopsis. Here, we investigated whether this regulatory system operates similarly in soybean roots and is likewise important for H. glycines infection. We found that a network involving nine MIR396 and 23 GRF genes is important for normal development of soybean roots and that GRF function is specified in the root apical meristem by miR396. All MIR396 genes are down-regulated in the syncytium during its formation phase while, specifically, 11 different GRF genes are up-regulated. The switch to the syncytium maintenance phase coincides with up-regulation of MIR396 and down-regulation of the 11 GRF genes specifically via post-transcriptional regulation by miR396. Furthermore, interference with the miR396-GRF6/8-13/15-17/19 regulatory network, through either overexpression or knockdown experiments, does not affect the number of H. glycines juveniles that enter the vascular cylinder to initiate syncytia, but specifically inhibits efficient H. glycines development to adult females. Therefore, homeostasis in the miR396-GRF6/8-13/15-17/19 regulatory network is essential for productive H. glycines infections.

Identification of copper (Cu) stress-responsive grapevine microRNAs and their target genes by high-throughput sequencing

MicroRNAs (miRNAs) are a class of single-stranded non-coding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play important roles in plant growth, development and stress responses. With more copper (Cu) and copper containing compounds used as bactericides and fungicides in plants, Cu stress has become one of the serious environmental problems that affect plant growth and development. In order to uncover the hidden response mechanisms of Cu stress, two small RNA libraries were constructed from Cu-treated and water-treated (Control) leaves of 'Summer Black' grapevine. Following high-throughput sequencing and filtering, a total of 158 known and 98 putative novel miRNAs were identified in the two libraries. Among these, 100 known and 47 novel miRNAs were identified as differentially expressed under Cu stress. Subsequently, the expression patterns of nine Cu-responsive miRNAs were validated by quantitative real-time PCR (qRT-PCR). There existed some consistency in expression levels of Cu-responsive miRNAs between high throughput sequencing and qRT-PCR assays. The targets prediction of miRNAs indicates that miRNA may regulate some transcription factors, including AP2, SBP, NAC, MYB and ARF during Cu stress. The target genes for two known and two novel miRNAs showed specific cleavage sites at the 10th and/or 11th nucleotide from the 5'-end of the miRNA corresponding to their miRNA complementary sequences. The findings will lay the foundation for exploring the role of the regulation of miRNAs in response to Cu stress and provide valuable gene information for breeding some Cu-tolerant grapevine cultivars.

miRNA Detection by Stem-Loop RT-qPCR in Studying microRNA Biogenesis and microRNA Responsiveness to Abiotic Stresses

This chapter is devoted to a PCR-based method for analyzing the expression level of mature miRNAs which utilizes the TaqMan^® technology. Stem-loop RT-qPCR requires preparation of separate cDNA templates for each analyzed miRNA as reverse transcription occurs in the presence of a miRNA-specific stem-loop reverse primer. In quantitative analysis, SYBR^® Green is not used but the more sensitive TaqMan^® probe that on 5' end contains a covalently attached fluorophore and on 3' quencher. When quencher and fluorophore are spatially separated due to nucleolytic DNA polymerase activity, the signal is released and quantified. This section provides a detailed and comprehensive protocol allowing for the successful analysis of mature miRNA levels in analyzed sample. Reverse transcription combined with classic real-time PCR as well as ddPCR™ (Droplet Digital™ PCR) will be presented.

Legume nodulation: The host controls the party

Global demand to increase food production and simultaneously reduce synthetic nitrogen fertilizer inputs in agriculture are underpinning the need to intensify the use of legume crops. The symbiotic relationship that legume plants establish with nitrogen-fixing rhizobia bacteria is central to their advantage. This plant-microbe interaction results in newly developed root organs, called nodules, where the rhizobia convert atmospheric nitrogen gas into forms of nitrogen the plant can use. However, the process of developing and maintaining nodules is resource intensive; hence, the plant tightly controls the number of nodules forming. A variety of molecular mechanisms are used to regulate nodule numbers under both favourable and stressful growing conditions, enabling the plant to conserve resources and optimize development in response to a range of circumstances. Using genetic and genomic approaches, many components acting in the regulation of nodulation have now been identified. Discovering and functionally characterizing these components can provide genetic targets and polymorphic markers that aid in the selection of superior legume cultivars and rhizobia strains that benefit agricultural sustainability and food security. This review addresses recent findings in nodulation control, presents detailed models of the molecular mechanisms driving these processes, and identifies gaps in these processes that are not yet fully explained.

Identification and characterisation of microRNAs and their target genes in phosphate-starved Nicotiana benthamiana by small RNA deep sequencing and 5'RACE analysis

BACKGROUND

Phosphorus is an important macronutrient that is severely lacking in soils. In plants, specific microRNAs (miRNAs) essential for nutrient management and the regulation of stress responses are responsible for the control of many phosphate starvation responses. Further understanding of conserved and species-specific microRNA species has potential implications for the development of crops tolerant to soils with low phosphate.

RESULTS

This study identified and characterised phosphate starvation-responsive miRNAs in the native Australian tobacco Nicotiana benthamiana. Small RNA libraries were constructed and sequenced from phosphate-starved plant leaves, stems and roots. Twenty-four conserved miRNA families and 36 species-specific miRNAs were identified. The majority of highly phosphate starvation-responsive miRNAs were highly conserved, comprising of members from the miR399, miR827, and miR2111 families. In addition, two miRNA-star species were identified to be phosphate starvation-responsive. A total of seven miRNA targets were confirmed using RLM-5'RACE to be cleaved by five miRNA families, including two confirmed cleavage targets for Nbe-miR399 species, one for Nbe-miR2111, and two for Nbe-miR398. A number of N. benthamiana-specific features for conserved miRNAs were identified, including species-specific miRNA targets predicted or confirmed for miR399, miR827, and miR398.

CONCLUSIONS

Our results give an insight into the phosphate starvation-responsive miRNAs of Nicotiana benthamiana, and indicate that the phosphate starvation response pathways in N. benthamiana contain both highly conserved and species-specific components.

The PILNCR1-miR399 Regulatory Module Is Important for Low Phosphate Tolerance in Maize

The regulation of adaptive responses to phosphorus (P) deficiency by the microRNA399 (miR399)/PHOSPHATE2 (PHO2) pathway has been well studied in Arabidopsis (Arabidopsis thaliana) but not in maize (Zea mays). Here, we show that miR399 transcripts are strongly induced in maize by phosphate (Pi) deficiency. Transgenic maize plants that overexpressed MIR399b accumulated excessive amounts of P in their shoots and displayed typical Pi-toxicity phenotypes. We reannotated ZmPHO2 with an additional 1,165 bp of the 5' untranslated region. miR399-guided posttranscriptional repression of ZmPHO2 was mainly observed in the P-efficient lines. We identified Pi-deficiency-induced long-noncoding RNA1 (PILNCR1) from our strand-specific RNA libraries. Transient expression assays in Nicotiana benthamiana and maize leaf protoplasts demonstrated that PILNCR1 inhibits ZmmiR399-guided cleavage of ZmPHO2 The abundance of PILNCR1 was significantly higher in P-inefficient lines than in P-efficient lines, which is consistent with the abundance of ZmmiR399 transcripts. These results indicate that the interaction between PILNCR1 and miR399 is important for tolerance to low Pi in maize.

Overexpression of miR169o, an Overlapping MicroRNA in Response to Both Nitrogen Limitation and Bacterial Infection, Promotes Nitrogen Use Efficiency and Susceptibility to Bacterial Blight in Rice

Limiting nitrogen (N) supply contributes to improved resistance to bacterial blight (BB) caused by Xanthomonas oryzae pv. oryzae (Xoo) in susceptible rice (Oryza sativa). To understand the regulatory roles of microRNAs (miRNAs) in this phenomenon, 63 differentially expressed overlapping miRNAs in response to Xoo infection and N limitation stress in rice were identified through deep RNA sequencing and stem-loop quantitative real-time PCR. Among these, miR169o was further assessed as a typical overlapping miRNA through the overexpression of the miR169o primary gene. Osa-miR169o-OX plants were taller, and had more biomass accumulation with significantly increased nitrate and total amino acid contents in roots than the wild type (WT). Transcript level assays showed that under different N supply conditions, miR169o oppositely regulated NRT2, and this is reduced under normal N supply conditions but remarkably induced under N-limiting stress. On the other hand, osa-miR169o-OX plants also displayed increased disease lesion lengths and reduced transcriptional levels of defense gene (PR1b, PR10a, PR10b and PAL) compared with the WT after inoculation with Xoo. In addition, miR169o impeded Xoo-mediated NRT transcription. Therefore, the overlapping miR169o contributes to increase N use efficiency and negatively regulates the resistance to BB in rice. Consistently, transient expression of NF-YA genes in rice protoplasts promoted the transcripts of PR genes and NRT2 genes, while it reduced the transcripts of NRT1 genes. Our results provide novel and additional insights into the co ordinated regulatory mechanisms of cross-talk between Xoo infection and N deficiency responses in rice.

Down-regulation of genes coding for core RNAi components and disease resistance proteins via corresponding microRNAs might be correlated with successful Soybean mosaic virus infection in soybean

Plants protect themselves from virus infections by several different defence mechanisms. RNA interference (RNAi) is one prominent antiviral mechanism, which requires the participation of AGO (Argonaute) and Dicer/DCL (Dicer-like) proteins. Effector-triggered immunity (ETI) is an antiviral mechanism mediated by resistance (R) genes, most of which encode nucleotide-binding site-leucine-rich repeat (NBS-LRR) family proteins. MicroRNAs (miRNAs) play important regulatory roles in plants, including the regulation of host defences. Soybean mosaic virus (SMV) is the most common virus in soybean and, in this work, we identified dozens of SMV-responsive miRNAs by microarray analysis in an SMV-susceptible soybean line. Amongst the up-regulated miRNAs, miR168a, miR403a, miR162b and miR1515a predictively regulate the expression of AGO1, AGO2, DCL1 and DCL2, respectively, and miR1507a, miR1507c and miR482a putatively regulate the expression of several NBS-LRR family disease resistance genes. The regulation of target gene expression by these seven miRNAs was validated by both transient expression assays and RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE) experiments. Transcript levels for AGO1, DCL1, DCL2 and five NBS-LRR family genes were repressed at different time points after SMV infection, whereas the corresponding miRNA levels were up-regulated at these same time points. Furthermore, inhibition of miR1507a, miR1507c, miR482a, miR168a and miR1515a by short tandem target mimic (STTM) technology compromised SMV infection efficiency in soybean. Our results imply that SMV can counteract soybean defence responses by the down-regulation of several RNAi pathway genes and NBS-LRR family resistance genes via the induction of the accumulation of their corresponding miRNA levels.

Functional Characterization of the Steroid Reductase Genes GmDET2a and GmDET2b form Glycine max

Brassinosteroids are important phytohormones for plant growth and development. In soybean (Glycine max), BR receptors have been identified, but the genes encoding BR biosynthesis-related enzymes remain poorly understood. Here, we found that the soybean genome encodes eight steroid reductases (GmDET2a to GmDET2h). Phylogenetic analysis grouped 105 steroid reductases from moss, fern and higher plants into five subgroups and indicated that the steroid reductase family has experienced purifying selection. GmDET2a and GmDET2b, homologs of the Arabidopsis thaliana steroid 5 α -reductase AtDET2, are proteins of 263 amino acids. Ectopic expression of GmDET2a and GmDET2b rescued the defects of the Atdet2-1 mutant in both darkness and light. Compared to the mutant, the hypocotyl length and plant height of the transgenic lines GmDET2a and GmDET2b increased significantly, in both darkness and light, and the transcript levels of the BR biosynthesis-related genes CPD, DWF4, BR6ox-1 and BR6ox-2 were downregulated in GmDET2aOX-23 and GmDET2bOX-16 lines compared to that in Atdet2-1. Quantitative real-time PCR revealed that GmDET2a and GmDET2b are ubiquitously expressed in all tested soybean organs, including roots, leaves and hypocotyls. Moreover, epibrassinosteroid negatively regulated GmDET2a and GmDET2b expression. Sulfate deficiency downregulated GmDET2a in leaves and GmDET2b in leaves and roots; by contrast, phosphate deficiency upregulated GmDET2b in roots and leaves. Taken together, our results revealed that GmDET2a and GmDET2b function as steroid reductases.