Boron-deficiency-responsive microRNAs and their targets in Citrus sinensis leaves

Yerba mate (Ilex paraguariensis) genome provides new insights into convergent evolution of caffeine biosynthesis

Yerba mate (YM, Ilex paraguariensis) is an economically important crop marketed for the elaboration of mate, the third-most widely consumed caffeine-containing infusion worldwide. Here, we report the first genome assembly of this species, which has a total length of 1.06 Gb and contains 53,390 protein-coding genes. Comparative analyses revealed that the large YM genome size is partly due to a whole-genome duplication (Ip-α) during the early evolutionary history of Ilex, in addition to the hexaploidization event (γ) shared by core eudicots. Characterization of the genome allowed us to clone the genes encoding methyltransferase enzymes that catalyse multiple reactions required for caffeine production. To our surprise, this species has converged upon a different biochemical pathway compared to that of coffee and tea. In order to gain insight into the structural basis for the convergent enzyme activities, we obtained a crystal structure for the terminal enzyme in the pathway that forms caffeine. The structure reveals that convergent solutions have evolved for substrate positioning because different amino acid residues facilitate a different substrate orientation such that efficient methylation occurs in the independently evolved enzymes in YM and coffee. While our results show phylogenomic constraint limits the genes coopted for convergence of caffeine biosynthesis, the X-ray diffraction data suggest structural constraints are minimal for the convergent evolution of individual reactions.

Genome-Wide Comprehensive Identification and In Silico Characterization of Lectin Receptor-Like Kinase Gene Family in Barley (Hordeum vulgare L.)

Lectin receptor-like kinases (LecRLKs) are a significant subgroup of the receptor-like kinases (RLKs) protein family. They play crucial roles in plant growth, development, immune responses, signal transduction, and stress tolerance. However, the genome-wide identification and characterization of LecRLK genes and their regulatory elements have not been explored in a major cereal crop, barley (Hordeum vulgare L.). Therefore, in this study, integrated bioinformatics tools were used to identify and characterize the LecRLK gene family in barley. Based on the phylogenetic tree and domain organization, a total of 113 LecRLK genes were identified in the barley genome (referred to as HvlecRLK) corresponding to the LecRLK genes of Arabidopsis thaliana. These putative HvlecRLK genes were classified into three groups: 62 G-type LecRLKs, 1 C-type LecRLK, and 50 L-type LecRLKs. They were unevenly distributed across eight chromosomes, including one unknown chromosome, and were predominantly located in the plasma membrane (G-type HvlecRLK (96.8%), C-type HvlecRLK (100%), and L-type HvlecRLK (98%)). An analysis of motif composition and exon-intron configuration revealed remarkable homogeneity with the members of AtlecRLK. Notably, most of the HvlecRLKs (27 G-type, 43 L-type) have no intron, suggesting their rapid functionality. The Ka/Ks and syntenic analysis demonstrated that HvlecRLK gene pairs evolved through purifying selection and gene duplication was the major factor for the expansion of the HvlecRLK gene family. Exploration of gene ontology (GO) enrichment indicated that the identified HvlecRLK genes are associated with various cellular processes, metabolic pathways, defense mechanisms, kinase activity, catalytic activity, ion binding, and other essential pathways. The regulatory network analysis identified 29 transcription factor families (TFFs), with seven major TFFs including bZIP, C2H2, ERF, MIKC_MADS, MYB, NAC, and WRKY participating in the regulation of HvlecRLK gene functions. Most notably, eight TFFs were found to be linked to the promoter region of both L-type HvleckRLK64 and HvleckRLK86. The promoter cis-acting regulatory element (CARE) analysis of barley identified a total of 75 CARE motifs responsive to light responsiveness (LR), tissue-specific (TS), hormone responsiveness (HR), and stress responsiveness (SR). The maximum number of CAREs was identified in HvleckRLK11 (25 for LR), HvleckRLK69 (17 for TS), and HvleckRLK80 (12 for HR). Additionally, HvleckRLK14, HvleckRLK16, HvleckRLK33, HvleckRLK50, HvleckRLK52, HvleckRLK56, and HvleckRLK110 were predicted to exhibit higher responses in stress conditions. In addition, 46 putative miRNAs were predicted to target 81 HvlecRLK genes and HvlecRLK13 was the most targeted gene by 8 different miRNAs. Protein-protein interaction analysis demonstrated higher functional similarities of 63 HvlecRLKs with 7 Arabidopsis STRING proteins. Our overall findings provide valuable information on the LecRLK gene family which might pave the way to advanced research on the functional mechanism of the candidate genes as well as to develop new barley cultivars in breeding programs.

