Identification and characterization of microRNAs in Phaseolus vulgaris by high-throughput sequencing

Genome-Wide Comparative Analysis of Five Amaranthaceae Species Reveals a Large Amount of Repeat Content

Amaranthus is a genus of C4 dicotyledonous herbaceous plant species that are widely distributed in Asia, Africa, Australia, and Europe and are used as grain, vegetables, forages, and ornamental plants. Amaranth species have gained significant attention nowadays as potential sources of nutritious food and industrial products. In this study, we performed a comparative genome analysis of five amaranth species, namely, Amaranthus hypochondriacus, Amaranthus tuberculatus, Amaranthus hybridus, Amaranthus palmeri, and Amaranthus cruentus. The estimated repeat content ranged from 54.49% to 63.26% and was not correlated with the genome sizes. Out of the predicted repeat classes, the majority of repetitive sequences were Long Terminal Repeat (LTR) elements, which account for about 13.91% to 24.89% of all amaranth genomes. Phylogenetic analysis based on 406 single-copy orthologous genes revealed that A. hypochondriacus is most closely linked to A. hybridus and distantly related to A. cruentus. However, dioecious amaranth species, such as A. tuberculatus and A. palmeri, which belong to the subgenera Amaranthus Acnida, have formed their distinct clade. The comparative analysis of genomic data of amaranth species will be useful to identify and characterize agronomically important genes and their mechanisms of action. This will facilitate genomics-based, evolutionary studies, and breeding strategies to design faster, more precise, and predictable crop improvement programs.

Biochemical and Phylogenetic Analysis of Italian Phaseolus vulgaris Cultivars as Sources of α-Amylase and α-Glucosidase Inhibitors

Phaseolus vulgaris α-amylase inhibitor (α-AI) is a protein that has recently gained commercial interest, as it inhibits mammalian α-amylase activity, reducing the absorption of dietary carbohydrates. Numerous studies have reported the efficacy of preparations based on this protein on the control of glycaemic peaks in type-2 diabetes patients and in overweight subjects. A positive influence on microbiota regulation has also been described. In this work, ten insufficiently studied Italian P. vulgaris cultivars were screened for α-amylase- and α-glucosidase-inhibiting activity, as well as for the absence of antinutritional compounds, such as phytohemagglutinin (PHA). All the cultivars presented α-glucosidase-inhibitor activity, while α-AI was missing in two of them. Only the Nieddone cultivar (ACC177) had no haemagglutination activity. In addition, the partial nucleotide sequence of the α-AI gene was identified with the degenerate hybrid oligonucleotide primer (CODEHOP) strategy to identify genetic variability, possibly linked to functional α-AI differences, expression of the α-AI gene, and phylogenetic relationships. Molecular studies showed that α-AI was expressed in all the cultivars, and a close similarity between the Pisu Grogu and Fasolu cultivars' α-AI and α-AI-4 isoform emerged from the comparison of the partially reconstructed primary structures. Moreover, mechanistic models revealed the interaction network that connects α-AI with the α-amylase enzyme characterized by two interaction hotspots (Asp38 and Tyr186), providing some insights for the analysis of the α-AI primary structure from the different cultivars, particularly regarding the structure-activity relationship. This study can broaden the knowledge about this class of proteins, fuelling the valorisation of Italian agronomic biodiversity through the development of commercial preparations from legume cultivars.

Iridoid glycosides from Picrorhiza genus endemic to the Himalayan region: phytochemistry, biosynthesis, pharmacological potential and biotechnological intercessions to boost production

Iridoid glycosides are monoterpenoids synthesized in several plant species known to exhibit a diverse range of pharmacological activities. They are used as important bioactive ingredients in many commercially available drug formulations and as lead compounds in pharmaceutical research. The genus Picrorhiza comprises two medicinally important herbs endemic to the Himalayan region viz. Picrorhiza kurrooa Royle and Picrorhiza scrophulariiflora Hong. The medicinal properties of these two species are mainly due to iridoid glycosides present in their root, rhizome, and leaves. Unregulated harvesting from the wild, habitat specificity, narrow distribution range, small population size and lack of organized cultivation led to the enrolling of these species in the endangered category by the International Union for Conservation of Nature and Natural Resources (IUCN). Therefore, there is a need for immediate biotechnological and molecular interventions. Such intercessions will open up new vistas for large-scale propagation, development of genomic/transcriptomic resources for understanding the biosynthetic pathway, the possibility of genetic/metabolic manipulations, and possible commercialization of iridoid glycosides. The current review article elucidates the phytochemistry and pharmacological importance of iridoid glycosides from the genus Picrorhiza. In addition, the role of biotechnological approaches and opportunities offered by next-generation sequencing technologies in overcoming challenges associated with the genetic engineering of these species are also discussed.

In silico Identification of miRNAs and Their Targets in Cluster Bean for Their Role in Development and Physiological Responses

Cluster bean popularly known as “guar” is a drought-tolerant, annual legume that has recently emerged as an economically important crop, owing to its high protein and gum content. The guar gum has wide range of applications in food, pharma, and mining industries. India is the leading exporter of various cluster bean-based products all across the globe. Non-coding RNAs (miRNAs) are involved in regulating the expression of the target genes leading to variations in the associated pathways or final protein concentrations. The understanding of miRNAs and their associated targets in cluster bean is yet to be used to its full potential. In the present study, cluster bean EST (Expressed Sequence Tags) database was exploited to identify the miRNA and their predicted targets associated with metabolic and biological processes especially response to diverse biotic and abiotic stimuli using in silico approach. Computational analysis based on cluster bean ESTs led to the identification of 57 miRNAs along with their targets. To the best of our knowledge, this is the first report on identification of miRNAs and their targets using ESTs in cluster bean. The miRNA related to gum metabolism was also identified. Most abundant miRNA families predicted in our study were miR156, miR172, and miR2606. The length of most of the mature miRNAs was found to be 21nt long and the range of minimal folding energy (MFE) was 5.8–177.3 (−kcal/mol) with an average value of 25.4 (−kcal/mol). The identification of cluster bean miRNAs and their targets is predicted to hasten the miRNA discovery, resulting in better knowledge of the role of miRNAs in cluster bean development, physiology, and stress responses.

Identification and expression analysis of ten novel small non-coding RNAs (sncRNAs) in cancer cells using a high-throughput sequencing approach

Non-coding RNAs are characterized as RNA molecules, which lack the capacity to encode protein structures and appear to include a level of internal signals. Moreover, they control various stages of gene expression, thus controlling the cell physiology and development. In this study, we implemented a high-throughput sequencing approach based on the primary semi-conductor technology and computational tools, in order to identity novel small non-coding RNAs. Fourteen human cancer cell lines were cultured, and RNA samples were enriched for small RNAs following semi-conductor next generation sequencing (NGS). Bioinformatics analysis of NGS data revealed the existence of several classes of ncRNAs using the miRDeep* and CPSS 2.0 software. To investigate the existence of the predicted non-coding RNA sequences in cDNA pools of cell lines, a developed qPCR-based assay was implemented. The structure of each novel small ncRNA was visualized, using the RNAfold algorithm. Our results support the existence of twenty (20) putative new small ncRNAs, ten (10) of which have had their expression experimentally validated and presented differential profiles in cancerous and normal cells. A deeper comprehension of the ncRNAs interactive network and its role in cancer can therefore be translated into a wide range of clinical applications. Despite this progress, further scientific research from different perspectives and in different fields is needed, so that the riddle of the human transcriptome can be solved.

