MicroRNAs in the shoot apical meristem of soybean

Flowering onset time is regulated by microRNA-mediated trehalose-6-phosphate signaling in Cajanus cajan L. under elevated CO2

CO2 levels are known to have an impact on plant development and physiology. In the current study, we have investigated the effect of elevated CO2 on flowering and its regulation through miRNA mediated sugar signaling. We also unraveled small RNA transcriptome of pigeonpea under ambient and elevated CO2 conditions and predicted the targets for crucial miRNAs through computational methods. The results have shown that the delayed flowering in pigeonpea under elevated CO2 was due to an imbalance in C:N stoichiometry and differential expression pattern of aging pathway genes, including SQUAMOSA PROMOTER BINDING PROTEIN-LIKE. Furthermore, qRT PCR analysis has revealed the role of miR156 and miR172 in mediating trehalose-6-phosphate dependent flowering regulation. The current study is crucial in understanding the responses of flowering patterns in a legume crop to elevated CO2 which showed a significant impact on its final yields. Also, these findings are crucial in devising effective crop improvement strategies for developing climate resilient crops, including pigeonpea.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01434-9.

Circular RNAs modulate the floral fate acquisition in soybean shoot apical meristem

BACKGROUND

Soybean (Glycine max), a major oilseed and protein source, requires a short-day photoperiod for floral induction. Though key transcription factors controlling flowering have been identified, the role of the non-coding genome is limited. Circular RNAs (circRNAs) recently emerged as a novel class of RNAs with critical regulatory functions. However, a study on circRNAs during the floral transition of a crop plant is lacking. We investigated the expression and potential function of circRNAs in floral fate acquisition by soybean shoot apical meristem in response to short-day treatment.

RESULTS

Using deep sequencing and in-silico analysis, we denoted 384 circRNAs, with 129 exhibiting short-day treatment-specific expression patterns. We also identified 38 circRNAs with predicted binding sites for miRNAs that could affect the expression of diverse downstream genes through the circRNA-miRNA-mRNA network. Notably, four different circRNAs with potential binding sites for an important microRNA module regulating developmental phase transition in plants, miR156 and miR172, were identified. We also identified circRNAs arising from hormonal signaling pathway genes, especially abscisic acid, and auxin, suggesting an intricate network leading to floral transition.

CONCLUSIONS

This study highlights the gene regulatory complexity during the vegetative to reproductive transition and paves the way to unlock floral transition in a crop plant.

Functional analysis of soybean miR156 and miR172 in tobacco highlights their role in plant morphology and floral transition

Soybean (Glycine max), a significant oilseed and protein source for humans and livestock feed, needs short day photoperiod for floral induction. Further, soybean has a paleopolyploid genome with multiple copies of flowering genes adding to the complexity of genetic regulation of flowering, and seed set, especially in investigating the role of the noncoding genome. microRNAs, a class of noncoding RNA, play a regulatory role in plant development. miR156 and miR172 are major components of the essential regulatory hub controlling juvenile and vegetative developments and initiation of reproductive phase change leading to flowering. These microRNAs have been originally isolated and studied from model plant, Arabidopsis. However, a study on soybean microRNAs is lacking. We investigated the temporal expression patterns of gma-miR156a and gma-miR172a and found inversely related - gma-miR156a expression was higher in the vegetative stage, and gma-miR172a expression was elevated under inductive flowering conditions. The functions of gma-miR156a and gma-miR172a were evaluated via heterologous expressions in transgenic tobacco plants (Nicotiana tabacum L.). The analysis of overexpression transgenic lines highlighted that gma-miR156a plays a role in juvenile development via repression of the SPL transcription factor family. In contrast, gma-miR172a plays a pivotal role in the reproductive development phase by down-regulating its target genes, AP2. In addition, ectopic expression of gma-miR156a and gma-miR172a affected plant morphology and physiology during plant growth. Collectively, our results suggest that gma-miR156a and gma-miR172a regulate multiple morpho-physiological traits that could be used to enhance crop yield under changing climate conditions.

What Do We Know about Barley miRNAs?

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

Transcriptome and miRNAome analysis reveals components regulating tissue differentiation of bamboo shoots

Primary thickening determines bamboo yield and wood property. However, little is known about the regulatory networks involved in this process. This study identified a total of 58,652 genes and 150 miRNAs via transcriptome and small RNA sequencing using the underground thickening shoot samples of wild-type (WT) Moso bamboo (Phyllostachys edulis) and a thick wall (TW) variant (P. edulis "Pachyloen") at five developmental stages (WTS1/TWS1-WTS5/TWS5). A total of 14,029 (65.17%) differentially expressed genes and 68 (45.33%) differentially expressed miRNAs were identified from the WT, TW, and WTTW groups. The first two groups were composed of four pairwise combinations, each between two successive stages (WTS2/TWS2_versus_WTS1/TWS1, WTS3/TWS3_versus_WTS2/TWS2, WTS4/TWS4_versus_WTS3/TWS3, and WTS5/TWS5_versus_WTS4/TWS4), and the WTTW group was composed of five combinations, each between two relative stages (TWS1-5_versus_WTS1-5). Additionally, among the phytohormones, zeatin showed more remarkable changes in concentrations than indole-3-acetic acid, gibberellic acid, and abscisic acid throughout the five stages in the WT and the TW groups. Moreover, 125 cleavage sites were identified for 387 miRNA-mRNA pairs via degradome sequencing (P < 0.05). The dual-luciferase reporter assay confirmed that 13 miRNAs bound to 12 targets. Fluorescence in situ hybridization localized miR166 and miR160 in the shoot apical meristem and the procambium of Moso bamboo shoots at the S1 stage. Thus, primary thickening is a complex process regulated by miRNA-gene-phytohormone networks, and the miRNAome and transcriptome dynamics regulate phenotypic plasticity. These findings provide insights into the molecular mechanisms underlying wood formation and properties and propose targets for bamboo breeding.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Biological Parts for Engineering Abiotic Stress Tolerance in Plants

It is vital to ramp up crop production dramatically by 2050 due to the increasing global population and demand for food. However, with the climate change projections showing that droughts and heatwaves becoming common in much of the globe, there is a severe threat of a sharp decline in crop yields. Thus, developing crop varieties with inbuilt genetic tolerance to environmental stresses is urgently needed. Selective breeding based on genetic diversity is not keeping up with the growing demand for food and feed. However, the emergence of contemporary plant genetic engineering, genome-editing, and synthetic biology offer precise tools for developing crops that can sustain productivity under stress conditions. Here, we summarize the systems biology-level understanding of regulatory pathways involved in perception, signalling, and protective processes activated in response to unfavourable environmental conditions. The potential role of noncoding RNAs in the regulation of abiotic stress responses has also been highlighted. Further, examples of imparting abiotic stress tolerance by genetic engineering are discussed. Additionally, we provide perspectives on the rational design of abiotic stress tolerance through synthetic biology and list various bioparts that can be used to design synthetic gene circuits whose stress-protective functions can be switched on/off in response to environmental cues.

Green Chemistry Production of Codlemone, the Sex Pheromone of the Codling Moth (Cydia pomonella), by Metabolic Engineering of the Oilseed Crop Camelina (Camelina sativa)

