Small RNA and Degradome Deep Sequencing Reveals the Roles of microRNAs in Seed Expansion in Peanut (Arachis hypogaea L.)

Omics-driven advances in the understanding of regulatory landscape of peanut seed development

Peanuts (Arachis hypogaea) are an essential oilseed crop known for their unique developmental process, characterized by aerial flowering followed by subterranean fruit development. This crop is polyploid, consisting of A and B subgenomes, which complicates its genetic analysis. The advent and progression of omics technologies-encompassing genomics, transcriptomics, proteomics, epigenomics, and metabolomics-have significantly advanced our understanding of peanut biology, particularly in the context of seed development and the regulation of seed-associated traits. Following the completion of the peanut reference genome, research has utilized omics data to elucidate the quantitative trait loci (QTL) associated with seed weight, oil content, protein content, fatty acid composition, sucrose content, and seed coat color as well as the regulatory mechanisms governing seed development. This review aims to summarize the advancements in peanut seed development regulation and trait analysis based on reference genome-guided omics studies. It provides an overview of the significant progress made in understanding the molecular basis of peanut seed development, offering insights into the complex genetic and epigenetic mechanisms that influence key agronomic traits. These studies highlight the significance of omics data in profoundly elucidating the regulatory mechanisms of peanut seed development. Furthermore, they lay a foundational basis for future research on trait-related functional genes, highlighting the pivotal role of comprehensive genomic analysis in advancing our understanding of plant biology.

The genus Arachis: an excellent resource for studies on differential gene expression for stress tolerance

Peanut Arachis hypogaea is a segmental allotetraploid in the section Arachis of the genus Arachis along with the Section Rhizomataceae. Section Arachis has several diploid species along with Arachis hypogaea and A. monticola. The section Rhizomataceae comprises polyploid species. Several species in the genus are highly tolerant to biotic and abiotic stresses and provide excellent sets of genotypes for studies on differential gene expression. Though there were several studies in this direction, more studies are needed to identify more and more gene combinations. Next generation RNA-seq based differential gene expression study is a powerful tool to identify the genes and regulatory pathways involved in stress tolerance. Transcriptomic and proteomic study of peanut plants under biotic stresses reveals a number of differentially expressed genes such as R genes (NBS-LRR, LRR-RLK, protein kinases, MAP kinases), pathogenesis related proteins (PR1, PR2, PR5, PR10) and defense related genes (defensin, F-box, glutathione S-transferase) that are the most consistently expressed genes throughout the studies reported so far. In most of the studies on biotic stress induction, the differentially expressed genes involved in the process with enriched pathways showed plant-pathogen interactions, phenylpropanoid biosynthesis, defense and signal transduction. Differential gene expression studies in response to abiotic stresses, reported the most commonly expressed genes are transcription factors (MYB, WRKY, NAC, bZIP, bHLH, AP2/ERF), LEA proteins, chitinase, aquaporins, F-box, cytochrome p450 and ROS scavenging enzymes. These differentially expressed genes are in enriched pathways of transcription regulation, starch and sucrose metabolism, signal transduction and biosynthesis of unsaturated fatty acids. These identified differentially expressed genes provide a better understanding of the resistance/tolerance mechanism, and the genes for manipulating biotic and abiotic stress tolerance in peanut and other crop plants. There are a number of differentially expressed genes during biotic and abiotic stresses were successfully characterized in peanut or model plants (tobacco or Arabidopsis) by genetic manipulation to develop stress tolerance plants, which have been detailed out in this review and more concerted studies are needed to identify more and more gene/gene combinations.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Integrated agronomic, physiological, microstructure, and whole-transcriptome analyses reveal the role of biomass accumulation and quality formation during Se biofortification in alfalfa