Genome-wide miRNAs profiles of pearl millet under contrasting high vapor pressure deficit reveal their functional roles in drought stress adaptations

Pearl millet (Pennisetum glaucum [L.] R. Br.) is an important crop capable of growing in harsh and marginal environments, with the highest degree of tolerance to drought and heat stresses among cereals. Diverse germplasm of pearl millet shows a significant phenotypic variation in response to abiotic stresses, making it a unique model to study the mechanisms responsible for stress mitigation. The present study focuses on identifying the physiological response of two pearl millet high-resolution cross (HRC) genotypes, ICMR 1122 and ICMR 1152, in response to low and high vapor pressure deficit (VPD). Under high VPD conditions, ICMR 1152 exhibited a lower transpiration rate (Tr), higher transpiration efficiency, and lower root sap exudation than ICMR 1122. Further, Pg-miRNAs expressed in the contrasting genotypes under low and high VPD conditions were identified by deep sequencing analysis. A total of 116 known and 61 novel Pg-miRNAs were identified from ICMR 1152, while 26 known and six novel Pg-miRNAs were identified from ICMR 1122 genotypes, respectively. While Pg-miR165, 168, 170, and 319 families exhibited significant differential expression under low and high VPD conditions in both genotypes, ICMR 1152 showed abundant expression of Pg-miR167, Pg-miR172, Pg-miR396 Pg-miR399, Pg-miR862, Pg-miR868, Pg-miR950, Pg-miR5054, and Pg-miR7527 indicating their direct and indirect role in root physiology and abiotic stress responses. Drought responsive Pg-miRNA targets showed upregulation in response to high VPD stress, further narrowing down the miRNAs involved in regulation of drought tolerance in pearl millet.

Epigenetics and its role in effecting agronomical traits

Climate-resilient crops with improved adaptation to the changing climate are urgently needed to feed the growing population. Hence, developing high-yielding crop varieties with better agronomic traits is one of the most critical issues in agricultural research. These are vital to enhancing yield as well as resistance to harsh conditions, both of which help farmers over time. The majority of agronomic traits are quantitative and are subject to intricate genetic control, thereby obstructing crop improvement. Plant epibreeding is the utilisation of epigenetic variation for crop development, and has a wide range of applications in the field of crop improvement. Epigenetics refers to changes in gene expression that are heritable and induced by methylation of DNA, post-translational modifications of histones or RNA interference rather than an alteration in the underlying sequence of DNA. The epigenetic modifications influence gene expression by changing the state of chromatin, which underpins plant growth and dictates phenotypic responsiveness for extrinsic and intrinsic inputs. Epigenetic modifications, in addition to DNA sequence variation, improve breeding by giving useful markers. Also, it takes epigenome diversity into account to predict plant performance and increase crop production. In this review, emphasis has been given for summarising the role of epigenetic changes in epibreeding for crop improvement.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

Barley Seeds miRNome Stability during Long-Term Storage and Aging

Seed aging is a complex biological process that has been attracting scientists’ attention for many years. High-throughput small RNA sequencing was applied to examine microRNAs contribution in barley seeds senescence. Unique samples of seeds that, despite having the same genetic makeup, differed in viability after over 45 years of storage in a dry state were investigated. In total, 61 known and 81 novel miRNA were identified in dry seeds. The highest level of expression was found in four conserved miRNA families, i.e., miR159, miR156, miR166, and miR168. However, the most astonishing result was the lack of significant differences in the level of almost all miRNAs in seed samples with significantly different viability. This result reveals that miRNAs in dry seeds are extremely stable. This is also the first identified RNA fraction that is not deteriorating along with the loss of seed viability. Moreover, the novel miRNA hvu-new41, with higher expression in seeds with the lowest viability as detected by RT-qPCR, has the potential to become an indicator of the decreasing viability of seeds during storage in a dry state.

Phaseolus vulgaris MIR1511 genotypic variations differentially regulate plant tolerance to aluminum toxicity

The common-bean (Phaseolus vulgaris), a widely consumed legume, originated in Mesoamerica and expanded to South America, resulting in the development of two geographically distinct gene pools. Poor soil condition, including metal toxicity, are often constraints to common-bean crop production. Several P. vulgaris miRNAs, including miR1511, respond to metal toxicity. The MIR1511 gene sequence from the two P. vulgaris model sequenced genotypes revealed that, as opposed to BAT93 (Mesoamerican), the G19833 (Andean) accession displays a 58-bp deletion, comprising the mature and star miR1511 sequences. Genotyping-By-Sequencing data analysis from 87 non-admixed Phaseolus genotypes, comprising different Phaseolus species and P. vulgaris populations, revealed that all the P. vulgaris Andean genotypes and part of the Mesoamerican (MW1) genotypes analyzed displayed a truncated MIR1511 gene. The geographic origin of genotypes with a complete versus truncated MIR1511 showed a distinct distribution. The P. vulgaris ALS3 (Aluminum Sensitive Protein 3) gene, known to be important for aluminum detoxification in several plants, was experimentally validated as the miR1511 target. Roots from BAT93 plants showed decreased miR1511 and increased ALS3 transcript levels at early stages under aluminum toxicity (AlT), while G19833 plants, lacking mature miR1511, showed higher and earlier ALS3 response. Root architecture analyses evidenced higher tolerance of G19833 plants to AlT. However, G19833 plants engineered for miR1511 overexpression showed lower ALS3 transcript level and increased sensitivity to AlT. Absence of miR1511 in Andean genotypes, resulting in a diminished ALS3 transcript degradation, appears to be an evolutionary advantage to high Al levels in soils with increased drought conditions.

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

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

MicroRNAs expression dynamics reveal post-transcriptional mechanisms regulating seed development in Phaseolus vulgaris L

The knowledge on post-transcriptional regulation mechanisms implicated in seed development (SD) is still limited, particularly in one of the most consumed grain legumes, Phaseolus vulgaris L. We explore for the first time the miRNA expression dynamics in P. vulgaris developing seeds. Seventy-two known and 39 new miRNAs were found expressed in P. vulgaris developing seeds. Most of the miRNAs identified were more abundant at 10 and 40 days after anthesis, suggesting that late embryogenesis/early filling and desiccation were SD stages in which miRNA action is more pronounced. Degradome analysis and target prediction identified targets for 77 expressed miRNAs. While several known miRNAs were predicted to target HD-ZIP, ARF, SPL, and NF-Y transcription factors families, most of the predicted targets for new miRNAs encode for functional proteins. MiRNAs-targets expression profiles evidenced that these miRNAs could tune distinct seed developmental stages. MiRNAs more accumulated at early SD stages were implicated in regulating the end of embryogenesis, postponing the seed maturation program, storage compound synthesis and allocation. MiRNAs more accumulated at late SD stages could be implicated in seed quiescence, desiccation tolerance, and longevity with still uncovered roles in germination. The miRNAs herein described represent novel P. vulgaris resources with potential application in future biotechnological approaches to modulate the expression of genes implicated in legume seed traits with impact in horticultural production systems.

Identification of Small RNAs During High Light Acclimation in Arabidopsis thaliana

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

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

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

Comprehensive Identification and Expression Profiling of Circular RNAs During Nodule Development in Phaseolus vulgaris

