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

From molecular basics to agronomic benefits: Insights into noncoding RNA-mediated gene regulation in plants

The development of plants is largely dependent on their growth environment. To better adapt to a particular habitat, plants have evolved various subtle regulatory mechanisms for altering gene expression. Non coding RNAs (ncRNAs) constitute a major portion of the transcriptomes of eukaryotes. Various ncRNAs have been recognized as important regulators of the expression of genes involved in essential biological processes throughout the whole life cycles of plants. In this review, we summarize the current understanding of the biogenesis and contributions of small nucle olar RNA (snoRNA)- and regulatory long non coding RNA (lncRNA)-mediated gene regulation in plant development and environmental responses. Many regulatory ncRNAs appear to be associated with increased yield, quality and disease resistance of various species and cultivars. These ncRNAs may potentially be used as genetic resources for improving agronomic traits and for molecular breeding. The challenges in understanding plant ncRNA biology and the possibilities to make better use of these valuable gene resources in the future are discussed in this review.

A R2R3 MYB Transcription Factor, TaMYB391, Is Positively Involved in Wheat Resistance to Puccinia striiformis f. sp. tritici

A biotrophic fungus, Puccinia striiformis f.sp. tritici (Pst), which causes stripe rust disease in wheat is the most yield-limiting factor in wheat production. Plants have complex defense mechanisms against invading pathogens. Hypersensitive response (HR), a kind of programmed cell death (PCD) at the infection site, is among these defense mechanisms. Transcription factors (TFs) play a crucial role in plant defense response against invading pathogens. Myeloblastosis (MYB) TFs are among the largest TFs families that are involved in response to both biotic and abiotic stresses. However, little is known about the mechanisms of MYB TFs during the interaction between wheat and the stripe rust fungus. Here, we identified an R2R3 MYB TF from wheat, designated as TaMYB391, and characterized its functional role during wheat-Pst interaction. Our data indicated that TaMYB391 is induced by Pst infection and exogenous application of salicylic acid (SA) and abscisic acid (ABA). TaMYB391 is localized in the nucleus of both wheat and Nicotiana benthamiana. Transient overexpression of TaMYB391 in N. benthamiana triggered HR-related PCD accompanied by increased electrolyte leakage, high accumulation of reactive oxygen species (ROS), and transcriptional accumulation of SA defense-related genes and HR-specific marker genes. Overexpression of TaMYB391 in wheat significantly enhanced wheat resistance to stripe rust fungus through the induction of pathogenesis-related (PR) genes, ROS accumulation and hypersensitive cell death. On the other hand, RNAi-mediated silencing of TaMYB391 decreased the resistance of wheat to Pst accompanied by enhanced growth of the pathogen. Together our findings demonstrate that TaMYB391 acts as a positive regulator of HR-associated cell death and positively contributes to the resistance of wheat to the stripe rust fungus by regulating certain PR genes, possibly through SA signaling pathways.

Prospects of genetics and breeding for low-phosphate tolerance: an integrated approach from soil to cell

Improving phosphorus (P) crop nutrition has emerged as a key factor toward achieving a more resilient and sustainable agriculture. P is an essential nutrient for plant development and reproduction, and phosphate (Pi)-based fertilizers represent one of the pillars that sustain food production systems. To meet the global food demand, the challenge for modern agriculture is to increase food production and improve food quality in a sustainable way by significantly optimizing Pi fertilizer use efficiency. The development of genetically improved crops with higher Pi uptake and Pi-use efficiency and higher adaptability to environments with low-Pi availability will play a crucial role toward this end. In this review, we summarize the current understanding of Pi nutrition and the regulation of Pi-starvation responses in plants, and provide new perspectives on how to harness the ample repertoire of genetic mechanisms behind these adaptive responses for crop improvement. We discuss on the potential of implementing more integrative, versatile, and effective strategies by incorporating systems biology approaches and tools such as genome editing and synthetic biology. These strategies will be invaluable for producing high-yielding crops that require reduced Pi fertilizer inputs and to develop a more sustainable global agriculture.

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.

Parallel analysis of RNA ends reveals global microRNA-mediated target RNA cleavage in maize

MicroRNAs (miRNAs) are endogenous 20-24-nucleotide non-coding RNAs that play important regulatory roles in many biological processes in eukaryotes. miRNAs modulate the expression of target genes at the post-transcriptional level by transcript cleavage or translational inhibition. The identification of miRNA target genes has been extensively investigated in Arabidopsis and rice, but an in-depth global analysis of miRNA-mediated target regulation is still lacking in maize. Here, we report a transcriptome-wide identification of miRNA targets by analyzing parallel analysis of RNA ends (PARE) datasets derived from nine different tissues at five developmental stages of the maize (Zea mays L.) B73 cultivar. In total, 246 targets corresponding to 60 miRNAs from 25 families were identified, including transcription factors and other genes. In addition, PARE analysis revealed that miRNAs guide specific target transcript cleavage in a tissue-preferential manner. Primary transcripts of MIR159c and MIR169e were found to be cleaved by mature miR159 and miR169, respectively, indicating a negative-feedback regulatory mechanism in miRNA biogenesis. Moreover, several miRNA-target gene pairs involved in seed germination were identified and experimentally validated. Our PARE analyses generated a wide and detailed miRNA-target interaction atlas, which provides a valuable resource for investigating the roles of miRNAs and their targets in maize.

Transcriptome Level Analysis of Genes of Exogenous Ethylene Applied under Phosphorus Stress in Chinese Fir

Chinese fir (Cunninghamia lanceolata (Lamb.) Hook) is a widely grown gymnosperm in China. Phosphorus (P) is an indispensable nutrient for the growth of Chinese fir. Inorganic phosphate (Pi) deficiency exists in soils of many Chinese fir planting area regions, and the trees themselves have limited efficiency in utilizing P from the soil. Ethylene is important in regulation responses to nutrient deficiencies. However, little is known about how ethylene signals participate in Pi stress in Chinese fir. A total of six different treatments were performed to reveal the transcript levels of Chinese fir under Pi, ethephon (an ethylene-releasing compound), and CoCl₂ (cobalt chloride, an ethylene biosynthesis inhibitor) treatments. We assembled a full-length reference transcriptome containing 22,243 unigenes as a reference for UMI RNA-seq (Digital RNA-seq). There were 586 Differentially Expressed Genes (DEGs) in the Pi starvation (NP) group, while DEGs from additional ethephon or CoCl₂ in NP were 708 and 292, respectively. Among the DEGs in each treatment, there were 83 TFs in these treatment groups. MYB (v-myb avian myeloblastosis viral oncogene homolog) family was the most abundant transcription factors (TFs). Three ERF (Ethylene response factor) family genes were identified when only ethylene content was imposed as a variable. Enrichment analysis indicated that the ascorbate and aldarate metabolism pathway plays a key role in resistance to Pi deficiency. This study provides insights for further elucidating the regulatory mechanism of Pi deficiency in Chinese fir.

Gypsy retrotransposon-derived maize lncRNA GARR2 modulates gibberellin response

