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

Identification of Black Cumin (Nigella sativa) MicroRNAs by Next-Generation Sequencing and Their Implications in Secondary Metabolite Biosynthesis

Secondary metabolites are bioactive compounds believed to contribute to the pharmacological properties of plants. MicroRNAs (miRNAs) are small non-coding RNA molecules involved in post-transcriptional regulation and are thought to play an important role in regulating secondary metabolism biosynthesis. Nevertheless, the extent of miRNA involvement in secondary metabolism remains minimal. Nigella sativa (black cumin/black seed) is a popular medicinal and culinary plant known for its pharmaceutical properties; however, its genomic information is scarce. In this study, next-generation sequencing (NGS) technology was employed to obtain the miRNA profile of N. sativa, and their involvement in secondary metabolite biosynthesis was explored. A total of 25,139,003 unique reads ranging from 16 to 40 nucleotides were attained, out of which 240 conserved and 34 novel miRNAs were identified. Moreover, 6083 potential target genes were recognized in this study. Several conserved and novel black cumin miRNAs were found to target enzymes involved in the terpenoid, diterpenoid, phenylpropanoid, carotenoid, flavonoid, steroid, and ubiquinone biosynthetic pathways, among others, for example, beta-carotene 3-hydroxylase, gibberellin 3 beta-dioxygenase, trimethyltridecatetraene synthase, carboxylic ester hydrolases, acetyl-CoA C-acetyltransferase, isoprene synthase, peroxidase, shikimate O-hydroxycinnamoyltransferase, etc. Furthermore, sequencing data were validated through qPCR by checking the relative expression of eleven randomly selected conserved and novel miRNAs (nsa-miR164d, nsa-miR166a, nsa-miR167b, nsa-miR171a, nsa-miR390b, nsa-miR396, nsa-miR159a, nsa-miRN1, nsa-miRN29, nsa-miRN32, and nsa-miRN34) and their expression patterns were found to be corroborated with the sequencing data. We anticipate that this work will assist in clarifying the implications of miRNAs in plant secondary metabolism and aid in the generation of artificial miRNA-based strategies to overproduce highly valuable secondary metabolites from N. sativa.

Advancements in multi-omics for nutraceutical enhancement and traits improvement in buckwheat

Buckwheat (Fagopyrum spp.) is a typical pseudocereal, valued for its extensive nutraceutical potential as well as its centuries-old cultivation. Tartary buckwheat and common buckwheat have been used globally and become well-known nutritious foods due to their high quantities of: proteins, flavonoids, and minerals. Moreover, its increasing demand makes it critical to improve nutraceutical, traits and yield. In this review, bioactive compounds accumulated in buckwheat were comprehensively evaluated according to their chemical structure, properties, and physiological function. Biosynthetic pathways of flavonoids, phenolic acids, and fagopyrin were methodically summarized, with the regulation of flavonoid biosynthesis. Although there are classic synthesis pathways presented in the previous research, the metabolic flow of how these certain compounds are being synthesized in buckwheat still remains uncovered. The functional genes involved in the biosynthesis of flavonols, stress response, and plant development were identified based on multi-omics research. Furthermore, it delves into the applications of multi-omics in improving buckwheat's agronomic traits, including: yield, nutritional content, stress resilience, and bioactive compounds biosynthesis. While pangenomics combined with other omics to mine elite genes, the regulatory network and mechanism of specific agronomic traits and biosynthetic of bioactive components, and developing a more efficient genetic transformation system for genetic engineering require further investigation for the execution of breeding designs aimed at enhancing desirable traits in buckwheat. This critical review will provide a comprehensive understanding of multi-omics for nutraceutical enhancement and traits improvement in buckwheat.

Combined miRNA and mRNA sequencing reveals the defensive strategies of resistant YHY15 rice against differentially virulent brown planthoppers

Introduction

The brown planthopper (BPH) poses a significant threat to rice production in Asia. The use of resistant rice varieties has been effective in managing this pest. However, the adaptability of BPH to resistant rice varieties has led to the emergence of virulent populations, such as biotype Y BPH. YHY15 rice, which carries the BPH resistance gene Bph15, exhibits notable resistance to biotype 1 BPH but is susceptible to biotype Y BPH. Limited information exists regarding how resistant rice plants defend against BPH populations with varying levels of virulence.