Symbiotic nitrogen fixation by legume nodules provides an abundant nitrogen source for plants, and understanding this process is key for developing green agriculture. Circular RNA (circRNA), a type of endogenous RNA produced by reverse splicing of mRNA precursors, plays important regulatory roles in plants at the transcriptional and post-transcriptional levels. However, the relationship between circRNAs and legume-rhizobium is unknown. Here, we performed comprehensive identification and expression profiling of circRNAs during nodulation in common bean (Phaseolus vulgaris) compared to uninoculated roots of corresponding ages by constructing circRNA-seq and mRNA-seq libraries. We identified 8,842 high-confident circRNAs, 3,448 of which were specifically produced during symbiosis, with the highest number at the nitrogen-fixing stage. Significantly, more circRNAs were derived from exons than from intergenic regions or introns in all samples. The lengths and GC contents of the circRNAs were similar in roots and nodules. However, circRNAs showed specific spatiotemporal expression patterns during nodule and root development. GO and other functional annotation of parental genes of differentially expressed circRNAs indicated their potential involvement in different biological processes. The expression of major circRNAs during symbiosis is independent of parental genes' expression to a certain degree, while expression of the remaining minor circRNAs showed positive correlation to parental genes. Functional annotation of the targeted mRNAs in the circRNA-miRNA-mRNA network showed that circRNAs may be involved in transmembrane transport and positive regulation of kinase activity during nodulation and nitrogen fixation as miRNA sponges. Our comprehensive analysis of the expression profile of circRNAs and their potential functions suggests that circRNAs may function as new post-transcriptional regulators in legume-rhizobium symbiosis.

Origin and Evolutionary Dynamics of the miR2119 and ADH1 Regulatory Module in Legumes

MicroRNAs are important regulators of gene expression in eukaryotes. Previously, we reported that in Phaseolus vulgaris, the precursor for miR2119 is located in the same gene as miR398a, conceiving a dicistronic MIR gene. Both miRNA precursors are transcribed and processed from a single transcript resulting in two mature microRNAs that regulate the mRNAs encoding ALCOHOL DEHYDROGENASE 1 (ADH1) and COPPER-ZINC SUPEROXIDE DISMUTASE 1 (CSD1). Genes for miR398 are distributed throughout the spermatophytes; however, miR2119 is only found in Leguminosae species, indicating its recent emergence. Here, we used public databases to explore the presence of the miR2119 sequence in several plant species. We found that miR2119 is present only in specific clades within the Papilionoideae subfamily, including important crops used for human consumption and forage. Within this subfamily, MIR2119 and MIR398a are found together as a single gene in the genomes of the Millettioids and Hologalegina. In contrast, in the Dalbergioids MIR2119 is located in a different locus from MIR398a, suggesting this as the ancestral genomic organization. To our knowledge, this is a unique example where two separate MIRNA genes have merged to generate a single polycistronic gene. Phylogenetic analysis of ADH1 gene sequences in the Papilionoideae subfamily revealed duplication events resulting in up to four ADH1 genes in certain species. Notably, the presence of MIR2119 correlates with the conservation of target sites in particular ADH1 genes in each clade. Our results suggest that post-transcriptional regulation of ADH1 genes by miR2119 has contributed to shaping the expansion and divergence of this gene family in the Papilionoideae. Future experimental work on ADH1 regulation by miR2119 in more legume species will help to further understand the evolutionary history of the ADH1 gene family and the relevance of miRNA regulation in this process.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana

BACKGROUND

Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes.

RESULT