Se-biofortified agricultural products receive considerable interest due to the worldwide severity of selenium (Se) deficiency. Alfalfa (Medicago sativa L.), the king of forage, has a large biomass, a high protein content, and a high level of adaptability, making it a good resource for Se biofortification. Analyses of agronomic, quality, physiological, and microstructure results indicated the mechanism of biomass increase and quality development in alfalfa during Se treatment. Se treatment effectively increased Se content, biomass accumulation, and protein levels in alfalfa. The enhancement of antioxidant capacity contributes to the maintenance of low levels of reactive oxygen species (ROS), which, in turn, serves to increase alfalfa's stress resistance and the stability of its intracellular environment. An increase in the rate of photosynthesis contributes to the accumulation of biomass in alfalfa. To conduct a more comprehensive investigation of the regulatory networks induced by Se treatment, the transcriptome sequencing of non-coding RNA (ncRNA) was employed to compare 100 mg/kg Se treatment and control groups. The analysis identified 1,414, 62, and 5 genes as DE-long non-coding RNAs (DE-lncRNA), DE-microRNAs (DE-miRNA), and DE-circular RNA (DE-circRNA), respectively. The function of miRNA-related regulatory networks during Se biofortification in alfalfa was investigated. Subsequent enrichment analysis revealed significant involvement of transcription factors, DNA replication and repair mechanisms, photosynthesis, carbohydrate metabolism, and protein processing. The antioxidant capacity and protein accumulation of alfalfa were regulated by the modulation of signal transduction, the glyoxalase pathway, proteostasis, and circRNA/lncRNA-related regulatory networks. The findings offer new perspectives on the regulatory mechanisms of Se in plant growth, biomass accumulation, and stress responses, and propose potential strategies for enhancing its utilization in the agricultural sector.

Characterization of PsmiR319 during flower development in early- and late-flowering tree peonies cultivars

The flowering period is the most important ornamental trait of tree peony, while industrial development of tree peony has been limited by short flowering period. miR319 plays an important regulatory role in plant flowering. In the current study, the expression characteristics and evolution of PsmiR319 in tree peony flowering was explored using 'Feng Dan' and 'Lian He', which are early-flowering and late-flowering varieties of tree peony, respectively. The structure, evolution, and target(s) of PsmiR319 were analyzed by bioinformatics. Evolution analysis showed that pre-PsmiR319 was distributed in 41 plant species, among which the length of the precursor sequence exhibited marked differences (between 52 and 308 bp). Pre-PsmiR319 of tree peony was located close to the corresponding sequences of Linum usitatissimum and Picea abies in the phylogenetic tree, and in addition, could form a typical hairpin structure including a mature body with a length of 20 bp located on the 3p arm and part of the loop sequence. The mature sequence of miR319 was highly conserved among different species. Target genes of PsmiR319 include MYB-related transcription factor in tree peony. Expression of PsmiR319, assayed by qRT-PCR, differed between 'Feng Dan' and 'Lian He' during different flower development periods. PsmiR319 and its target gene showed a negative expression regulation relationship during the periods of CE (color exposure), BS (blooming stage), IF (initial flowering), and HO (half opening) in the early-flowering 'Feng Dan', and the same in FB (Full blooming) periods of late-flowering 'Lian He'. Findings from this study provide a reference for further investigation into the mechanism of miR319 in the development of different varieties of tree peony.

Degradome Assisted Plant MicroRNA Prediction Under Alternative Annotation Criteria

Current microRNA (miRNA) prediction methods are generally based on annotation criteria that tend to miss potential functional miRNAs. Recently, new miRNA annotation criteria have been proposed that could lead to improvements in miRNA prediction methods in plants. Here, we investigate the effect of the new criteria on miRNA prediction in Arabidopsis thaliana and present a new degradome assisted functional miRNA prediction approach. We investigated the effect by applying the new criteria, and a more permissive criteria on miRNA prediction using existing miRNA prediction tools. We also developed an approach to miRNA prediction that is assisted by the functional information extracted from the analysis of degradome sequencing. We demonstrate the improved performance of degradome assisted miRNA prediction compared to unassisted prediction and evaluate the approach using miRNA differential expression analysis. We observe how the miRNA predictions fit under the different criteria and show a potential novel miRNA that has been missed within Arabidopsis thaliana. Additionally, we introduce a freely available software 'PAREfirst' that employs the degradome assisted approach. The study shows that some miRNAs could be missed due to the stringency of the former annotation criteria, and combining a degradome assisted approach with more permissive miRNA criteria can expand confident miRNA predictions.

Integration of Small RNA and Degradome Sequencing Reveals the Regulatory Network of Al-Induced Programmed Cell Death in Peanut