Long non-coding RNAs (lncRNAs) mediate diverse biological events mainly through the modulation of transcriptional hierarchy. The phytohormone gibberellin (GA) is essential for various aspects of plant growth and development. However, the roles of lncRNAs in the regulation of the GA response remain largely unknown. Through sequencing multiple strand-specific and ribosomal-depleted RNA libraries, we delineated the landscape of lncRNAs in maize (Zea mays). Out of identified lncRNAs, 445 GIBBERELLIN-RESPONSIVE lncRNAs (GARRs) were differentially expressed upon GA application. By the intersection of GARRs from normal-height and dwarf plants from an advanced backcross population, four shared GARRs (GARR1 to GARR4) were identified. Out of these four shared GARRs, GARR2 was derived from a Gypsy LTR retrotransposon. GA-responsive element P-boxes were identified upstream of GARR2. GARR2-edited lines exhibited a GA-induced phenotype. Editing of GARR2 resulted in changes in the transcriptional abundance of GA pathway components and endogenous GA contents. Besides GA, GARR2 affected the primary auxin response. An RNA pull-down assay revealed the HECT ubiquitin-protein ligase family member ZmUPL1 as a potential interaction target of GARR2. GARR2 influenced the abundance of ZmUPL1 in the GA response. Our study uncovers lncRNA players involved in the modulation of the GA response and guides the development of plant height ideotype driven by knowledge of the phytohormone GA.

Analysis of miRNAs responsive to long-term calcium deficiency in tef (Eragrostis tef (Zucc.) Trotter)

MicroRNAs (miRNAs) play an important role in growth, development, stress resilience, and epigenetic modifications of plants. However, the effect of calcium (Ca²⁺) deficiency on miRNA expression in the orphan crop tef (Eragrostis tef) remains unknown. In this study, we analyzed expression of miRNAs in roots and shoots of tef in response to Ca²⁺ treatment. miRNA-seq followed by bioinformatic analysis allowed us to identify a large number of small RNAs (sRNAs) ranging from 17 to 35 nt in length. A total of 1380 miRNAs were identified in tef experiencing long-term Ca²⁺ deficiency while 1495 miRNAs were detected in control plants. Among the miRNAs identified in this study, 161 miRNAs were similar with those previously characterized in other plant species and 348 miRNAs were novel, while the remaining miRNAs were uncharacterized. Putative target genes and their functions were predicted for all the known and novel miRNAs that we identified. Based on gene ontology (GO) analysis, the predicted target genes are known to have various biological and molecular functions including calcium uptake and transport. Pairwise comparison of differentially expressed miRNAs revealed that some miRNAs were specifically enriched in roots or shoots of low Ca²⁺-treated plants. Further characterization of the miRNAs and their targets identified in this study may help in understanding Ca²⁺ deficiency responses in tef and related orphan crops.

The Characters of Non-Coding RNAs and Their Biological Roles in Plant Development and Abiotic Stress Response

Plant growth and development are greatly affected by the environment. Many genes have been identified to be involved in regulating plant development and adaption of abiotic stress. Apart from protein-coding genes, more and more evidence indicates that non-coding RNAs (ncRNAs), including small RNAs and long ncRNAs (lncRNAs), can target plant developmental and stress-responsive mRNAs, regulatory genes, DNA regulatory regions, and proteins to regulate the transcription of various genes at the transcriptional, posttranscriptional, and epigenetic level. Currently, the molecular regulatory mechanisms of sRNAs and lncRNAs controlling plant development and abiotic response are being deeply explored. In this review, we summarize the recent research progress of small RNAs and lncRNAs in plants, focusing on the signal factors, expression characters, targets functions, and interplay network of ncRNAs and their targets in plant development and abiotic stress responses. The complex molecular regulatory pathways among small RNAs, lncRNAs, and targets in plants are also discussed. Understanding molecular mechanisms and functional implications of ncRNAs in various abiotic stress responses and development will benefit us in regard to the use of ncRNAs as potential character-determining factors in molecular plant breeding.

Significance of miRNA in enhancement of flavonoid biosynthesis

Flavonoid metabolism shows very strong plasticity in plant development and coping with the changing environment. Flavonoid biosynthesis is regulated by many metabolic pathways, including transcriptional regulation, post-transcriptional control, post-translational regulationand epigenetic regulation. miRNA is a form of endogenous noncoding single-strand small molecule RNA that primarily regulates the expression of target genes horizontally after transcription through splicing and translational suppression. It also plays an important role in regulating plant growth and development, secondary metabolism and biotic and abiotic stress. miRNA can regulate the formation of flavonoids by acting on structural genes or indirectly by using an MBW transcription complex comprising MYB-bHLH-WD40. This study summarizes the biosynthesis and mechanisms of miRNA, and provides a summary of the mechanisms of miRNAs involved in production of flavonoids, in order to elucidate the biosynthesis pathway and complex regulatory network of plant flavonoids. We aim to provide new insights into improving the content of flavonoid active ingredients in plants.

GreeNC 2.0: a comprehensive database of plant long non-coding RNAs

The Green Non-Coding Database (GreeNC) is one of the reference databases for the study of plant long non-coding RNAs (lncRNAs). Here we present our most recent update where 16 species have been updated, while 78 species have been added, resulting in the annotation of more than 495 000 lncRNAs. Moreover, sequence clustering was applied providing information about sequence conservation and gene families. The current version of the database is available at: http://greenc.sequentiabiotech.com/wiki2/Main_Page.

Integrated Analysis of miRNAs Associated With Sugarcane Responses to Low-Potassium Stress

Sugarcane is among the most important global crops and a key bioenergy source. Sugarcane production is restricted by limited levels of available soil potassium (K⁺). The ability of plants to respond to stressors can be regulated by a range of microRNAs (miRNAs). However, there have been few studies regarding the roles of miRNAs in the regulation of sugarcane responses to K⁺-deficiency. To understand how these non-coding RNAs may influence sugarcane responses to low-K⁺ stress, we conducted expression profiling of miRNAs in sugarcane roots under low-K⁺ conditions via high-throughput sequencing. This approach led to the identification of 324 and 42 known and novel miRNAs, respectively, of which 36 were found to be differentially expressed miRNAs (DEMs) under low-K⁺ conditions. These results also suggested that miR156-x/z and miR171-x are involved in these responses as potential regulators of lateral root formation and the ethylene signaling pathway, respectively. A total of 705 putative targets of these DEMs were further identified through bioinformatics predictions and degradome analyses, and GO and KEGG enrichment analyses revealed these target mRNAs to be enriched for catalytic activity, binding functions, metabolic processes, plant hormone signal transduction, and mitogen-activated protein kinase (MAPK) signaling. In summary, these data provide an overview of the roles of miRNAs in the regulation of sugarcane response to low-K⁺ conditions.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

tae-miR399-UBC24 Module Enhances Freezing Tolerance in Winter Wheat via a CBF Signaling Pathway

Although the regulation of Pi homeostasis by miR399 has been studied in various plant species, its underlying molecular mechanism in response to freezing stress is still poorly understood. In this work, we found that the expression of tae-miR399 and its target gene TaUBC24 in the tillering nodes of the strong cold-resistant winter wheat cultivar Dongnongdongmai1 (Dn1) was not only significantly altered after severe winters but also responsive to short-term freezing stress. TaUBC24 physically interacted with TaICE1. Enhanced freezing tolerance was observed for tae-miR399-overexpressing Arabidopsis lines. Under freezing stress, overexpression of tae-miR399 ultimately decreased the expression of AtUBC24, inhibiting the degradation of AtICE1, which increased the expression of genes involved in the CBF signaling pathway and starch metabolism and promoted the activities of antioxidant enzymes. These results will improve our understanding of the molecular mechanism through which the miR399-UBC24 module plays a cardinal role in regulating plant freezing stress tolerance through mediation of downstream pathways.

Identification of SPX family genes in the maize genome and their expression under different phosphate regimes