We subjected A. thaliana plants to cold acclimation conditions (4 °C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cis- and trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation.

CONCLUSION

In A. thaliana, miRNAs and sRNAs derived from cis- and trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.

Gene Expression in Nitrogen-Fixing Symbiotic Nodule Cells in Medicago truncatula and Other Nodulating Plants

Root nodules formed by plants of the nitrogen-fixing clade (NFC) are symbiotic organs that function in the maintenance and metabolic integration of large populations of nitrogen-fixing bacteria. These organs feature unique characteristics and processes, including their tissue organization, the presence of specific infection structures called infection threads, endocytotic uptake of bacteria, symbiotic cells carrying thousands of intracellular bacteria without signs of immune responses, and the integration of symbiont and host metabolism. The early stages of nodulation are governed by a few well-defined functions, which together constitute the common symbiosis-signaling pathway (CSSP). The CSSP activates a set of transcription factors (TFs) that orchestrate nodule organogenesis and infection. The later stages of nodule development require the activation of hundreds to thousands of genes, mostly expressed in symbiotic cells. Many of these genes are only active in symbiotic cells, reflecting the unique nature of nodules as plant structures. Although how the nodule-specific transcriptome is activated and connected to early CSSP-signaling is poorly understood, candidate TFs have been identified using transcriptomic approaches, and the importance of epigenetic and chromatin-based regulation has been demonstrated. We discuss how gene regulation analyses have advanced our understanding of nodule organogenesis, the functioning of symbiotic cells, and the evolution of symbiosis in the NFC.

Profiling miRNA expression in photo-thermo-sensitive male genic sterility line (PTGMS) PA64S under high and low temperature

Peiai64S (PA64S) is a photo-thermo-sensitive genic male sterile line (PTGMS), with wide application in hybrid seed production in rice (Oryza sativa L.). Micro-RNAs are 21-24 nt, endogenously expressed small RNAs that have been characterized in various developmental stages of rice, but none have been studied with respect to the regulation of TGMS in rice. Here, we employed high-throughput sequencing to identify expression profiles of miRNAs in the anthers of PA64S at high (PA64S-H) and low temperature (PA64S-L). Two small RNA libraries from PA64S-H and PA64-L anthers were sequenced, and 263 known and 321 novel candidate miRNAs were identified. Based on the number of sequencing reads, a total of 133 known miRNAs were found to be differentially expressed between PA64S-H and PA64S-L. Target prediction showed that the target genes encode MYB and TCP transcription factors, and bHLH proteins. These target genes are related to pollen development and male sterility, suggesting that miRNA/targets may play roles in regulating TGMS in rice. Further, starch and sucrose metabolism pathways, sphingolipid metabolism, arginine and proline metabolism, and plant hormone signal transduction pathways were enriched by KEGG pathway annotation. These findings contribute to our understanding of the role of miRNAs during anther development and TGMS occurrence in rice.

The role of microRNAs in the healing of diabetic ulcers

MicroRNAs (miRNAs) are small protected molecules with a length of 18 to 25 nucleotides. Many studies have recently been conducted on miRNAs, illustrating their role in regulating many biological, physiological, and pathological activities, such as maintaining cellular signalling and regulating cellular pathways. The main role of miRNAs is to regulate the expression of genes after translation, which can lead to the destruction or suppression of translation by binding to mRNAs. As any change in the regulation of miRNAs is associated with several physiological abnormalities, such as type 2 diabetes and its complications, these molecules can be used for therapeutic purposes or as biomarkers for the diagnosis of diseases such as diabetes and its complications. In this review article, we will discuss important findings about the miRNAs and the role of these molecules in different phases of the wound-healing process of chronic wounds, especially diabetic ulcer.

Genome-wide association analysis of the lipid and fatty acid metabolism regulatory network in the mesocarp of oil palm (Elaeis guineensis Jacq.) based on small noncoding RNA sequencing

Oil palm (Elaeis guineensis Jacq.) is the highest oil-yielding crop in the plant kingdom and accumulates 90% of palm oil in the mesocarp. However, the regulatory mechanisms of lipid and fatty acid (FA) metabolism in oil palm are just beginning to be understood, and more studies are needed, especially in the understanding of small noncoding RNA (ncRNA) and mRNA. Based on the deep sequencing of small noncoding RNAs and the degradome in five developmental mesocarp stages, 452 microRNAs (miRNAs), including 170 conserved known-miRNAs (kn-miRNAs) and 282 novel-miRNA (nov-miRNAs), were identified. After predicting the targets of those miRNAs to 37 FA synthesis-related genes, we found that 22 kn-miRNAs and 14 nov-miRNAs might be involved in FA metabolism pathways. Among them, eg-miR156c, eg-miR397, eg-miR444b and nov-miR129 regulated FA synthesis in plastids and the transport of FA-ACP from plastids to the endoplasmic reticulum by targeting acetyl-CoA carboxylase 1 (ACC1), long-chain acyl-CoA synthetase 9 (LACS9), LACS4 and enoyl-ACP reductase (ENR), respectively. Nov-miR138 and nov-miR59 targeted glycerol-3-phosphate acyltransferase (GPAT), and nov-miR274 targeted phosphatidate phosphatase 1 (PAP1). Both target genes are involved in triacylglycerol synthesis in the endoplasmic reticulum. Eg-miR156e and eg-miR156j played pivotal roles by targeting β-ketoacyl-CoA synthase 12 (KCS12), and nov-miR201 targets very-long-chain enoyl-CoA reductase (ECR). Several miRNAs were also predicted to indirectly regulate FA synthesis and lipid metabolism through the squamosa promoter-binding protein-like gene (SPL), NAC and MYB transcription factors. As a whole, indications of a complex and extensive miRNA-mRNA regulatory network associated with FA metabolism in the mesocarp of the oil palm is presented. The results help to broaden the knowledge of potential mechanisms that might be regulated by miRNAs through modulation of the expression of FA-related target gene metabolism in the oil palm.

Identification and characterization of differentially expressed Phaseolus vulgaris miRNAs and their targets during mungbean yellow mosaic India virus infection reveals new insight into Phaseolus-MYMIV interaction

Phaseolus vulgaris is an economically important legume in tropical and subtropical regions of Asia, Africa, Latin-America and parts of USA and Europe. However, its production gets severely affected by mungbean yellow mosaic India virus (MYMIV). We aim to identify and characterize differentially expressed miRNAs during MYMIV-infection in P. vulgaris. A total of 422 miRNAs are identified of which 292 are expressed in both MYMIV-treated and mock-treated samples, 109 are expressed only in MYMIV-treated and 21 are expressed only in mock-treated samples. Selected up- and down-regulated miRNAs are validated by RT-qPCR. 3367 target ORFs are identified for 270 miRNAs. Selected targets are validated by 5' RLM-RACE. Differentially expressed miRNAs regulate transcription factors and are involved in improving stress tolerance to MYMIV. These findings will provide an insight into the role of miRNAs during MYMIV infection in P. vulgaris in particular and during any biotic stress conditions in Leguminosae family in general.

Identification and expression analysis of microRNAs during ovule development in rice (Oryza sativa) by deep sequencing

MicroRNA (miRNA) expression profiles during rice ovule development revealed the possible miRNA-mediated regulation between ovule sporophytic tissue and female gametophyte and the involvement of miRNAs in programmed cell death. MiRNAs are 20-24-nucleotide small RNAs that play key roles in the regulation of many growth and developmental processes in plants. Rice ovule development comprises a series of biological events, which are regulated by complex molecular mechanisms. To gain insight into miRNA-mediated regulation of rice ovule development, Illumina sequencing was used to examine the expression of miRNAs from the megaspore mother cell meiosis stage to the fertilized ovule stage. Based on the sequencing data, 486 known and 204 novel miRNAs were identified during rice ovule development. Moreover, 56, 65 and 11 differentially expressed miRNAs between adjacent developmental stages were identified. By analyzing transcriptome and degradome data, we identified 41, 65 and 12 coherent target genes for the differentially expressed miRNAs in ovule development. We found that changes in the expression of plant hormone-related miRNAs may play important roles in embryo sac development, providing evidence for cross-talk communication between sporophytic tissue and the female gametophyte. Additionally, we revealed that miRNAs may be involved in programmed cell death after fertilization. Finally, we constructed miRNA-mediated regulatory networks that are active during rice ovule development.

microRNAs participate in gene expression regulation and phytohormone cross-talk in barley embryo during seed development and germination

BACKGROUND

Small RNA and degradome sequencing have identified a large number of miRNA-target pairs in plant seeds. However, detailed spatial and temporal studies of miRNA-mediated regulation, which can reflect links between seed development and germination are still lacking.

RESULTS

In this study, we extended our investigation on miRNAs-involved gene regulation by a combined analysis of seed maturation and germination in barley. Through bioinformatics analysis of small RNA sequencing data, a total of 1324 known miRNA families and 448 novel miRNA candidates were identified. Of those, 16 known miRNAs with 40 target genes, and three novel miRNAs with four target genes were confirmed based on degradome sequencing data. Conserved miRNA families such as miR156, miR168, miR166, miR167, and miR894 were highly expressed in embryos of developing and germinating seeds. A barley-specific miRNA, miR5071, which was predicted to target an OsMLA10-like gene, accumulated at a high level, suggesting its involvement in defence response during these two developmental stages. Based on target prediction and Kyoto Encyclopedia of Genes and Genomes analysis of putative targets, nine highly expressed miRNAs were found to be related to phytohormone signalling and hormone cross-talk. Northern blot and qRT-PCR analysis showed that these miRNAs displayed differential expression patterns during seed development and germination, indicating their different roles in hormone signalling pathways. In addition, we showed that miR393 affected seed development through targeting two genes encoding the auxin receptors TIR1/AFBs in barley, as over-expression of miR393 led to an increased length-width ratio of seeds, whereas target mimic (MIM393)-mediated inhibition of its activity decreased the 1000-grain weight of seeds. Furthermore, the expression of auxin-responsive genes, abscisic acid- and gibberellic acid-related genes was altered in miR393 misexpression lines during germination and early seedling growth.

CONCLUSIONS

Our work indicates that miRNA-target pairs participate in gene expression regulation and hormone interaction in barley embryo and provides evidence that miR393-mediated auxin response regulation affects grain development and influences gibberellic acid and abscisic acid homeostasis during germination.

MicroRNAs associated with drought response in the pulse crop common bean (Phaseolus vulgaris L.)

Drought stress significantly reduces common bean yields. Recently, some drought-related miRNAs were found in various plants. However, reports of miRNAs involved in drought stress in common bean are limited. Here, we obtained four sRNA samples from drought-tolerant and -sensitive cultivars of common bean that experienced with or without drought treatment. A total of 49 novel miRNAs and 120 known miRNAs were detected. Under drought treatment, 9 and 7 known miRNAs were down and up-regulated, respectively, and 5 and 3 of the novel miRNAs were increased and decreased, respectively. Among these miRNAs, four miRNAs shared the same pattern of expression between Long 22-0579 and Naihua. Target genes of these miRNAs included transcription factors, protein kinases, and nuclear transcription factors. Finally, we verified all of the differentially expressed miRNAs by RT-qPCR, and we identified 16 miRNAs that are potentially associated with the drought stress response. These miRNAs and target genes will be useful in future basic studies and in applied studies investigating how miRNA regulation can be used to enhance drought resistance in plant species.

The legume miR1514a modulates a NAC transcription factor transcript to trigger phasiRNA formation in response to drought

Recent studies have identified microRNAs as post-transcriptional regulators involved in stress responses. miR1514a is a legume microRNA that is induced in response to drought stress in Phaseolus vulgaris (common bean) and shows differential accumulation levels in roots during water deficit in two cultivars with different drought tolerance phenotypes. A recent degradome analysis revealed that miR1514a targets the transcripts of two NAC transcription factors (TFs), Phvul.010g121000 and Phvul.010g120700. Furthermore, expression studies and small RNA-seq data indicate that only Phvul.010g120700 generates phasiRNAs, which also accumulate under water deficit conditions. To confirm these results, we over-expressed miR1514a in transgenic hairy roots, and observed a reduced accumulation of Phvul.010g120700 and an increase in NAC-derived phasiRNAs; inhibition of miR1514a activity resulted in the opposite effect. Moreover, we determined that a NAC-derived phasiRNA associates with ARGONAUTE 1 (AGO1), suggesting that it is functional. In addition, a transcriptome analysis of transgenic hairy roots with reduced miR1514a levels revealed several differentially expressed transcripts, mainly involved in metabolism and stress responses, suggesting they are regulated by the NAC TF and/or by phasiRNAs. This work therefore demonstrates the participation of miR1514 in the regulation of a NAC transcription factor transcript through phasiRNA production during the plant response to water deficit.

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.

Small RNAs Derived from the T-DNA of Agrobacterium rhizogenes in Hairy Roots of Phaseolus vulgaris

Agrobacterium rhizogenes is a pathogenic bacteria that causes hairy root disease by transferring bacterial DNA into the plant genome. It is an essential tool for industry and research due to its capacity to produce genetically modified roots and whole organisms. Here, we identified and characterized small RNAs generated from the transfer DNA (T-DNA) of A. rhizogenes in hairy roots of common bean (Phaseolus vulgaris). Distinct abundant A. rhizogenes T-DNA-derived small RNAs (ArT-sRNAs) belonging to several oncogenes were detected in hairy roots using high-throughput sequencing. The most abundant and diverse species of ArT-sRNAs were those of 21- and 22-nucleotides in length. Many T-DNA encoded genes constituted phasiRNA producing loci (PHAS loci). Interestingly, degradome analysis revealed that ArT-sRNAs potentially target genes of P. vulgaris. In addition, we detected low levels of ArT-sRNAs in the A. rhizogenes-induced calli generated at the wound site before hairy root emergence. These results suggest that RNA silencing targets several genes from T-DNA of A. rhizogenes in hairy roots of common bean. Therefore, the role of RNA silencing observed in this study has implications in our understanding and usage of this unique plant-bacteria interaction.

Gene Expression Characteristics and Regulation Mechanisms of Superoxide Dismutase and Its Physiological Roles in Plants under Stress

Superoxide dismutases (SODs) are key enzymes functioning as the first line of antioxidant defense by virtue of the ability to convert highly reactive superoxide radicals to hydrogen peroxide and molecular oxygen. SOD plays a central role in protecting plants against the toxic effects of reactive oxygen species generated during normal cellular metabolic activity or as a result of various environmental stresses. Our review focuses on the characteristics of expression of SOD genes, the mechanisms regulating expression of SOD genes at transcriptional, posttranscriptional, and translation levels, and their functional role(s) during development and in response to biotic or abiotic stresses. We propose two important research directions: studying SOD at the genome-wide or proteome-wide level, and improving plant stress tolerances by selecting varieties using transgenic technology.

Comparison of Small RNA Profiles of Glycine max and Glycine soja at Early Developmental Stages

Small RNAs, including microRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs; from PHAS loci), play key roles in plant development. Cultivated soybean, Glycine max, contributes a great deal to food production, but, compared to its wild kin, Glycine soja, it may lose some genetic information during domestication. In this work, we analyzed the sRNA profiles of different tissues in both cultivated (C08) and wild soybeans (W05) at three stages of development. A total of 443 known miRNAs and 15 novel miRNAs showed varying abundances between different samples, but the miRNA profiles were generally similar in both accessions. Based on a sliding window analysis workflow that we developed, 50 PHAS loci generating 55 21-nucleotide phasiRNAs were identified in C08, and 46 phasiRNAs from 41 PHAS loci were identified in W05. In germinated seedlings, phasiRNAs were more abundant in C08 than in W05. Disease resistant TIR-NB-LRR genes constitute a very large family of PHAS loci. PhasiRNAs were also generated from several loci that encode for NAC transcription factors, Dicer-like 2 (DCL2), Pentatricopeptide Repeat (PPR), and Auxin Signaling F-box 3 (AFB3) proteins. To investigate the possible involvement of miRNAs in initiating the PHAS-phasiRNA pathway, miRNA target predictions were performed and 17 C08 miRNAs and 15 W05 miRNAs were predicted to trigger phasiRNAs biogenesis. In summary, we provide a comprehensive description of the sRNA profiles of wild versus cultivated soybeans, and discuss the possible roles of sRNAs during soybean germination.