Synthetic pheromones have been used for pest control over several decades. The conventional synthesis of di-unsaturated pheromone compounds is usually complex and costly. Camelina (Camelina sativa) has emerged as an ideal, non-food biotech oilseed platform for production of oils with modified fatty acid compositions. We used Camelina as a plant factory to produce mono- and di-unsaturated C₁₂ chain length moth sex pheromone precursors, (E)-9-dodecenoic acid and (E,E)-8,10-dodecadienoic acid, by introducing a fatty acyl-ACP thioesterase FatB gene UcTE from California bay laurel (Umbellularia californica) and a bifunctional ∆9 desaturase gene Cpo_CPRQ from the codling moth, Cydia pomonella. Different transgene combinations were investigated for increasing pheromone precursor yield. The most productive Camelina line was engineered with a vector that contained one copy of UcTE and the viral suppressor protein encoding P19 transgenes and three copies of Cpo_CPRQ transgene. The T₂ generation of this line produced 9.4% of (E)-9-dodecenoic acid and 5.5% of (E,E)-8,10-dodecadienoic acid of the total fatty acids, and seeds were selected to advance top-performing lines to homozygosity. In the T₄ generation, production levels of (E)-9-dodecenoic acid and (E,E)-8,10-dodecadienoic acid remained stable. The diene acid together with other seed fatty acids were converted into corresponding alcohols, and the bioactivity of the plant-derived codlemone was confirmed by GC-EAD and a flight tunnel assay. Trapping in orchards and home gardens confirmed significant and specific attraction of C. pomonella males to the plant-derived codlemone.

MicroRNA Omics Analysis of Camellia sinesis Pollen Tubes in Response to Low-Temperature and Nitric Oxide

Nitric oxide (NO) as a momentous signal molecule participates in plant reproductive development and responds to various abiotic stresses. Here, the inhibitory effects of the NO-dominated signal network on the pollen tube growth of Camellia sinensis under low temperature (LT) were studied by microRNA (miRNA) omics analysis. The results showed that 77 and 71 differentially expressed miRNAs (DEMs) were induced by LT and NO treatment, respectively. Gene ontology (GO) analysis showed that DEM target genes related to microtubules and actin were enriched uniquely under LT treatment, while DEM target genes related to redox process were enriched uniquely under NO treatment. In addition, the target genes of miRNA co-regulated by LT and NO are only located on the cell membrane and cell wall, and most of them are enriched in metal ion binding and/or transport and cell wall organization. Furthermore, DEM and its target genes related to metal ion binding/transport, redox process, actin, cell wall organization and carbohydrate metabolism were identified and quantified by functional analysis and qRT-PCR. In conclusion, miRNA omics analysis provides a complex signal network regulated by NO-mediated miRNA, which changes cell structure and component distribution by adjusting Ca²⁺ gradient, thus affecting the polar growth of the C. sinensis pollen tube tip under LT.

miRNA expression profiling and zeatin dynamic changes in a new model system of in vivo indirect regeneration of tomato

Callus formation and adventitious shoot differentiation could be observed on the cut surface of completely decapitated tomato plants. We propose that this process can be used as a model system to investigate the mechanisms that regulate indirect regeneration of higher plants without the addition of exogenous hormones. This study analyzed the patterns of trans-zeatin and miRNA expression during in vivo regeneration of tomato. Analysis of trans-zeatin revealed that the hormone cytokinin played an important role in in vivo regeneration of tomato. Among 183 miRNAs and 1168 predicted target genes sequences identified, 93 miRNAs and 505 potential targets were selected based on differential expression levels for further characterization. Expression patterns of six miRNAs, including sly-miR166, sly-miR167, sly-miR396, sly-miR397, novel 156, and novel 128, were further validated by qRT-PCR. We speculate that sly-miR156, sly-miR160, sly-miR166, and sly-miR397 play major roles in callus formation of tomato during in vivo regeneration by regulating cytokinin, IAA, and laccase levels. Overall, our microRNA sequence and target analyses of callus formation during in vivo regeneration of tomato provide novel insights into the regulation of regeneration in higher plants.

Uncloaking lncRNA-meditated gene expression as a potential regulator of CMS in cotton (Gossypium hirsutum L.)

Cytoplasmic male sterility is a well-proven mechanism for cotton hybrid production. Long non-coding RNAs belong to a class of transcriptional regulators that function in multiple biological processes. The cDNA libraries from the flower buds of the cotton CGMS, it's restorer (Rf) and maintainer lines were sequenced using high throughput NGS technique. A total of 1531 lncRNAs showed significant differential expression patterns between these three lines. Functional analysis of the co-expression network of lncRNA-mRNA using gene ontology vouchsafes that, lncRNAs play a crucial role in cytoplasmic male sterility and fertility restoration through pollen development, INO80 complex, development of anther wall tapetum, chromatin remodeling, and histone modification. Additionally, 94 lncRNAs were identified as putative precursors of 49 miRNAs. qRT-PCR affirms the concordance of expression pattern to RNA-seq data. These findings divulge the lncRNA driven miRNA-mediated regulation of gene expression profiling superintended for a better understanding of the CMS mechanisms of cotton.

MicroRNAs and their targeted genes associated with phase changes of stem explants during tissue culture of tea plant

Elucidation of the molecular mechanism related to the dedifferentiation and redifferentiation during tissue culture will be useful for optimizing regeneration system of tea plant. In this study, an integrated sRNAome and transcriptome analyses were carried out during phase changes of the stem explant culture. Among 198 miRNAs and 8001 predicted target genes, 178 differentially expressed miRNAs and 4264 potential targets were screened out from explants, primary calli, as well as regenerated roots and shoots. According to KEGG analysis of the potential targets, pathway of "aminoacyl-tRNA biosynthesis", "proteasome" and "glutathione metabolism" was of great significance during the dedifferentiation, and pathway of "porphyrin and chlorophyll metabolism", "mRNA surveillance pathway", "nucleotide excision repair" was indispensable for redifferentiation of the calli. Expression pattern of 12 miRNAs, including csn-micR390e, csn-miR156b-5p, csn-miR157d-5p, csn-miR156, csn-miR166a-3p, csn-miR166e, csn-miR167d, csn-miR393c-3p, csn-miR394, csn-miR396a-3p, csn-miR396 and csn-miR396e-3p, was validated by qRT-PCR among 57 differentially expressed phase-specific miRNAs. Validation also confirmed that regulatory module of csn-miR167d/ERF3, csn-miR156/SPB1, csn-miR166a-3p/ATHB15, csn-miR396/AIP15A, csn-miR157d-5p/GST and csn-miR393c-3p/ATG18b might play important roles in regulating the phase changes during tissue culture of stem explants.

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

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

Identification and Analysis of microRNAs in the SAM and Leaves of Populus tomentosa

The shoot apical meristem (SAM) is a crucial tissue located at the tops of plants which can continually grow and differentiate to develop into all aboveground parts. SAM development is controlled by a series of complicated molecular regulation networks, among which microRNAs (miRNAs) and their target genes play key roles. However, little is known about these miRNAs in woody plants. In this study, we used small RNA (sRNA) sequencing to build four libraries derived from shoot tips and mature leaf tissues of Populus tomentosa, and identified 99 known miRNA families. In addition, 193 known miRNAs, including phytohormone-, developmental-, and cellular process-related miRNAs, showed significant differential expression. Interestingly, quantitative real-time reverse transcription polymerase chain reaction (PCR) analysis of miR172, miR164, and miR393 expression showed marked changes in expression patterns during the development of shoot tips. The target genes of these miRNAs were involved in the regulation of hormone responses and stem cell function. In particular, the miR172 target APETALA2 (AP2), involved in the maintenance of stem cells in the shoot apex, was expressed specifically during the initial active stage of development. These findings provide new insights into the regulatory mechanisms of miRNAs involved in SAM development and differentiation in tree species.