Peanut is one of the most important oil crops in the world. In China, the peanut is highly produced in its southern part, in which the arable land is dominated by acid soil. At present, miRNAs have been identified in stress response, but their roles and mechanisms are not clear, and no miRNA studies have been found related to aluminum (Al)-induced programmed cell death (PCD). In the present study, transcriptomics, sRNAs, and degradome analysis in the root tips of two peanut cultivars ZH2 (Al-sensitive, S) and 99-1507 (Al-tolerant, T) were carried out. Here, we generated a comprehensive resource focused on identifying key regulatory miRNA-target circuits that regulate PCD under Al stress. Through deep sequencing, 2284 miRNAs were identified and 147 miRNAs were differentially expressed under Al stress. Furthermore, 19237 target genes of 749 miRNAs were validated by degradome sequencing. GO and KEGG analyses of differential miRNA targets showed that the pathways of synthesis and degradation of ketone bodies, citrate cycle (TCA cycle), and peroxisome were responded to Al stress. The combined analysis of the degradome data sets revealed 89 miRNA-mRNA interactions that may regulate PCD under Al stress. Ubiquitination may be involved in Al-induced PCD in peanut. The regulatory networks were constructed based on the differentially expressed miRNAs and their targets related to PCD. Our results will provide a useful platform to research on PCD induced by Al and new insights into the genetic engineering for plant stress response.

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.

Integrated microRNA and transcriptome profiling reveal key miRNA-mRNA interaction pairs associated with seed development in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat seed development is an extremely complex process involving many gene regulatory pathways. MicroRNAs (miRNAs) have been identified as the important negative regulators of gene expression and performed crucial regulatory roles in various plant biological processes. However, whether miRNAs participate in Tartary buckwheat seed development remains unexplored.

RESULTS

In this study, we first identified 26 miRNA biosynthesis genes in the Tartary buckwheat genome and described their phylogeny and expression profiling. Then we performed small RNA (sRNA) sequencing for Tartary buckwheat seeds at three developmental stages to identify the miRNAs associated with seed development. In total, 230 miRNAs, including 101 conserved and 129 novel miRNAs, were first identified in Tartary buckwheat, and 3268 target genes were successfully predicted. Among these miRNAs, 76 exhibited differential expression during seed development, and 1534 target genes which correspond to 74 differentially expressed miRNAs (DEMs) were identified. Based on integrated analysis of DEMs and their targets expression, 65 miRNA-mRNA interaction pairs (25 DEMs corresponding to 65 target genes) were identified that exhibited significantly opposite expression during Tartary buckwheat seed development, and 6 of the miRNA-mRNA pairs were further verified by quantitative real-time polymerase chain reaction (qRT-PCR) and ligase-mediated rapid amplification of 5' cDNA ends (5'-RLM-RACE). Functional annotation of the 65 target mRNAs showed that 56 miRNA-mRNA interaction pairs major involved in cell differentiation and proliferation, cell elongation, hormones response, organogenesis, embryo and endosperm development, seed size, mineral elements transport, and flavonoid biosynthesis, which indicated that they are the key miRNA-mRNA pairs for Tartary buckwheat seed development.

CONCLUSIONS

Our findings provided insights for the first time into miRNA-mediated regulatory pathways in Tartary buckwheat seed development and suggested that miRNAs play important role in Tartary buckwheat seed development. These findings will be help to study the roles and regulatory mechanism of miRNAs in Tartary buckwheat seed development.

MicroRNA-mediated regulation of agronomically important seed traits: a treasure trove with shades of grey!

Seed development is an intricate process with multiple levels of regulation. MicroRNAs (miRNAs) have emerged as one of the crucial components of molecular networks underlying agronomically important seed traits in diverse plant species. In fact, loss of function of the genes regulating miRNA biogenesis also exhibits defects in seed development. A total of 21 different miRNAs have experimentally been shown to regulate seed size, nutritional content, vigor, and shattering, and have been reviewed here. The mechanism details of the associated regulatory cascades mediated through transcriptional regulators, phytohormones, basic metabolic machinery, and secondary siRNAs are elaborated. Co-localization of miRNAs and their target regions with seed-related QTLs provides new avenues for engineering these traits using conventional breeding programs or biotechnological interventions. While global analysis of miRNAs using small RNA sequencing studies are expanding the repertoire of candidate miRNAs, recent revelations on their inheritance, transport, and mechanism of action would be instrumental in designing better strategies for optimizing agronomically relevant seed traits.

Genome-wide profiling of miRNAs during seed development reveals their functional relevance in seed size/weight determination in chickpea