Identification and characterization of microRNAs expressed in the African malaria vector Anopheles funestus life stages using high throughput sequencing

BACKGROUND

Over the past several years, thousands of microRNAs (miRNAs) have been identified in the genomes of various insects through cloning and sequencing or even by computational prediction. However, the number of miRNAs identified in anopheline species is low and little is known about their role. The mosquito Anopheles funestus is one of the dominant malaria vectors in Africa, which infects and kills millions of people every year. Therefore, small RNA molecules isolated from the four life stages (eggs, larvae, pupae and unfed adult females) of An. funestus were sequenced using next generation sequencing technology.

RESULTS

High throughput sequencing of four replicates in combination with computational analysis identified 107 mature miRNA sequences expressed in the An. funestus mosquito. These include 20 novel miRNAs without sequence identity in any organism and eight miRNAs not previously reported in the Anopheles genus but are known in non-anopheles mosquitoes. Finally, the changes in the expression of miRNAs during the mosquito development were determined and the analysis showed that many miRNAs have stage-specific expression, and are co-transcribed and co-regulated during development.

CONCLUSIONS

This study presents the first direct experimental evidence of miRNAs in An. funestus and the first profiling study of miRNA associated with the maturation in this mosquito. Overall, the results indicate that miRNAs play important roles during the growth and development. Silencing such molecules in a specific life stage could decrease the vector population and therefore interrupt malaria transmission.

Identification and characterization of microRNAs in maize endosperm response to exogenous sucrose using small RNA sequencing

Sucrose acts as a signaling molecule for genes critical to starch biosynthesis in maize endosperm. Previously, we showed that sucrose could regulate starch biosynthesis in maize via transcription factors. To better understand the complex regulation of starch biosynthesis, the 10days after pollination endosperm from Zea mays L. B73 inbred line was collected and treated with sucrose for small RNA sequencing. The sequencing results revealed that 24 known miRNAs and 190 novel miRNAs were significantly differentially expressed in response to sucrose. In addition, most of target mRNAs were characterized as transcription factors, mainly including, MYB, ARF, NAC, AP2/ERF, WRKY, and GRAS, which play important roles in starch biosynthesis and seed development in maize endosperm. The expression profiles of miR398a/b and miR159b/j/k followed opposite expression trends to their target genes when analyzed by qPCR. In conclusion, these results show that sucrose regulates the expression of starch synthetic genes through miRNAs.