Methods

In this study, we integrated miRNA and mRNA expression profiling analyses to study the differential responses of YHY15 rice to both avirulent (biotype 1) and virulent (biotype Y) BPH.

Results

YHY15 rice demonstrated a rapid response to biotype Y BPH infestation, with significant transcriptional changes occurring within 6 hours. The biotype Y-responsive genes were notably enriched in photosynthetic processes. Accordingly, biotype Y BPH infestation induced more intense transcriptional responses, affecting miRNA expression, defenserelated metabolic pathways, phytohormone signaling, and multiple transcription factors. Additionally, callose deposition was enhanced in biotype Y BPH-infested rice seedlings.

Discussion

These findings provide comprehensive insights into the defense mechanisms of resistant rice plants against virulent BPH, and may potentially guide the development of insect-resistant rice varieties.

Integrated analysis of miRNA, transcriptome, and degradome sequencing provides new insights into lipid metabolism in perilla seed

MicroRNAs (miRNA) are small noncoding RNAs that play a crucial as molecular regulators in lipid metabolism in various oil crops. Perilla (Perilla frutescens) is a specific oil crop known for its high alpha-linolenic acid (C18:3n3, ALA) content (>65 %) in their seed oils. In view of the regulatory mechanism of miRNAs in perilla remains unclear, we conducted miRNAs and transcriptome sequencing in two cultivars with distinct lipid compositions. A total of 525 unique miRNAs, including 142 differentially expressed miRNAs was identified in perilla seeds. The 318 miRNAs targeted 7,761 genes. Furthermore, we identified 112 regulated miRNAs and their 610 target genes involved in lipid metabolism. MiR159b and miR167a as the core nodes to regulate the expression of genes in oil biosynthesis (e.g., KAS, FATB, GPAT, FAD, DGK, LPAAT) and key regulatory TFs (e.g., MYB, ARF, DOF, SPL, NAC, TCP, and bHLH). The 1,219 miRNA-mRNA regulation modules were confirmed through degradome sequencing. Notably, pf-miR159b-MYBs and pf-miR167a-ARFs regulation modules were confirmed. They exhibited significantly different expression levels in two cultivars and believed to play important roles in oil biosynthesis in perilla seeds. This provides valuable insights into the functional analysis of miRNA-regulated lipid metabolism in perilla seeds.

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.

Transcriptome and metabolome analyses reveal the key genes related to grain size of big grain mutant in Tartary Buckwheat (Fagopyrum tartaricum)

Grain size with high heritability and stability is an important selection target during Tartary buckwheat breeding. However, the mechanisms that regulate Tartary buckwheat grain development are unknown. We generated transcriptome and metabolome sequencing from 10 and 15 days past anthesis (DPA) grains of big grain mutant (bg1) and WT, and identified 4108 differentially expressed genes (DEGs) including 93 significantly up-regulated differential genes and 85 significantly down-regulated genes in both stages, simultaneously. Meanwhile, we identified DEGs involved in ubiquitin-proteasome pathway, HAI-KU (IKU) pathway, mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone (auxin, brassinosteroids and cytokinins) transduction pathway and five transcription factor families, including APETALA (AP2), GROWTH-REGULATING FACTORS (GRF), AUXIN RESPONSE FACTOR (ARF), WRKY and MYB. Weighted gene co-expression network analysis (WGCNA) was performed and obtained 9 core DEGs. Conjoint analyses of transcriptome and metabolome sequencing screened out 394 DEGs. Using a combined comprehensive analysis, we identified 24 potential candidate genes that encode E3 ubiquitin-protein ligase HIP1, EMBRYO-DEFECTIVE (EMB) protein, receptor-like protein kinase FERONIA (FER), kinesin-4 protein SRG1, and so on, which may be associated with the big-grain mutant bg1. Finally, a quantitative real-time Polymerase Chain Reaction (qRT-PCR) assay was conducted to validate the identified DEGs. Our results provide additional knowledge for identification and functions of causal candidate genes responsible for the variation in grain size and will be an invaluable resource for the genetic dissection of Tartary buckwheat high-yield molecular breeding.

Comparative physiological, transcriptomic, and WGCNA analyses reveal the key genes and regulatory pathways associated with drought tolerance in Tartary buckwheat