Genome-wide identification and profiling of microRNAs in Paulownia tomentosa cambial tissues in response to seasonal changes

MicroRNAs (miRNAs), a group of endogenous small non-coding RNAs, have been shown to play essential roles in the regulation of gene expression at the post-transcriptional level. Although Paulownia tomentosa is an ecologically and economically important timber species due to its rapid growth, few efforts have focused on small RNAs (sRNAs) in the cambial tissues during winter and summer transition. In the present study, we identified 33 known miRNA families and 29 novel miRNAs which include 20 putative novel miRNAs* in P. tomentosa cambial tissues during winter and summer transition. Through differential expression analysis, we showed that 15 known miRNAs and 8 novel miRNAs were preferentially abundant in certain stage of cambial tissues. Based on the P. tomentosa mRNA transcriptome database, 1667 and 78 potential targets were predicted for 29 known and 20 novel miRNAs, respectively and the predicted targets are mostly transcription factors and functional genes. The targets of these miRNAs were enriched in "metabolic process" and "transcription regulation" by using Gene Ontology enrichment analysis. In addition, KEGG pathway analyses revealed the involvement of miRNAs in starch and sucrose metabolism and plant-pathogen interaction metabolism pathways. Noticeably, qRT-PCR expression analysis demonstrated that 9 miRNAs and their targets were existed a negative correlation in P. tomentosa cambial tissues. This study is the first to examine known and novel miRNAs and their potential targets in P. tomentosa cambial tissues during winter and summer transition and identify several candidate genes potentially regulating cambial phase transition, and thus provide a framework for further understanding of miRNAs functions in the regulation of cambial phase transition and wood formation in trees.

The MicroRNA319d/TCP10 Node Regulates the Common Bean - Rhizobia Nitrogen-Fixing Symbiosis

Micro-RNAs from legume plants are emerging as relevant regulators of the rhizobia nitrogen-fixing symbiosis. In this work we functionally characterized the role of the node conformed by micro-RNA319 (miR319) - TEOSINTE BRANCHED/CYCLOIDEA/PCF (TCP) transcription factor in the common bean (Phaseolus vulgaris) - Rhizobium tropici symbiosis. The miR319d, one of nine miR319 isoforms from common bean, was highly expressed in root and nodules from inoculated plants as compared to roots from fertilized plants. The miR319d targets TCP10 (Phvul.005G067950), identified by degradome analysis, whose expression showed a negative correlation with miR319d expression. The phenotypic analysis of R. tropici-inoculated composite plants with transgenic roots/nodules overexpressing or silencing the function of miR319d demonstrated the relevant role of the miR319d/TCP10 node in the common bean rhizobia symbiosis. Increased miR319d resulted in reduced root length/width ratio, increased rhizobial infection evidenced by more deformed root hairs and infection threads, and decreased nodule formation and nitrogenase activity per plant. In addition, these plants with lower TCP10 levels showed decreased expression level of the jasmonic acid (JA) biosynthetic gene: LOX2. The transcription of LOX2 by TCPs has been demonstrated for Arabidopsis and in several plants LOX2 level and JA content have been associate with TCP levels. On this basis, we propose that in roots/nodules of inoculated common bean plants TCP10 could be the transcriptional regulator of LOX2 and the miR319d/TCP10 node could affect nodulation through JA signaling. However, given the complexity of nodulation, the participation of other signaling pathways in the phenotypes observed cannot be ruled out.

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

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

High-throughput sequencing analysis revealed the regulation patterns of small RNAs on the development of A. comosus var. bracteatus leaves

Studies of the molecular mechanisms involved in the formation of the albino leaf cells are important for understanding the development of chimera leaves in Ananas comosus var. bracteatus. In this study, we identified a total of 163 novel miRNAs involved in the development of complete white (CWh) and complete green (CGr) leaves using high-throughput sequencing method. The potential miRNA target genes were predicted and annotated using the NR, Swiss-Prot, GO, COG, KEGG, KOG and Pfam databases. The main biological processes regulated by miRNAs were revealed. The miRNAs which potentially play important roles in the development of the leaves and the albino of the CWh leaf cells were selected and their expression patterns were analyzed. The expression levels of nine miRNAs and their potential target genes were studied using qRT-PCR. These results will help to elucidate the functional and regulatory roles of miRNAs in the formation of the albino cells and the development of the leaves of A. comosus var. bracteatus. These data may also be helpful as a resource for studies of small RNA in other leaf color chimeric plant species.

Whole Mount in situ Localization of miRNAs and mRNAs During Somatic Embryogenesis in Arabidopsis

Somatic embryogenesis (SE) results from the transition of differentiated plant somatic cells into embryogenic cells that requires the extensive reprogramming of the somatic cell transcriptome. Commonly, the SE-involved genes are identified by analyzing the heterogeneous population of explant cells and thus, it is necessary to validate the expression of the candidate genes in the cells that are competent for embryogenic transition. Here, we optimized and implemented the whole mount in situ hybridization (WISH) method (Bleckmann and Dresselhaus, 2016; Dastidar et al., 2016) in order to analyze the spatiotemporal localization of miRNAs (miR156, miR166, miR390, miR167) and mRNAs such as WOX5 and PHABULOSA-target of miR165/166 during the SE that is induced in Arabidopsis explants. This study presents a detailed step-by-step description of the WISH procedure in which DIG-labeled LNA and RNA probes were used to detect miRNAs and mRNAs, respectively. The usefulness of the WISH in the functional analysis of the SE-involved regulatory pathways is demonstrated and the advantages of this method are highlighted: (i) the ability to analyze intact non-sectioned plant tissue; (ii) the specificity of transcript detection; (iii) the detection of miRNA; and (iv) a semi-quantitative assessment of the RNA abundance.

Genome wide identification of microRNAs involved in fatty acid and lipid metabolism of Brassica napus by small RNA and degradome sequencing

Rapeseed (Brassica napus) is an important cash crop considered as the third largest oil crop worldwide. Rapeseed oil contains various saturation or unsaturation fatty acids, these fatty acids, whose could incorporation with TAG form into lipids stored in seeds play various roles in the metabolic activity. The different fatty acids in B. napus seeds determine oil quality, define if the oil is edible or must be used as industrial material. miRNAs are kind of non-coding sRNAs that could regulate gene expressions through post-transcriptional modification to their target transcripts playing important roles in plant metabolic activities. We employed high-throughput sequencing to identify the miRNAs and their target transcripts involved in fatty acids and lipids metabolism in different development of B. napus seeds. As a result, we identified 826 miRNA sequences, including 523 conserved and 303 newly miRNAs. From the degradome sequencing, we found 589 mRNA could be targeted by 236 miRNAs, it includes 49 novel miRNAs and 187 conserved miRNAs. The miRNA-target couple suggests that bna-5p-163957_18, bna-5p-396192_7, miR9563a-p3, miR9563b-p5, miR838-p3, miR156e-p3, miR159c and miR1134 could target PDP, LACS9, MFPA, ADSL1, ACO32, C0401, GDL73, PlCD6, OLEO3 and WSD1. These target transcripts are involving in acetyl-CoA generate and carbon chain desaturase, regulating the levels of very long chain fatty acids, β-oxidation and lipids transport and metabolism process. At the same, we employed the q-PCR to valid the expression of miRNAs and their target transcripts that involve in fatty acid and lipid metabolism, the result suggested that the miRNA and their transcript expression are negative correlation, which in accord with the expression of miRNA and its target transcript. The study findings suggest that the identified miRNA may play important role in the fatty acids and lipids metabolism in seeds of B. napus.