Many studies have revealed that SPX (SYG1/Pho81/XPR1) family genes play a key role in signal transduction related to phosphorus (P) deficiency in plants. Here, we identified 33 SPX gene family members in maize through genome-wide analysis and classified them into 4 subfamilies according to SPX structural characteristics (SPX, SPX-MFS, SPX-EXS and SPX-RING). The promoter regions of ZmSPXs are rich in biotic/abiotic-related stress elements. The quantitative real-time PCR analysis of 33 ZmSPXs revealed that all members except for ZmSPX3 of the SPX subfamily were significantly induced under P-deficient conditions, especially ZmSPX4.1 and ZmSPX4.2, which showed strong responses to low P stress and exhibited remarkably different expression patterns in low Pi sensitive and insensitive cultivars of maize. These results suggested that the SPX subfamily might play pivotal roles in P stress sensing and response. Experimental observations of subcellular localization in maize protoplasts indicated the following results, implying multiple roles in cell metabolism: ZmSPX2, ZmSPX5 and ZmSPX6 localized in the nucleus; ZmSPX1 and ZmSPX3 localized in the nucleus and cytoplasm; and ZmSPX4.2 localized in the chloroplast. A Y2H assay suggested that ZmPHR1 could interact with ZmSPX3, ZmSPX4.2, ZmSPX5, and ZmSPX6, indicating the involvement of these proteins in the P stress response in a ZmPHR1-mediated manner.

Mining plant endogenous target mimics from miRNA-lncRNA interactions based on dual-path parallel ensemble pruning method

The interactions between microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) play important roles in biological activities. Specially, lncRNAs as endogenous target mimics (eTMs) can bind miRNAs to regulate the expressions of target messenger RNAs (mRNAs). A growing number of studies focus on animals, but the studies on plants are scarce and many functions of plant eTMs are unknown. This study proposes a novel ensemble pruning protocol for predicting plant miRNA-lncRNA interactions at first. It adaptively prunes the base models based on dual-path parallel ensemble method to meet the challenge of cross-species prediction. Then potential eTMs are mined from predicted results. The expression levels of RNAs are identified through biological experiment to construct the lncRNA-miRNA-mRNA regulatory network, and the functions of potential eTMs are inferred through enrichment analysis. Experiment results show that the proposed protocol outperforms existing methods and state-of-the-art predictors on various plant species. A total of 17 potential eTMs are verified by biological experiment to involve in 22 regulations, and 14 potential eTMs are inferred by Gene Ontology enrichment analysis to involve in 63 functions, which is significant for further research.

Long noncoding RNA lncRNA354 functions as a competing endogenous RNA of miR160b to regulate ARF genes in response to salt stress in upland cotton

Long non-coding RNAs (lncRNAs) play important roles in response to biotic and abiotic stress through acting as competing endogenous RNAs (ceRNAs) to decoy mature miRNAs. However, whether this mechanism is involved in cotton salt stress response remains unknown. We report the characterization of an endogenous lncRNA, lncRNA354, whose expression was reduced in salt-treated cotton and was localized at the nucleus and cytoplasm. Using endogenous target mimic (eTM) analysis, we predicted that lncRNA354 had a potential binding site for miR160b. Transient expression in tobacco demonstrated that lncRNA354 was a miR160b eTM and attenuated miR160b suppression of its target genes, including auxin response factors (ARFs). Silencing or overexpressing lncRNA354 affected the expression of miR160b and target ARFs. Silencing lncRNA354 and targets GhARF17/18 resulted in taller cotton plants and enhanced the resistant to salt stress. Overexpression of lncRNA354 and targets GhARF17/18 in Arabidopsis led to dwarf plants, decreased root dry weight and reduced salt tolerance. Our results show that the lncRNA354-miR160b effect on GhARF17/18 expression may modulate auxin signalling and thus affect growth. These results also shed new light on a mechanism of lncRNA-associated responses to salt stress.

Identification and characterization of circRNAs in maize seedlings under deficient nitrogen

Here, deep sequencing results of the maize transcriptome in leaves and roots were compared under high-nitrogen (HN) and low-nitrogen (LN) conditions to identify differentially expressed circRNAs (DECs). Circular RNAs (circRNAs) are covalently closed non-coding RNA with widely regulatory potency that has been identified in animals and plants. However, the understanding of circRNAs involved in responsive nitrogen deficiency remains to be elucidated. A total of 24 and 22 DECs were obtained from the leaves and roots, respectively. Ten circRNAs were validated by divergent and convergent primers, and 6 DECs showed the same expression tendency validated by reverse transcriptase-quantitative PCR. Integrating the identified differentially expressed miRNAs, 34 circRNAs could act as miRNA decoys, which might play important roles in multiple biological processes, including organonitrogen compound biosynthesis and regulation of the metabolic process. A total of 51 circRNA-parent genes located in the genome-wide association study identified loci were assessed between HN and LN conditions and were associated with root growth and development. In summary, our results provide valuable information regarding further study of maize circRNAs under nitrogen deficiency and provide new insights into screening of candidate genes as well as the improvement of maize regarding nitrogen deficiency resistance. CircRNA-miRNA-mRNA co-expression networks were constructed to explore the circRNAs that participated in biological development and nitrogen metabolism.

Genetic Dissection of Phosphorus Use Efficiency in a Maize Association Population under Two P Levels in the Field

Phosphorus (P) deficiency is an important challenge the world faces while having to increase crop yields. It is therefore necessary to select maize (Zea may L.) genotypes with high phosphorus use efficiency (PUE). Here, we extensively analyzed the biomass, grain yield, and PUE-related traits of 359 maize inbred lines grown under both low-P and normal-P conditions. A significant decrease in grain yield per plant and biomass, an increase in PUE under low-P condition, as well as significant correlations between the two treatments were observed. In a genome-wide association study, 49, 53, and 48 candidate genes were identified for eleven traits under low-P, normal-P conditions, and in low-P tolerance index (phenotype under low-P divided by phenotype under normal-P condition) datasets, respectively. Several gene ontology pathways were enriched for the genes identified under low-P condition. In addition, seven key genes related to phosphate transporter or stress response were molecularly characterized. Further analyses uncovered the favorable haplotype for several core genes, which is less prevalent in modern lines but often enriched in a specific subpopulation. Collectively, our research provides progress in the genetic dissection and molecular characterization of PUE in maize.

Integrated analysis of miRNA-mRNA regulatory networks of potato (Solanum tuberosum L.) in response to cadmium stress

Cadmium (Cd) stress is a ubiquitous abiotic stress affecting plant growth worldwide and negatively impacting crop yield and food safety. Potato is the most important non-grain crop globally, but there is limited research available on the response of this crop to Cd stress. This study explored the coping mechanism for Cd stress in potato through analyses of miRNA and mRNA. Tissue culture seedlings (20-day-old) of potato variety 'Atlantic' were cultured for up to 48 h in liquid medium containing 5 mmol/L CdCl2, and phenotypic, physiological, and transcriptomic changes were observed at specific times. With the extension of Cd stress time, the potato leaves gradually wilted and curled, and root salicylic acid (SA), glutathione (GSH), and lignin contents and peroxidase (POD) activity increased, while indole-3-acetic acid (IAA) and zeatin (ZT) contents decreased. Using miRNA-seq, 161 existing miRNAs, 383 known miRNAs, and 7361 novel miRNAs were identified, and, 18 miRNAs were differentially expressed in response to Cd stress. Based on mRNA-seq, 7340 differentially expressed mRNAs (DEGs) were found. Through mRNA-miRNA integrated analysis, miRNA-target gene pairs consisting of 23 DEGs and 33 miRNAs were identified. Furthermore, "glutathione metabolism" "plant hormone signal transduction" and "phenylpropanoid biosynthesis" were established as crucial pathways in the Cd stress response of potato. Novel miRNAs novel-m3483-5p and novel-m2893-5p participate in these pathways through targeted regulation of cinnamic alcohol dehydrogenase (CAD; PG0005359) and alanine aminotransferase (POP; PG0024281), respectively. This study provides information that will help elucidate the complex mechanism of the Cd stress response in potato. Moreover, candidate miRNAs and mRNAs could yield new strategies for the development of Cd-tolerant potato breeding.

