Transcriptome Analysis Reveals Key Seed-Development Genes in Common Buckwheat (Fagopyrum esculentum)

Buckwheat OMICS: present status and future prospects

Buckwheat (Fagopyrum spp.) is an underutilized resilient crop of North Western Himalayas belonging to the family Polygonaceae and is a source of essential nutrients and therapeutics. Common Buckwheat and Tatary Buckwheat are the two main cultivated species used as food. It is the only grain crop possessing rutin, an important metabolite with high nutraceutical potential. Due to its inherent tolerance to various biotic and abiotic stresses and a short life cycle, Buckwheat has been proposed as a model crop plant. Nutritional security is one of the major concerns, breeding for a nutrient-dense crop such as Buckwheat will provide a sustainable solution. Efforts toward improving Buckwheat for nutrition and yield are limited due to the lack of available: genetic resources, genomics, transcriptomics and metabolomics. In order to harness the agricultural importance of Buckwheat, an integrated breeding and OMICS platforms needs to be established that can pave the way for a better understanding of crop biology and developing commercial varieties. This, coupled with the availability of the genome sequences of both Buckwheat species in the public domain, should facilitate the identification of alleles/QTLs and candidate genes. There is a need to further our understanding of the molecular basis of the genetic regulation that controls various economically important traits. The present review focuses on: the food and nutritional importance of Buckwheat, its various omics resources, utilization of omics approaches in understanding Buckwheat biology and, finally, how an integrated platform of breeding and omics will help in developing commercially high yielding nutrient rich cultivars in Buckwheat.

Chromatographic authentication of botanical origin: Herbaceous pollen profiling with HPLC, HPTLC and GC-MS analysis

This study describes a robust chromatographic authentication methodology for herbaceous pollen, employing gas chromatography-mass spectrometry (GC-MS), high-performance liquid chromatography (HPLC) and high-performance thin liquid chromatography (HPTLC) protocols. The comprehensive profiling of organic compounds not only distinguishes between different botanical sources but also establishes a reliable framework for quality control and assessment of herbaceous pollen authenticity. Traces of quercetin were detectable using HPTLC in Chaenomeles japonica, and the composition of the mobile phase led to distinct phenolic acid tracks in the extracts of free phenolic compounds. In Lonicera nummulariifolia, prominent chlorogenic acid signal and traces of 3,4-dihydroxybenzoic acid were identified, along with the presence of vanillic, trans-ferulic, p-coumaric and p-hydroxybenzoic and sinapic as phenolic acid standards. The HPLC chromatogram identified six peaks representing bioactive phenolic compounds such as gallic acid measuring 5.89 ± 0.56 mg g-1, hydroxybenzoic acid 2.39 ± 0.78 mg g-1 and caffeic acid 2.83 ± 0.11 mg g-1. The combined use of GC-MS, HPTLC and HPLC techniques provides a powerful and reliable means of authenticating the botanical origin of herbaceous pollen, offering valuable insights for quality control and ensuring the accuracy of botanical source identification.

Transcriptome Landscape Analyses of the Regulatory Network for Zygotic Embryo Development in Paeonia ostii

Paeonia ostii is a worldwide ornamental flower and an emerging oil crop. Zyotic embryogenesis is a critical process during seed development, and it can provide a basis for improving the efficiency of somatic embryogenesis (SE). In this study, transcriptome sequencing of embryo development was performed to investigate gene expression profiling in P. ostii and identified Differentially expressed genes (DEGs) related to transcription factors, plant hormones, and antioxidant enzymes. The results indicated that IAA (Indole-3-acetic acid), GA (Gibberellin), BR (Brassinosteroid) and ETH (Ethylene) were beneficial to early embryonic morphogenesis, while CTK (Cytokinin) and ABA (Abscisic Acid) promoted embryo morphogenesis and maturation. The antioxidant enzymes' activity was the highest in early embryos and an important participant in embryo formation. The high expression of the genes encoding fatty acid desaturase was beneficial to fast oil accumulation. Representative DEGs were selected and validated using qRT-PCR. Protein-protein interaction network (PPI) was predicted, and six central node proteins, including AUX1, PIN1, ARF6, LAX3, ABCB19, PIF3, and PIF4, were screened. Our results provided new insights into the formation of embryo development and even somatic embryo development in tree peonies.