The Class II Trehalose 6-phosphate Synthase Gene PvTPS9 Modulates Trehalose Metabolism in Phaseolus vulgaris Nodules

Legumes form symbioses with rhizobia, producing nitrogen-fixing nodules on the roots of the plant host. The network of plant signaling pathways affecting carbon metabolism may determine the final number of nodules. The trehalose biosynthetic pathway regulates carbon metabolism and plays a fundamental role in plant growth and development, as well as in plant-microbe interactions. The expression of genes for trehalose synthesis during nodule development suggests that this metabolite may play a role in legume-rhizobia symbiosis. In this work, PvTPS9, which encodes a Class II trehalose-6-phosphate synthase (TPS) of common bean (Phaseolus vulgaris), was silenced by RNA interference in transgenic nodules. The silencing of PvTPS9 in root nodules resulted in a reduction of 85% (± 1%) of its transcript, which correlated with a 30% decrease in trehalose contents of transgenic nodules and in untransformed leaves. Composite transgenic plants with PvTPS9 silenced in the roots showed no changes in nodule number and nitrogen fixation, but a severe reduction in plant biomass and altered transcript profiles of all Class II TPS genes. Our data suggest that PvTPS9 plays a key role in modulating trehalose metabolism in the symbiotic nodule and, therefore, in the whole plant.

The First Report of miRNAs from a Thysanopteran Insect, Thrips palmi Karny Using High-Throughput Sequencing

Thrips palmi Karny (Thysanoptera: Thripidae) is the sole vector of Watermelon bud necrosis tospovirus, where the crop loss has been estimated to be around USD 50 million annually. Chemical insecticides are of limited use in the management of T. palmi due to the thigmokinetic behaviour and development of high levels of resistance to insecticides. There is an urgent need to find out an effective futuristic management strategy, where the small RNAs especially microRNAs hold great promise as a key player in the growth and development. miRNAs are a class of short non-coding RNAs involved in regulation of gene expression either by mRNA cleavage or by translational repression. We identified and characterized a total of 77 miRNAs from T. palmi using high-throughput deep sequencing. Functional classifications of the targets for these miRNAs revealed that majority of them are involved in the regulation of transcription and translation, nucleotide binding and signal transduction. We have also validated few of these miRNAs employing stem-loop RT-PCR, qRT-PCR and Northern blot. The present study not only provides an in-depth understanding of the biological and physiological roles of miRNAs in governing gene expression but may also lead as an invaluable tool for the management of thysanopteran insects in the future.

Deep-sequence profiling of miRNAs and their target prediction in Monotropa hypopitys

Myco-heterotroph Monotropa hypopitys is a widely spread perennial herb used to study symbiotic interactions and physiological mechanisms underlying the development of non-photosynthetic plant. Here, we performed, for the first time, transcriptome-wide characterization of M. hypopitys miRNA profile using high throughput Illumina sequencing. As a result of small RNA library sequencing and bioinformatic analysis, we identified 55 members belonging to 40 families of known miRNAs and 17 putative novel miRNAs unique for M. hypopitys. Computational screening revealed 206 potential mRNA targets for known miRNAs and 31 potential mRNA targets for novel miRNAs. The predicted target genes were described in Gene Ontology terms and were found to be involved in a broad range of metabolic and regulatory pathways. The identification of novel M. hypopitys-specific miRNAs, some with few target genes and low abundances, suggests their recent evolutionary origin and participation in highly specialized regulatory mechanisms fundamental for non-photosynthetic biology of M. hypopitys. This global analysis of miRNAs and their potential targets in M. hypopitys provides a framework for further investigation of miRNA role in the evolution and establishment of non-photosynthetic myco-heterotrophs.

Identification of miRNAs and their targets through high-throughput sequencing and degradome analysis in male and female Asparagus officinalis

BACKGROUND

MicroRNAs (miRNAs), a class of non-coding small RNAs (sRNAs), regulate various biological processes. Although miRNAs have been identified and characterized in several plant species, miRNAs in Asparagus officinalis have not been reported. As a dioecious plant with homomorphic sex chromosomes, asparagus is regarded as an important model system for studying mechanisms of plant sex determination.

RESULTS

Two independent sRNA libraries from male and female asparagus plants were sequenced with Illumina sequencing, thereby generating 4.13 and 5.88 million final clean reads, respectively. Both libraries predominantly contained 24-nt sRNAs, followed by 21-nt sRNAs. Further analysis identified 154 conserved miRNAs, which belong to 26 families, and 39 novel miRNA candidates seemed to be specific to asparagus. Comparative profiling revealed that 63 miRNAs exhibited significant differential expression between male and female plants, which was confirmed by real-time quantitative PCR analysis. Among them, 37 miRNAs were significantly up-regulated in the female library, whereas the others were preferentially expressed in the male library. Furthermore, 40 target mRNAs representing 44 conserved and seven novel miRNAs were identified in asparagus through high-throughput degradome sequencing. Functional annotation showed that these target mRNAs were involved in a wide range of developmental and metabolic processes.

CONCLUSIONS

We identified a large set of conserved and specific miRNAs and compared their expression levels between male and female asparagus plants. Several asparagus miRNAs, which belong to the miR159, miR167, and miR172 families involved in reproductive organ development, were differentially expressed between male and female plants, as well as during flower development. Consistently, several predicted targets of asparagus miRNAs were associated with floral organ development. These findings suggest the potential roles of miRNAs in sex determination and reproductive developmental processes in asparagus.

Genome-wide identification and in silico characterisation of microRNAs, their targets and processing pathway genes in Phaseolus vulgaris L

Common bean (Phaseolus vulgaris L., Fabaceae) is a globally important staple crop, which is an important source of calories, protein and essential micronutrients. At the genomic level little is known regarding the small non-coding RNAs within the common bean genome. One of the most important classes of such small non-coding RNAs is microRNAs (miRNAs), which control mRNA and protein expression levels in many eukaryotes. Computational methods have been applied to identify putative miRNAs in the genomes of different organisms. In this study, our objective was to comprehensively identify and characterise miRNAs from the genome and transcriptome of P. vulgaris, including both mature and precursor miRNA forms. We also sought to identify the putative proteins involved in miRNA processing and the likely target genes of common bean miRNAs. We identified 221 mature miRNAs and 136 precursor miRNAs distributed across 52 different miRNA families in the P. vulgaris genome. Amongst these, we distinguished 129 novel mature miRNAs and 123 miRNA precursors belonging to 24 different miRNA families. We also identified 31 proteins predicted to participate in the miRNA-processing pathway in P. vulgaris. Finally, we also identified 483 predicted miRNA targets, including many which corroborate results from other species, suggesting that miRNA regulatory systems are evolutionarily conserved and important for plant development. Our results expand the study of miRNAs and their target genes in common bean, and provide new opportunities to understand their roles in the biology of this important staple crop.

Identification of microRNAs and their mRNA targets during soybean nodule development: functional analysis of the role of miR393j-3p in soybean nodulation