Integrated transcriptome and miRNA sequencing approaches provide insights into salt tolerance in allotriploid Populus cathayana

Some salt-stress responsive DEGs, mainly involved in ion transmembrane transport, hormone regulation, antioxidant system, osmotic regulation, and some miRNA jointly regulated the salt response process in allotriploid Populus cathayana. The molecular mechanism of plant polyploid stress resistance has been a hot topic in biological research. In this study, Populus diploids and first division restitution (FDR) and second division restitution (SDR) triploids were selected as research materials. All materials were treated with 70 mM NaCl solutions for 30 days in the same pot environment. We observed the growth state of triploids and diploids and determined the ratio of potassium and sodium ions, peroxidase (POD) activity, proline content, and ABA and jasmonic acid (JA) hormone content in leaves in the same culture environment with the same concentration of NaCl solution treatment. In addition, RNA-seq technology was used to study the differential expression of mRNA and miRNA. The results showed that triploid Populus grew well and the K⁺ content and the K⁺/Na⁺ ratio in the salt treatment were significantly lower than those in the control. The contents of ABA, JA, POD, and proline were increased compared with contents in diploid under salt stress. The salt-stress responsive DEGs were mainly involved in ion transport, cell homeostasis, the MAPK signaling pathway, peroxisome, citric acid cycle, and other salt response and growth pathways. The transcription factors mainly included NAC, MYB, MYB_related and AP2/ERF. Moreover, the differentially expressed miRNAs involved 32 families, including 743 miRNAs related to predicted target genes, among which 22 miRNAs were significantly correlated with salt-stress response genes and related to the regulation of hormones, ion transport, reactive oxygen species (ROS) and other biological processes. Our results provided insights into the physiological and molecular aspects for further research into the response mechanisms of allotriploid Populus cathayana to salt stress. This study provided valuable information for the salt tolerance mechanism of allopolyploids.

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.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting soybean genotypes subjected to phosphate starvation

BACKGROUND

Phosphorus (P) is essential for plant growth and development, and low-phosphorus (LP) stress is a major factor limiting the growth and yield of soybean. Long noncoding RNAs (lncRNAs) have recently been reported to be key regulators in the responses of plants to stress conditions, but the mechanism through which LP stress mediates the biogenesis of lncRNAs in soybean remains unclear.

RESULTS

In this study, to explore the response mechanisms of lncRNAs to LP stress, we used the roots of two representative soybean genotypes that present opposite responses to P deficiency, namely, a P-sensitive genotype (Bogao) and a P-tolerant genotype (NN94156), for the construction of RNA sequencing (RNA-seq) libraries. In total, 4,166 novel lncRNAs, including 525 differentially expressed (DE) lncRNAs, were identified from the two genotypes at different P levels. GO and KEGG analyses indicated that numerous DE lncRNAs might be involved in diverse biological processes related to phosphate, such as lipid metabolic processes, catalytic activity, cell membrane formation, signal transduction, and nitrogen fixation. Moreover, lncRNA-mRNA-miRNA and lncRNA-mRNA networks were constructed, and the results identified several promising lncRNAs that might be highly valuable for further analysis of the mechanism underlying the response of soybean to LP stress.

CONCLUSIONS

These results revealed that LP stress can significantly alter the genome-wide profiles of lncRNAs, particularly those of the P-sensitive genotype Bogao. Our findings increase the understanding of and provide new insights into the function of lncRNAs in the responses of soybean to P stress.

From "Dark Matter" to "Star": Insight Into the Regulation Mechanisms of Plant Functional Long Non-Coding RNAs

Long non-coding RNAs (lncRNAs) play a vital role in a variety of biological functions in plant growth and development. In this study, we provided an overview of the molecular mechanisms of lncRNAs in interacting with other biomolecules with an emphasis on those lncRNAs validated only by low-throughput experiments. LncRNAs function through playing multiple roles, including sponger for sequestering RNA or DNA, guider or decoy for recruiting or hijacking transcription factors or peptides, and scaffold for binding with chromatin modification complexes, as well as precursor of microRNAs or small interfering RNAs. These regulatory roles have been validated in several plant species with a comprehensive list of 73 lncRNA-molecule interaction pairs in 16 plant species found so far, suggesting their commonality in the plant kingdom. Such initial findings of a small number of functional plant lncRNAs represent the beginning of what is to come as lncRNAs with unknown functions were found in orders of magnitude more than proteins.

Low nitrogen availability inhibits the phosphorus starvation response in maize (Zea mays ssp. mays L.)

BACKGROUND

Nitrogen (N) and phosphorus (P) are macronutrients essential for crop growth and productivity. In cultivated fields, N and P levels are rarely sufficient, contributing to the gap between realized and potential production. Fertilizer application increases nutrient availability, but is not available to all farmers, nor are current rates of application sustainable or environmentally desirable. Transcriptomic studies of cereal crops have revealed dramatic responses to either low N or low P single stress treatments. In the field, however, levels of both N and P may be suboptimal. The interaction between N and P starvation responses remains to be fully characterized.

RESULTS

We characterized growth and root and leaf transcriptomes of young maize plants under nutrient replete, low N, low P or combined low NP conditions. We identified 1555 genes to respond to our nutrient treatments, in one or both tissues. A large group of genes, including many classical P starvation response genes, were regulated antagonistically between low N and P conditions. An additional experiment over a range of N availability indicated that a mild reduction in N levels was sufficient to repress the low P induction of P starvation genes. Although expression of P transporter genes was repressed under low N or low NP, we confirmed earlier reports of P hyper accumulation under N limitation.

CONCLUSIONS

Transcriptional responses to low N or P were distinct, with few genes responding in a similar way to the two single stress treatments. In combined NP stress, the low N response dominated, and the P starvation response was largely suppressed. A mild reduction in N availability was sufficient to repress the induction of P starvation associated genes. We conclude that activation of the transcriptional response to P starvation in maize is contingent on N availability.

Interplay between miRNAs and lncRNAs: Mode of action and biological roles in plant development and stress adaptation

Plants employ sophisticated mechanisms to control developmental processes and to cope with environmental changes at transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) and long noncoding RNAs (lncRNAs), two classes of endogenous noncoding RNAs, are key regulators of gene expression in plants. Recent studies have identified the interplay between miRNAs and lncRNAs as a novel regulatory layer of gene expression in plants. On one hand, miRNAs target lncRNAs for the production of phased small interfering RNAs (phasiRNAs). On the other hand, lncRNAs serve as origin of miRNAs or regulate the accumulation or activity of miRNAs at transcription and post-transcriptional levels. Theses lncRNA-miRNA interplays are crucial for plant development, physiology and responses to biotic and abiotic stresses. In this review, we summarize recent advances in the biological roles, interaction mechanisms and computational predication methods of the interplay between miRNAs and lncRNAs in plants.

Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the MiR156, MiR5488 and MiR399 are Involved in the Regulation of Male Sterility in PTGMS Rice