Genome-wide analysis of long intergenic non-coding RNAs in chickpea and their potential role in flower development

Non-coding RNAs constitute a major portion of the transcriptome in most of eukaryotes. Long non-coding transcripts originating from the DNA segment present between the protein coding genes are termed as long intergenic non-coding RNAs (lincRNAs). Several evidences suggest the role of lincRNAs in regulation of various biological processes. In this study, we identified a total of 2248 lincRNAs in chickpea using RNA-seq data from eight successive stages of flower development and three vegetative tissues via an optimized pipeline. Different characteristic features of lincRNAs were studied and compared with those of predicted mRNAs in chickpea. Further, we utilized a method using network propagation algorithm to reveal the putative function of lincRNAs in plants. In total, at least 79% of the identified chickpea lincRNAs were assigned with a putative function. A comprehensive expression profiling revealed differential expression patterns and tissue specificity of lincRNAs in different stages of flower development in chickpea. In addition, potential lincRNAs-miRNA interactions were explored for the predicted lincRNAs in chickpea. These findings will pave the way for understanding the role of lincRNAs in the regulatory mechanism underlying flower development in chickpea and other legumes.

Identification and Comparative Analysis of microRNA in Wheat (Triticum aestivum L.) Callus Derived from Mature and Immature Embryos during In vitro Culture

Feasible and efficient tissue culture plays an important role in plant genetic engineering. Wheat (Triticum aestivum L.) immature embryos (IMEs) are preferred for tissue culture to mature embryos (MEs) because IMEs easily generate embryogenic callus, producing large number of plants. The molecular mechanisms of regulation and the biological pathways involved in embryogenic callus formation in wheat remain unclear. Here, microRNAs (miRNAs) potentially involved in embryogenic callus formation and somatic embryogenesis were identified through deep sequencing of small RNAs (sRNAs) and analyzed with bioinformatics tools. Six sRNA libraries derived from calli of IMEs and MEs after 3, 6, or 15 d of culture (DC) were constructed and sequenced. A total of 85 known miRNAs were identified, of which 30, 33, and 18 were differentially expressed (P < 0.05) between the IME and ME libraries at 3, 6, and 15 DC, respectively. Additionally, 171 novel and 41 candidate miRNAs were also identified, of the novel miRNA, 69, 67, and 37 were differentially expressed (P < 0.05) between the two types of libraries at 3, 6, and 15 DC, respectively. The expression patterns of eight known and eight novel miRNAs were validated using quantitative real-time polymerase chain reaction. Gene ontology annotation of differentially expressed miRNA targets provided information regarding the underlying molecular functions, biological processes, and cellular components involved in embryogenic callus development. Functional miRNAs, such as miR156, miR164, miR1432, miR398, and miR397, differentially expressed in IMEs and MEs might be related to embryogenic callus formation and somatic embryogenesis. This study suggests that miRNA plays an important role in embryogenic callus formation and somatic embryogenesis in wheat, and our data provide a useful resource for further research.

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.

Identification of MicroRNAs in Response to Different Day Lengths in Soybean Using High-Throughput Sequencing and qRT-PCR

MicroRNAs (miRNAs) are short, non-coding single-strand RNA molecules that play important roles in plant growth, development and stress responses. Flowering time affects the seed yield and quality of soybean. However, the miRNAs involved in the regulation of flowering time in soybean have not been reported until recently. Here, high-throughput sequencing and qRT-PCR were used to identify miRNAs involved in soybean photoperiodic pathways. The first trifoliate leaves of soybean that receive the signal of light treatment were used to construct six libraries (0, 8, and 16 h under short-day (SD) treatment and 0, 8, and 16 h under long-day (LD) treatment). The libraries were sequenced using Illumina Solexa. A total of 318 known plant miRNAs belonging to 163 miRNA families and 81 novel predicted miRNAs were identified. Among these, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h, showed differences in abundance between LD and SD treatments. Furthermore, the results of GO and KEGG analyses indicated that most of the miRNA targets were transcription factors. Seven miRNAs at 0 h, 23 miRNAs (including four novel predicted miRNAs) at 8 h, 16 miRNAs (including one novel predicted miRNA) at 16 h and miRNA targets were selected for qRT-PCR analysis to assess the accuracy of the sequencing and target prediction. The results indicated that the expression patterns of the selected miRNAs and miRNA targets showed no differences between the qRT-PCR and sequencing results. In addition, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h responded to day length changes in soybean, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h. These results provided an important molecular basis to understand the regulation of flowering time through photoperiodic pathways in soybean.

High-throughput sequencing of small RNAs and anatomical characteristics associated with leaf development in celery

MicroRNAs (miRNAs) exhibit diverse and important roles in plant growth, development, and stress responses and regulate gene expression at the post-transcriptional level. Knowledge about the diversity of miRNAs and their roles in leaf development in celery remains unknown. To elucidate the roles of miRNAs in celery leaf development, we identified leaf development-related miRNAs through high-throughput sequencing. Small RNA libraries were constructed using leaves from three stages (10, 20, and 30 cm) of celery cv.'Ventura' and then subjected to high-throughput sequencing and bioinformatics analysis. At Stage 1, Stage 2, and Stage 3 of 'Ventura', a total of 333, 329, and 344 conserved miRNAs (belonging to 35, 35, and 32 families, respectively) were identified. A total of 131 miRNAs were identified as novel in 'Ventura'. Potential miRNA target genes were predicted and annotated using the eggNOG, GO, and KEGG databases to explore gene functions. The abundance of five conserved miRNAs and their corresponding potential target genes were validated. Expression profiles of novel potential miRNAs were also detected. Anatomical characteristics of the leaf blades and petioles at three leaf stages were further analyzed. This study contributes to our understanding on the functions and molecular regulatory mechanisms of miRNAs in celery leaf development.

The functional scope of plant microRNA-mediated silencing

Deep sequencing has identified a complex set of plant miRNAs that potentially regulates many target genes of high complementarity. Furthermore, the discovery that many plant miRNAs work through a translational repression mechanism, along with the identification of noncanonical targets, has encouraged bioinformatic searches with less stringent parameters, identifying an even wider range of potential targets. Together, these findings suggest that any given plant miRNA family may regulate a highly diverse set of mRNAs. Here we present evolutionary, genetic, and mechanistic evidence that opposes this idea but instead suggests that families of sequence-related miRNAs regulate very few functionally related targets. We propose that complexities beyond complementarity impact plant miRNA target recognition, possibly explaining the current disparity between bioinformatic prediction and functional evidence.

Identification and characterization of Wilt and salt stress-responsive microRNAs in chickpea through high-throughput sequencing