Plant microRNAs (miRNAs) play important regulatory roles in a number of developmental processes. The present work investigated the roles of miRNAs during nodule development in the crop legume soybean (Glycine max). Fifteen soybean small RNA libraries were sequenced from different stages of nodule development, including young nodules, mature nodules and senescent nodules. In order to identify the regulatory targets of the miRNAs, five parallel analysis of RNA ends (PARE) libraries were also sequenced from the same stages of nodule development. Sequencing identified 284 miRNAs, including 178 novel soybean miRNAs. Analysis of miRNA abundance identified 139 miRNAs whose expression was significantly regulated during nodule development, including 12 miRNAs whose expression changed > 10-fold. Analysis of the PARE libraries identified 533 miRNA targets, including three nodulation-related genes and eight nodule-specific genes. miR393j-3p was selected for detailed analysis as its expression was significantly regulated during nodule formation, and it targeted a nodulin gene, Early Nodulin 93 (ENOD93). Strong, ectopic expression of miR393j-3p, as well as RNAi silencing of ENOD93 expression, significantly reduced nodule formation. The data indicate that miR393j-3p regulation of ENOD93 mRNA abundance is a key control point for soybean nodule formation.

Responses of symbiotic nitrogen-fixing common bean to aluminum toxicity and delineation of nodule responsive microRNAs

Aluminum (Al) toxicity is widespread in acidic soils where the common bean (Phaseolus vulgaris), the most important legume for human consumption, is produced and it is a limiting factor for crop production and symbiotic nitrogen fixation. We characterized the nodule responses of common bean plants inoculated with Rhizobioum tropici CIAT899 and the root responses of nitrate-fertilized plants exposed to excess Al in low pH, for long or short periods. A 43-50% reduction in nitrogenase activity indicates that Al toxicity (Alt) highly affected nitrogen fixation in common bean. Bean roots and nodules showed characteristic symptoms for Alt. In mature nodules Al accumulation and lipoperoxidation were observed in the infected zone, while callose deposition and cell death occurred mainly in the nodule cortex. Regulatory mechanisms of plant responses to metal toxicity involve microRNAs (miRNAs) along other regulators. Using a miRNA-macroarray hybridization approach we identified 28 (14 up-regulated) Alt nodule-responsive miRNAs. We validated (quantitative reverse transcriptase-PCR) the expression of eight nodule responsive miRNAs in roots and in nodules exposed to high Al for long or short periods. The inverse correlation between the target and miRNA expression ratio (stress:control) was observed in every case. Generally, miRNAs showed a higher earlier response in roots than in nodules. Some of the common bean Alt-responsive miRNAs identified has also been reported as differentially expressed in other plant species subjected to similar stress condition. miRNA/target nodes analyzed in this work are known to be involved in relevant signaling pathways, thus we propose that the participation of miR164/NAC1 (NAM/ATAF/CUC transcription factor) and miR393/TIR1 (TRANSPORT INHIBITOR RESPONSE 1-like protein) in auxin and of miR170/SCL (SCARECROW-like protein transcription factor) in gibberellin signaling is relevant for common bean response/adaptation to Al stress. Our data provide a foundation for evaluating the individual roles of miRNAs in the response of common bean nodules to Alt.

Computational prediction of miRNAs and their targets in Phaseolus vulgaris using simple sequence repeat signatures

BACKGROUND

MicroRNAs (miRNAs) are endogenous, noncoding, short RNAs directly involved in regulating gene expression at the post-transcriptional level. In spite of immense importance, limited information of P. vulgaris miRNAs and their expression patterns prompted us to identify new miRNAs in P. vulgaris by computational methods. Besides conventional approaches, we have used the simple sequence repeat (SSR) signatures as one of the prediction parameter. Moreover, for all other parameters including normalized Shannon entropy, normalized base pairing index and normalized base-pair distance, instead of taking a fixed cut-off value, we have used 99% probability range derived from the available data.

RESULTS

We have identified 208 mature miRNAs in P. vulgaris belonging to 118 families, of which 201 are novel. 97 of the predicted miRNAs in P. vulgaris were validated with the sequencing data obtained from the small RNA sequencing of P. vulgaris. Randomly selected predicted miRNAs were also validated using qRT-PCR. A total of 1305 target sequences were identified for 130 predicted miRNAs. Using 80% sequence identity cut-off, proteins coded by 563 targets were identified. The computational method developed in this study was also validated by predicting 229 miRNAs of A. thaliana and 462 miRNAs of G. max, of which 213 for A. thaliana and 397 for G. max are existing in miRBase 20.

CONCLUSIONS

There is no universal SSR that is conserved among all precursors of Viridiplantae, but conserved SSR exists within a miRNA family and is used as a signature in our prediction method. Prediction of known miRNAs of A. thaliana and G. max validates the accuracy of our method. Our findings will contribute to the present knowledge of miRNAs and their targets in P. vulgaris. This computational method can be applied to any species of Viridiplantae for the successful prediction of miRNAs and their targets.

Identification and characterization of microRNAs at different flowering developmental stages in moso bamboo (Phyllostachys edulis) by high-throughput sequencing

Researching moso bamboo flowering has been difficult because of its unknown flowering interval and the rarity of florescent samples. To identify microRNAs (miRNAs) and study their expression patterns during the flower developmental process of moso bamboo, small RNAs from non-flowering leaves and four flower developmental periods were sequenced using Illumina technology. In total, 409 known miRNAs and 492 differentially expressed novel miRNAs were identified in moso bamboo. Of the known miRNAs that were differentially expressed between non-flowering and flowering samples, 64 were predicted to have a total of 308 targets. Among the miRNAs, seven known and five novel miRNAs were selected, as were four of their target genes, and their expression profiles were validated using qRT-PCR. The results indicated that the miRNA expression levels were negatively correlated with those of their targets. The research comprehensively revealed that the differentially expressed miRNAs and their targets participated in diverse biological pathways and played significant regulatory roles in moso bamboo flowering. The data provide a significant resource for understanding the molecular mechanisms in moso bamboo flowering and senescence, and serve as the primary foundation for further studies on metabolic regulatory networks that involve miRNAs.

Genome-wide identification of the Phaseolus vulgaris sRNAome using small RNA and degradome sequencing

BACKGROUND

MiRNAs and phasiRNAs are negative regulators of gene expression. These small RNAs have been extensively studied in plant model species but only 10 mature microRNAs are present in miRBase version 21, the most used miRNA database, and no phasiRNAs have been identified for the model legume Phaseolus vulgaris. Thanks to the recent availability of the first version of the common bean genome, degradome data and small RNA libraries, we are able to present here a catalog of the microRNAs and phasiRNAs for this organism and, particularly, we suggest new protagonists in the symbiotic nodulation events.

RESULTS

We identified a set of 185 mature miRNAs, including 121 previously unpublished sequences, encoded by 307 precursors and distributed in 98 families. Degradome data allowed us to identify a total of 181 targets for these miRNAs. We reveal two regulatory networks involving conserved miRNAs: those known to play crucial roles in the establishment of nodules, and novel miRNAs present only in common bean, suggesting a specific role for these sequences. In addition, we identified 125 loci that potentially produce phased small RNAs, with 47 of them having all the characteristics of being triggered by a total of 31 miRNAs, including 14 new miRNAs identified in this study.

CONCLUSIONS

We provide here a set of new small RNAs that contribute to the broader knowledge of the sRNAome of Phaseolus vulgaris. Thanks to the identification of the miRNA targets from degradome analysis and the construction of regulatory networks between the mature microRNAs, we present here the probable functional regulation associated with the sRNAome and, particularly, in N2-fixing symbiotic nodules.

The micro-RNA72c-APETALA2-1 node as a key regulator of the common bean-Rhizobium etli nitrogen fixation symbiosis