Differential Expression of miRNAs Involved in Response to Candidatus Liberibacter asiaticus Infection in Mexican Lime at Early and Late Stages of Huanglongbing Disease

Huanglongbing (HLB) is one of the most destructive diseases threatening citriculture worldwide. This disease has been associated with α-proteobacteria species, namely Candidatus Liberibacter. Due to the unculturable nature of the causal agent, it has been difficult to mitigate the disease, and nowadays a cure is not available. MicroRNAs (miRNAs) are key regulators of gene expression, playing an essential role in abiotic and biotic stress in plants including antibacterial responses. However, knowledge derived from non-model systems including Candidatus Liberibacter asiaticus (CLas)-citrus pathosystem remains largely unknown. In this study, small RNA profiles from Mexican lime (Citrus aurantifolia) plants infected with CLas at asymptomatic and symptomatic stages were generated by sRNA-Seq, and miRNAs were obtained with ShortStack software. A total of 46 miRNAs, including 29 known miRNAs and 17 novel miRNAs, were identified in Mexican lime. Among them, six miRNAs were deregulated in the asymptomatic stage, highlighting the up regulation of two new miRNAs. Meanwhile, eight miRNAs were differentially expressed in the symptomatic stage of the disease. The target genes of miRNAs were related to protein modification, transcription factors, and enzyme-coding genes. Our results provide new insights into miRNA-mediated regulation in C. aurantifolia in response to CLas infection. This information will be useful to understand molecular mechanisms behind the defense and pathogenesis of HLB.

Characterization of plant laccase genes and their functions

Laccase is a copper-containing polyphenol oxidase found in different organisms. The multigene family that encodes laccases is widely distributed in plant genomes. Plant laccases oxidize monolignols to produce lignin which is important for plant growth and stress responses. Industrial applications of fungal and bacterial laccases are extensively explored and addressed. Recently many studies have focused on the significance of plant laccase, particularly in crop yield, and its functions in different environmental conditions. This review summarizes the transcriptional and posttranscriptional regulation of plant laccase genes and their functions in plant growth and development. It especially describes the responses of laccase genes to various stresses and their contributions to plant biotic and abiotic stress resistance. In-depth explanations and scientific advances will serve as foundations for research into plant laccase genes' function, mechanism, and possible applications.

TIR1/AFB proteins: Active players in abiotic and biotic stress signaling

The TIR1/AFB family of proteins is a group of functionally diverse auxin receptors that are only found in plants. TIR1/AFB family members are characterized by a conserved N-terminal F-box domain followed by 18 leucine-rich repeats. In the past few decades, extensive research has been conducted on the role of these proteins in regulating plant development, metabolism, and responses to abiotic and biotic stress. In this review, we focus on TIR1/AFB proteins that play crucial roles in plant responses to diverse abiotic and biotic stress. We highlight studies that have shed light on the mechanisms by which TIR1/AFB proteins are regulated at the transcriptional and post-transcriptional as well as the downstream in abiotic or biotic stress pathways regulated by the TIR1/AFB family.

The miR393-Target Module Regulates Plant Development and Responses to Biotic and Abiotic Stresses

MicroRNAs (miRNAs), a class of endogenous small RNAs, are broadly involved in plant development, morphogenesis and responses to various environmental stresses, through manipulating the cleavage, translational expression, or DNA methylation of target mRNAs. miR393 is a conserved miRNA family present in many plants, which mainly targets genes encoding the transport inhibitor response1 (TIR1)/auxin signaling F-box (AFB) auxin receptors, and thus greatly affects the auxin signal perception, Aux/IAA degradation, and related gene expression. This review introduces the advances made on the miR393/target module regulating plant development and the plant’s responses to biotic and abiotic stresses. This module is valuable for genetic manipulation of optimized conditions for crop growth and development and would also be helpful in improving crop yield through molecular breeding.

MicroRNA Mediated Plant Responses to Nutrient Stress