Chickpea (Cicer arietinum) is the second most widely grown legume worldwide and is the most important pulse crop in the Indian subcontinent. Chickpea productivity is adversely affected by a large number of biotic and abiotic stresses. MicroRNAs (miRNAs) have been implicated in the regulation of plant responses to several biotic and abiotic stresses. This study is the first attempt to identify chickpea miRNAs that are associated with biotic and abiotic stresses. The wilt infection that is caused by the fungus Fusarium oxysporum f.sp. ciceris is one of the major diseases severely affecting chickpea yields. Of late, increasing soil salinization has become a major problem in realizing these potential yields. Three chickpea libraries using fungal-infected, salt-treated and untreated seedlings were constructed and sequenced using next-generation sequencing technology. A total of 12,135,571 unique reads were obtained. In addition to 122 conserved miRNAs belonging to 25 different families, 59 novel miRNAs along with their star sequences were identified. Four legume-specific miRNAs, including miR5213, miR5232, miR2111 and miR2118, were found in all of the libraries. Poly(A)-based qRT-PCR (Quantitative real-time PCR) was used to validate eleven conserved and five novel miRNAs. miR530 was highly up regulated in response to fungal infection, which targets genes encoding zinc knuckle- and microtubule-associated proteins. Many miRNAs responded in a similar fashion under both biotic and abiotic stresses, indicating the existence of cross talk between the pathways that are involved in regulating these stresses. The potential target genes for the conserved and novel miRNAs were predicted based on sequence homologies. miR166 targets a HD-ZIPIII transcription factor and was validated by 5′ RLM-RACE. This study has identified several conserved and novel miRNAs in the chickpea that are associated with gene regulation following exposure to wilt and salt stress.

Identification of Cold-Responsive miRNAs and Their Target Genes in Nitrogen-Fixing Nodules of Soybean

As a warm climate species, soybean is highly sensitive to chilling temperatures. Exposure to chilling temperatures causes a significant reduction in the nitrogen fixation rate in soybean plants and subsequent yield loss. However, the molecular basis for the sensitivity of soybean to chilling is poorly understood. In this study, we identified cold-responsive miRNAs in nitrogen-fixing nodules of soybean. Upon chilling, the expression of gma-miR397a, gma-miR166u and gma-miR171p was greatly upregulated, whereas the expression of gma-miR169c, gma-miR159b, gma-miR319a/b and gma-miR5559 was significantly decreased. The target genes of these miRNAs were predicted and validated using 5' complementary DNA ends (5'-RACE) experiments, and qPCR analysis identified putative genes targeted by the cold-responsive miRNAs in response to chilling temperatures. Taken together, our results reveal that miRNAs may be involved in the protective mechanism against chilling injury in mature nodules of soybean.

Unique and conserved features of floral evocation in legumes

Legumes, with their unique ability to fix atmospheric nitrogen, play a vital role in ensuring future food security and mitigating the effects of climate change because they use less fossil energy and produce less greenhouse gases compared with N-fertilized systems. Grain legumes are second only to cereal crops as a source of human and animal food, and they contribute approximately one third of the protein consumed by the human population. The productivity of seed crops, such as grain legumes, is dependent on flowering. Despite the genetic variation and importance of flowering in legume production, studies of the molecular pathways that control flowering in legumes are limited. Recent advances in genomics have revealed that legume flowering pathways are divergent from those of such model species as Arabidopsis thaliana. Here, we discuss the current understanding of flowering time regulation in legumes and highlight the unique and conserved features of floral evocation in legumes.

The 2HA line of Medicago truncatula has characteristics of an epigenetic mutant that is weakly ethylene insensitive

BACKGROUND

The Medicago truncatula 2HA seed line is highly embryogenic while the parental line Jemalong rarely produces embryos. The 2HA line was developed from one of the rare Jemalong regenerates and this method for obtaining a highly regenerable genotype in M. truncatula is readily reproducible suggesting an epigenetic mechanism. Microarray transcriptomic analysis showed down regulation of an ETHYLENE INSENSITIVE 3-like gene in 2HA callus which provided an approach to investigating epigenetic regulation of genes related to ethylene signalling and the 2HA phenotype. Ethylene is involved in many developmental processes including somatic embryogenesis (SE) and is associated with stress responses.

RESULTS

Microarray transcriptomic analysis showed a significant number of up-regulated transcripts in 2HA tissue culture, including nodule and embryo specific genes and transposon-like genes, while only a few genes were down-regulated, including an EIN3-like gene we called MtEIL1. This reduced expression was associated with ethylene insensitivity of 2HA plants that was further investigated. The weak ethylene insensitivity affected root and nodule development. Sequencing of MtEIL1 found no difference between 2HA and wild-type plants. DNA methylation analysis of MtEIL1 revealed significant difference between 2HA and wild-type plants. Tiling arrays demonstrated an elevated level of miRNA in 2HA plants that hybridised to the antisense strand of the MtEIL1 gene. AFLP-like methylation profiling revealed more differences in DNA methylation between 2HA and wild-type. Segregation analysis demonstrated the recessive nature of the eil1 phenotype and the dominant nature of the SE trait.

CONCLUSIONS

We have demonstrated that EIL1 of Medicago truncatula (MtEIL1) is epigenetically silenced in the 2HA seed line. The possible cause is an elevated level of miRNA that targets its 3'UTR and is also associated with DNA methylation of MtEIL1. Down regulation of MtEIL1 makes it possible to form nodules in the presence of ethylene and affects root growth under normal conditions. Segregation analysis showed no association between MtEIL1 expression and SE in culture but the role and mechanism of ethylene signalling in the process of plant regeneration through SE requires further investigation. The work also suggests that epigenetic changes to a particular gene induced in culture can be fixed in regenerated plants.

Analysis of promoters of microRNAs from a Glycine max degradome library

OBJECTIVE

MicroRNAs (miRNAs) are genome-encoded, small non-coding RNAs that play important functions in development, biotic and abiotic stress responses, and other processes. Our aim was to explore the regulation of miRNA expression.

METHODS

We used bioinformatics methods to predict the core promoters of 440 miRNAs identified from a soybean (Glycine max) degradome library and to analyze cis-acting elements for 369 miRNAs.

RESULTS

The prediction results showed that 83.86% of the 440 miRNAs contained promoters in their upstream sequences, and 8.64% (38 loci) in their downstream sequences. The distributions of two core promoter elements, TATA-boxes and transcription start sites (TSSs), were similar. The cis-acting elements were examined to provide clues to the function and regulation of spatiotemporal expression of the miRNAs. Analyses of miRNA cis-elements and targets indicated a potential auxin response factor (ARF)- and gibberellin response factor (GARF)-mediated negative feedback loop for miRNA expression.

CONCLUSIONS

The features of miRNAs from a Glycine max degradome library obtained here provide insights into the transcription regulation and functions of miRNAs in soybean.

Stress-induced early flowering is mediated by miR169 in Arabidopsis thaliana

Plants interact with their environment and they often flower earlier under stress conditions, but how such stress-induced flowering is regulated remains poorly understood. Here evidence is presented that the miR169 family plays a key role in stress-induced flowering in plants. The microRNA (miRNA) miR169 family members are up-regulated in Arabidopsis, maize, and soybean under abiotic stresses. Overexpression of miR169d in Arabidopsis results in early flowering, and overexpression of the miR169d target gene, AtNF-YA2, especially a miR169d-resistant version of AtNF-YA2, results in late flowering. The results suggest that the miR169 family regulates stress-induced flowering by repressing the AtNF-YA transcription factor, which in turn reduces the expression of FLOWERING LOCUS C (FLC), allowing for the expression of FLC target genes such as FLOWERING LOCUS T (FT) and LEAFY (LFY) to promote flowering. It was shown that the expression of genes or miRNAs involved in the other flowering pathways, namely the photoperiod (CO), ambient temperature (SVP), ageing (miR156), and gibberelin (SOC1) pathways, was not affected in miR169d-overexpressing plants, suggesting that stress-induced early flowering is a novel signalling pathway mediated by miR169.