Drought stress is one of the major abiotic stress factors that affect plant growth and crop productivity. Tartary buckwheat is a nutritionally balanced and flavonoid-rich pseudocereal crop and also has strong adaptability to different adverse environments including drought. However, little is known about its drought tolerance mechanism. In this study, we performed comparative physiological and transcriptomic analyses of two contrasting drought-resistant Tartary buckwheat genotypes under nature drought treatment in the reproductive stage. Under drought stress, the drought-tolerant genotype XZSN had significantly higher contents of relative water, proline, and soluble sugar, as well as lower relative electrolyte leakage in the leaves than the drought-susceptible LK3. A total of 5,058 (2,165 upregulated and 2,893 downregulated) and 5,182 (2,358 upregulated and 2,824 downregulated) potential drought-responsive genes were identified in XZSN and LK3 by transcriptome sequencing analysis, respectively. Among the potential drought-responsive genes of XZSN, 1,206 and 1,274 genes were identified to be potential positive and negative contributors for XZSN having higher drought resistance ability than LK3. Furthermore, 851 out of 1,206 positive drought-resistant genes were further identified to be the core drought-resistant genes of XZSN based on WGCNA analysis, and most of them were induced earlier and quicker by drought stress than those in LK3. Functional annotation of the 851 core drought-resistant genes found that a large number of stress-responsive genes were involved in TFs, abscisic acid (ABA) biosynthesis, signal transduction and response, non-ABA signal molecule biosynthesis, water holding, oxygen species scavenging, osmotic adjustment, cell damage prevention, and so on. Transcriptional regulatory network analyses identified the potential regulators of these drought-resistant functional genes and found that the HD-ZIP and MYB TFs might be the key downstream TFs of drought resistance in Tartary buckwheat. Taken together, these results indicated that the XZSN genotype was more drought-tolerant than the LK3 genotype as evidenced by triggering the rapid and dramatic transcriptional reprogramming of drought-resistant genes to reduce water loss, prevent cell damage, and so on. This research expands our current understanding of the drought tolerance mechanisms of Tartary buckwheat and provides important information for its further drought resistance research and variety breeding.

Genome-wide identification of miRNAs and targets associated with cell wall biosynthesis: Differential roles of dlo-miR397a and dlo-miR408-3p during early somatic embryogenesis in longan

The dynamic alterations in cell wall (CW) biosynthesis play an essential role in physiological isolation during the plant somatic embryogenesis (SE). However, the mechanisms underlying the functions of cell wall-associated miRNAs (CW-miRNA) remain poorly understood in plant SE. Here, we have identified 36 distinct candidate miRNAs associated with CW biosynthesis from longan third-generation genome as well as miRNA transcriptome, and modified RLM-RACE validated four distinct miRNA, which specifically targeted four CW-related genes. More importantly, we found that the dlo-miR397a-antagomir significantly enhanced DlLAC7 expression and improved laccase activity. Interestingly, inhibition of dlo-miR397a increased CW lignin deposition and promoted the tightening of protodermal cell by miRNA-mimic technology during early SE. Moreover, overexpression of dlo-miR408-3p (dlo-miR408-3p-agomir) markedly decreased DlLAC12 expression. dlo-miR408-3p-agomir activated rapid cell division, thus promoting the globular embryo (GE) development, which might be due to high DNA synthesis activity in protoepidermal cells, rather than affecting lignin synthesis. The subcellular location also indicated that both DlLAC7 and DlLAC12 proteins were primarily localized in CW and regulated CW biosynthesis. Overall, our findings provided new insight on the molecular regulatory networks comprising various miRNAs associated with cell wall, and established that dlo-miR397a and dlo-miR408-3p played differential roles during early SE in longan. The findings also shed some light on the potential role of miRNA target DlLAC regulating in vivo embryonic development of plant.

Cytological, physiological and transcriptomic analysis of variegated Leaves in Primulina pungentisepala offspring

BACKGROUND

Primulina pungentisepala is suitable for use as a potted plant because of its beautiful leaf variegation, which is significantly different in its selfed offspring. However, the mechanism of P. pungentisepala leaf variegation is unclear. In this study, two types of offspring showing the greatest differences were compared in terms of leaf structure, chlorophyll contents, chlorophyll fluorescence parameters and transcriptomes to provide a reference for studying the molecular mechanism of structural leaf variegation.

