Integrated microRNA and transcriptome profiling reveals a miRNA-mediated regulatory network of embryo abortion under calcium deficiency in peanut (Arachis hypogaea L.)

Characterization of tae-miR156(s) and their response to abiotic stress in wheat (Triticum aestivum L.)

The microRNA156 (miR156) has been widely studied in plants, however, the characterization of the miR156 family of genes in wheat and their expression patterns under abiotic stress are not completely clear. In this study, a total of 20 miR156 family members, referred to as tae-miR156a to tae-miR156t, were identified in wheat with their loci mapped to various chromosomes. These members were divided into five subgroups: miR156a/b/c/d/e/f, miR156g/h/i, miR156j/k, miR156l/m/n/o/p/q, and miR156r/s/t. They were highly conserved during evolution. The prediction of cis-elements in the tae-MIR156(s) promoter region revealed that the tae-MIR156(s) had diverse cis-acting elements; of these, 15 tae-MIR156(s) and 6 tae-MIR156(s) were found to be drought-responsive elements and cold-responsive elements, respectively. And the prediction target genes of tae-miR156(s) are mainly SPL transcription factor genes. Expression analysis based on quantitative real-time polymerase chain reaction (qRT‒PCR) showed that miR156(s) have different expression levels in the various wheat tissues, and the subgroups' response to abiotic stress varied. Among them, miR156g/h/i were strongly induced in the root of cold and heat stress, and miR156a/b/c/d/e/f were significantly increased in roots after drought stress, whereas miR156r/s/t were highly inhibited in leaves and roots after salt stress. These findings imply that tae-miR156(s) are involved in stress response in wheat, and they provide new fundamental knowledge for further analysis of the function of miR156 and its regulatory mechanism in response to abiotic stress.

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

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.

Functional Characterization of Core and Unique Calcite-Dissolving Bacteria Communities from Peanut Fields

Calcium deficiency is a leading cause of reduced peanut (Arachis hypogaea) seed quality and has been linked to increased disease susceptibility, specifically to soilborne fungal pathogens. Sufficient calcium at flowering time is critical to ensure proper pod development. Calcite-dissolving bacteria (CDB) isolated from farming fields can dissolve calcite (CaCO3) on plates and increase soluble calcium levels in soil. However, the phylogenetic diversity and geographic distribution of CDB is unclear. Here, we surveyed soil samples from 15 peanut-producing fields in three regions in southern Georgia, representing distinct soil compositions. We isolated CDB through differentiating media and identified 52 CDB strains. CDB abundance was not associated with any of the soil characteristics we evaluated. Three core genera, represented by 43 strains, were found in all three regions. Paenibacillus was the most common CDB found in all regions, making up 30 of the 52 identified strains. Six genera, represented by eight strains, are unique to one region. Members of the core and unique communities showed comparable solubilization indexes on plates. We conclude that a diversified phylogenetic population of CDB is present in Georgia peanut fields. Despite the phylogenetic diversity, as a population, they exhibit comparable functions in solubilizing calcite on plates.

Comparative genome-wide characterization of salt responsive micro RNA and their targets through integrated small RNA and de novo transcriptome profiling in sugarcane and its wild relative Erianthus arundinaceus

Soil salinity and saline irrigation water are major constraints in sugarcane affecting the production of cane and sugar yield. To understand the salinity induced responses and to identify novel genomic resources, integrated de novo transcriptome and small RNA sequencing in sugarcane wild relative, Erianthus arundinaceus salt tolerant accession IND 99-907 and salt-sensitive sugarcane genotype Co 97010 were performed. A total of 362 known miRNAs belonging to 62 families and 353 miRNAs belonging to 63 families were abundant in IND 99-907 and Co 97010 respectively. The miRNA families such as miR156, miR160, miR166, miR167, miR169, miR171, miR395, miR399, miR437 and miR5568 were the most abundant with more than ten members in both genotypes. The differential expression analysis of miRNA reveals that 221 known miRNAs belonging to 48 families and 130 known miRNAs belonging to 42 families were differentially expressed in IND 99-907 and Co 97010 respectively. A total of 12,693 and 7982 miRNA targets against the monoploid mosaic genome and a total of 15,031 and 12,152 miRNA targets against the de novo transcriptome were identified for differentially expressed known miRNAs of IND 99-907 and Co 97010 respectively. The gene ontology (GO) enrichment analysis of the miRNA targets revealed that 24, 12 and 14 enriched GO terms (FDR < 0.05) for biological process, molecular function and cellular component respectively. These miRNAs have many targets that associated in regulation of biotic and abiotic stresses. Thus, the genomic resources generated through this study are useful for sugarcane crop improvement through biotechnological and advanced breeding approaches.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03867-7.

Genome-Wide Identification and Characterization of CDPK Gene Family in Cultivated Peanut (Arachis hypogaea L.) Reveal Their Potential Roles in Response to Ca Deficiency