MicroRNAs (miRNAs) are non-coding small RNAs that regulate gene expression at transcriptional and post-transcriptional levels. The role of miRNAs in seed development and seed size/weight determination is poorly understood in legumes. In this study, we profiled miRNAs at seven successive stages of seed development in a small-seeded and a large-seeded chickpea cultivar via small RNA sequencing. In total, 113 known and 243 novel miRNAs were identified. Gene ontology analysis revealed the enrichment of seed/reproductive/post-embryonic development and signaling pathways processes among the miRNA target genes. A large fraction of the target genes exhibited antagonistic correlation with miRNA expression. The sets of co-expressed miRNAs showing differential expression between the two cultivars were recognized. Known transcription factor (TF) encoding genes involved in seed size/weight determination, including SPL, GRF, MYB, ARF, HAIKU1, SHB1, KLUH/CYP78A5, and E2Fb along with novel genes were found to be targeted by the predicted miRNAs. Differential expression analysis revealed higher transcript levels of members of SPL and REVOLUTA TF families and lower expression of their corresponding miRNAs in the large-seeded cultivar. At least 19 miRNAs known to be involved in seed development or differentially expressed between small-seeded and large-seeded cultivars at late-embryogenesis and/or mid-maturation stages were located within known quantitative trait loci (QTLs) associated with seed size/weight determination. Moreover, 41 target genes of these miRNAs were also located within these QTLs. Altogether, we revealed important roles of miRNAs in seed development and identified candidate miRNAs and their target genes that have functional relevance in determining seed size/weight in chickpea.

Identification and molecular characterization of miRNAs and their target genes associated with seed development through small RNA sequencing in chickpea

Multiple studies have attempted to dissect the molecular mechanism underlying seed development in chickpea (Cicer arietinum L.). These studies highlight the need to focus on the role of miRNAs in regulating storage protein accumulation in seeds. Therefore, a total of 8,856,691 short-read sequences were generated from a small RNA library of developing chickpea seeds and were analyzed using miRDeep-P to identify 74 known and 26 novel miRNA sequences. Known miRNAs were classified into 22 miRNA families with miRNA156 family being most abundant. Of the 26 putative novel miRNAs identified, only 22 could be experimentally validated using stem loop end point PCR. Differential expression analyses led to the identification of known as well as novel miRNAs that could regulate various stages of chickpea seed development. In silico target prediction revealed several important target genes and transcription factors like SPL, mediator of RNA Polymerase II transcription subunit 12, aspartic proteinase and NACs, which were further validated by real-time PCR analysis. A comparative expression analysis in chickpea genotypes with contrasting seed protein content revealed one known (Car-miR156h) and two novel miRNA (Car-novmiR7 and Car-novmiR23) candidates to be highly expressed in the LPC (low protein content) chickpea genotypes, targets of which are known to regulate seed storage protein accumulation. Therefore, this study provides a useful resource in the form of miRNA and their targets which can be further utilized to understand and manipulate various regulatory mechanisms involved in seed development with the overall aim of improving yield and nutrition attributes in chickpea.

Integrated Analysis of MicroRNA and Target Genes in Brachypodium distachyon Infected by Magnaporthe oryzae by Small RNA and Degradome Sequencing

Rice blast caused by Magnaporthe oryzae is one of the most important diseases that seriously threaten rice production. Brachypodium distachyon is a grass species closely related to grain crops, such as rice, barley, and wheat, and has become a new model plant of Gramineae. In this study, 15 small RNA samples were sequenced to examine the dynamic changes in microRNA (miRNA) expression in B. distachyon infected by M. oryzae at 0, 24, and 48 h after inoculation. We identified 432 conserved miRNAs and 288 predicted candidate miRNAs in B. distachyon. Additionally, there were 7 and 19 differentially expressed miRNAs at 24 and 48 h post-inoculation, respectively. Furthermore, using degradome sequencing, we identified 2,126 genes as targets for 308 miRNAs; using quantitative real-time PCR (qRT-PCR), we validated five miRNA/target regulatory units involved in B. distachyon-M. oryzae interactions. Moreover, using co-transformation technology, we demonstrated that BdNAC21 was negatively regulated by miR164c. This study provides a new approach for identifying resistance genes in B. distachyon by mining the miRNA regulatory network of host-pathogen interactions.

Integration of Small RNA and Transcriptome Sequencing Reveal the Roles of miR395 and ATP Sulfurylase in Developing Seeds of Chinese Kale