Genome-wide identification of non-coding RNAs interacted with microRNAs in soybean

A wide range of RNA species interacting with microRNAs (miRNAs) form a complex gene regulation network and play vital roles in diverse biological processes. In this study, we performed a genome-wide identification of endogenous target mimics (eTMs) for miRNAs and phased-siRNA-producing loci (PHAS) in soybean with a focus on those involved in lipid metabolism. The results showed that a large number of eTMs and PHAS genes could be found in soybean. Additionally, we found that lipid metabolism related genes were potentially regulated by 28 miRNAs, and nine of them were potentially further regulated by a number of eTMs with expression evidence. Thirty-three miRNAs were found to trigger production of phasiRNAs from 49 PHAS genes, which were able to target lipid metabolism related genes. Degradome data supported miRNA- and/or phasiRNA-mediated cleavage of genes involved in lipid metabolism. Most eTMs for miRNAs involved in lipid metabolism and phasiRNAs targeting lipid metabolism related genes showed a tissue-specific expression pattern. Our bioinformatical evidences suggested that lipid metabolism in soybean is potentially regulated by a complex non-coding network, including miRNAs, eTMs, and phasiRNAs, and the results extended our knowledge on functions of non-coding RNAs.

An RNA-seq transcriptome analysis of histone modifiers and RNA silencing genes in soybean during floral initiation process

Epigenetics has been recognised to play vital roles in many plant developmental processes, including floral initiation through the epigenetic regulation of gene expression. The histone modifying proteins that mediate these modifications involve the SET domain-containing histone methyltransferases, JmjC domain-containing demethylase, acetylases and deacetylases. In addition, RNA interference (RNAi)-associated genes are also involved in epigenetic regulation via RNA-directed DNA methylation and post-transcriptional gene silencing. Soybean, a major crop legume, requires a short day to induce flowering. How histone modifications regulate the plant response to external cues that initiate flowering is still largely unknown. Here, we used RNA-seq to address the dynamics of transcripts that are potentially involved in the epigenetic programming and RNAi mediated gene silencing during the floral initiation of soybean. Soybean is a paleopolyploid that has been subjected to at least two rounds of whole genome duplication events. We report that the expanded genomic repertoire of histone modifiers and RNA silencing genes in soybean includes 14 histone acetyltransferases, 24 histone deacetylases, 47 histone methyltransferases, 15 protein arginine methyltransferases, 24 JmjC domain-containing demethylases and 47 RNAi-associated genes. To investigate the role of these histone modifiers and RNA silencing genes during floral initiation, we compared the transcriptional dynamics of the leaf and shoot apical meristem at different time points after a short-day treatment. Our data reveal that the extensive activation of genes that are usually involved in the epigenetic programming and RNAi gene silencing in the soybean shoot apical meristem are reprogrammed for floral development following an exposure to inductive conditions.

Identification of microRNAs in Caragana intermedia by high-throughput sequencing and expression analysis of 12 microRNAs and their targets under salt stress

142 miRNAs were identified and 38 miRNA targets were predicted, 4 of which were validated, in C. intermedia . The expression of 12 miRNAs in salt-stressed leaves was assessed by qRT-PCR. MicroRNAs (miRNAs) are endogenous small RNAs that play important roles in various biological and metabolic processes in plants. Caragana intermedia is an important ecological and economic tree species prominent in the desert environment of west and northwest China. To date, no investigation into C. intermedia miRNAs has been reported. In this study, high-throughput sequencing of small RNAs and analysis of transcriptome data were performed to identify both conserved and novel miRNAs, and also their target mRNA genes in C. intermedia. Based on sequence similarity and hairpin structure prediction, 132 putative conserved miRNAs (12 of which were confirmed to form hairpin precursors) belonging to 31 known miRNA families were identified. Ten novel miRNAs (including the miRNA* sequences of three novel miRNAs) were also discovered. Furthermore, 36 potential target genes of 17 known miRNA families and 2 potential target genes of 1 novel miRNA were predicted; 4 of these were validated by 5' RACE. The expression of 12 miRNAs was validated in different tissues, and these and five target mRNAs were assessed by qRT-PCR after salt treatment. The expression levels of seven miRNAs (cin-miR157a, cin-miR159a, cin-miR165a, cin-miR167b, cin-miR172b, cin-miR390a and cin-miR396a) were upregulated, while cin-miR398a expression was downregulated after salt treatment. The targets of cin-miR157a, cin-miR165a, cin-miR172b and cin-miR396a were downregulated and showed an approximately negative correlation with their corresponding miRNAs under salt treatment. These results would help further understanding of miRNA regulation in response to abiotic stress in C. intermedia.

Novel members of the AGAMOUS LIKE 6 subfamily of MIKCC-type MADS-box genes in soybean

BACKGROUND

The classical (C) MIKC-type MADS-box transcription factors comprise one gene family that plays diverse roles in the flowering process ranging from floral initiation to the development of floral organs. Despite their importance in regulating developmental processes that impact crop yield, they remain largely unexplored in the major legume oilseed crop, soybean.

RESULTS

We identified 57 MIKC(c)-type transcription factors from soybean and determined the in silico gene expression profiles of the soybean MIKC(c)-type genes across different tissues. Our study implicates three MIKC(c)-type transcription factors as novel members of the AGAMOUS LIKE 6 (AGL6) subfamily of the MIKC(C)-type MADS-box genes, and we named this sister clade PsMADS3. While similar genes were identified in other legume species, poplar and grape, no such gene is represented in Arabidopsis thaliana or rice. RT-PCR analysis on these three soybean PsMADS3 genes during early floral initiation processes revealed their temporal expression similar to that of APETALA1, a gene known to function as a floral meristem identity gene. However, RNA in situ hybridisation showed that their spatial expression patterns are markedly different from those of APETALA1.

CONCLUSION

Legume flower development system differs from that in the model plant, Arabidopsis. There is an overlap in the initiation of different floral whorls in soybean, and inflorescent meristems can revert to leaf production depending on the environmental conditions. MIKC(C)-type MADS-box genes have been shown to play key regulatory roles in different stages of flower development. We identified members of the PsMADS3 sub-clade in legumes that show differential spatial expression during floral initiation, indicating their potential novel roles in the floral initiation process. The results from this study will contribute to a better understanding of legume-specific floral developmental processes.

The dynamics of soybean leaf and shoot apical meristem transcriptome undergoing floral initiation process