This study identified 45 calcium-dependent protein kinase (CDPK) genes in cultivated peanut (Arachis hypogaea L.), which are integral in plant growth, development, and stress responses. These genes, classified into four subgroups based on phylogenetic relationships, are unevenly distributed across all twenty peanut chromosomes. The analysis of the genetic structure of AhCDPKs revealed significant similarity within subgroups, with their expansion primarily driven by whole-genome duplications. The upstream promoter sequences of AhCDPK genes contained 46 cis-acting regulatory elements, associated with various plant responses. Additionally, 13 microRNAs were identified that target 21 AhCDPK genes, suggesting potential post-transcriptional regulation. AhCDPK proteins interacted with respiratory burst oxidase homologs, suggesting their involvement in redox signaling. Gene ontology and KEGG enrichment analyses affirmed AhCDPK genes' roles in calcium ion binding, protein kinase activity, and environmental adaptation. RNA-seq data revealed diverse expression patterns under different stress conditions. Importantly, 26 AhCDPK genes were significantly induced when exposed to Ca deficiency during the pod stage. During the seedling stage, four AhCDPKs (AhCDPK2/-25/-28/-45) in roots peaked after three hours, suggesting early signaling roles in pod Ca nutrition. These findings provide insights into the roles of CDPK genes in plant development and stress responses, offering potential candidates for predicting calcium levels in peanut seeds.

miRNAs for crop improvement

Climate change significantly impacts crop production by inducing several abiotic and biotic stresses. The increasing world population, and their food and industrial demands require focused efforts to improve crop plants to ensure sustainable food production. Among various modern biotechnological tools, microRNAs (miRNAs) are one of the fascinating tools available for crop improvement. miRNAs belong to a class of small non-coding RNAs playing crucial roles in numerous biological processes. miRNAs regulate gene expression by post-transcriptional target mRNA degradation or by translation repression. Plant miRNAs have essential roles in plant development and various biotic and abiotic stress tolerance. In this review, we provide propelling evidence from previous studies conducted around miRNAs and provide a one-stop review of progress made for breeding stress-smart future crop plants. Specifically, we provide a summary of reported miRNAs and their target genes for improvement of plant growth and development, and abiotic and biotic stress tolerance. We also highlight miRNA-mediated engineering for crop improvement and sequence-based technologies available for the identification of miRNAs associated with stress tolerance and plant developmental events.

In-silico identification and characterization of O-methyltransferase gene family in peanut (Arachis hypogaea L.) reveals their putative roles in development and stress tolerance

Cultivated peanut (Arachis hypogaea) is a leading protein and oil-providing crop and food source in many countries. At the same time, it is affected by a number of biotic and abiotic stresses. O-methyltransferases (OMTs) play important roles in secondary metabolism, biotic and abiotic stress tolerance. However, the OMT genes have not been comprehensively analyzed in peanut. In this study, we performed a genome-wide investigation of A. hypogaea OMT genes (AhOMTs). Gene structure, motifs distribution, phylogenetic history, genome collinearity and duplication of AhOMTs were studied in detail. Promoter cis-elements, protein-protein interactions, and micro-RNAs targeting AhOMTs were also predicted. We also comprehensively studied their expression in different tissues and under different stresses. We identified 116 OMT genes in the genome of cultivated peanut. Phylogenetically, AhOMTs were divided into three groups. Tandem and segmental duplication events played a role in the evolution of AhOMTs, and purifying selection pressure drove the duplication process. AhOMT promoters were enriched in several key cis-elements involved in growth and development, hormones, light, and defense-related activities. Micro-RNAs from 12 different families targeted 35 AhOMTs. GO enrichment analysis indicated that AhOMTs are highly enriched in transferase and catalytic activities, cellular metabolic and biosynthesis processes. Transcriptome datasets revealed that AhOMTs possessed varying expression levels in different tissues and under hormones, water, and temperature stress. Expression profiling based on qRT-PCR results also supported the transcriptome results. This study provides the theoretical basis for further work on the biological roles of AhOMT genes for developmental and stress responses.

Genome-Wide Investigation of Apyrase (APY) Genes in Peanut (Arachis hypogaea L.) and Functional Characterization of a Pod-Abundant Expression Promoter AhAPY2-1p

Peanut (Arachis hypogaea L.) is an important food and feed crop worldwide and is affected by various biotic and abiotic stresses. The cellular ATP levels decrease significantly during stress as ATP molecules move to extracellular spaces, resulting in increased ROS production and cell apoptosis. Apyrases (APYs) are the nucleoside phosphatase (NPTs) superfamily members and play an important role in regulating cellular ATP levels under stress. We identified 17 APY homologs in A. hypogaea (AhAPYs), and their phylogenetic relationships, conserved motifs, putative miRNAs targeting different AhAPYs, cis-regulatory elements, etc., were studied in detail. The transcriptome expression data were used to observe the expression patterns in different tissues and under stress conditions. We found that the AhAPY2-1 gene showed abundant expression in the pericarp. As the pericarp is a key defense organ against environmental stress and promoters are the key elements regulating gene expression, we functionally characterized the AhAPY2-1 promoter for its possible use in future breeding programs. The functional characterization of AhAPY2-1P in transgenic Arabidopsis plants showed that it effectively regulated GUS gene expression in the pericarp. GUS expression was also detected in flowers of transgenic Arabidopsis plants. Overall, these results strongly suggest that APYs are an important future research subject for peanut and other crops, and AhPAY2-1P can be used to drive the resistance-related genes in a pericarp-specific manner to enhance the defensive abilities of the pericarp.