Seed development is closely related to plant production and reproduction, and MicroRNAs (miRNA) is widely involved in plant development including seed development. Chinese kale, as a Brassicaceae vegetable, mainly depends on seed for proper reproduction. In the present study, Chinese kale seed and silique at different stages were selected to establish small RNA (sRNA) libraries including silique wall sRNA libraries at torpedo-embryo stage (PC), silique wall sRNA libraries at cotyledonary-embryo stage (PD), seed sRNA libraries at torpedo-embryo stage (SC), and seed sRNA libraries at cotyledonary-embryo stage (SD). The results showed that miRNA expressed differentially in the seeds and corresponding siliques at different stages. To further clarify the functional mode of miRNA in the process of seed development, Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis was performed on target genes of the differentially expressed miRNAs, and these target genes were mainly enriched in plant hormone signal transduction, primary and secondary metabolic pathways. After joint analysis with the transcriptome change of the corresponding period, miR156-SPL10/SPL11, miR395-APS3, and miR397-LAC2/LAC11 modules were identified to be directly involved in the development of Chinese kale seeds. What's more, modified 5'RLM-RACE and Agrobacteria-mediated Chinese kale transient transformation suggest miR395b_2 is involved in sulfur metabolism during seed development by regulating its target gene APS3.

Alternative splicing profiling provides insights into the molecular mechanisms of peanut peg development

BACKGROUND

The cultivated peanut (Arachis hypogaea) is one of the most important oilseed crops worldwide, and the generation of pegs and formation of subterranean pods are essential processes in peanut reproductive development. However, little information has been reported about alternative splicing (AS) in peanut peg formation and development.

RESULTS

Herein, we presented a comprehensive full-length (FL) transcriptome profiling of AS isoforms during peanut peg and early pod development. We identified 1448, 1102, 832, and 902 specific spliced transcripts in aerial pegs, subterranean pegs, subterranean unswollen pegs, and early swelling pods, respectively. A total of 184 spliced transcripts related to gravity stimulation, light and mechanical response, hormone mediated signaling pathways, and calcium-dependent proteins were identified as possibly involved in peanut peg development. For aerial pegs, spliced transcripts we got were mainly involved in gravity stimulation and cell wall morphogenetic processes. The genes undergoing AS in subterranean peg were possibly involved in gravity stimulation, cell wall morphogenetic processes, and abiotic response. For subterranean unswollen pegs, spliced transcripts were predominantly related to the embryo development and root formation. The genes undergoing splice in early swelling pods were mainly related to ovule development, root hair cells enlargement, root apex division, and seed germination.

CONCLUSION

This study provides evidence that multiple genes are related to gravity stimulation, light and mechanical response, hormone mediated signaling pathways, and calcium-dependent proteins undergoing AS express development-specific spliced isoforms or exhibit an obvious isoform switch during the peanut peg development. AS isoforms in subterranean pegs and pods provides valuable sources to further understand post-transcriptional regulatory mechanisms of AS in the generation of pegs and formation of subterranean pods.

Analysis of small RNA populations generated in peanut leaves after exogenous application of dsRNA and dsDNA targeting aflatoxin synthesis genes

Previously, we have shown that RNA interference (RNAi) can prevent aflatoxin accumulation in transformed peanuts. To explore aflatoxin control by exogenous delivery of double-strand RNA (dsRNA) it is necessary to understand the generation of small RNA (sRNA) populations. We sequenced 12 duplicate sRNA libraries of in-vitro-grown peanut plants, 24 and 48 h after exogenous application of five gene fragments (RNAi-5x) related to aflatoxin biosynthesis in Aspergillus flavus. RNAi-5x was applied either as double-stranded RNA (dsRNA) or RNAi plasmid DNA (dsDNA). Small interfering RNAs (siRNAs) derived from RNAi-5x were significantly more abundant at 48 h than at 24 h, and the majority mapped to the fragment of aflatoxin efflux-pump gene. RNAi-5x-specific siRNAs were significantly, three to fivefold, more abundant in dsDNA than dsRNA treatments. Further examination of known micro RNAs related to disease-resistance, showed significant down-regulation of miR399 and up-regulation of miR482 in leaves treated with dsDNA compared to the control. These results show that sRNA sequencing is useful to compare exogenous RNAi delivery methods on peanut plants, and to analyze the efficacy of molecular constructs to generate siRNAs against specific gene targets. This work lays the foundation for non-transgenic delivery of RNAi in controlling aflatoxins in peanut.