To complete their life cycles, plants require several minerals that are found in soil. Plant growth and development can be affected by nutrient shortages or high nutrient availability. Several adaptations and evolutionary changes have enabled plants to cope with inappropriate growth conditions and low or high nutrient levels. MicroRNAs (miRNAs) have been recognized for transcript cleavage and translational reduction, and can be used for post-transcriptional regulation. Aside from regulating plant growth and development, miRNAs play a crucial role in regulating plant’s adaptations to adverse environmental conditions. Additionally, miRNAs are involved in plants’ sensory functions, nutrient uptake, long-distance root transport, and physiological functions related to nutrients. It may be possible to develop crops that can be cultivated in soils that are either deficient in nutrients or have extreme nutrient supplies by understanding how plant miRNAs are associated with nutrient stress. In this review, an overview is presented regarding recent advances in the understanding of plants’ responses to nitrogen, phosphorus, potassium, sulfur, copper, iron, boron, magnesium, manganese, zinc, and calcium deficiencies via miRNA regulation. We conclude with future research directions emphasizing the modification of crops for improving future food security.

Comprehensive Identification and Profiling of miRNAs Involved in Terpenoid Synthesis of Gleditsia sinensis Lam

Gleditsia sinensis Lam. is a tree with worldwide distribution and important economic and medicinal values; its pods contain terpenoids including gleditsioside, thiamine, and brassinosteroids. However, thus far, there are few studies on the terpenoid regulation of G. sinensis at the molecular level. microRNA (miRNA) is a class of small RNAs with conserved and crucial roles in the regulation of diverse biological processes during plant growth and development. To identify the miRNAs of G. sinensis and evaluate their involvement in terpenoid synthesis, this investigation quantified the content changes in saponins in pods at three developmental stages: May (pod-setting stage), July (elongation stage), and September (browning stage), and then we performed genome-wide miRNA profiles during the three development stages of the G. sinensis pods. A total of 351 conserved miRNAs belonging to 216 families were identified, among which 36 conserved miRNAs exist specifically in legumes. Through target analysis, 708 unigenes were predicted to be candidate targets of 37 differentially expressed miRNAs. The targets of miR838-3p and miR2093-5p were involved in the derived branches of monoterpenes and gleditsioside, in brassinosteroid biosynthesis (BRB), and in indole alkaloid biosynthesis (IAB). Intriguingly, the targets of miR829-3p.1 were predicted to take part in thiamine biosynthesis, and the targets of miR4414b and miR5037a were involved in the main process of cytokinin synthesis. The corresponding targets participated in BRB, IAB, and terpenoid backbone biosynthesis, which were enriched significantly, suggesting that miR2093-5p, miR4414b, miR5037a, miR829-3p.1, and miR838-3p play indispensable roles in the regulation of triterpenoid saponin and monoterpenoid biosynthesis. To date, this is the first report of miRNA identification in G. sinensis and miRNA expression profiles at different developmental stages of G. sinensis pods, which provides a basis for further uncovering the molecular regulation of terpenoid synthesis in G. sinensis and new insights into the role of miRNAs in legumes.

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.

Citrus Physiological and Molecular Response to Boron Stresses

Since the essentiality of boron (B) to plant growth was reported nearly one century ago, the implication of B in physiological performance, productivity and quality of agricultural products, and the morphogenesis of apical meristem in plants has widely been studied. B stresses (B deficiency and toxicity), which lead to atrophy of canopy and deterioration of Citrus fruits, have long been discovered in citrus orchards. This paper reviews the research progress of B stresses on Citrus growth, photosynthesis, light use efficiency, nutrient absorption, organic acid metabolism, sugar metabolism and relocation, and antioxidant system. Moreover, the beneficial effects of B on plant stress tolerance and further research in this area were also discussed.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Next-Generation Sequencing Identification and Characterization of MicroRNAs in Dwarfed Citrus Trees Infected With Citrus Dwarfing Viroid in High-Density Plantings

Citrus dwarfing viroid (CDVd) induces stunting on sweet orange trees [Citrus sinensis (L.) Osbeck], propagated on trifoliate orange rootstock [Citrus trifoliata (L.), syn. Poncirus trifoliata (L.) Raf.]. MicroRNAs (miRNAs) are a class of non-coding small RNAs (sRNAs) that play important roles in the regulation of tree gene expression. To identify miRNAs in dwarfed citrus trees, grown in high-density plantings, and their response to CDVd infection, sRNA next-generation sequencing was performed on CDVd-infected and non-infected controls. A total of 1,290 and 628 miRNAs were identified in stem and root tissues, respectively, and among those, 60 were conserved in each of these two tissue types. Three conserved miRNAs (csi-miR479, csi-miR171b, and csi-miR156) were significantly downregulated (adjusted p-value < 0.05) in the stems of CDVd-infected trees compared to the non-infected controls. The three stem downregulated miRNAs are known to be involved in various physiological and developmental processes some of which may be related to the characteristic dwarfed phenotype displayed by CDVd-infected C. sinensis on C. trifoliata rootstock field trees. Only one miRNA (csi-miR535) was significantly downregulated in CDVd-infected roots and it was predicted to target genes controlling a wide range of cellular functions. Reverse transcription quantitative polymerase chain reaction analysis performed on selected miRNA targets validated the negative correlation between the expression levels of these targets and their corresponding miRNAs in CDVd-infected trees. Our results indicate that CDVd-responsive plant miRNAs play a role in regulating important citrus growth and developmental processes that may participate in the cellular changes leading to the observed citrus dwarf phenotype.