Transcriptome sequencing and gene expression analysis revealed early ovule abortion of Paeonia ludlowii

BACKGROUND

Paeonia ludlowii (Stern & G. Taylor D.Y. Hong) belongs to the peony group of the genus Paeonia in the Paeoniaceae family and is now classified as a "critically endangered species" in China. Reproduction is important for this species, and its low fruiting rate has become a critical factor limiting both the expansion of its wild population and its domestic cultivation.

RESULTS

In this study, we investigated possible causes of the low fruiting rate and ovule abortion in Paeonia ludlowii. We clarified the characteristics of ovule abortion and the specific time of abortion in Paeonia ludlowii, and used transcriptome sequencing to investigate the mechanism of abortion of ovules in Paeonia ludlowii.

CONCLUSIONS

In this paper, the ovule abortion characteristics of Paeonia ludlowii were systematically studied for the first time and provide a theoretical basis for the optimal breeding and future cultivation of Paeonia ludlowii.

Molecular cloning and functional characterization of the promoter of a novel Aspergillus flavus inducible gene (AhOMT1) from peanut

Peanut is an important oil and food legume crop grown in more than one hundred countries, but the yield and quality are often impaired by different pathogens and diseases, especially aflatoxins jeopardizing human health and causing global concerns. For better management of aflatoxin contamination, we report the cloning and characterization of a novel A. flavus inducible promoter of the O-methyltransferase gene (AhOMT1) from peanut. The AhOMT1 gene was identified as the highest inducible gene by A. flavus infection through genome-wide microarray analysis and verified by qRT-PCR analysis. AhOMT1 gene was studied in detail, and its promoter, fussed with the GUS gene, was introduced into Arabidopsis to generate homozygous transgenic lines. Expression of GUS gene was studied in transgenic plants under the infection of A. flavus. The analysis of AhOMT1 gene characterized by in silico assay, RNAseq, and qRT-PCR revealed minute expression in different organs and tissues with trace or no response to low temperature, drought, hormones, Ca2+, and bacterial stresses, but highly induced by A. flavus infection. It contains four exons encoding 297 aa predicted to transfer the methyl group of S-adenosyl-L-methionine (SAM). The promoter contains different cis-elements responsible for its expression characteristics. Functional characterization of AhOMT1P in transgenic Arabidopsis plants demonstrated highly inducible behavior only under A. flavus infection. The transgenic plants did not show GUS expression in any tissue(s) without inoculation of A. flavus spores. However, GUS activity increased significantly after inoculation of A. flavus and maintained a high level of expression after 48 hours of infection. These results provided a novel way for future management of peanut aflatoxins contamination through driving resistance genes in A. flavus inducible manner.

Genome-wide identification of germin-like proteins in peanut (Arachis hypogea L.) and expression analysis under different abiotic stresses

Peanut is an important food and feed crop, providing oil and protein nutrients. Germins and germin-like proteins (GLPs) are ubiquitously present in plants playing numerous roles in defense, growth and development, and different signaling pathways. However, the GLP members have not been comprehensively studied in peanut at the genome-wide scale. We carried out a genome-wide identification of the GLP genes in peanut genome. GLP members were identified comprehensively, and gene structure, genomic positions, motifs/domains distribution patterns, and phylogenetic history were studied in detail. Promoter Cis-elements, gene duplication, collinearity, miRNAs, protein-protein interactions, and expression were determined. A total of 84 GLPs (AhGLPs ) were found in the genome of cultivated peanut. These GLP genes were clustered into six groups. Segmental duplication events played a key role in the evolution of AhGLPs, and purifying selection pressure was underlying the duplication process. Most AhGLPs possessed a well-maintained gene structure and motif organization within the same group. The promoter regions of AhGLPs contained several key cis-elements responsive to 'phytohormones', 'growth and development', defense, and 'light induction'. Seven microRNAs (miRNAs) from six families were found targeting 25 AhGLPs. Gene Ontology (GO) enrichment analysis showed that AhGLPs are highly enriched in nutrient reservoir activity, aleurone grain, external encapsulating structure, multicellular organismal reproductive process, and response to acid chemicals, indicating their important biological roles. AhGLP14, AhGLP38, AhGLP54, and AhGLP76 were expressed in most tissues, while AhGLP26, AhGLP29, and AhGLP62 showed abundant expression in the pericarp. AhGLP7, AhGLP20, and AhGLP21, etc., showed specifically high expression in embryo, while AhGLP12, AhGLP18, AhGLP40, AhGLP78, and AhGLP82 were highly expressed under different hormones, water, and temperature stress. The qRT-PCR results were in accordance with the transcriptome expression data. In short, these findings provided a foundation for future functional investigations on the AhGLPs for peanut breeding programs.