A photoperiod- and thermo-sensitive genic male sterile (PTGMS) line is the basic material for two-hybrid rice and is an important genetic breeding resource. Peiai64S (PA64S) is an important germplasm resource of PTGMS rice, and it has been applied to two-line hybrid rice systems in China. Pollen fertility in PA64S is regulated by the temperature and photoperiod, but the mechanism of the fertility transition is unclear. In this study, we obtained the male fertile plant PA64S(F) and the male sterile plant PA64S(S) by controlling different temperatures under long light conditions and used the male fertile and sterile plants to investigate the role of microRNAs (miRNAs) in regulating male fertility in rice. We performed the small RNA library sequencing of anthers from PA64S(S) and PA64S(F). A total of 196 miRNAs were identified-166 known miRNAs among 27 miRNA families and 30 novel miRNAs. In the transcriptome analysis, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially expressed genes revealed significant enrichment in the synthesis and metabolism of fatty acids and some secondary metabolism pathways such as fatty acid metabolism and phenylalanine metabolism. With a comprehensive analysis of miRNA, transcriptome, and degradome sequencing, we identified that 13 pairs of miRNA/target genes regulated male fertility in rice by responding to temperature change, among which the miR156, miR5488, and miR399 affect the male fertility of PA64S by influencing SPLs, the lignin synthesis of anther walls, and the flavonoid metabolism pathway. The results provide a new understanding of PTGMS rice, which will help us better understand the potential regulatory mechanisms of male sterility in the future.

Plant Long Noncoding RNAs: New Players in the Field of Post-Transcriptional Regulations

The first reference to the "C-value paradox" reported an apparent imbalance between organismal genome size and morphological complexity. Since then, next-generation sequencing has revolutionized genomic research and revealed that eukaryotic transcriptomes contain a large fraction of non-protein-coding components. Eukaryotic genomes are pervasively transcribed and noncoding regions give rise to a plethora of noncoding RNAs with undeniable biological functions. Among them, long noncoding RNAs (lncRNAs) seem to represent a new layer of gene expression regulation, participating in a wide range of molecular mechanisms at the transcriptional and post-transcriptional levels. In addition to their role in epigenetic regulation, plant lncRNAs have been associated with the degradation of complementary RNAs, the regulation of alternative splicing, protein sub-cellular localization, the promotion of translation and protein post-translational modifications. In this review, we report and integrate numerous and complex mechanisms through which long noncoding transcripts regulate post-transcriptional gene expression in plants.

Transcriptomic and genome-wide association study reveal long noncoding RNAs responding to nitrogen deficiency in maize

BACKGROUND

Long noncoding RNAs (lncRNAs) play important roles in essential biological processes. However, our understanding of lncRNAs as competing endogenous RNAs (ceRNAs) and their responses to nitrogen stress is still limited.

RESULTS

Here, we surveyed the lncRNAs and miRNAs in maize inbred line P178 leaves and roots at the seedling stage under high-nitrogen (HN) and low-nitrogen (LN) conditions using lncRNA-Seq and small RNA-Seq. A total of 894 differentially expressed lncRNAs and 38 different miRNAs were identified. Co-expression analysis found that two lncRNAs and four lncRNA-targets could competitively combine with ZmmiR159 and ZmmiR164, respectively. To dissect the genetic regulatory by which lncRNAs might enable adaptation to limited nitrogen availability, an association mapping panel containing a high-density single-nucleotide polymorphism (SNP) array (56,110 SNPs) combined with variable LN tolerant-related phenotypes obtained from hydroponics was used for a genome-wide association study (GWAS). By combining GWAS and RNA-Seq, 170 differently expressed lncRNAs within the range of significant markers were screened. Moreover, 40 consistently LN-responsive genes including those involved in glutamine biosynthesis and nitrogen acquisition in root were identified. Transient expression assays in Nicotiana benthamiana demonstrated that LNC_002923 could inhabit ZmmiR159-guided cleavage of Zm00001d015521.

CONCLUSIONS

These lncRNAs containing trait-associated significant SNPs could consider to be related to root development and nutrient utilization. Taken together, the results of our study can provide new insights into the potential regulatory roles of lncRNAs in response to LN stress, and give valuable information for further screening of candidates as well as the improvement of maize resistance to LN stress.

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

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.

Identification and Functional Characterization of Plant MiRNA Under Salt Stress Shed Light on Salinity Resistance Improvement Through MiRNA Manipulation in Crops

Salinity, as a major environmental stressor, limits plant growth, development, and crop yield remarkably. However, plants evolve their own defense systems in response to salt stress. Recently, microRNA (miRNA) has been broadly studied and considered to be an important regulator of the plant salt-stress response at the post-transcription level. In this review, we have summarized the recent research progress on the identification, functional characterization, and regulatory mechanism of miRNA involved in salt stress, have discussed the emerging manipulation of miRNA to improve crop salt resistance, and have provided future direction for plant miRNA study under salt stress, suggesting that the salinity resistance of crops could be improved by the manipulation of microRNA.

Association mapping and genomic selection for sorghum adaptation to tropical soils of Brazil in a sorghum multiparental random mating population

A multiparental random mating population used in sorghum breeding is amenable for the detection of QTLs related to tropical soil adaptation, fine mapping of underlying genes and genomic selection approaches. Tropical soils where low phosphorus (P) and aluminum (Al) toxicity limit sorghum [Sorghum bicolor (L.) Moench] production are widespread in the developing world. We report on BRP13R, a multiparental random mating population (MP-RMP), which is commonly used in sorghum recurrent selection targeting tropical soil adaptation. Recombination dissipated much of BRP13R's likely original population structure and average linkage disequilibrium (LD) persisted up to 2.5 Mb, establishing BRP13R as a middle ground between biparental populations and sorghum association panels. Genome-wide association mapping (GWAS) identified conserved QTL from previous studies, such as for root morphology and grain yield under low-P, and indicated the importance of dominance in the genetic architecture of grain yield. By overlapping consensus QTL regions, we mapped two candidate P efficiency genes to a ~ 5 Mb region on chromosomes 6 (ALMT) and 9 (PHO2). Remarkably, we find that only 200 progeny genotyped with ~ 45,000 markers in BRP13R can lead to GWAS-based positional cloning of naturally rare, subpopulation-specific alleles, such as for SbMATE-conditioned Al tolerance. Genomic selection was found to be useful in such MP-RMP, particularly if markers in LD with major genes are fitted as fixed effects into GBLUP models accommodating dominance. Shifts in allele frequencies in progeny contrasting for grain yield indicated that intermediate to minor-effect genes on P efficiency, such as SbPSTOL1 genes, can be employed in pre-breeding via allele mining in the base population. Therefore, MP-RMPs such as BRP13R emerge as multipurpose resources for efficient gene discovery and deployment for breeding sorghum cultivars adapted to tropical soils.

Long Non-Coding RNAs, the Dark Matter: An Emerging Regulatory Component in Plants

Long non-coding RNAs (lncRNAs) are pervasive transcripts of longer than 200 nucleotides and indiscernible coding potential. lncRNAs are implicated as key regulatory molecules in various fundamental biological processes at transcriptional, post-transcriptional, and epigenetic levels. Advances in computational and experimental approaches have identified numerous lncRNAs in plants. lncRNAs have been found to act as prime mediators in plant growth, development, and tolerance to stresses. This review summarizes the current research status of lncRNAs in planta, their classification based on genomic context, their mechanism of action, and specific bioinformatics tools and resources for their identification and characterization. Our overarching goal is to summarize recent progress on understanding the regulatory role of lncRNAs in plant developmental processes such as flowering time, reproductive growth, and abiotic stresses. We also review the role of lncRNA in nutrient stress and the ability to improve biotic stress tolerance in plants. Given the pivotal role of lncRNAs in various biological processes, their functional characterization in agriculturally essential crop plants is crucial for bridging the gap between phenotype and genotype.

Independent recurrent evolution of MICRORNA genes converging onto similar non-canonical organisation across green plant lineages is driven by local and segmental duplication events in species, family and lineages