RESULTS

Air spaces were found between water storage tissue, and the palisade tissue cells were spherical in the white type. The content of chlorophyll a and total chlorophyll (chlorophyll a + b) was significantly lower in the white type, but there were no significant differences in the content of chlorophyll b, chlorophyll a/b or chlorophyll fluorescence parameters between the white and green types. We performed transcriptomic sequencing to identify differentially expressed genes (DEGs) involved in cell division and differentiation, chlorophyll metabolism and photosynthesis. Among these genes, the expression of the cell division- and differentiation-related leucine-rich repeat receptor-like kinases (LRR-RLKs), xyloglucan endotransglycosylase/hydrolase (XET/H), pectinesterase (PE), expansin (EXP), cellulose synthase-like (CSL), VARIEGATED 3 (VAR3), and ZAT10 genes were downregulated in the white type, which might have promoted the development air spaces and variant palisade cells. Chlorophyll biosynthesis-related hydroxymethylbilane synthase (HEMC) and the H subunit of magnesium chelatase (CHLH) were downregulated, while chlorophyll degradation-related chlorophyllase-2 (CHL2) was upregulated in the white type, which might have led to lower chlorophyll accumulation.

CONCLUSION

Leaf variegation in P. pungentisepala was caused by a combination of mechanisms involving structural variegation and low chlorophyll levels. Our research provides significant insights into the molecular mechanisms of structural leaf variegation.

Identification of Nitrogen Starvation-Responsive miRNAs to Reveal the miRNA-Mediated Regulatory Network in Betula luminifera

Because of the immobility, plants encounter a series of stresses, such as varied nutrient concentrations in soil, which regulate plant growth, development, and phase transitions. Nitrogen (N) is one of the most limiting factors for plants, which was exemplified by the fact that low nitrogen (LN) has a great adverse effect on plant growth and development. In the present study, we explored the potential role of microRNAs (miRNAs) in response to LN stress in Betula luminifera. We identified 198 miRNAs using sRNA sequencing, including 155 known and 43 novel miRNAs. Among them, 98 known miRNAs and 31 novel miRNAs were differentially expressed after 0.5 h or 24 h of LN stress. Based on degradome data, 122 differential expressed miRNAs (DEmiRNAs) including 102 known miRNAs and 20 novel miRNAs targeted 203 genes, comprising 321 miRNA–target pairs. A big proportion of target genes were transcription factors and functional proteins, and most of the Gene Ontology terms were enriched in biological processes; moreover, one Kyoto Encyclopedia of Genes and Genomes term “ascorbate and aldarate metabolism” was significantly enriched. The expression patterns of six miRNAs and their corresponding target genes under LN stress were monitored. According to the potential function for targets of DEmiRNAs, a proposed regulatory network mediated by miRNA–target pairs under LN stress in B. luminifera was constructed. Taken together, these findings provide useful information to elucidate miRNA functions and establish a framework for exploring N signaling networks mediated by miRNAs in B. luminifera. It may provide new insights into the genetic engineering of the high use efficiency of N in forestry trees.

In Vitro Effects of Tartary Buckwheat-Derived Nanovesicles on Gut Microbiota

Evidence suggests that plant-derived nanovesicles may play a significant role in human health. Tartary buckwheat has several physiological activities; however, its underlying health-promoting mechanism remains unclear. In this study, first, Tartary buckwheat-derived nanovesicles (TBDNs) were collected, their structures were analyzed, and microRNA sequencing was performed. Next, target prediction and functional verification were conducted. Finally, the effects of TBDNs on gut microbiota and short-chain fatty acid levels were evaluated. The average size of TBDNs was 141.8 nm diameter. Through the sequencing analyses, 129 microRNAs, including 11 novel microRNAs were identified. Target gene prediction showed that some microRNAs could target functional genes in Escherichia coli and Lactobacillus rhamnosus-related physiological processes. TBDNs significantly promoted the growth of E. coli and L. rhamnosus, enhanced the diversity of fecal microorganisms and increased the short-chain fatty acid levels. These findings provided a new nutritional perspective for Tartary buckwheat and were conducive to promote the development and utilization of Tartary buckwheat.