Integrated microRNA and transcriptome profiling reveals the regulatory network of embryo abortion in jujube

Hybridization is an important approach to the production of new varieties with exceptional traits. Although the kernel rate of wild jujube (Ziziphus jujuba Mill. var. spinosa Hu.) is generally high, that of cultivated jujube (Z. jujuba Mill.) is low, greatly hampering the jujube breeding process. However, the mechanism by which this trait changed during jujube domestication remains unclear. Here, we explored the potential regulatory network that governs jujube embryo abortion using correlation analysis of population traits, artificial pollination, sugar content measurements and multi-omics analysis. The results showed that embryo abortion was an important reason for the low kernel rate of cultivated jujube, and kernel rate was negatively correlated with edible rate. Twenty-one days after pollination was a critical period for embryo abortion. At this time, the sugar content of cultivated 'Junzao' kernels decreased significantly compared with that of the pulp, but sugar content remained relatively stable in kernels of wild 'Suanzao'. A total of 1142 differentially expressed genes targeted by 93 microRNAs (miRNAs) were identified by transcriptome, miRNA and degradome sequencing, and may be involved in the regulation of embryo abortion during kernel development. Among them, DELLA protein, TCP14 and bHLH93 transcription factors have been shown to participate in the regulation of embryonic development. Our findings suggest that carbohydrate flow between different tissues of cultivated jujube exhibits a bias toward the pulp at 21 days after pollination, thereby restricting the process of kernel development. This information enhances our understanding of the embryo abortion process and reveals miRNA-target gene pairs that may be useful for molecular-assisted breeding.

Roles of non-coding RNAs in the hormonal and nutritional regulation in nodulation and nitrogen fixation

Symbiotic nitrogen fixation is an important component in the nitrogen cycle and is a potential solution for sustainable agriculture. It is the result of the interactions between the plant host, mostly restricted to legume species, and the rhizobial symbiont. From the first encounter between the host and the symbiont to eventual successful nitrogen fixation, there are delicate processes involved, such as nodule organogenesis, rhizobial infection thread progression, differentiation of the bacteroid, deregulation of the host defense systems, and reallocation of resources. All these processes are tightly regulated at different levels. Recent evidence revealed that non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), participate in these processes by controlling the transcription and translation of effector genes. In general, ncRNAs are functional transcripts without translation potential and are important gene regulators. MiRNAs, negative gene regulators, bind to the target mRNAs and repress protein production by causing the cleavage of mRNA and translational silencing. LncRNAs affect the formation of chromosomal loops, DNA methylation, histone modification, and alternative splicing to modulate gene expression. Both lncRNAs and circRNAs could serve as target mimics of miRNA to inhibit miRNA functions. In this review, we summarized and discussed the current understanding of the roles of ncRNAs in legume nodulation and nitrogen fixation in the root nodule, mainly focusing on their regulation of hormone signal transduction, the autoregulation of nodulation (AON) pathway and nutrient homeostasis in nodules. Unraveling the mediation of legume nodulation by ncRNAs will give us new insights into designing higher-performance leguminous crops for sustainable agriculture.

Genome-Wide Characterization of Ascorbate Peroxidase Gene Family in Peanut (Arachis hypogea L.) Revealed Their Crucial Role in Growth and Multiple Stress Tolerance

Ascorbate peroxidase (APX), an important antioxidant enzyme, plays a significant role in ROS scavenging by catalyzing the decrease of hydrogen peroxide under various environmental stresses. Nevertheless, information about the APX gene family and their evolutionary and functional attributes in peanut (Arachis hypogea L.) was not reported. Therefore, a comprehensive genome-wide study was performed to discover the APX genes in cultivated peanut genome. This study identified 166 AhAPX genes in the peanut genome, classified into 11 main groups. The gene duplication analysis showed that AhAPX genes had experienced segmental duplications and purifying selection pressure. Gene structure and motif investigation indicated that most of the AhAPX genes exhibited a comparatively well-preserved exon-intron pattern and motif configuration contained by the identical group. We discovered five phytohormones-, six abiotic stress-, and five growth and development-related cis-elements in the promoter regions of AhAPX. Fourteen putative ah-miRNAs from 12 families were identified, targeting 33 AhAPX genes. Furthermore, we identified 3,257 transcription factors from 38 families (including AP2, ARF, B3, bHLH, bZIP, ERF, MYB, NAC, WRKY, etc.) in 162 AhAPX genes. Gene ontology and KEGG enrichment analysis confirm the role of AhAPX genes in oxidoreductase activity, catalytic activity, cell junction, cellular response to stimulus and detoxification, biosynthesis of metabolites, and phenylpropanoid metabolism. Based on transcriptome datasets, some genes such as AhAPX4/7/17/77/82/86/130/133 and AhAPX160 showed significantly higher expression in diverse tissues/organs, i.e., flower, leaf, stem, roots, peg, testa, and cotyledon. Likewise, only a few genes, including AhAPX4/17/19/55/59/82/101/102/137 and AhAPX140, were significantly upregulated under abiotic (drought and cold), and phytohormones (ethylene, abscisic acid, paclobutrazol, brassinolide, and salicylic acid) treatments. qRT-PCR-based expression profiling presented the parallel expression trends as generated from transcriptome datasets. Our discoveries gave new visions into the evolution of APX genes and provided a base for further functional examinations of the AhAPX genes in peanut breeding programs.

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

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