The relationship between evolutionary history, organisation and transcriptional regulation of genes are intrinsically linked. These have been well studied in canonically organised protein-coding genes but not of MIRNAs. In the present study, we investigated the non-canonical arrangement of MIRNAs across taxonomic boundaries from algae to angiosperms employing a combination of genome organization, phylogeny and synteny. We retrieved the complete dataset of MIRNA from twenty-five species to identify and classify based on organisational patterns. The median size of cluster was between 2-5 kb and between 1-20 % of all MIRNAs are organized in head-to-head (with bidirectional promoter), head-to-tail (tandem), and overlapping manner. Although majority of the clusters are composed of MIRNA homologs, 25% of all clusters comprises of non-homologous genes with a potential of generating functional and regulatory complexity. A comparison of phylogeny and organizational patterns revealed that multiple independent events, some of which are species-specific, and some ancient, in different lineages, are responsible for non-canonical organization. Detailed investigation of MIR395 family across the plants revealed a complex origin of non-canonical arrangement through ancient and recent, segmental and local duplications; analysis of MIR399 family revealed major expansion occurred prior to monocot-dicot split, with few lineage-specific events. Evolution of "convergent" organization pattern of non-canonical arrangement originating from independent loci through recurrent event highlights our poor understanding of evolutionary process of MIRNA genes. The present investigation thus paves way for comparative functional genomics to understand the role of non-canonical organization on transcriptional regulation and regulatory diversity in MIRNA gene families.

Interactions and links among the noncoding RNAs in plants under stresses

The complex interplay among sRNAs, lncRNAs and circRNAs has been implicated in plants under biotic and abiotic stresses. Here, we review current advances in our understanding of ncRNA interactions and links, which have considerable potential for improving the agronomic traits and the environmental adaptability of plants. Plants can respond to biotic or abiotic stresses. To cope with various conditions, numerous intricate molecular regulatory mechanisms have evolved in plants. Noncoding RNAs (ncRNAs) can be divided into small noncoding RNAs (sRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs). Emerging evidence has demonstrated that interplay among the ncRNAs acts as a novel layer in the regulatory mechanisms, which has attracted substantial interest. Links between sRNAs can affect plant immune responses and development in synergistic or antagonistic manners. Additionally, multiple interactions between lncRNAs and sRNAs are involved in crop breeding, disease resistance and high tolerance to environmental stresses. Here, we summarize current knowledge of the interactions and links among the ncRNAs in plant responses to stresses and the methods for identifying ncRNA interactions. Furthermore, challenges and prospects for further progress in elucidating ncRNA interactions and links are highlighted.

MicroRNA transcriptomic analysis of the sixth leaf of maize (Zea mays L.) revealed a regulatory mechanism of jointing stage heterosis

BACKGROUND

Zhengdan 958 (Zheng 58 × Chang 7-2), a commercial hybrid that is produced in a large area in China, is the result of the successful use of the heterotic pattern of Reid × Tang-SPT. The jointing stage of maize is the key period from vegetative to reproductive growth, which determines development at later stages and heterosis to a certain degree. MicroRNAs (miRNAs) play vital roles in the regulation of plant development, but how they function in the sixth leaf at the six-leaf (V6) stage to influence jointing stage heterosis is still unclear.

RESULT

Our objective was to study miRNAs in four hybrid combinations developed in accordance with the Reid × Tang-SPT pattern, Zhengdan 958, Anyu 5 (Ye 478 × Chang 7-2), Ye 478 × Huangzaosi, Zheng 58 × Huangzaosi, and their parental inbred lines to explore the mechanism related to heterosis. A total of 234 miRNAs were identified in the sixth leaf at the V6 stage, and 85 miRNAs were differentially expressed between the hybrid combinations and their parental inbred lines. Most of the differentially expressed miRNAs were non-additively expressed, which indicates that miRNAs may participate in heterosis at the jointing stage. miR164, miR1432 and miR528 families were repressed in the four hybrid combinations, and some miRNAs, such as miR156, miR399, and miR395 families, exhibited different expression trends in different hybrid combinations, which may result in varying effects on the heterosis regulatory mechanism.

CONCLUSIONS

The potential targets of the identified miRNAs are related to photosynthesis, the response to plant hormones, and nutrient use. Different hybrid combinations employ different mature miRNAs of the same miRNA family and exhibit different expression trends that may result in enhanced or repressed gene expression to regulate heterosis. Taken together, our results reveal a miRNA-mediated network that plays a key role in jointing stage heterosis via posttranscriptional regulation.

Root Adaptation via Common Genetic Factors Conditioning Tolerance to Multiple Stresses for Crops Cultivated on Acidic Tropical Soils

Crop tolerance to multiple abiotic stresses has long been pursued as a Holy Grail in plant breeding efforts that target crop adaptation to tropical soils. On tropical, acidic soils, aluminum (Al) toxicity, low phosphorus (P) availability and drought stress are the major limitations to yield stability. Molecular breeding based on a small suite of pleiotropic genes, particularly those with moderate to major phenotypic effects, could help circumvent the need for complex breeding designs and large population sizes aimed at selecting transgressive progeny accumulating favorable alleles controlling polygenic traits. The underlying question is twofold: do common tolerance mechanisms to Al toxicity, P deficiency and drought exist? And if they do, will they be useful in a plant breeding program that targets stress-prone environments. The selective environments in tropical regions are such that multiple, co-existing regulatory networks may drive the fixation of either distinctly different or a smaller number of pleiotropic abiotic stress tolerance genes. Recent studies suggest that genes contributing to crop adaptation to acidic soils, such as the major Arabidopsis Al tolerance protein, AtALMT1, which encodes an aluminum-activated root malate transporter, may influence both Al tolerance and P acquisition via changes in root system morphology and architecture. However, trans-acting elements such as transcription factors (TFs) may be the best option for pleiotropic control of multiple abiotic stress genes, due to their small and often multiple binding sequences in the genome. One such example is the C2H2-type zinc finger, AtSTOP1, which is a transcriptional regulator of a number of Arabidopsis Al tolerance genes, including AtMATE and AtALMT1, and has been shown to activate AtALMT1, not only in response to Al but also low soil P. The large WRKY family of transcription factors are also known to affect a broad spectrum of phenotypes, some of which are related to acidic soil abiotic stress responses. Hence, we focus here on signaling proteins such as TFs and protein kinases to identify, from the literature, evidence for unifying regulatory networks controlling Al tolerance, P efficiency and, also possibly drought tolerance. Particular emphasis will be given to modification of root system morphology and architecture, which could be an important physiological "hub" leading to crop adaptation to multiple soil-based abiotic stress factors.

The miR399-CsUBC24 Module Regulates Reproductive Development and Male Fertility in Citrus

MicroRNA399 (miR399) regulates phosphate homeostasis in plants by down-regulating the expression of PHOSPHATE2 (PHO2, or UBC24 encoding the ubiquitin-conjugating E2 enzyme). We previously identified CsmiR399a.1 in a small RNA sequencing screen of a male-sterile somatic cytoplasmic hybrid (or cybrid) of pummelo (Citrus grandis). Here, we report that miR399 affects reproductive development and male fertility in citrus. Down-regulation of CsmiR399a.1 using a short tandem target mimic (STTM) led to abnormal floral development, inhibition of anther dehiscence, and decreased pollen fertility. When grown in inorganic phosphate (Pi)-sufficient conditions, CsmiR399a.1-STTM plants had lower total phosphorus content in their leaves than the wild type and showed typical symptoms of Pi deficiency. In CsmiR399a.1-STTM plants, the expression of genes involved in starch metabolism and Pi homeostasis was significantly different than in the wild type. Thus, we conclude that miR399-STTM mimicked Pi deficiency, disturbed starch metabolism, and was responsible for pollen grain collapse in the transgenic lines. We identified CsUBC24, a citrus homolog of Arabidopsis (Arabidopsis thaliana) AtUBC24 (PHO2), as a target of CsmiR399a.1 that physically interacts with the floral development regulators SEPALLATA family (CsSEP1.1, CsSEP1.2, and CsSEP3) and the anther dehiscence regulator INDUCER OF CBF EXPRESSION1 (CsICE1). We hypothesize that CsUBC24 downregulates the CsSEPs, which disrupts the floral meristem identity regulatory network and leads to developmental abnormalities in flowers. By interacting with CsICE1, CsUBC24 disturbs stomate function on the anther surface, which inhibits anther dehiscence. These findings indicate that a miR399-based mechanism influences both reproductive development and male fertility in citrus.