Integrated small RNA and Degradome sequencing provide insights into salt tolerance in sesame (Sesamum indicum L.)

BACKGROUND

MicroRNAs (miRNAs) exhibit important regulatory roles in the response to abiotic stresses by post-transcriptionally regulating the target gene expression in plants. However, their functions in sesame response to salt stress are poorly known. To dissect the complex mechanisms underlying salt stress response in sesame, miRNAs and their targets were identified from two contrasting sesame genotypes by a combined analysis of small RNAs and degradome sequencing.

RESULTS

A total of 351 previously known and 91 novel miRNAs were identified from 18 sesame libraries. Comparison of miRNA expressions between salt-treated and control groups revealed that 116 miRNAs were involved in salt stress response. Using degradome sequencing, potential target genes for some miRNAs were also identified. The combined analysis of all the differentially expressed miRNAs and their targets identified miRNA-mRNA regulatory networks and 21 miRNA-mRNA interaction pairs that exhibited contrasting expressions in sesame under salt stress.

CONCLUSIONS

This comprehensive integrated analysis may provide new insights into the genetic regulation mechanism of miRNAs underlying the adaptation of sesame to salt stress.

Genome-wide identification and analysis of long noncoding RNAs (lncRNAs) during seed development in peanut (Arachis hypogaea L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) have several known functions involving various biological regulatory processes in plant. However, the possible roles of lncRNAs during peanut seed development have not been fully explored.

RESULTS

In this study, two peanut recombinant inbred lines (RIL₈) that differ in seed size were used to investigate comprehensive lncRNA profiles derived from the seed development at 15 and 35 days after flowering (DAF). We identified a total of 9388 known and 4037 novel lncRNAs, from which 1437 were differentially expressed lncRNAs. Interestingly, the expression patterns of a number of lncRNAs can be very different between two closely related inbred lines and these lncRNAs were expressed predominantly in only one RIL at 35 DAF. Some differentially expressed lncRNAs were found related to putative cis-acting target genes and predicted to be involved in transcription, transport, cell division, and plant hormone biosynthesis. The expression patterns of several representative lncRNAs and 12 protein-coding genes were validated by qPCR. Same expression pattern was observed between most lncRNAs and their target genes. 11 lncRNAs, XR_001593099.1, MSTRG.18462.1, MSTRG.34915.1, MSTRG.41848.1, MSTRG.22884.1, MSTRG.12404.1, MSTRG.26719.1, MSTRG.35761.1, MSTRG.20033.1, MSTRG.13500.1, and MSTRG.9304.1 and their cis-acting target genes may play key roles in peanut seed development.

CONCLUSIONS

These results provided new information on lncRNA-mediated regulatory roles in peanut seed development, contributing to the comprehensive understanding of the molecular mechanisms involved in peanut seed development.

Identification of microRNAs and their Endonucleolytic Cleavaged target mRNAs in colorectal cancer

Background

Colorectal cancer (CRC) ranks the third among the most common malignancies globally. It is well known that microRNAs (miRNAs) play vital roles in destabilizing mRNAs and repressing their translations in this disease. However, the mechanism of miRNA-induced mRNA cleavage remains to be investigated.

Method

In this study, high-throughput small RNA (sRNA) sequencing was utilized to identify and profile miRNAs from six pairs of colorectal cancer tissues (CTs) and adjacent tissues (CNs). Degradome sequencing (DS) was employed to detect the cleaved target genes. The Database for Annotation, Visualization and Integrated Discovery (DAVID) software was used for GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis.

Results

In total, 1278 known miRNAs (clustered into 337 families) and 131 novel miRNAs were characterized in the CT and CN libraries, respectively. Of those, 420 known and eight novel miRNAs were defined as differentially expressed miRNAs (DEmiRNAs) by comparing the expression levels observed in the CT and CN libraries. Furthermore, through DS, 9685 and 202 potential target transcripts were characterized as target genes for 268 known and 33 novel miRNAs, respectively. It was further predicted that a total of 264 targeted genes for the 85 DEmiRNAs are involved in proteoglycans in cancer and the AMP-activated protein kinase signaling pathway. After systemic analysis of prognosis-related miRNA targets in those cancer-related signal pathways, we found that two targets ezrin (EZR) and hematopoietic cell-specific Lyn substrate 1 (HCLS1) had the potential prognostic characteristics with CRC regarding over survival (OS) or recurrence.