Transcriptome and Small RNA Sequencing Reveal the Mechanisms Regulating Harvest Index in Brassica napus

Harvest index (HI), the ratio of harvested seed weight to total aboveground biomass weight, is an economically critical value reflecting the convergence of complex agronomic traits. HI values in rapeseed (Brassica napus) remain much lower than in other major crops, and the underlying regulatory network is largely unknown. In this study, we performed mRNA and small RNA sequencing to reveal the mechanisms shaping HI in B. napus during the seed-filling stage. A total of 8,410 differentially expressed genes (DEGs) between high-HI and low-HI accessions in four tissues (silique pericarp, seed, leaves, and stem) were identified. Combining with co-expression network, 72 gene modules were identified, and a key gene BnaSTY46 was found to participate in retarded establishment of photosynthetic capacity to influence HI. Further research found that the genes involved in circadian rhythms and response to stimulus may play important roles in HI and that their transcript levels were modulated by differentially expressed microRNAs (DEMs), and we identified 903 microRNAs (miRNAs), including 46 known miRNAs and 857 novel miRNAs. Furthermore, transporter activity-related genes were critical to enhancing HI in good cultivation environments. Of 903 miRNAs, we found that the bna-miR396-Bna.A06SRp34a/Bna.A01EMB3119 pair may control the seed development and the accumulation of storage compounds, thus contributing to higher HI. Our findings uncovered the underlying complex regulatory network behind HI and offer potential approaches to rapeseed improvement.

Genome-wide Identification of Metal Tolerance Protein Genes in Peanut: Differential Expression in the Root of Two Contrasting Cultivars Under Metal Stresses

Metal tolerance proteins (MTP) are Me²⁺/H⁺(K⁺) antiporters that play important roles in the transport of divalent cations in plants. However, their functions in peanut are unknown. In the present study, a total of 24 AhMTP genes were identified in peanut, which were divided into seven groups belonging to three substrate-specific clusters (Zn-CDFs, Zn/Fe-CDFs, and Mn-CDFs). All AhMTP genes underwent whole genome or segmental gene duplication events except AhMTP12. Most AhMTP members within the same subfamily or group generally have similar gene and protein structural characteristics. However, some genes, such as AhMTP1.3, AhMTP2.4, and AhMTP12, showed wide divergences. Most of AhMTP genes preferentially expressed in reproductive tissues, suggesting that these genes might play roles in metal transport during the pod and seed development stages. Excess metal exposure induced expressions for most of AhMTP genes in peanut roots depending on cultivars. By contrast, AhMTP genes in the root of Fenghua 1 were more sensitive to excess Fe, Cd, and Zn exposure than that of Silihong. Stepwise linear regression analysis showed that the percentage of Fe in shoots significantly and positively correlated with the expression of AhMTP4.1, AhMTP9.1, and AhMTPC4.1, but negatively correlated with that of AhMTPC2.1 and AhMTP12. The expression of AhMTP1.1 showed a significant and negative correlation with the percentage of Mn in shoots. The percentage of Zn in shoots was significantly and positively correlated with the expression of AhMTP2.1 but was negatively correlated with that of AhMTPC2.1. The differential responses of AhMTP genes to metal exposure might be, at least partially, responsible for the different metal translocation from roots to shoots between Fenghua 1 and Silihong.

Integrated analysis of the lncRNA/circRNA-miRNA-mRNA expression profiles reveals novel insights into potential mechanisms in response to root-knot nematodes in peanut

Background

Peanut is the most essential oil and food crop globally due to its high oil and protein content. Root-knot nematode infects peanut roots, causing poor development and severely limiting peanut yields worldwide. The discovery of peanut genome identified a considerable number of genetic loci controlling the peanut root-knot nematode; however, the molecular mechanism of root-knot nematode remains unknown.

Results

The heterogeneous response to root-knot nematode stress in peanut roots was identified using whole-transcriptome RNA-seq. A total of 430 mRNAs, 111 miRNAs, 4453 lncRNAs, and 123 circRNAs were found to have differential expression between infected and non-infected peanuts. The expression profiles of the lncRNA/circRNA-miRNA-mRNA network were developed to understand the potential pathways that lead to root-knot nematodes in peanut roots. During root-knot nematodes stress, a total of 10 lncRNAs, 4 circRNAs, 5 miRNAs, and 13 mRNAs can create competing endogenous RNA and participate in the oxidation–reduction process as well as other biological metabolism processes in peanuts. The findings will highlight the role of peanut ceRNAs in response to root-knot nematodes.

Conclusion

The GO classification and KEGG pathway enrichment study of core regulatory networks revealed that ceRNAs are involved in oxidation–reduction, peroxidase activity, lignin synthesis in the xylem, and flavonoid synthesis. Overall, these findings may help researchers better understand the role of non-coding RNAs in response to root-knot nematodes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08470-3.