Effects of phosphorus deficiency on the absorption of mineral nutrients, photosynthetic system performance and antioxidant metabolism in Citrus grandis

Phosphorus (P) is an essential macronutrient for plant growth, development and production. However, little is known about the effects of P deficiency on nutrient absorption, photosynthetic apparatus performance and antioxidant metabolism in citrus. Seedlings of 'sour pummelo' (Citrus grandis) were irrigated with a nutrient solution containing 0.2 mM (Control) or 0 mM (P deficiency) KH2PO4 until saturated every other day for 16 weeks. P deficiency significantly decreased the dry weight (DW) of leaves and stems, and increased the root/shoot ratio in C. grandis but did not affect the DW of roots. The decreased DW of leaves and stems might be induced by the decreased chlorophyll (Chl) contents and CO2 assimilation in P deficient seedlings. P deficiency heterogeneously affected the nutrient contents of leaves, stems and roots. The analysis of Chl a fluorescence transients showed that P deficiency impaired electron transport from the donor side of photosystem II (PSII) to the end acceptor side of PSI, which showed a greater impact on the performance of the donor side of PSII than that of the acceptor side of PSII and photosystem I (PSI). P deficiency increased the contents of ascorbate (ASC), H2O2 and malondialdehyde (MDA) as well as the activities of superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APX), dehydroascorbate reductase (DHAR) and glutathione reductase (GR) in leaves. In contrast, P deficiency increased the ASC content, reduced the glutathione (GSH) content and the activities of SOD, CAT, APX and monodehydroascorbate reductase (MDHAR), but did not increase H2O2 production, anthocyanins and MDA content in roots. Taking these results together, we conclude that P deficiency affects nutrient absorption and lowers photosynthetic performance, leading to ROS production, which might be a crucial cause of the inhibited growth of C. grandis.

Integration of small RNA, degradome and proteome sequencing in Oryza sativa reveals a delayed senescence network in tetraploid rice seed

Seed of rice is an important strategic resource for ensuring the security of China's staple food. Seed deterioration as a result of senescence is a major problem during seed storage, which can cause major economic losses. Screening among accessions in rice germplasm resources for traits such as slow senescence and increased seed longevity during storage is, therefore, of great significance. However, studies on delayed senescence in rice have been based mostly on diploid rice seed to date. Despite better tolerance have been verified by the artificial aging treatment for polyploid rice seed, the delayed senescence properties and delayed senescence related regulatory mechanisms of polyploid rice seed are rarely reported, due to the lack of polyploid rice materials with high seed set. High-throughput sequencing was applied to systematically investigate variations in small RNAs, the degradome, and the proteome between tetraploid and diploid rice seeds. Degradome sequencing analysis of microRNAs showed that expression of miR-164d, which regulates genes encoding antioxidant enzymes, was changed significantly, resulting in decreased miRNA-mediated cleavage of target genes in tetraploid rice. Comparisons of the expression levels of small RNAs (sRNAs) in the tetraploid and diploid libraries revealed that 12 sRNAs changed significantly, consistent with the findings from degradome sequencing. Furthermore, proteomics also showed that antioxidant enzymes were up-regulated in tetraploid rice seeds, relative to diploids.

Knock-down of δ-aminolevulinic acid dehydratase via virus-induced gene silencing alters the microRNA biogenesis and causes stress-related reactions in citrus plants