Conclusion

In total, we found that endonucleolytic miRNA-directed mRNA cleavage occurs in CRC. A number of potential genes targeted by CRC-related miRNAs were identified and some may have the potential as prognosis markers of CRC. The present findings may lead to an improved better appreciation of the novel interaction mode between miRNAs and target genes in CRC.

Peg Biology: Deciphering the Molecular Regulations Involved During Peanut Peg Development

Peanut or groundnut is one of the most important legume crops with high protein and oil content. The high nutritional qualities of peanut and its multiple usage have made it an indispensable component of our daily life, in both confectionary and therapeutic food industries. Given the socio-economic significance of peanut, understanding its developmental biology is important in providing a molecular framework to support breeding activities. In peanut, the formation and directional growth of a specialized reproductive organ called a peg, or gynophore, is especially relevant in genetic improvement. Several studies have indicated that peanut yield can be improved by improving reproductive traits including peg development. Therefore, we aim to identify unifying principles for the genetic control, underpinning molecular and physiological basis of peg development for devising appropriate strategy for peg improvement. This review discusses the current understanding of the molecular aspects of peanut peg development citing several studies explaining the key mechanisms. Deciphering and integrating recent transcriptomic, proteomic, and miRNA-regulomic studies provide a new perspective for understanding the regulatory events of peg development that participate in pod formation and thus control yield.

Genome-Wide Analysis of the Growth-Regulating Factor Family in Peanut (Arachis hypogaea L.)

Growth-regulating factors (GRFs) are plant-specific transcription factors that perform important functions in plant growth and development. Herein, we identified and characterised 24 AhGRF genes in peanut (Arachis hypogaea). AhGRF family genes were divided into six classes with OLQ and WRC domains. Transcriptome expression profile showed that more AhGRF genes, such as AhGRF5a gene, were at higher expression during pod development in Arachis monticola than cultivated species, especially at the pod rapid-expansion stage. AhGRF5a and AhGRF5b genes expressed at higher levels in pods than roots, leaves and stems tissues, existing in the difference between Arachis monticola and H8107. Exogenous GA3 application can activate AhGRF5a and AhGRF5b genes and H8107 line showed more positive response than Arachis monticola species. These results imply that these two AhGRF genes may be active during the peanut pod development.

Genome-wide identification of circular RNAs in peanut (Arachis hypogaea L.)

Background

Circular RNAs (circRNAs), a class of widely expressed endogenous regulatory RNAs, are involved in diverse physiological and developmental processes in eukaryotic cells. However, there have been no related studies on the number of circRNAs and their overall characteristics including circRNA abundance and expression profiles in peanut, which is one of the most important edible oil seed crops in the world.

Results

We performed a genome-wide identification of circular RNAs using ribosomal-depleted RNA-sequencing from the seeds of two peanut eighth-generation recombinant inbred lines (RIL8): ‘RIL 8106’ (a medium-pod variety) and ‘RIL 8107’ (a super-pod variety), at 15 and 35 days after flowering (DAF), respectively. A total of 347 circRNA candidates were detected by two computational pipelines: CIRCexplorer and CIRI, with at least two supporting junction reads. All these circRNAs were generated from exons of annotated genes, and widespread on the 20 peanut chromosomes. The expression profiles revealed that circRNAs were differentially expressed between two stages and between two lines. GO enrichment analysis of the host genes produced differentially-expressed circRNAs suggested that circRNAs are involved in seed development and regulation of seed size. Fifteen circRNAs were experimentally analyzed by qRT-PCR with divergent primers, and six circRNAs were resistant to digestion with RNase R exonuclease, and the back-splicing sites were further validated by Sanger DNA sequencing.

Conclusions

We present the first systematical investigation of the genomic characteristics and expression profiles of circRNAs in peanut. The results revealed that circRNAs are abundant and widespread in peanut, and the differentially-expressed circRNAs between two lines suggested that they might play regulatory roles in peanut seeds development.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6020-7) contains supplementary material, which is available to authorized users.

Overexpression of MiR482c in Tomato Induces Enhanced Susceptibility to Late Blight