Identification of Key Gene Networks and Deciphering Transcriptional Regulators Associated With Peanut Embryo Abortion Mediated by Calcium Deficiency

Peanut embryo development is easily affected by a variety of nutrient elements in the soil, especially the calcium level. Peanut produces abortive embryos in calcium-deficient soil, but underlying mechanism remains unclear. Thus, identifying key transcriptional regulators and their associated regulatory networks promises to contribute to a better understanding of this process. In this study, cellular biology and gene expression analyses were performed to investigate peanut embryo development with the aim to discern the global architecture of gene regulatory networks underlying peanut embryo abortion under calcium deficiency conditions. The endomembrane systems tended to disintegrate, impairing cell growth and starch, protein and lipid body accumulation, resulting in aborted seeds. RNA-seq analysis showed that the gene expression profile in peanut embryos was significantly changed under calcium deficiency. Further analysis indicated that multiple signal pathways were involved in the peanut embryo abortion. Differential expressed genes (DEGs) related to cytoplasmic free Ca²⁺ were significantly altered. DEGs in plant hormone signaling pathways tended to be associated with increased IAA and ethylene but with decreased ABA, gibberellin, cytokinin, and brassinosteroid levels. Certain vital genes, including apoptosis-inducing factor, WRKYs and ethylene-responsive transcription factors, were up-regulated, while key regulators of embryo development, such as TCP4, WRI1, FUS3, ABI3, and GLK1 were down-regulated. Weighted gene co-expression network analysis (WGCNA) identified 16 significant modules associated with the plant hormone signaling, MAPK signaling, ubiquitin mediated proteolysis, reserve substance biosynthesis and metabolism pathways to decipher regulatory network. The most significant module was darkolivegreen2 and FUS3 (AH06G23930) had the highest connectivity among this module. Importantly, key transcription factors involved in embryogenesis or ovule development including TCP4, GLK1, ABI3, bHLH115, MYC2, etc., were also present in this module and down regulated under calcium deficiency. This study presents the first global view of the gene regulatory network involved in peanut embryo abortion under calcium deficiency conditions and lays foundation for improving peanut tolerances to calcium deficiency by a targeted manipulation of molecular breeding.

MicroRNA Mediated Plant Responses to Nutrient Stress

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

miR-222-3p is involved in neural tube closure by directly targeting Ddit4 in RA induced NTDs mouse model

Previously our results showed miR-222-3p was significantly downregulated in retinoic acid-induced neural tube defect (NTD) mouse model through transcriptome. Down-regulation of miR-222-3p may be a causative biomarker in NTDs. In this study, RNA was extracted from mouse embryos at E8.5, E9.5 and E10.5, and the expression level of miR-222-3p was measured by quantitative real-time PCR analysis. The preliminary mechanism of miR-222-3p in NTDs involved in cell proliferation, apoptosis and migration was investigated in mouse HT-22 cell line. The expression of miR-222-3p was significantly decreased at E8.5, E9.5 and E10.5 developed in mouse embryos which were consistent with our transcriptome sequencing. Suppression of miR-222-3p in HT-22 cells resulted in the inhibition of cell proliferation and migration, cell cycle and apoptosis. Moreover, DNA damage transcript 4 (Ddit4) was identified as a direct and functional target of miR-222-3p. miR-222-3p is negatively regulated by Ddit4. The mutation of binding site of Ddit4 3'UTR abrogated the responsiveness of luciferase reporters to miR-222-3p and showed that Ddit4 expression partially attenuated the function of miR-222-3p. We preliminatively confirmed that low expression of miR-222-3p has reduced the expression of β-catenin, TCF4 and other related genes in the Wnt/β-catenin signaling pathway.Collectively, these results demonstrated that miR-222-3p regulates the Wnt/β-catenin signaling pathway through Ddit4 inhibition in HT-22 cells, resulted in cell proliferation and apoptosis imbalance, and thus led to neural tube defects.

The resistance of peanut to soil-borne pathogens improved by rhizosphere probiotics under calcium treatment

Background

Peanut (Arachis hypogaea L.) is an important oil and economic crop. Calcium modulates plants in response to abiotic stresses and improves plant resistance to pathogens. Enrichment of beneficial microorganisms in the rhizosphere is associated with plant disease resistance and soil development. The purpose of this study was to analyze the differences in peanut rhizosphere microbial community structure between the calcium treatment and the control during two growth stages and to explain why calcium application could improve the resistance of peanuts to soil-borne pathogens.

Results

The 16S rDNA amplicon sequencing of rhizosphere microbiome showed that calcium application significantly enriched Serratia marcescens and other three dominant strains at the seedling stage. At the pod filling stage, ten dominant stains such as Sphingomonas changbaiensis and Novosphingobium panipatense were enriched by calcium. Serratia marcescens aseptic fermentation filtrate was mixed with PDA medium and inoculated with the main soil-borne pathogens in the seedling stage, which could inhibit the growth of Fusarium solani and Aspergillus flavus. The aseptic fermentation filtrate of Novosphingobium panipatense was mixed with PDA medium and inoculated with the main soil-borne pathogens in the pod filling stage, which could inhibit the growth of Sclerotium rolfsii and Leptosphaerulina arachidicola.