The δ-aminolevulinic acid (δ-ALA) is an intermediate in the biosynthetic pathway of tetrapyrroles. Tetrapyrroles play vital roles in many biological processes such as photosynthesis, respiration, and light-sensing. ALA-dehydratase (ALAD) combines two molecules of δ-ALA to form porphobilinogen. In citrus, the silencing of ALAD caused discrete yellow spots and necrosis in leaves and stems. Additionally, it caused rapid death in developing new shoots. Herein, we hypothesize that the accumulation of δ-ALA results in severe stress and reduced meristem development. For that reason, we investigated the dynamic changes in the expression profiles of 23 microRNA (miRNA) identified through small RNA sequencing, from CTV-tALAD plants in comparison with healthy C. macrophylla and C. macrophylla infiltrated with CTV-wt. Furthermore, we reported the effect of ALAD silencing on the total phenolics, H₂O_2, and reactive oxygen species (ROS) levels, to examine the possibilities of miRNAs involving the regulation of these pathways. Our results showed that the total phenolics content, H₂O₂, and O₂^- levels were increased in CTV-tALAD plants. Moreover, 63 conserved miRNA members belonging to 23 different miRNA families were differentially expressed in CTV-tALAD plants compared to controls. The identified miRNAs are implicated in auxin biosynthesis and signaling, axillary shoot meristem formation and leaf morphology, starch metabolism, and oxidative stress. Collectively, our findings suggested that ALAD silencing initiates stress on citrus plants. As a result, CTV-tALAD plants exhibit reduced metabolic rate, growth, and development in order to cope with the stress that resulted from the accumulation of δ-ALA. This cascade of events led to leaf, stem, and meristem necrosis and failure of new shoot development.

Integrated analyses of miRNAome and transcriptome reveal zinc deficiency responses in rice seedlings

BACKGROUND

Zinc (Zn) deficiency is one of the most widespread soil constraints affecting rice productivity, but the molecular mechanisms underlying the regulation of Zn deficiency response is still limited. Here, we aim to understand the molecular mechanisms of Zn deficiency response by integrating the analyses of the global miRNA and mRNA expression profiles under Zn deficiency and resupply in rice seedlings by integrating Illumina's high-throughput small RNA sequencing and transcriptome sequencing.

RESULTS

The transcriptome sequencing identified 360 genes that were differentially expressed in the shoots and roots of Zn-deficient rice seedlings, and 97 of them were recovered after Zn resupply. A total of 68 miRNAs were identified to be differentially expressed under Zn deficiency and/or Zn resupply. The integrated analyses of miRNAome and transcriptome data showed that 12 differentially expressed genes are the potential target genes of 10 Zn-responsive miRNAs such as miR171g-5p, miR397b-5p, miR398a-5p and miR528-5p. Some miRNA genes and differentially expressed genes were selected for validation by quantitative RT-PCR, and their expressions were similar to that of the sequencing results.

CONCLUSION

These results provide insights into miRNA-mediated regulatory pathways in Zn deficiency response, and provide candidate genes for genetic improvement of Zn deficiency tolerance in rice.

Boron mitigates citrus root injuries by regulating intracellular pH and reactive oxygen species to resist H+-toxicity

Boron (B)-deficiency and H⁺-toxicity are important limiting factors for plants growth in acid soils. High B supply may reduce H⁺-toxicity-induced inhibition of growth in citrus. Trifoliate orange rootstock seedlings were irrigated with nutrient solution containing either 0 μM or 10 μM H₃BO₃ at two pH levels (pH4 (H⁺-toxicity) and pH6 (normal)). The results showed that H⁺-toxicity without B severely hampered main root elongation. Simultaneously, oxidative damage caused by H⁺-toxicity led to severe damage to the apical structure of root such as root crown abscission. However, B application promoted the root length, root cell viability and reduced cell wall (CW) thickness of root tips under H⁺-toxicity. Additionally, B application reduced the H⁺-toxicity-induced reactive oxygen species (ROS) accumulation in roots as characterized by lower fluorescence intensity of H₂O₂ and O₂^- staining. Moreover, ³¹P-NMR (³¹P nuclear magnetic resonance) spectra revealed B application regulated the pH of vacuoles and cytoplasm in root tips by reducing phosphoenolpyruvate carboxykinase (PEPCase) activity while enhancing NADP malic enzyme (NADP-ME) activity during H⁺-toxicity. Collectively, our results demonstrate that B supply alleviates H⁺-toxicity and promotes root growth by reducing ROS accumulation, attenuating intracellular acidic microenvironment to ensure normal chemical reactions in root tip cells.

Tissue-specific transcriptional regulation of seven heavy metal stress-responsive miRNAs and their putative targets in nickel indicator castor bean (R. communis L.) plants