Differential Gene Expression Responding to Low Phosphate Stress in Leaves and Roots of Maize by cDNA-SRAP

Phosphate (Pi) deficiency in soil can have severe impacts on the growth, development, and production of maize worldwide. In this study, a cDNA-sequence-related amplified polymorphism (cDNA-SRAP) transcript profiling technique was used to evaluate the gene expression in leaves and roots of maize under Pi stress for seven days. A total of 2494 differentially expressed fragments (DEFs) were identified in response to Pi starvation with 1202 and 1292 DEFs in leaves and roots, respectively, using a total of 60 primer pairs in the cDNA-SRAP analysis. These DEFs were categorized into 13 differential gene expression patterns. Results of sequencing and functional analysis showed that 63 DEFs (33 in leaves and 30 in roots) were annotated to a total of 54 genes involved in diverse groups of biological pathways, including metabolism, photosynthesis, signal transduction, transcription, transport, cellular processes, genetic information, and organismal system. This study demonstrated that (1) the cDNA-SRAP transcriptomic profiling technique is a powerful method to analyze differential gene expression in maize showing advantageous features among several transcriptomic methods; (2) maize undergoes a complex adaptive process in response to low Pi stress; and (3) a total of seven differentially expressed genes were identified in response to low Pi stress in leaves or roots of maize and could be used in the genetic modification of maize.

A long non-coding apple RNA, MSTRG.85814.11, acts as a transcriptional enhancer of SAUR32 and contributes to the Fe-deficiency response

Iron (Fe) is an essential plant nutrient and its deficiency typically limits plant growth. Long non-coding (lnc) RNAs are involved in adaptive responses to nutrient stress; however, it is not known whether they function in the regulation of the canonical Fe-deficiency response. The expression of Malus domestica (apple) lncRNA MSTRG.85814 is induced by Fe deficiency, as identified by high-throughput strand-specific RNA-seq analysis of an apple homograft system. MSTRG.85814 has a complex structure, with 13 predicted RNA sequence variants, four of which are upregulated in the roots of plants experiencing Fe deficiency. We found that one MSTRG.85814 splice variant (MSTRG.85814.11) positively modulated its cis target mRNA derived from the small auxin upregulated gene SAUR32. This in turn promoted the expression of SAUR32 and caused an increase in the expression of a plasma membrane proton ATPase, AHA10. Using a pH imaging technique, a significant decrease in the apoplastic pH was observed to occur in the root tips of MSTRG.85814.11 or SAUR32-overexpressing apple plants. Thus MSTRG.85814.11 was shown to positively promote SAUR32 expression, which then activated proton extrusion involved in the Fe-deficiency response. These results reveal a mechanism by which lncRNA promotes environmental Fe-deficiency stress adaption.

Epigenomic landscape and epigenetic regulation in maize

Epigenetic regulation has been implicated in the control of multiple agronomic traits in maize. Here, we review current advances in our understanding of epigenetic regulation, which has great potential for improving agronomic traits and the environmental adaptability of crops. Epigenetic regulation plays vital role in the control of complex agronomic traits. Epigenetic variation could contribute to phenotypic diversity and can be used to improve the quality and productivity of crops. Maize (Zea mays L.), one of the most widely cultivated crops for human food, animal feed, and ethanol biofuel, is a model plant for genetic studies. Recent advances in high-throughput sequencing technology have made possible the study of epigenetic regulation in maize on a genome-wide scale. In this review, we discuss recent epigenetic studies in maize many achieved by Chinese research groups. These studies have explored the roles of DNA methylation, posttranslational modifications of histones, chromatin remodeling, and noncoding RNAs in the regulation of gene expression in plant development and environment response. We also provide our future prospects for manipulating epigenetic regulation to improve crops.

MiR399d and epigenetic modification comodulate anthocyanin accumulation in Malus leaves suffering from phosphorus deficiency

Inorganic phosphorus (Pi) deficiency induces anthocyanin accumulation in the leaves of some plant species; however, the molecular mechanisms underlying this phenomenon have not been well characterized. Here, we showed that microRNA399d (miR399d), high-affinity Pi transporter McPHT1;4, and McMYB10 are strongly induced in Malus leaves suffering from Pi deficiency. By culturing explants of transiently transformed plants in MS medium under conditions of Pi sufficiency and Pi deficiency, miR399d and McPHT1;4 were shown to play essential roles in the response to Pi deficiency and to play positive roles in the regulation of anthocyanin biosynthesis. Silencing of McHDA6 expression and treatment with the inhibitor trichostatin A suggested that the low expression of McHDA6 simultaneously reduced the transcription of McMET1 and decreased the methylation level of the McMYB10 promoter; however, the expression of McMYB10 and anthocyanin content were increased. Bimolecular fluorescence complementation and yeast two-hybrid assays revealed that McHDA6 binds directly to McMET1 through its BAH2 and DNMT1-RFD domains. Based on the results of our study, we propose a mechanism for the molecular regulation of anthocyanin biosynthesis, namely, the miR399d and epigenetic modification comodulation model, to explain the phenomenon in which leaves turn red under conditions of Pi deficiency.

Distinct non-coding RNAs confer root-dependent sense transgene-induced post-transcriptional gene silencing and nitrogen-dependent post-transcriptional regulation to AtAMT1;1 transcripts in Arabidopsis roots

High-affinity ammonium uptake in roots mediate by AMT1-type ammonium transporters, which are tightly controlled at multiple regulatory levels for adapting various nitrogen availability. For Arabidopsis AtAMT1;1 gene, in addition to the transcriptional and post-translational controls, an organ-dependent and N-dependent post-transcriptional regulation was suggested as an additional regulatory step for fine tuning ammonium uptake, but the underlying mechanisms remain to be elucidated. Here, we showed that degradation of AtAMT1;1 transcript in roots of Pro35s:AtAMT1;1-transformed atamt1;1-1 Arabidopsis plants resulted from RDR6-dependent sense transgene-induced post-transcriptional gene silencing (S-PTGS). The siRNAs for S-PTGS may derive from the aberrant RNA, of which the production was co-determined by sequence feature and excessive expression of AtAMT1;1. Switching to the expression of AtAMT1;1 driven by ProAtUBQ10 or of AtAMT1;1 mutated at two siRNA-targeted hotspots reduced AtAMT1;1-specific siRNAs and overcame S-PTGS in roots. In roots of these lines, however, the steady-state transcript levels of AtAMT1;1 still significantly decreased under conditions of N-sufficiency compared with N-deficiency, confirming a N-dependent post-transcriptional regulatory manner. A crucial role of the 207-bp 3'-end sequence of AtAMT1;1 was further demonstrated by N-dependent accumulation of chimeric-AtAMT1;1 transcript in T-DNA insertion lines and of GFP-tagged chimeric-AtAMT1;1 transcript in transgenic lines. A novel non-coding RNA (ncRNA), which was highly abundant in N-sufficient roots, may target the above-identified 3'-end region for the degrading AtAMT1;1 transcript. This degradation could be prevented by a mutation on the AtAMT1;1 transcript at a potential cleavage site (+1458). These results suggested two distinct mechanisms of regulating AtAMT1;1 mRNA turnover by ncRNA for strictly control of ammonium uptake in roots.