Conclusions

Calcium application increases the resistance of peanuts to soil-borne pathogens by enriching them with specific dominant bacteria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02355-3.

Investigations on Regulation of MicroRNAs in Rice Reveal [Ca2+]cyt Signal Transduction Regulated MicroRNAs

MicroRNAs (miRNAs) are critical components of the multidimensional regulatory networks in eukaryotic systems. Given their diverse spectrum of function, it is apparent that the transcription, processing, and activity of the miRNAs themselves, is very dynamically regulated. One of the most important and universally implicated signaling molecules is [Ca2+]cyt. It is known to regulate a plethora of developmental and metabolic processes in both plants and animals; however, its impact on the regulation of miRNA expression is relatively less explored. The current study employed a combination of internal and external calcium channel inhibitors to establishing that [Ca2+]cyt signatures actively regulate miRNA expression in rice. Involvement of [Ca2+]cyt in the regulation of miRNA expression was further confirmed by treatment with calcimycin, the calcium ionophore. Modulation of the cytosolic calcium levels was also found to regulate the drought-responsive expression as well as ABA-mediated response of miRNA genes in rice seedlings. The study further establishes the role of calmodulins and Calmodulin-binding Transcription Activators (CAMTAs) as important components of the signal transduction schema that regulates miRNA expression. Yeast one-hybrid assay established that OsCAMTA4 & 6 are involved in the transcriptional regulation of miR156a and miR167h. Thus, the study was able to establish that [Ca2+]cyt is actively involved in regulating the expression of miRNA genes both under control and stress conditions.

Noncoding-RNA-Mediated Regulation in Response to Macronutrient Stress in Plants

Macronutrient elements including nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), and sulfur (S) are required in relatively large and steady amounts for plant growth and development. Deficient or excessive supply of macronutrients from external environments may trigger a series of plant responses at phenotypic and molecular levels during the entire life cycle. Among the intertwined molecular networks underlying plant responses to macronutrient stress, noncoding RNAs (ncRNAs), mainly microRNAs (miRNAs) and long ncRNAs (lncRNAs), may serve as pivotal regulators for the coordination between nutrient supply and plant demand, while the responsive ncRNA-target module and the interactive mechanism vary among elements and species. Towards a comprehensive identification and functional characterization of nutrient-responsive ncRNAs and their downstream molecules, high-throughput sequencing has produced massive omics data for comparative expression profiling as a first step. In this review, we highlight the recent findings of ncRNA-mediated regulation in response to macronutrient stress, with special emphasis on the large-scale sequencing efforts for screening out candidate nutrient-responsive ncRNAs in plants, and discuss potential improvements in theoretical study to provide better guidance for crop breeding practices.

Genome-Wide Identification and Characterization of HSP90-RAR1-SGT1-Complex Members From Arachis Genomes and Their Responses to Biotic and Abiotic Stresses

The molecular chaperone complex HSP90-RAR1-SGT1 (HRS) plays important roles in both biotic and abiotic stress responses in plants. A previous study showed that wild peanut Arachis diogoi SGT1 (AdSGT1) could enhance disease resistance in transgenic tobacco and peanut. However, no systematic analysis of the HRS complex in Arachis has been conducted to date. In this study, a comprehensive analysis of the HRS complex were performed in Arachis. Nineteen HSP90, two RAR1 and six SGT1 genes were identified from the allotetraploid peanut Arachis hypogaea, a number close to the sum of those from the two wild diploid peanut species Arachis duranensis and Arachis ipaensis. According to phylogenetic and chromosomal location analyses, thirteen orthologous gene pairs from Arachis were identified, all of which except AhHSP90-A8, AhHSP90-B9, AdHSP90-9, and AiHSP90-9 were localized on the syntenic locus, and they shared similar exon-intron structures, conserved motifs and expression patterns. Phylogenetic analysis showed that HSP90 and RAR1 from dicot and monocot plants diverged into different clusters throughout their evolution. Chromosomal location analysis indicated that AdSGT1 (the orthologous gene of AhSGT1-B3 in this study) might provide resistance to leaf late spot disease dependent on the orthologous genes of AhHSP90-B10 and AhRAR1-B in the wild peanut A. diogoi. Several HRS genes exhibited tissue-specific expression patterns, which may reflect the sites where they perform functions. By exploring published RNA-seq data, we found that several HSP90 genes play major roles in both biotic and abiotic stress responses, especially salt and drought responses. Autoactivation assays showed that AhSGT1-B1 could not be used as bait for yeast two-hybrid (Y2H) library screening. AhRAR1 and AhSGT1 could strongly interact with each other and interact with AhHSP90-B8. The present study represents the first systematic analysis of HRS complex genes in Arachis and provides valuable information for functional analyses of HRS complex genes. This study also offers potential stress-resistant genes for peanut improvement.

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

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