R.communis L. has high capability to accumulate nickel which is a trace nutrient for higher plants and also an environmental contaminant causes toxicity related symptoms at higher concentrations. MicroRNAs (miRNAs) are known to be important modulators of responses against heavy metal stress for detoxification of the metal. In this study, we experimentally measured and validated the transcript levels of the seven heavy metal stress response-related miRNAs and the expression levels of target genes in both leaf and root tissues of R. communis L. subjected to three different concentrations of nickel stress via qRT-PCR quantification. The results demonstrated differential regulations of heavy metal stress-responsive miRNAs and their putative targets in both tissues in same stress treatments. This dynamic regulation suggest that regulatory processes differ between the tissues under nickel stress. Our data suggest that, miR838 was the most responsive to the Ni2+ stress. miR398 target gene Cu-Zn/SOD was found to be up-regulated in both root and leaf tissues. The relations between TCP and expression levels of miR159 and miR319 were also found statistically significant exclusive to leaf tissues. In leaf tissue, changes in miR395 level and its putative target genes, sulphate transporter and sulphate adenyltransferase gene were found in relation whereas, only expression level of sulphate transporter represented a statistically significant relation in root tissue. The sharp decrease in transcript levels of 2r3 myb gene at lower nickel dose suggest to investigate the role of r2r3 myb and the all MYB family members in primary and secondary metabolisms against nickel stress.

Comparative Study of Withanolide Biosynthesis-Related miRNAs in Root and Leaf Tissues of Withania somnifera

Withania somnifera, popularly known as Indian ginseng, is one of the most important medicinal plants. The plant is well studied in terms of its pharmaceutical activities and genes involved in biosynthetic pathways. However, not much is known about the regulatory mechanism of genes responsible for the production of secondary metabolites. The idea was to identify miRNA transcriptome responsible for the regulation of withanolide biosynthesis, specifically of root and leaf tissues individually. The transcriptome data of in vitro culture of root and leaf tissues of the plant was considered for miRNA identification. A total of 24 and 39 miRNA families were identified in root and leaf tissues, respectively. Out of these, 15 and 27 miRNA families have shown their involvement in different biological functions in root and leaf tissues, respectively. We report here, specific miRNAs and their corresponding target genes for corresponding root and leaf tissues. The target genes have also been analyzed for their role in withanolide metabolism. Endogenous root-miR5140, root-miR159, leaf-miR477, and leaf-miR530 were reported for regulation of withanolide biosynthesis.

In Silico Analysis of Small RNAs Suggest Roles for Novel and Conserved miRNAs in the Formation of Epigenetic Memory in Somatic Embryos of Norway Spruce

Epigenetic memory in Norway spruce affects the timing of bud burst and bud set, vitally important adaptive traits for this long-lived forest species. Epigenetic memory is established in response to the temperature conditions during embryogenesis. Somatic embryogenesis at different epitype inducing (EpI) temperatures closely mimics the natural processes of epigenetic memory formation in seeds, giving rise to epigenetically different clonal plants in a reproducible and predictable manner, with respect to altered bud phenology. MicroRNAs (miRNAs) and other small non-coding RNAs (sRNAs) play an essential role in the regulation of plant gene expression and may affect this epigenetic mechanism. We used NGS sequencing and computational in silico methods to identify and profile conserved and novel miRNAs among small RNAs in embryogenic tissues of Norway spruce at three EpI temperatures (18, 23 and 28°C). We detected three predominant classes of sRNAs related to a length of 24 nt, followed by a 21–22 nt class and a third 31 nt class of sRNAs. More than 2100 different miRNAs within the prevailing length 21–22 nt were identified. Profiling these putative miRNAs allowed identification of 1053 highly expressed miRNAs, including 523 conserved and 530 novels. 654 of these miRNAs were found to be differentially expressed (DEM) depending on EpI temperature. For most DEMs, we defined their putative mRNA targets. The targets represented mostly by transcripts of multiple-repeats proteins, like TIR, NBS-LRR, PPR and TPR repeat, Clathrin/VPS proteins, Myb-like, AP2, etc. Notably, 124 DE miRNAs targeted 203 differentially expressed epigenetic regulators. Developing Norway spruce embryos possess a more complex sRNA structure than that reported for somatic tissues. A variety of the predicted miRNAs showed distinct EpI temperature dependent expression patterns. These putative EpI miRNAs target spruce genes with a wide range of functions, including genes known to be involved in epigenetic regulation, which in turn could provide a feedback process leading to the formation of epigenetic marks. We suggest that TIR, NBS and LRR domain containing proteins could fulfill more general functions for signal transduction from external environmental stimuli and conversion them into molecular response. Fine-tuning of the miRNA production likely participates in both developmental regulation and epigenetic memory formation in Norway spruce.

MicroRNA-mediated responses to long-term magnesium-deficiency in Citrus sinensis roots revealed by Illumina sequencing