Tomato is the highest-value fruit/vegetable crop worldwide. However, the quality and yield of tomatoes are severely affected by late blight. MicroRNA482s (miR482s) are involved in the plant's immune system. In this study, miR482c was transiently and stably overexpressed in tomatoes in transgenic plants to explore its mechanism in tomato resistance against late blight. Transgenic tomato plants with transiently overexpressed miR482c displayed a larger lesion area than the control plants upon infection. Furthermore, compared with wild-type (WT) tomato plants, the transgenic tomato plants stably overexpressing miR482c displayed a decreased expression of target genes accompanied by lower peroxidase (POD), superoxide dismutase (SOD), and phenylalanine ammonia-lyase (PAL) activity activities and higher malondialdehyde (MDA) content, thereby leading to a decline in reactive oxygen species (ROS) scavenging ability and aggravating the damage of lipid peroxidation product accumulation on the cell membrane, eventually enhancing plant susceptibility. This finding indicates that miR482c may act as a negative regulator in tomato resistance by regulating nucleotide binding sites and leucine-rich repeat (NBS-LRR) expression levels and ROS levels.

Unravelling the MicroRNA-Mediated Gene Regulation in Developing Pongamia Seeds by High-Throughput Small RNA Profiling

Pongamia (Millettia pinnata syn. Pongamia pinnata) is a multipurpose biofuel tree which can withstand a variety of abiotic stresses. Commercial applications of Pongamia trees may substantially benefit from improvements in their oil-seed productivity, which is governed by complex regulatory mechanisms underlying seed development. MicroRNAs (miRNAs) are important molecular regulators of plant development, while relatively little is known about their roles in seed development, especially for woody plants. In this study, we identified 236 conserved miRNAs within 49 families and 143 novel miRNAs via deep sequencing of Pongamia seeds sampled at three developmental phases. For these miRNAs, 1327 target genes were computationally predicted. Furthermore, 115 differentially expressed miRNAs (DEmiRs) between successive developmental phases were sorted out. The DEmiR-targeted genes were preferentially enriched in the functional categories associated with DNA damage repair and photosynthesis. The combined analyses of expression profiles for DEmiRs and functional annotations for their target genes revealed the involvements of both conserved and novel miRNA-target modules in Pongamia seed development. Quantitative Real-Time PCR validated the expression changes of 15 DEmiRs as well as the opposite expression changes of six targets. These results provide valuable miRNA candidates for further functional characterization and breeding practice in Pongamia and other oilseed plants.

Conserved miR396b-GRF Regulation Is Involved in Abiotic Stress Responses in Pitaya (Hylocereus polyrhizus)

MicroRNA396 (miR396) is a conserved microRNA family that targets growth-regulating factors (GRFs), which play significant roles in plant growth and stress responses. Available evidence justifies the idea that miR396-targeted GRFs have important functions in many plant species; however, no genome-wide analysis of the pitaya (Hylocereus polyrhizus) miR396 gene has yet been reported. Further, its biological functions remain elusive. To uncover the regulatory roles of miR396 and its targets, the hairpin sequence of pitaya miR396b and the open reading frame (ORF) of its target, HpGRF6, were isolated from pitaya. Phylogenetic analysis showed that the precursor miR396b (MIR396b) gene of plants might be clustered into three major groups, and, generally, a more recent evolutionary relationship in the intra-family has been demonstrated. The sequence analysis indicated that the binding site of hpo-miR396b in HpGRF6 is located at the conserved motif which codes the conserved "RSRKPVE" amino acid in the Trp-Arg-Cys (WRC) region. In addition, degradome sequencing analysis confirmed that four GRFs (GRF1, c56908.graph_c0; GRF4, c52862.graph_c0; GRF6, c39378.graph_c0 and GRF9, c54658.graph_c0) are hpo-miR396b targets that are regulated by specific cleavage at the binding site between the 10th and 11th nucleotides from the 5' terminus of hpo-miR396b. Furthermore, quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that hpo-miR396b is down-regulated when confronted with drought stress (15% polyethylene glycol, PEG), and its expression fluctuates under other abiotic stresses, i.e., low temperature (4 ± 1 °C), high temperature (42 ± 1 °C), NaCl (100 mM), and abscisic acid (ABA; 0.38 mM). Conversely, the expression of HpGRF6 showed the opposite trend to exposure to these abiotic stresses. Taken together, hpo-miR396b plays a regulatory role in the control of HpGRF6, which might influence the abiotic stress response of pitaya. This is the first documentation of this role in pitaya and improves the understanding of the molecular mechanisms underlying the tolerance to drought stress in this fruit.