Micro-RNAs are recognized as important posttranscriptional regulators in plants. The relevance of micro-RNAs as regulators of the legume-rhizobia nitrogen-fixing symbiosis is emerging. The objective of this work was to functionally characterize the role of micro-RNA172 (miR172) and its conserved target APETALA2 (AP2) transcription factor in the common bean (Phaseolus vulgaris)-Rhizobium etli symbiosis. Our expression analysis revealed that mature miR172c increased upon rhizobial infection and continued increasing during nodule development, reaching its maximum in mature nodules and decaying in senescent nodules. The expression of AP2-1 target showed a negative correlation with miR172c expression. A drastic decrease in miR172c and high AP2-1 mRNA levels were observed in ineffective nodules. Phenotypic analysis of composite bean plants with transgenic roots overexpressing miR172c or a mutated AP2-1 insensitive to miR172c cleavage demonstrated the pivotal regulatory role of the miR172 node in the common bean-rhizobia symbiosis. Increased miR172 resulted in improved root growth, increased rhizobial infection, increased expression of early nodulation and autoregulation of nodulation genes, and improved nodulation and nitrogen fixation. In addition, these plants showed decreased sensitivity to nitrate inhibition of nodulation. Through transcriptome analysis, we identified 114 common bean genes that coexpressed with AP2-1 and proposed these as being targets for transcriptional activation by AP2-1. Several of these genes are related to nodule senescence, and we propose that they have to be silenced, through miR172c-induced AP2-1 cleavage, in active mature nodules. Our work sets the basis for exploring the miR172-mediated improvement of symbiotic nitrogen fixation in common bean, the most important grain legume for human consumption.

Transcriptome-Wide Identification of miRNAs and Their Targets from Typha angustifolia by RNA-Seq and Their Response to Cadmium Stress

MicroRNAs (miRNAs) play important roles in plant responses to environmental stress. In this work, we used high-throughput sequencing to analyze transcriptome and small RNAs (sRNAs) in Typha angustifolia under cadmium (Cd) stress. 57,608,230 raw reads were obtained from deep sequencing of a pooled cDNA library. Sequence assembly and analysis yielded 102,473 unigenes. We subsequently sequenced two sRNA libraries from T. angustifolia with or without Cd exposure respectively. Based on transcriptome data of T. angustifolia, we catalogued and analyzed the sRNAs, resulting in the identification of 114 conserved miRNAs and 41 novel candidate miRNAs in both small RNA libraries. In silico analysis revealed 764 targets for 89 conserved miRNAs and 21 novel miRNAs. Statistical analysis on sequencing reads abundance and experimental validation revealed that 4 conserved and 6 novel miRNAs showed specific expression. Combined with function of target genes, these results suggested that miRNAs might play a role in plant Cd stress response. This study provided the first transcriptome-based analysis of miRNAs and their targets responsive to Cd stress in T. angustifolia, which provide a framework for further analysis of miRNAs and their role in regulating plant responses to Cd stress.

Lipopolysaccharide perception leads to dynamic alterations in the microtranscriptome of Arabidopsis thaliana cells and leaf tissues

BACKGROUND

MicroRNAs (miRNAs) are non-coding RNA molecules which have recently emerged as important gene regulators in plants and their gene expression analysis is becoming increasingly important. miRNAs regulate gene expression at the post-transcriptional level by translational repression or target degradation of specific mRNAs and gene silencing. In order to profile the microtranscriptome of Arabidopsis thaliana leaf and callus tissues in response to bacterial lipopolysaccharide (LPS), small RNA libraries were constructed at 0 and 3 h post induction with LPS and sequenced by Illumina sequencing technology.

RESULTS

Differential regulation of subset of miRNAs in response to LPS treament was observed. Small RNA reads were mapped to the miRNA database and 358 miRNAs belonging to 49 miRNA families in the callus tissues and 272 miRNAs belonging to 40 miRNA families in the leaf tissues were identified. Moreover, target genes for all the identified miRNAs families in the leaf tissues and 44 of the 49 miRNAs families in the callus tissues were predicted. The sequencing analysis showed that in both callus and leaf tissues, various stress regulated-miRNAs were differentially expressed and real time PCR validated the expression profile of miR156, miR158, miR159, miR169, miR393, miR398, miR399 and miR408 along with their target genes.

CONCLUSION

A. thaliana callus and leaf callus tissues respond to LPS as a microbe-associated molecular pattern molecule through dynamic changes to the microtranscriptome associated with differential transcriptional regulation in support of immunity and basal resistance.

Identification and characterization of microRNAs from tree peony (Paeonia ostii) and their response to copper stress

MicroRNAs (miRNAs) are a class of non-coding, small RNAs recognized as important regulators of gene expression. Although plant miRNAs have been extensively studied in model systems, less is known in other plants with limited genome sequence data, including Paeonia ostii. In this work, we used high-throughput sequencing to identify conserved and nonconserved miRNAs and other short RNAs in Paeonia ostii under control and copper stressed condition. 102 previously known plant miRNAs were identified and classified into 89 families according to their gene sequence identity. Some miRNAs were highly conserved in the plant kingdom suggesting that these miRNA play important and conserved roles. Combined our transcriptome sequencing data of Paeonia ostii under same conditions, 34 novel potential miRNAs were identified. The potential targets of the identified known and novel miRNAs were also predicted based on sequence homology search. Comparing the two libraries, it was observed that 12 conserved miRNAs and 18 novel miRNAs showed significantly changes in response to copper stress. Some of the new identified potential miRNAs might be involved in Paeonia ostii-specific regulating mechanisms under copper stress. These results provide a framework for further analysis of miRNAs and their role in regulating Paeonia ostii response to copper stress.

Comparative expression profiling of miRNAs between the cytoplasmic male sterile line MeixiangA and its maintainer line MeixiangB during rice anther development

miRNAs are involved in the pollen development during the CMS occurrence in rice. miRNAs are 20-24 nt endogenously expressed small RNAs that play key roles in the regulation of many growth and developmental processes in plants. The knowledge on cytoplasmic male sterility (CMS) regulation by miRNAs in rice is rather limited. In this study, Illumina sequencing was employed to examine the expression profiles of rice anther miRNAs from the CMS line MeixiangA (MxA) and its maintainer line MeixiangB (MxB). A total of 518 known miRNAs and 144 novel miRNAs were identified during rice anther development. Based on the number of sequencing reads, a total of 24 miRNAs were discovered to be differentially expressed between MxA and MxB, and the results were partially validated by qRT-PCR. Among these, 16 miRNAs were decreased and 8 miRNAs were increased in MxA compared with MxB. Target prediction showed that they target genes encoding EF-hand family proteins, F-box domain-containing proteins, MYB transcription factors, PPR-containing proteins and transposons. The expression patterns for targets of osa-miR528, osa-miR5793, osa-miR1432, osa-miR159, osa-miR812d, osa-miR2118c, osa-miR172d and osa-miR5498 were selectively examined, and the results showed that there was a negative correlation on the expression patterns between miRNAs and their targets. These targets have previously been reported to be related with pollen development and male sterility, suggesting that miRNAs might act as regulators of CMS occurrence in rice anthers. Furthermore, miRNA editing events were observed. The U → C and U → A editing phenomenon was validated by molecular cloning and sequencing. These findings contribute to our understanding of the roles of miRNAs during anther development and CMS occurrence in rice.