Background

Magnesium (Mg)-deficiency occurs most frequently in strongly acidic, sandy soils. Citrus are grown mainly on acidic and strong acidic soils. Mg-deficiency causes poor fruit quality and low fruit yield in some Citrus orchards. For the first time, we investigated Mg-deficiency-responsive miRNAs in ‘Xuegan’ (Citrus sinensis) roots using Illumina sequencing in order to obtain some miRNAs presumably responsible for Citrus Mg-deficiency tolerance.

Results

We obtained 101 (69) miRNAs with increased (decreased) expression from Mg-starved roots. Our results suggested that the adaptation of Citrus roots to Mg-deficiency was related to the several aspects: (a) inhibiting root respiration and related gene expression via inducing miR158 and miR2919; (b) enhancing antioxidant system by down-regulating related miRNAs (miR780, miR6190, miR1044, miR5261 and miR1151) and the adaptation to low-phosphorus (miR6190); (c) activating transport-related genes by altering the expression of miR6190, miR6485, miR1044, miR5029 and miR3437; (d) elevating protein ubiquitination due to decreased expression levels of miR1044, miR5261, miR1151 and miR5029; (e) maintaining root growth by regulating miR5261, miR6485 and miR158 expression; and (f) triggering DNA repair (transcription regulation) by regulating miR5176 and miR6485 (miR6028, miR6190, miR6485, miR5621, miR160 and miR7708) expression. Mg-deficiency-responsive miRNAs involved in root signal transduction also had functions in Citrus Mg-deficiency tolerance.

Conclusions

We obtained several novel Mg-deficiency-responsive miRNAs (i.e., miR5261, miR158, miR6190, miR6485, miR1151 and miR1044) possibly contributing to Mg-deficiency tolerance. These results revealed some novel clues on the miRNA-mediated adaptation to nutrient deficiencies in higher plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3999-5) contains supplementary material, which is available to authorized users.

Identification of protoplast-isolation responsive microRNAs in Citrus reticulata Blanco by high-throughput sequencing

Protoplast isolation is a stress-inducing process, during which a variety of physiological and molecular alterations take place. Such stress response affects the expression of totipotency of cultured protoplasts. MicroRNAs (miRNAs) play important roles in plant growth, development and stress responses. However, the underlying mechanism of miRNAs involved in the protoplast totipotency remains unclear. In this study, high-throughput sequencing technology was used to sequence two populations of small RNA from calli and callus-derived protoplasts in Citrus reticulata Blanco. A total of 67 known miRNAs from 35 families and 277 novel miRNAs were identified. Among these miRNAs, 18 known miRNAs and 64 novel miRNAs were identified by differentially expressed miRNAs (DEMs) analysis. The expression patterns of the eight DEMs were verified by qRT-PCR. Target prediction showed most targets of the miRNAs were transcription factors. The expression levels of half targets showed a negative correlation to those of the miRNAs. Furthermore, the physiological analysis showed high levels of antioxidant activities in isolated protoplasts. In short, our results indicated that miRNAs may play important roles in protoplast-isolation response.

Systems biology study of transcriptional and post-transcriptional co-regulatory network sheds light on key regulators involved in important biological processes in Citrus sinensis

Transcriptional and post-transcriptional regulators including transcription regulator, transcription factor and miRNA are the main endogenous molecular elements which control complex cellular mechanisms such as development, growth and response to biotic and abiotic stresses in a coordinated manner in plants. Utilizing the most recent information on such relationships in a plant species, obtained from high-throughput experimental technologies and advanced computational tools, we can reconstruct its co-regulatory network which consequently sheds light on key regulators involved in its important biological processes. In this article, combined systems biology approaches such as mining the literatures, various databases and network reconstruction, analysis, and visualization tools were employed to infer and visualize the coregulatory relationships between miRNAs and transcriptional regulators in Citrus sinensis. Using computationally and experimentally verified miRNA-target interactions and constructed co-expression networks on array-based data, 10 coregulatory networks and 10 corresponding subgraphs include FFL motifs were obtained for 10 distinct tissues/conditions. Then PPI subnetworks were extracted for transcripts/genes included in mentioned subgraphs in order to the functional analysis of extracted coregulatory circuits. These proposed coregulatory connections shed light on precisely identifying C. sinensis metabolic pathways key switches, which are demanded for ultimate goals such as genome editing.