Flowering process governs seed set and thus affects agricultural productivity. Soybean, a major legume crop, requires short-day photoperiod conditions for flowering. While leaf-derived signal(s) are essential for the photoperiod-induced floral initiation process at the shoot apical meristem, molecular events associated with early floral transition stages in either leaves or shoot apical meristems are not well understood. To provide novel insights into the molecular basis of floral initiation, RNA-Seq was used to characterize the soybean transcriptome of leaf and micro-dissected shoot apical meristem at different time points after short-day treatment. Shoot apical meristem expressed a higher number of transcripts in comparison to that of leaf highlighting greater diversity and abundance of transcripts expressed in the shoot apical meristem. A total of 2951 shoot apical meristem and 13,609 leaf sequences with significant profile changes during the time course examined were identified. Most changes in mRNA level occurred after 1short-day treatment. Transcripts involved in mediating responses to stimulus including hormones or in various metabolic processes represent the top enriched GO functional category for the SAM and leaf dataset, respectively. Transcripts associated with protein degradation were also significantly changing in leaf and SAM implicating their involvement in triggering the developmental switch. RNA-Seq analysis of shoot apical meristem and leaf from soybean undergoing floral transition reveal major reprogramming events in leaves and the SAM that point toward hormones gibberellins (GA) and cytokinin as key regulators in the production of systemic flowering signal(s) in leaves. These hormones may form part of the systemic signals in addition to the established florigen, FLOWERING LOCUS T (FT). Further, evidence is emerging that the conversion of shoot apical meristem to inflorescence meristem is linked with the interplay of auxin, cytokinin and GA creating a low cytokinin and high GA environment.

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

BACKGROUND

Phosphorus (P) plays important roles in plant growth and development. MicroRNAs involved in P signaling have been identified in Arabidopsis and rice, but P-responsive microRNAs and their targets in soybean leaves and roots are poorly understood.

RESULTS

Using high-throughput sequencing-by-synthesis (SBS) technology, we sequenced four small RNA libraries from leaves and roots grown under phosphate (Pi)-sufficient (+Pi) and Pi-depleted (-Pi) conditions, respectively, and one RNA degradome library from Pi-depleted roots at the genome-wide level. Each library generated ~21.45-28.63 million short sequences, resulting in ~20.56-27.08 million clean reads. From those sequences, a total of 126 miRNAs, with 154 gene targets were computationally predicted. This included 92 new miRNA candidates with 20-23 nucleotides that were perfectly matched to the Glycine max genome 1.0, 70 of which belong to 21 miRNA families and the remaining 22 miRNA unassigned into any existing miRNA family in miRBase 18.0. Under both +Pi and -Pi conditions, 112 of 126 total miRNAs (89%) were expressed in both leaves and roots. Under +Pi conditions, 12 leaf- and 2 root-specific miRNAs were detected; while under -Pi conditions, 10 leaf- and 4 root-specific miRNAs were identified. Collectively, 25 miRNAs were induced and 11 miRNAs were repressed by Pi starvation in soybean. Then, stem-loop real-time PCR confirmed expression of four selected P-responsive miRNAs, and RLM-5' RACE confirmed that a PHO2 and GmPT5, a kelch-domain containing protein, and a Myb transcription factor, respectively are targets of miR399, miR2111, and miR159e-3p. Finally, P-responsive cis-elements in the promoter regions of soybean miRNA genes were analyzed at the genome-wide scale.

CONCLUSIONS

Leaf- and root-specific miRNAs, and P-responsive miRNAs in soybean were identified genome-wide. A total of 154 target genes of miRNAs were predicted via degradome sequencing and computational analyses. The targets of miR399, miR2111, and miR159e-3p were confirmed. Taken together, our study implies the important roles of miRNAs in P signaling and provides clues for deciphering the functions for microRNA/target modules in soybean.

Identification and Dynamic Regulation of microRNAs Involved in Salt Stress Responses in Functional Soybean Nodules by High-Throughput Sequencing

Both symbiosis between legumes and rhizobia and nitrogen fixation in functional nodules are dramatically affected by salt stress. Better understanding of the molecular mechanisms that regulate the salt tolerance of functional nodules is essential for genetic improvement of nitrogen fixation efficiency. microRNAs (miRNAs) have been implicated in stress responses in many plants and in symbiotic nitrogen fixation (SNF) in soybean. However, the dynamic regulation of miRNAs in functioning nodules during salt stress response remains unknown. We performed deep sequencing of miRNAs to understand the miRNA expression profile in normal or salt stressed-soybean mature nodules. We identified 110 known miRNAs belonging to 61 miRNA families and 128 novel miRNAs belonging to 64 miRNA families. Among them, 104 miRNAs were dramatically differentially expressed (>2-fold or detected only in one library) during salt stress. qRT-PCR analysis of eight miRNAs confirmed that these miRNAs were dynamically regulated in response to salt stress in functional soybean nodules. These data significantly increase the number of miRNAs known to be expressed in soybean nodules, and revealed for the first time a dynamic regulation of miRNAs during salt stress in functional nodules. The findings suggest great potential for miRNAs in functional soybean nodules during salt stress.

Small RNA pathways and diversity in model legumes: lessons from genomics

Small non-coding RNAs (smRNA) participate in the regulation of development, cell differentiation, adaptation to environmental constraints and defense responses in plants. They negatively regulate gene expression by degrading specific mRNA targets, repressing their translation or modifying chromatin conformation through homologous interaction with target loci. MicroRNAs (miRNA) and short-interfering RNAs (siRNA) are generated from long double stranded RNA (dsRNA) that are cleaved into 20-24-nucleotide dsRNAs by RNase III proteins called DICERs (DCL). One strand of the duplex is then loaded onto effective complexes containing different ARGONAUTE (AGO) proteins. In this review, we explored smRNA diversity in model legumes and compiled available data from miRBAse, the miRNA database, and from 22 reports of smRNA deep sequencing or miRNA identification genome-wide in three legumes: Medicago truncatula, soybean (Glycine max) and Lotus japonicus. In addition to conserved miRNAs present in other plant species, 229, 179, and 35 novel miRNA families were identified respectively in these 3 legumes, among which several seems legume-specific. New potential functions of several miRNAs in the legume-specific nodulation process are discussed. Furthermore, a new category of siRNA, the phased siRNAs, which seems to mainly regulate disease-resistance genes, was recently discovered in legumes. Despite that the genome sequence of model legumes are not yet fully completed, further analysis was performed by database mining of gene families and protein characteristics of DCLs and AGOs in these genomes. Although most components of the smRNA pathways are conserved, identifiable homologs of key smRNA players from non-legumes, like AGO10 or DCL4, could not yet be detected in M. truncatula available genomic and expressed sequence (EST) databases. In contrast to Arabidopsis, an important gene diversification was observed in the three legume models (for DCL2, AGO4, AGO2, and AGO10) or specifically in soybean for DCL1 and DCL4. Functional significance of these variant isoforms may reflect peculiarities of smRNA biogenesis and functions in legumes.

Identification of wild soybean miRNAs and their target genes responsive to aluminum stress

BACKGROUND

MicroRNAs (miRNAs) play important regulatory roles in development and stress response in plants. Wild soybean (Glycine soja) has undergone long-term natural selection and may have evolved special mechanisms to survive stress conditions as a result. However, little information about miRNAs especially miRNAs responsive to aluminum (Al) stress is available in wild soybean.

RESULTS

Two small RNA libraries and two degradome libraries were constructed from the roots of Al-treated and Al-free G. soja seedlings. For miRNA identification, a total of 7,287,655 and 7,035,914 clean reads in Al-treated and Al-free small RNAs libraries, respectively, were generated, and 97 known miRNAs and 31 novel miRNAs were identified. In addition, 49 p3 or p5 strands of known miRNAs were found. Among all the identified miRNAs, the expressions of 30 miRNAs were responsive to Al stress. Through degradome sequencing, 86 genes were identified as targets of the known miRNAs and five genes were found to be the targets of the novel miRNAs obtained in this study. Gene ontology (GO) annotations of target transcripts indicated that 52 target genes cleaved by conserved miRNA families might play roles in the regulation of transcription. Additionally, some genes, such as those for the auxin response factor (ARF), domain-containing disease resistance protein (NB-ARC), leucine-rich repeat and toll/interleukin-1 receptor-like protein (LRR-TIR) domain protein, cation transporting ATPase, Myb transcription factors, and the no apical meristem (NAM) protein, that are known to be responsive to stress, were found to be cleaved under Al stress conditions.