LncRNA39026 Enhances Tomato Resistance to Phytophthora infestans by Decoying miR168a and Inducing PR Gene Expression

Our previous study has indicated that a long noncoding RNA (lncRNA), lncRNA39026, can be responsive to Phytophthora infestans infection. However, the function and regulation mechanism of lncRNA39026 during tomato resistance to P. infestans are unknown. In this study, an lncRNA39026 sequence was cloned from tomato Zaofen No. 2, and this sequence contained an endogenous target mimicry for miR168a, which might suppress the expression of miR168a. LncRNA39026 was strongly downregulated at 3 h in the tomato plants infected with P. infestans, and its expression level displayed a negative correlation with the expression level of miR168a and a positive correlation with the expression levels of SlAGO1 genes (target gene of miR168a) upon P. infestans infection. Tomato plants in which lncRNA39026 was overexpressed displayed enhanced resistance to P. infestans, decreased level of miR168a, and increased level of SlAGO1, whereas the resistance was impaired, level of miR168a was increased, and level of SlAGO1 was decreased after lncRNA39026 silencing. In addition, lncRNA39026 could also induce the expression of pathogenesis-related (PR) genes, as shown by increased and decreased expression levels of PR genes in tomato plants with overexpressed and silenced lncRNA39026, respectively. The result demonstrated that lncRNA39026 might function to decoy miR168a and affect the expression of PR genes in tomato plants, increasing resistance to disease.

Research Progress on Plant Long Non-Coding RNA

Non-coding RNAs (ncRNAs) that were once considered “dark matter” or “transcriptional noise” in genomes are research hotspots in the field of epigenetics. The most well-known microRNAs (miRNAs) are a class of short non-coding, small molecular weight RNAs with lengths of 20–24 nucleotides that are highly conserved throughout evolution. Through complementary pairing with the bases of target sites, target gene transcripts are cleaved and degraded, or translation is inhibited, thus regulating the growth and development of organisms. Unlike miRNAs, which have been studied thoroughly, long non-coding RNAs (lncRNAs) are a group of poorly conserved RNA molecules with a sequence length of more than 200 nucleotides and no protein encoding capability; they interact with large molecules, such as DNA, RNA, and proteins, and regulate protein modification, chromatin remodeling, protein functional activity, and RNA metabolism in vivo through cis- or trans-activation at the transcriptional, post-transcriptional, and epigenetic levels. Research on plant lncRNAs is just beginning and has gradually emerged in the field of plant molecular biology. Currently, some studies have revealed that lncRNAs are extensively involved in plant growth and development and stress response processes by mediating the transmission and expression of genetic information. This paper systematically introduces lncRNA and its regulatory mechanisms, reviews the current status and progress of lncRNA research in plants, summarizes the main techniques and strategies of lncRNA research in recent years, and discusses existing problems and prospects, in order to provide ideas for further exploration and verification of the specific evolution of plant lncRNAs and their biological functions.

Transcriptome-wide identification of novel circular RNAs in soybean in response to low-phosphorus stress

Low-phosphorus (LP) stress is a major factor limiting the growth and yield of soybean. Circular RNAs (circRNAs) are novel noncoding RNAs that play a crucial role in plant responses to abiotic stress. However, how LP stress mediates the biogenesis of circRNAs in soybean remains unclear. Here, to explore the response mechanisms of circRNAs to LP stress, the roots of two representative soybean genotypes with different P-use efficiency, Bogao (a LP-sensitive genotype) and Nannong 94156 (a LP-tolerant genotype), were used for the construction of RNA sequencing (RNA-seq) libraries and circRNA identification. In total, 371 novel circRNA candidates, including 120 significantly differentially expressed (DE) circRNAs, were identified across different P levels and genotypes. More DE circRNAs were significantly regulated by LP stress in Bogao than in NN94156, suggesting that the tolerant genotype was less affected by LP stress than the sensitive genotype was; in other words, NN94156 may have a better ability to maintain P homeostasis under LP stress. Moreover, a positive correlation was observed between the expression patterns of P stress-induced circRNAs and their circRNA-host genes. Gene Ontology (GO) enrichment analysis of these circRNA-host genes and microRNA (miRNA)-targeted genes indicated that these DE circRNAs were involved mainly in defense responses, ADP binding, nucleoside binding, organic substance catabolic processes, oxidoreductase activity, and signal transduction. Together, our results revealed that LP stress can significantly alter the genome-wide profiles of circRNAs and indicated that the regulation of circRNAs was both genotype and environment specific in response to LP stress. LP-induced circRNAs might provide a rich resource for LP-responsive circRNA candidates for future studies.

Analysis of the genetic architecture of maize kernel size traits by combined linkage and association mapping

Kernel size-related traits are the most direct traits correlating with grain yield. The genetic basis of three kernel traits of maize, kernel length (KL), kernel width (KW) and kernel thickness (KT), was investigated in an association panel and a biparental population. A total of 21 single nucleotide polymorphisms (SNPs) were detected to be most significantly (P < 2.25 × 10^-6 ) associated with these three traits in the association panel under four environments. Furthermore, 50 quantitative trait loci (QTL) controlling these traits were detected in seven environments in the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, of which eight were repetitively identified in at least three environments. Combining the two mapping populations revealed that 56 SNPs (P < 1 × 10^-3 ) fell within 18 of the QTL confidence intervals. According to the top significant SNPs, stable-effect SNPs and the co-localized SNPs by association analysis and linkage mapping, a total of 73 candidate genes were identified, regulating seed development. Additionally, seven miRNAs were found to situate within the linkage disequilibrium (LD) regions of the co-localized SNPs, of which zma-miR164e was demonstrated to cleave the mRNAs of Arabidopsis CUC1, CUC2 and NAC6 in vitro. Overexpression of zma-miR164e resulted in the down-regulation of these genes above and the failure of seed formation in Arabidopsis pods, with the increased branch number. These findings provide insights into the mechanism of seed development and the improvement of molecular marker-assisted selection (MAS) for high-yield breeding in maize.

Dynamic architecture and regulatory implications of the miRNA network underlying the response to stress in melon

miRNAs are small RNAs that regulate mRNAs at both transcriptional and posttranscriptional level. In plants, miRNAs are involved in the regulation of different processes including development and stress-response. Elucidating how stress-responsive miRNAs are regulated is key to understand the global response to stress but also to develop efficient biotechnological tools that could help to cope with stress. Here, we describe a computational approach based on sRNA sequencing, transcript quantification and degradome data to analyse the accumulation, function and structural organization of melon miRNAs reactivated under seven biotic and abiotic stress conditions at two and four days post-treatment. Our pipeline allowed us to identify fourteen stress-responsive miRNAs (including evolutionary conserved such as miR156, miR166, miR172, miR319, miR398, miR399, miR894 and miR408) at both analysed times. According to our analysis miRNAs were categorized in three groups showing a broad-, intermediate- or narrow- response range. miRNAs reactive to a broad range of environmental cues appear as central components in the stress-response network. The strictly coordinated response of miR398 and miR408 (broad response-range) to the seven stress treatments during the period analysed here reinforces this notion. Although both, the amplitude and diversity of the miRNA-related response to stress changes during the exposition time, the architecture of the miRNA-network is conserved. This organization of miRNA response to stress is also conserved in rice and soybean supporting the conservation of miRNA-network organization in other crops. Overall, our work sheds light into how miRNA networks in plants organize and function during stress.