Functional Roles of microRNAs in Agronomically Important Plants-Potential as Targets for Crop Improvement and Protection

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have recently emerged as important regulators of gene expression, mainly through cleavage and/or translation inhibition of the target mRNAs during or after transcription. miRNAs play important roles by regulating a multitude of biological processes in plants which include maintenance of genome integrity, development, metabolism, and adaptive responses toward environmental stresses. The increasing population of the world and their food demands requires focused efforts for the improvement of crop plants to ensure sustainable food production. Manipulation of mRNA transcript abundance via miRNA control provides a unique strategy for modulating differential plant gene expression and miRNAs are thus emerging as the next generation targets for genetic engineering for improvement of the agronomic properties of crops. However, a deeper understanding of its potential and the mechanisms involved will facilitate the design of suitable strategies to obtain the desirable traits with minimum trade-offs in the modified crops. In this regard, this review highlights the diverse roles of conserved and newly identified miRNAs in various food and industrial crops and recent advances made in the uses of miRNAs to improve plants of agronomically importance so as to significantly enhance crop yields and increase tolerance to various environmental stress agents of biotic-or abiotic origin.

MicroRNA Regulatory Mechanisms on Citrus sinensis leaves to Magnesium-Deficiency

Magnesium (Mg)-deficiency, which affects crop productivity and quality, widespreadly exists in many agricultural crops, including citrus. However, very limited data are available on Mg-deficiency-responsive microRNAs (miRNAs) in higher plants. Using Illumina sequencing, we isolated 75 (73 known and 2 novel) up- and 71 (64 known and 7 novel) down-regulated miRNAs from Mg-deficient Citrus sinensis leaves. In addition to the remarkable metabolic flexibility as indicated by the great alteration of miRNA expression, the adaptive responses of leaf miRNAs to Mg-deficiency might also involve the following several aspects: (a) up-regulating stress-related genes by down-regulating miR164, miR7812, miR5742, miR3946, and miR5158; (b) enhancing cell transport due to decreased expression of miR3946 and miR5158 and increased expression of miR395, miR1077, miR1160, and miR8019; (c) activating lipid metabolism-related genes by repressing miR158, miR5256, and miR3946; (d) inducing cell wall-related gene expansin 8A by repressing miR779; and (e) down-regulating the expression of genes involved in the maintenance of S, K and Cu by up-regulating miR395 and miR6426. To conclude, we isolated some new known miRNAs (i.e., miR7812, miR8019, miR6218, miR1533, miR6426, miR5256, miR5742, miR5561, miR5158, and miR5818) responsive to nutrient deficiencies and found some candidate miRNAs that might contribute to Mg-deficiency tolerance. Therefore, our results not only provide novel information about the responses of plant to Mg-deficiency, but also are useful for obtaining the key miRNAs for plant Mg-deficiency tolerance.

Identification of miRNAs Affecting the Establishment of Brassica Alboglabra Seedling

MicroRNAs (miRNAs) are important for plant development including seed formation, dormancy, and germination, as well as seedling establishment. The Brassica vegetable seedling establishment stage influences the development of high quality seedlings, but also affects the nutrient content of sprouts. Chinese kale (Brassica alboglabra) seedlings at different growth stages were used to construct two small-RNA (sRNA) libraries. We comprehensively analyzed the miRNAs in 2- and 9-day-old seedlings. An average of 11,722,490 clean reads were generated after removing low-quality reads and adapter contaminants. The results revealed that 37.65 and 26.69% of the sRNAs in 2- and 9-day-old seedlings, respectively, were 24 nt long. In total, 254 known mature miRNA sequences from 228 miRNA families and 343 novel miRNAs were identified. Of these miRNAs, 224 were differentially expressed between the two analyzed libraries. The most abundant miRNAs identified by sequence homology were miR156, miR167, and miR157, each with more than 100,000 sequenced reads. Compared with the expression levels in 2-day-old seedlings, MiR8154 and miR390 were the most up- and down-regulated miRNAs respectively in 9-day-old seedlings. Gene ontology enrichment analysis of the differentially expressed-miRNA target genes affecting biological processes revealed that most genes were in the "regulation of transcription" category. Additionally, the expression patterns of some miRNAs and target genes were validated by quantitative real-time polymerase chain reaction. We determined that development-associated miRNAs (e.g., bal-miR156/157/159/166/167/172/396), were highly-expressed during seedling-establishment stage, as were stress-related (bal-miR408) and metabolism-related (bal-miR826) miRNAs. Combined with the low level of targets SPL9 and AP2, it was concluded that miR156-SPL9 and miR172-AP modules play key roles during the B. alboglabra seedling establishment stage.