CONCLUSIONS

A number of miRNAs and their targets were detected in wild soybean. Some of them that were responsive to biotic and abiotic stresses were regulated by Al stress. These findings provide valuable information to understand the function of miRNAs in Al tolerance.

Advances in genomics for flatfish aquaculture

Fish aquaculture is considered to be one of the most sustainable sources of protein for humans. Many different species are cultured worldwide, but among them, marine flatfishes comprise a group of teleosts of high commercial interest because of their highly prized white flesh. However, the aquaculture of these fishes is seriously hampered by the scarce knowledge on their biology. In recent years, various experimental 'omics' approaches have been applied to farmed flatfishes to increment the genomic resources available. These tools are beginning to identify genetic markers associated with traits of commercial interest, and to unravel the molecular basis of different physiological processes. This article summarizes recent advances in flatfish genomics research in Europe. We focus on the new generation sequencing technologies, which can produce a massive amount of DNA sequencing data, and discuss their potentials and applications for de novo genome sequencing and transcriptome analysis. The relevance of these methods in nutrigenomics and foodomics approaches for the production of healthy animals, as well as high quality and safety products for the consumer, is also briefly discussed.

Expression and tissue-specific localization of nitrate-responsive miRNAs in roots of maize seedlings

Nitrogen availability seriously affects crop productivity and environment. The knowledge of post-transcriptional regulation of plant response to nutrients is important to improve nitrogen use efficiency of crop. This research was aimed at understanding the role of miRNAs in the molecular control of plant response to nitrate. The expression profiles of six mature miRNAs were deeply studied by quantitative real time polymerase chain reaction and in situ hybridization (ISH). To this aim, a novel optimized protocol was set up for the use of digoxygenin-labelled Zip Nucleic Acid-modified oligonucleotides as probes for ISH. Significant differences in miRNAs' transcripts accumulation were evidenced between nitrate-supplied and nitrate-depleted roots. Real-time PCR analyses and in situ detection of miRNA confirmed the array data and allowed us to evidence distinct miRNAs spatio-temporal expression patterns in maize roots. Our results suggest that a prolonged nitrate depletion may induce post-transcriptionally the expression of target genes by repressing the transcription of specific miRNAs. In particular, the repression of the transcription of miR528a/b, miR528a*/b*, miR169i/j/k, miR169i*/j*/k*, miR166j/k/n and miR408/b upon nitrate shortage could represent a crucial step integrating nitrate signals into developmental changes in maize roots.

Genome organization and characteristics of soybean microRNAs

BACKGROUND

microRNAs (miRNAs) are key regulators of gene expression and play important roles in many aspects of plant biology. The role(s) of miRNAs in nitrogen-fixing root nodules of leguminous plants such as soybean is not well understood. We examined a library of small RNAs from Bradyrhizobium japonicum-inoculated soybean roots and identified novel miRNAs. In order to enhance our understanding of miRNA evolution, diversification and function, we classified all known soybean miRNAs based on their phylogenetic conservation (conserved, legume- and soybean-specific miRNAs) and examined their genome organization, family characteristics and target diversity. We predicted targets of these miRNAs and experimentally validated several of them. We also examined organ-specific expression of selected miRNAs and their targets.

RESULTS

We identified 120 previously unknown miRNA genes from soybean including 5 novel miRNA families. In the soybean genome, genes encoding miRNAs are primarily intergenic and a small percentage were intragenic or less than 1000 bp from a protein-coding gene, suggesting potential co-regulation between the miRNA and its parent gene. Difference in number and orientation of tandemly duplicated miRNA genes between orthologous genomic loci indicated continuous evolution and diversification. Conserved miRNA families are often larger in size and produce less diverse mature miRNAs than legume- and soybean-specific families. In addition, the majority of conserved and legume-specific miRNA families produce 21 nt long mature miRNAs with distinct nucleotide distribution and regulate a more conserved set of target mRNAs compared to soybean-specific families. A set of nodule-specific target mRNAs and their cognate regulatory miRNAs had inverse expression between root and nodule tissues suggesting that spatial restriction of target gene transcripts by miRNAs might govern nodule-specific gene expression in soybean.

CONCLUSIONS

Genome organization of soybean miRNAs suggests that they are actively evolving. Distinct family characteristics of soybean miRNAs suggest continuous diversification of function. Inverse organ-specific expression between selected miRNAs and their targets in the roots and nodules, suggested a potential role for these miRNAs in regulating nodule development.

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

Background

MicroRNAs (miRNAs) are endogenously encoded small RNAs that post-transcriptionally regulate gene expression. MiRNAs play essential roles in almost all plant biological processes. Currently, few miRNAs have been identified in the model food legume Phaseolus vulgaris (common bean). Recent advances in next generation sequencing technologies have allowed the identification of conserved and novel miRNAs in many plant species. Here, we used Illumina's sequencing by synthesis (SBS) technology to identify and characterize the miRNA population of Phaseolus vulgaris.

Results

Small RNA libraries were generated from roots, flowers, leaves, and seedlings of P. vulgaris. Based on similarity to previously reported plant miRNAs,114 miRNAs belonging to 33 conserved miRNA families were identified. Stem-loop precursors and target gene sequences for several conserved common bean miRNAs were determined from publicly available databases. Less conserved miRNA families and species-specific common bean miRNA isoforms were also characterized. Moreover, novel miRNAs based on the small RNAs were found and their potential precursors were predicted. In addition, new target candidates for novel and conserved miRNAs were proposed. Finally, we studied organ-specific miRNA family expression levels through miRNA read frequencies.

Conclusions

This work represents the first massive-scale RNA sequencing study performed in Phaseolus vulgaris to identify and characterize its miRNA population. It significantly increases the number of miRNAs, precursors, and targets identified in this agronomically important species. The miRNA expression analysis provides a foundation for understanding common bean miRNA organ-specific expression patterns. The present study offers an expanded picture of P. vulgaris miRNAs in relation to those of other legumes.

In situ localization of small RNAs in plants by using LNA probes

Small RNAs have crucial roles in numerous aspects of plant biology. Despite our current understanding of their biogenesis and mechanisms of action, the biological function of small RNAs, particularly miRNAs, remains largely unknown. To decipher small RNA function, knowledge about their spatiotemporal patterns of expression is essential. Here we report an in situ hybridization method for the precise localization of small RNAs in plants by using locked nucleic acid (LNA) oligonucleotide probes. This method has been adapted from protocols used to detect messenger RNAs in formaldehyde-fixed and paraffin-embedded tissue sections, but it includes essential optimizations in key prehybridization, hybridization and posthybridization steps. Most importantly, optimization of probe concentration and hybridization temperature is required for each unique LNA probe. We present the detailed protocol starting from sectioned tissues, and we include troubleshooting tips and recommended controls. This method has been used successfully in several plant species and can be completed within 2-6 d.