miR169 and PmRGL2 synergistically regulate the NF-Y complex to activate dormancy release in Japanese apricot (Prunus mume Sieb. et Zucc.)

This study is the first to demonstrate that GA₄-induced dormancy release is associated with the NF-Y complex, which interacts with gibberellin inhibitor RGL2 in Japanese apricot. Seasonal dormancy is not only vital for the survival in cold winter but also affects flowering of temperate fruit trees and the dormancy release depends on the accumulation of the cold temperatures (Chilling requirement-CR). To understand the mechanism of dormancy release in deciduous fruit crops, we compared miRNA sequencing data during the transition stage from paradormancy to dormancy release in the Japanese apricot and found that the miR169 family showed significant differentially up-regulated expression during dormancy induction and was down-regulated during the dormancy release periods. The 5' RACE assay and RT-qPCR validated its target gene NUCLEAR FACTOR-Y subunit A (NF-YA), which exhibited the opposite expression pattern. Further study showed that exogenous GA₄ could inhibit the expression of the gibberellic acid (GA) signal transduction suppressor PmRGL2 (RGA-LIKE 2) and promote the expression of NF-Y. Moreover, the interaction between the NF-Y family and GA inhibitor PmRGL2 was verified by the yeast-two-hybrid (Y2H) system and a bimolecular fluorescence complementarity (BiFC) experiment. These results suggest that synergistic regulation of the NF-Y and PmRGL2 complex leads to the activation of dormancy release induced by GA₄. These findings will help to elucidate the functional and regulatory roles of miR169 and NF-Y complex in seasonal bud dormancy induced by GA in Japanese apricot and provide new insights for the discovery of dormancy release mechanisms in woody plants.

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

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

Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the Water-Deficit and Heat Stress Response Network in Durum Wheat

Water-deficit and heat stress negatively impact crop production. Mechanisms underlying the response of durum wheat to such stresses are not well understood. With the new durum wheat genome assembly, we conducted the first multi-omics analysis with next-generation sequencing, providing a comprehensive description of the durum wheat small RNAome (sRNAome), mRNA transcriptome, and degradome. Single and combined water-deficit and heat stress were applied to stress-tolerant and -sensitive Australian genotypes to study their response at multiple time-points during reproduction. Analysis of 120 sRNA libraries identified 523 microRNAs (miRNAs), of which 55 were novel. Differentially expressed miRNAs (DEMs) were identified that had significantly altered expression subject to stress type, genotype, and time-point. Transcriptome sequencing identified 49,436 genes, with differentially expressed genes (DEGs) linked to processes associated with hormone homeostasis, photosynthesis, and signaling. With the first durum wheat degradome report, over 100,000 transcript target sites were characterized, and new miRNA-mRNA regulatory pairs were discovered. Integrated omics analysis identified key miRNA-mRNA modules (particularly, novel pairs of miRNAs and transcription factors) with antagonistic regulatory patterns subject to different stresses. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis revealed significant roles in plant growth and stress adaptation. Our research provides novel and fundamental knowledge, at the whole-genome level, for transcriptional and post-transcriptional stress regulation in durum wheat.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Integration of sRNA, degradome, transcriptome analysis and functional investigation reveals gma-miR398c negatively regulates drought tolerance via GmCSDs and GmCCS in transgenic Arabidopsis and soybean

BACKGROUND

Drought conditions adversely affect soybean growth, resulting in severe yield losses worldwide. Increasing experimental evidence indicates miRNAs are important post-transcriptional regulators of gene expression. However, the drought-responsive molecular mechanism underlying miRNA-mRNA interactions remains largely uncharacterized in soybean. Meanwhile, the miRNA-regulated drought response pathways based on multi-omics approaches remain elusive.

RESULTS

We combined sRNA, transcriptome and degradome sequencing to elucidate the complex regulatory mechanism mediating soybean drought resistance. One-thousand transcripts from 384 target genes of 365 miRNAs, which were enriched in the peroxisome, were validated by degradome-seq. An integrated analysis showed 42 miRNA-target pairs exhibited inversely related expression profiles. Among these pairs, a strong induction of gma-miR398c as a major gene negatively regulates multiple peroxisome-related genes (GmCSD1a/b, GmCSD2a/b/c and GmCCS). Meanwhile, we detected that alternative splicing of GmCSD1a/b might affect soybean drought tolerance by bypassing gma-miR398c regulation. Overexpressing gma-miR398c in Arabidopsis thaliana L. resulted in decreased percentage germination, increased leaf water loss, and reduced survival under water deficiency, which displayed sensitivity to drought during seed germination and seedling growth. Furthermore, overexpressing gma-miR398c in soybean decreased GmCSD1a/b, GmCSD2a/b/c and GmCCS expression, which weakened the ability to scavenge O2.-, resulting in increased relative electrolyte leakage and stomatal opening compared with knockout miR398c and wild-type soybean under drought conditions.

CONCLUSION

The study indicates that gma-miR398c negatively regulates soybean drought tolerance, and provides novel insights useful for breeding programs to improve drought resistance by CRISPR technology.