Genomic Sequence of Canadian Chenopodium berlandieri: A North American Wild Relative of Quinoa

Chenopodium berlandieri (pitseed goosefoot) is a widespread native North American plant, which was cultivated and consumed by indigenous peoples prior to the arrival of European colonists. Chenopodium berlandieri is closely related to, and freely hybridizes with the domesticated South American food crop C. quinoa. As such it is a potential source of wild germplasm for breeding with C. quinoa, for improved quinoa production in North America. The C. berlandieri genome sequence could also be a useful source of information for improving quinoa adaptation. To this end, we first optimized barcode markers in two chloroplast genes, rbcL and matK. Together these markers can distinguish C. berlandieri from the morphologically similar Eurasian invasive C. album (lamb's quarters). Second, we performed whole genome sequencing and preliminary assembly of a C. berlandieri accession collected in Manitoba, Canada. Our assembly, while fragmented, is consistent with the expected allotetraploid structure containing diploid Chenopodium sub-genomes A and B. The genome of our accession is highly homozygous, with only one variant site per 3-4000 bases in non-repetitive sequences. This is consistent with predominant self-fertilization. As previously reported for the genome of a partly domesticated Mexican accession of C. berlandieri, our genome assembly is similar to that of C. quinoa. Somewhat unexpectedly, the genome of our accession had almost as many variant sites when compared to the Mexican C. berlandieri, as compared to C. quinoa. Despite the overall similarity of our genome sequence to that of C. quinoa, there are differences in genes known to be involved in the domestication or genetics of other food crops. In one example, our genome assembly appears to lack one functional copy of the SOS1 (salt overly sensitive 1) gene. SOS1 is involved in soil salinity tolerance, and by extension may be relevant to the adaptation of C. berlandieri to the wet climate of the Canadian region where it was collected. Our genome assembly will be a useful tool for the improved cultivation of quinoa in North America.

Comparative proteomic analyses of Tartary buckwheat (Fagopyrum tataricum) seeds at three stages of development

Tartary buckwheat is among the valuable crops, utilized as both food and Chinese herbal medicine. To uncover the accumulation dynamics of the main nutrients and their regulatory mechanism of Tartary buckwheat seeds, microscopic observations and nutrient analysis were conducted which suggested that starch, proteins as well as flavonoid gradually accumulated among seed development. Comparative proteomic analysis of rice Tartary buckwheat at three different developmental stages was performed. A total of 78 protein spots showed differential expression with 74 of them being successfully identified by MALDI-TOF/TOF MS. Among them, granule bound starch synthase (GBSS1) might be the critical enzyme that determines starch biosynthesis, while 11 S seed storage protein and vicilin seemed to be the main globulin and affect seed storage protein accumulation in Tartary buckwheat seeds. Two enzymes, flavanone 3-hydroxylase (F3H) and anthocyanidin reductase (ANR), involved in the flavonoid biosynthesis pathway were identified. Further analysis on the expression profiles of flavonoid biosynthetic genes revealed that F3H might be the key enzyme that promote flavonoid accumulation. This study provides insights into the mechanism of nutrition accumulation at the protein level in Tartary buckwheat seeds and may facilitate in the breeding and enhancement of Tartary buckwheat germplasm.

The spatiotemporal regulations of epicatechin biosynthesis under normal flowering and the continuous inflorescence removal treatment in Fagopyrum dibotrys

BACKGROUND

Flowering is a critical physiological change that interferes with not only biomass yield but also secondary metabolism, such as the biosynthesis of flavonoids, in rhizome/root plants. The continuous inflorescence removal (CIR) treatment is frequently conducted to weaken this effect. Fagopyrum dibotrys (D.Don) H.Hara (Golden buckwheat) is a kind of rhizome medicinal plant rich in flavonoids and is widely used for the treatment of lung diseases. The CIR treatment is usually conducted in F. dibotrys because of its excessive reproductive growth. To uncover the molecular mechanisms, comprehensive analysis was performed using metabolome and transcriptome data obtained from normally bloomed and the CIR treated plants.

RESULTS

Metabolome results demonstrated that in the rhizomes of F. dibotrys, its bioactive compound called epicatechin has higher amount than most of the detected precursors. Compared with the normally bloomed plants, the level of epicatechin in the rhizomes of the CIR group increased by 25% at the withering stage. Based on 96 samples of the control and the CIR groups at 4 flowering stages for 4 tissues, RNA-Seq results revealed a 3 ~ 5 times upregulations of all the key enzyme genes involved in the biosynthesis of epicatechin in both time (from the bud stage to the withering stage) and spatial dimensions (from the top of branch to rhizome) under the CIR treatment compared to normal flowering. Integrated analysis of LC-MS/MS and transcriptome revealed the key roles of several key enzyme genes besides anthocyanidin reductase (ANR). A total of 93 transcription factors were identified to co-expressed with the genes in epicatechin biosynthetic pathway. The flowering activator SQUAMOSA promoter-binding protein like (SPLs) exhibited opposite spatiotemporal expression patterns to that of the epicatechin pathway genes; SPL3 could significantly co-express with all the key enzyme genes rather than the flowering repressor DELLA. Weighted gene co-expression network analysis (WGCNA) further confirmed the correlations among chalcone synthases (CHSs), chalcone isomerases (CHIs), ANRs, SPLs and other transcription factors.

CONCLUSIONS

SPL3 might dominantly mediate the effect of normal flowering and the CIR treatment on the biosynthesis of epicatechin in rhizomes mainly through the negative regulations of its key enzyme genes including CHS, CHI and ANR.

Transcriptome Analysis Reveals Key Pathways and Candidate Genes Controlling Seed Development and Size in Ricebean (Vigna umbellata)

Ricebean (Vigna umbellata) is a lesser known pulse with well-recognized potential. Recently, it has emerged as a legume with endowed nutritional potential because of high concentration of quality protein and other vital nutrients in its seeds. However, the genes and pathways involved in regulating seed development and size are not understood in this crop. In our study, we analyzed the transcriptome of two genotypes with contrasting grain size (IC426787: large seeded and IC552985: small seeded) at two different time points, namely, 5 and 10 days post-anthesis (DPA). The bold seeded genotype across the time points (B5_B10) revealed 6,928 differentially expressed genes (DEGs), whereas the small seeded genotype across the time point (S5_S10) contributed to 14,544 DEGs. We have also identified several candidate genes for seed development–related traits like seed size and 100-seed weight. On the basis of similarity search and domain analysis, some candidate genes (PHO1, cytokinin dehydrogenase, A-type cytokinin, and ARR response negative regulator) related to 100-seed weight and seed size showed downregulation in the small seeded genotype. The MapMan and KEGG analysis confirmed that auxin and cytokinin pathways varied in both the contrasting genotypes and can therefore be the regulators of the seed size and other seed development–related traits in ricebeans. A total of 51 genes encoding SCF^TIR1/AFB, Aux/IAA, ARFs, E3 ubiquitin transferase enzyme, and 26S proteasome showing distinct expression dynamics in bold and small genotypes were also identified. We have also validated randomly selected SSR markers in eight accessions of the Vigna species (V. umbellata: 6; Vigna radiata: 1; and Vigna mungo: 1). Cross-species transferability pattern of ricebean–derived SSR markers was higher in V. radiata (73.08%) than V. mungo (50%). To the best of our knowledge, this is the first transcriptomic study conducted in this crop to understand the molecular basis of any trait. It would provide us a comprehensive understanding of the complex transcriptome dynamics during the seed development and gene regulatory mechanism of the seed size determination in ricebeans.

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.

Reducing Flower Competition for Assimilates by Half Results in Higher Yield of Fagopyrum esculentum

Despite abundant flowering throughout the season, common buckwheat develops a very low number of kernels probably due to competition for assimilates. We hypothesized that plants with a shorter flowering period may give a higher seed yield. To verify the hypothesis, we studied nutrient stress in vitro and in planta and analyzed different embryological and yield parameters, including hormone profile in the flowers. In vitro cultivated flowers on media with strongly reduced nutrient content demonstrated a drastic increase in degenerated embryo sacs. In in planta experiments, where 50% or 75% of flowers or all lateral ramifications were removed, the reduction of the flower competition by half turned out to be the most promising treatment for improving yield. This treatment increased the frequency of properly developed embryo sacs, the average number of mature seeds per plant, and their mass. Strong seed compensation under 50% inflorescence removal could result from increased production of salicylic and jasmonic acid that both favor more effective pollinator attraction. Plants in single-shoot cultivation finished their vegetation earlier, and they demonstrated greater single seed mass per plant than in control. This result suggests that plants of common buckwheat with shorter blooming period could deliver higher seed yield.

Comparative cellular, physiological and transcriptome analyses reveal the potential easy dehulling mechanism of rice-tartary buckwheat (Fagopyrum Tararicum)

BACKGROUND

Tartary buckwheat has gained popularity in the food marketplace due to its abundant nutrients and high bioactive flavonoid content. However, its difficult dehulling process has severely restricted its food processing industry development. Rice-tartary buckwheat, a rare local variety, is very easily dehulled, but the cellular, physiological and molecular mechanisms responsible for this easy dehulling remains largely unclear.

RESULTS

In this study, we integrated analyses of the comparative cellular, physiological, transcriptome, and gene coexpression network to insight into the reason that rice-tartary buckwheat is easy to dehull. Compared to normal tartary buckwheat, rice-tartary buckwheat has significantly brittler and thinner hull, and thinner cell wall in hull sclerenchyma cells. Furthermore, the cellulose, hemicellulose, and lignin contents of rice-tartary buckwheat hull were significantly lower than those in all or part of the tested normal tartary buckwheat cultivars, respectively, and the significant difference in cellulose and hemicellulose contents between rice-tartary buckwheat and normal tartary buckwheat began at 10 days after pollination (DAP). Comparative transcriptome analysis identified a total of 9250 differentially expressed genes (DEGs) between the rice- and normal-tartary buckwheat hulls at four different development stages. Weighted gene coexpression network analysis (WGCNA) of all DEGs identified a key module associated with the formation of the hull difference between rice- and normal-tartary buckwheat. In this specific module, many secondary cell wall (SCW) biosynthesis regulatory and structural genes, which involved in cellulose and hemicellulose biosynthesis, were identified as hub genes and displayed coexpression. These identified hub genes of SCW biosynthesis were significantly lower expression in rice-tartary buckwheat hull than in normal tartary buckwheat at the early hull development stages. Among them, the expression of 17 SCW biosynthesis relative-hub genes were further verified by quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

Our results showed that the lower expression of SCW biosynthesis regulatory and structural genes in rice-tartary buckwheat hull in the early development stages contributes to its easy dehulling by reducing the content of cell wall chemical components, which further effects the cell wall thickness of hull sclerenchyma cells, and hull thickness and mechanical strength.

Transcriptome Analysis Provides Insights into Grain Filling in Foxtail Millet (Setaria italica L.)

Grain filling is an importantly developmental process which is associated with the yield and quality of foxtail millet (Setaria italic L.). However, the molecular mechanisms of grain filling are rarely reported in foxtail millet. In our study, RNA-seq was performed to investigate the transcriptional dynamics and identify the key genes involved in grain filling in foxtail millet at five different developmental stages. A total of 11,399 differentially expressed genes (DEGs), including 902 transcription factors (TFs), were identified. Certain important genes involved in grain filling were discovered through a function annotation and temporal expression patterns analysis. These genes included genes associated with starch biosynthesis, cell-wall invertases, hormone signal transduction, and polyamine metabolism pathways. The expression levels of seven randomly selected DEGs were validated by a quantitative real-time polymerase chain reaction (qRT-PCR). This study provides the first insight into the changes in the gene expression of grain filling at different developmental stages in foxtail millet. These results could help understand the complex molecular mechanisms of the panicle formation in foxtail millet and other cereal crops.