Global transcriptome analysis and identification of genes involved in nutrients accumulation during seed development of rice tartary buckwheat (Fagopyrum Tararicum)

Understanding the amylose biosynthesis and regulation mechanisms in Tartary buckwheat by the endosperm transcriptome

Starch serves as a crucial energy source for both plants and humans, predominantly synthesized and stored in endosperms, tubers, rhizomes, and cotyledons. Given the significant role of amylose in determining the quality of starchy crops, optimizing its content has become a key objective in current crop breeding efforts. Tartary buckwheat, a dicotyledonous plant, notably accumulates high levels of amylose in its endosperm, surpassing common cereals like rice and maize. However, the mechanisms underlying amylose accumulation, distribution, and regulation in Tartary buckwheat remain unclear. Here, amylose content was determined across various tissues and organs of Tartary buckwheat, identifying with the endosperm as the primary site for its biosynthesis and accumulation. RNA sequencing analysis of endosperms from different developmental stages identified 35 genes potentially involved in starch biosynthesis, with 13 genes showing high endosperm-specific expression, suggesting crucial roles in starch biosynthesis. Additionally, the transcription factor FtNF-YB2, which was specifically highly expressed in the endosperm, was discovered to enhance amylose synthesis. Moreover, promoters with potential endosperm-specific activity were identified, advancing our understanding of amylose regulation. Additionally, this study also demonstrates that brassinosteroids (BR) positively influence amylose biosynthesis in Tartary buckwheat endosperm. These findings provide essential insights into the mechanisms of understanding amylose biosynthesis, accumulation and regulation in Tartary buckwheat, offering significant implications for future breeding strategies.

QTL Mapping and Candidate Gene Analysis for Starch-Related Traits in Tartary Buckwheat (Fagopyrum tataricum (L.) Gaertn)

Starch is the main component that determines the yield and quality of Tartary buckwheat. As a quantitative trait, using quantitative trait locus (QTL) mapping to excavate genes associated with starch-related traits is crucial for understanding the genetic mechanisms involved in starch synthesis and molecular breeding of Tartary buckwheat varieties with high-quality starch. Employing a recombinant inbred line population as research material, this study used QTL mapping to investigate the amylose, amylopectin, and total starch contents across four distinct environments. The results identified a total of 20 QTLs spanning six chromosomes, which explained 4.07% to 14.41% of the phenotypic variation. One major QTL cluster containing three stable QTLs governing both amylose and amylopectin content, qClu-4-1, was identified and located in the physical interval of 39.85-43.34 Mbp on chromosome Ft4. Within this cluster, we predicted 239 candidate genes and analyzed their SNP/InDel mutations, expression patterns, and enriched KEGG pathways. Ultimately, five key candidate genes, namely FtPinG0004897100.01, FtPinG0002636200.01, FtPinG0009329200.01, FtPinG0007371600.01, and FtPinG0005109900.01, were highlighted, which are potentially involved in starch synthesis and regulation, paving the way for further investigative studies. This study, for the first time, utilized QTL mapping to detect major QTLs controlling amylose, amylopectin, and total starch contents in Tartary buckwheat. The QTLs and candidate genes would provide valuable insights into the genetic mechanisms underlying starch synthesis and improving starch-related traits of Tartary buckwheat.

Plant Hormone Pathway Is Involved in Regulating the Embryo Development Mechanism of the Hydrangea macrophylla Hybrid

The research is aimed to elucidate the role of plant hormones in regulating the development of hybrid embryos in Hydrangea macrophylla. Fruits from the intraspecific cross of H. macrophylla 'Otaksa' × 'Coerulea' were selected at the globular, heart, and torpedo stages of embryo development. Transcriptome sequencing and differential gene expression analysis were conducted. The results showed that fruit growth followed a single "S-shaped growth curve, with globular, heart, and torpedo embryos appearing at 30, 40, and 50 d post-pollination, respectively, and the embryo maintaining the torpedo shape from 60 to 90 d. A total of 12,933 genes was quantified across the three developmental stages, with 3359, 3803, and 3106 DEGs in the S1_vs_S2, S1_vs_S3, and S2_vs_S3 comparisons, respectively. Among these, 133 genes related to plant hormone biosynthesis and metabolism were differentially expressed, regulating the synthesis and metabolism of eight types of plant hormones, including cytokinin, auxin, gibberellin, abscisic acid, and jasmonic acid. The pathways with the most differentially expressed genes were cytokinin, auxin, and gibberellin, suggesting these hormones may play crucial roles in embryo development. In the cytokinin pathway, CKX (Hma1.2p1_0579F.1_g182670.gene, Hma1.2p1_1194F.1_g265700.gene, and NewGene_12164) genes were highly expressed during the globular embryo stage, promoting rapid cell division in the embryo. In the auxin pathway, YUC (Hma1.2p1_0271F.1_g109005.gene and Hma1.2p1_0271F.1_g109020.gene) genes were progressively up-regulated during embryo growth; the early response factor AUX/IAA (Hma1.2p1_0760F.1_g214260.gene) was down-regulated, while the later transcriptional activator ARF (NewGene_21460, NewGene_21461, and Hma1.2p1_0209F.1_g089090.gene) was up-regulated, sustaining auxin synthesis and possibly preventing the embryo from transitioning to maturity. In the gibberellin pathway, GA3ox (Hma1.2p1_0129F.1_g060100.gene) expression peaked during the heart embryo stage and then declined, while the negative regulator GA2ox (Hma1.2p1_0020F.1_g013915.gene) showed the opposite trend; and the gibberellin signaling repressor DELLA (Hma1.2p1_1054F.1_g252590.gene) increased over time, potentially inhibiting embryo development and maintaining the torpedo shape until fruit maturity. These findings preliminarily uncover the factors affecting the development of hybrid H. macrophylla embryos, laying a foundation for further research into the regulatory mechanisms of H. macrophylla hybrid embryo development.

New strategies to address world food security and elimination of malnutrition: future role of coarse cereals in human health

As we face increasing challenges of world food security and malnutrition, coarse cereals are coming into favor as an important supplement to human staple foods due to their high nutritional value. In addition, their functional components, such as flavonoids and polyphenols, make them an important food source for healthy diets. However, we lack a systematic understanding of the importance of coarse cereals for world food security and nutritional goals. This review summarizes the worldwide cultivation and distribution of coarse cereals, indicating that the global area for coarse cereal cultivation is steadily increasing. This paper also focuses on the special adaptive mechanisms of coarse cereals to drought and discusses the strategies to improve coarse cereal crop yields from the perspective of agricultural production systems. The future possibilities, challenges, and opportunities for coarse cereal production are summarized in the face of food security challenges, and new ideas for world coarse cereal production are suggested.

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.

Understanding the Potential Gene Regulatory Network of Starch Biosynthesis in Tartary Buckwheat by RNA-Seq

Starch is a major component of crop grains, and its content affects food quality and taste. Tartary buckwheat is a traditional pseudo-cereal used in food as well as medicine. Starch content, granule morphology, and physicochemical properties have been extensively studied in Tartary buckwheat. However, the complex regulatory network related to its starch biosynthesis needs to be elucidated. Here, we performed RNA-seq analyses using seven Tartary buckwheat varieties differing in starch content and combined the RNA-seq data with starch content by weighted correlation network analysis (WGCNA). As a result, 10,873 differentially expressed genes (DEGs) were identified and were functionally clustered to six hierarchical clusters. Fifteen starch biosynthesis genes had higher expression level in seeds. Four trait-specific modules and 3131 hub genes were identified by WGCNA, with the lightcyan and brown modules positively correlated with starch-related traits. Furthermore, two potential gene regulatory networks were proposed, including the co-expression of FtNAC70, FtPUL, and FtGBSS1-3 in the lightcyan module and FtbHLH5, C3H, FtBE2, FtISA3, FtSS3-5, and FtSS1 in the brown. All the above genes were preferentially expressed in seeds, further suggesting their role in seed starch biosynthesis. These results provide crucial guidance for further research on starch biosynthesis and its regulatory network in Tartary buckwheat.

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.

Integrative metabolomic and transcriptomic analyses reveal the mechanisms of Tibetan hulless barley grain coloration

Cereal grains accumulate anthocyanin during developmental process. The anthocyanin content increases at grain filling stages to develop grain coloration in cereals. However, anthocyanin biosynthesis responsible for grain coloring and its regulatory mechanisms controlled by structural and functional genes remain unclear. Therefore, this study aimed to explore the global map of metabolic changes linked to grain coloration of Tibetan hulless barley (qingke) using an integrative metabolome and transcriptome approach. Grains from three colored qingke cultivars at different developmental stages were considered for molecular and metabolic investigations. A total of 120 differentially accumulated metabolites (DAMs) and 8,327 differentially expressed genes (DEGs) were filtered. DEGs were mainly enriched in the phenylpropanoid and flavonoid pathways. The transcript levels of anthocyanin biosynthesis genes (PAL, C4H, 4CL, CHS, FLS, F3H, F3'H, DFR, ANS, GT, OMT, and MAT) significantly upregulate in colored qingke compared to the non-colored variety. During grain development and maturation, the strong correlation of HvMYC2 expression with anthocyanin contents and anthocyanin biosynthesis genes suggested it as a critical gene in anthocyanin accumulation. Further results confirmed that HvMYC2 could be activated by HvMYB and be a positive regulator of UV-B and cold tolerance in qingke. In addition, verification based on enzymatic assays indicated that six key modifier enzymes could catalyze glycosylation, malonylation, and methylation of anthocyanins, thereby dissecting the major anthocyanin modification pathway in colored qingke. Overall, our study provides global insight into anthocyanin accumulation and the mechanism underlying grain coloration in qingke.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Transcriptome and Coexpression Network Analyses Provide In-Sights into the Molecular Mechanisms of Hydrogen Cyanide Synthesis during Seed Development in Common Vetch (Vicia sativa L.)

The common vetch (Vicia sativa L.) seed is an ideal plant-based protein food for humans, but its edible value is mainly limited by the presence of cyanogenic glycosides that hydrolyze to produce toxic hydrogen cyanide (HCN), and the genes that regulate HCN synthesis in common vetch are unknown. In this study, seeds from common vetch at 5, 10, 15, 20, 25, 30, and 35 days after anthesis were sampled, and the seven stages were further divided into five developmental stages, S1, S2, S3, S4, and S5, based on morphological and transcriptome analyses. A total of 16,403 differentially expressed genes were identified in the five developmental stages. The HCN contents of seeds in these five stages were determined by alkaline titration, and weighted gene coexpression network analysis was used to explain the molecular regulatory mechanism of HCN synthesis in common vetch seeds. Eighteen key regulatory genes for HCN synthesis were identified, including the VsGT2, VsGT17 and CYP71A genes, as well as the VsGT1 gene family. VsGT1, VsGT2, VsGT17 and CYP71A jointly promoted HCN synthesis, from 5 to 25 days after anthesis, with VsGT1-1, VsGT1-4, VsGT1-11 and VsGT1-14 playing major roles. The HCN synthesis was mainly regulated by VsGT1, from 25 to 35 days after anthesis. As the expression level of VsGT1 decreased, the HCN content no longer increased. In-depth elucidation of seed HCN synthesis lays the foundations for breeding common vetch with low HCN content.

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.

Transcriptome analysis revealed gene regulatory network involved in PEG-induced drought stress in Tartary buckwheat (Fagopyrum Tararicum)

Tartary buckwheat is a nutritious pseudo-cereal crop that is resistant to abiotic stresses, such as drought. However, the buckwheat’s mechanisms for responding to drought stress remains unknown. We investigated the changes in physiology and gene expression under drought stress, which was simulated by treatment with polyethylene glycol (PEG). Five physiological indexes, namely MDA content, H₂O₂ content, CAT activity, SOD activity, and POD activity, were measured over time after 20% PEG treatment. All indexes showed dramatic changes in response to drought stress. A total of 1,190 differentially expressed genes (DEGs) were identified using RNA-seq and the most predominant were related to a number of stress-response genes and late embryogenesis abundant (LEA) proteins. DEGs were gathered into six clusters and were found to be involved in the ABA biosynthesis and signal pathway based on hierarchical clustering and GO and KEGG pathway enrichment. Transcription factors, such as NAC and bZIP, also took part in the response to drought stress. We determined an ABA-dependent and ABA-independent pathway in the regulation of drought stress in Tartary buckwheat. To the best of our knowledge, this is the first transcriptome analysis of drought stress in Tartary buckwheat, and our results provide a comprehensive gene regulatory network of this crop in response to drought stress.

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.

Mapping QTLs for 1000-grain weight and genes controlling hull type using SNP marker in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum), an important pseudocereal crop, has high economic value due to its nutritional and medicinal properties. However, dehulling of Tartary buckwheat is difficult owing to its thick and tough hull, which has greatly limited the development of the Tartary buckwheat processing industry. The construction of high-resolution genetic maps serves as a basis for identifying quantitative trait loci (QTLs) and qualitative trait genes for agronomic traits. In this study, a recombinant inbred lines (XJ-RILs) population derived from a cross between the easily dehulled Rice-Tartary type and Tartary buckwheat type was genotyped using restriction site-associated DNA (RAD) sequencing to construct a high-density SNP genetic map. Furthermore, QTLs for 1000-grain weight (TGW) and genes controlling hull type were mapped in multiple environments.

RESULTS

In total, 4151 bin markers comprising 122,185 SNPs were used to construct the genetic linkage map. The map consisted of 8 linkage groups and covered 1444.15 cM, with an average distance of 0.35 cM between adjacent bin markers. Nine QTLs for TGW were detected and distributed on four loci on chromosome 1 and 4. A major locus detected in all three trials was mapped in 38.2-39.8 cM region on chromosome 1, with an LOD score of 18.1-37.0, and explained for 23.6-47.5% of the phenotypic variation. The genes controlling hull type were mapped to chromosome 1 between marker Block330 and Block331, which was closely followed by the major locus for TGW. The expression levels of the seven candidate genes controlling hull type present in the region between Block330 and Block336 was low during grain development, and no significant difference was observed between the parental lines. Six non-synonymous coding SNPs were found between the two parents in the region.

CONCLUSIONS

We constructed a high-density SNP genetic map for the first time in Tartary buckwheat. The mapped major loci controlling TGW and hull type will be valuable for gene cloning and revealing the mechanism underlying grain development and easy dehulling, and marker-assisted selection in Tartary buckwheat.

Healthy and Resilient Cereals and Pseudo-Cereals for Marginal Agriculture: Molecular Advances for Improving Nutrient Bioavailability

With the ever-increasing world population, an extra 1.5 billion mouths need to be fed by 2050 with continuously dwindling arable land. Hence, it is imperative that extra food come from the marginal lands that are expected to be unsuitable for growing major staple crops under the adverse climate change scenario. Crop diversity provides right alternatives for marginal environments to improve food, feed, and nutritional security. Well-adapted and climate-resilient crops will be the best fit for such a scenario to produce seed and biomass. The minor millets are known for their high nutritional profile and better resilience for several abiotic stresses that make them the suitable crops for arid and salt-affected soils and poor-quality waters. Finger millet (Eleucine coracana) and foxtail millet (Setaria italica), also considered as orphan crops, are highly tolerant grass crop species that grow well in marginal and degraded lands of Africa and Asia with better nutritional profile. Another category of grains, called pseudo-cereals, is considered as rich foods because of their protein quality and content, high mineral content, and healthy and balance food quality. Quinoa (Chenopodium quinoa), amaranth (Amaranthus sp.), and buckwheat (Fagopyrum esculentum) fall under this category. Nevertheless, both minor millets and pseudo-cereals are morphologically different, although similar for micronutrient bioavailability, and their grains are gluten-free. The cultivation of these millets can make dry lands productive and ensure future food as well as nutritional security. Although the natural nutrient profile of these crop plant species is remarkably good, little development has occurred in advances in molecular genetics and breeding efforts to improve the bioavailability of nutrients. Recent advances in NGS have enabled the genome and transcriptome sequencing of these millets and pseudo-cereals for the faster development of molecular markers and application in molecular breeding. Genomic information on finger millet (1,196 Mb with 85,243 genes); S. italica, a model small millet (well-annotated draft genome of 420 Mb with 38,801 protein-coding genes); amaranth (466 Mb genome and 23,059 protein-coding genes); buckwheat (genome size of 1.12 Gb with 35,816 annotated genes); and quinoa (genome size of 1.5 Gb containing 54,438 protein-coding genes) could pave the way for the genetic improvement of these grains. These genomic resources are an important first step toward genetic improvement of these crops. This review highlights the current advances and available resources on genomics to improve nutrient bioavailability in these five suitable crops for the sustained healthy livelihood.

Metabolite Profiling and Transcriptome Analyses Provide Insights into the Flavonoid Biosynthesis in the Developing Seed of Tartary Buckwheat (Fagopyrum tataricum)

Tartary buckwheat (Fagopyrum tataricum) seeds are rich in flavonoids. However, the detailed flavonoid compositions and the molecular basis of flavonoid biosynthesis in tartary buckwheat seeds remain largely unclear. Here, we performed a combined metabolite profiling and transcriptome analysis to identify flavonoid compositions and characterize genes involved in flavonoid biosynthesis in the developing tartary buckwheat seeds. In total, 234 flavonoids, including 10 isoflavones, were identified. Of these, 80 flavonoids were significantly differential accumulation during seed development. Transcriptome analysis indicated that most structural genes and some potential regulatory genes of flavonoid biosynthesis were significantly differentially expressed in the course of seed development. Correlation analysis between transcriptome and metabolite profiling shown that the expression patterns of some differentially expressed structural genes and regulatory genes were more consistent with the changes in flavonoids profiles during seed development and promoted one SG7 subgroup R2R3-MYB transcription factors (FtPinG0009153900.01) was identified as the key regulatory gene of flavonoid biosynthesis. These findings provide valuable information for understanding the mechanism of flavonoid biosynthesis in tartary buckwheat seeds and the further development of tartary buckwheat health products.

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

Seed development is an essential and complex process, which is involved in seed size change and various nutrients accumulation, and determines crop yield and quality. Common buckwheat (Fagopyrum esculentum Moench) is a widely cultivated minor crop with excellent economic and nutritional value in temperate zones. However, little is known about the molecular mechanisms of seed development in common buckwheat (Fagopyrum esculentum). In this study, we performed RNA-Seq to investigate the transcriptional dynamics and identify the key genes involved in common buckwheat seed development at three different developmental stages. A total of 4619 differentially expressed genes (DEGs) were identified. Based on the results of Gene Ontology (GO) and KEGG analysis of DEGs, many key genes involved in the seed development, including the Ca²⁺ signal transduction pathway, the hormone signal transduction pathways, transcription factors (TFs), and starch biosynthesis-related genes, were identified. More importantly, 18 DEGs were identified as the key candidate genes for seed size through homologous query using the known seed size-related genes from different seed plants. Furthermore, 15 DEGs from these identified as the key genes of seed development were selected to confirm the validity of the data by using quantitative real-time PCR (qRT-PCR), and the results show high consistency with the RNA-Seq results. Taken together, our results revealed the underlying molecular mechanisms of common buckwheat seed development and could provide valuable information for further studies, especially for common buckwheat seed improvement.

Revisiting the versatile buckwheat: reinvigorating genetic gains through integrated breeding and genomics approach

Emerging insights in buckwheat molecular genetics allow the integration of genomics driven breeding to revive this ancient crop of immense nutraceutical potential from Asia. Out of several thousand known edible plant species, only four crops-rice, wheat, maize and potato provide the largest proportion of daily nutrition to billions of people. While these crops are the primary supplier of carbohydrates, they lack essential amino acids and minerals for a balanced nutrition. The overdependence on only few crops makes the future cropping systems vulnerable to the predicted climate change. Diversifying food resources through incorporation of orphan or minor crops in modern cropping systems is one potential strategy to improve the nutritional security and mitigate the hostile weather patterns. One such crop is buckwheat, which can contribute to the agricultural sustainability as it grows in a wide range of environments, requires relatively low inputs and possess balanced amino acid and micronutrient profiles. Additionally, gluten-free nature of protein and nutraceutical properties of secondary metabolites make the crop a healthy alternative of wheat-based diet in developed countries. Despite enormous potential, efforts for the genetic improvement of buckwheat are considerably lagged behind the conventional cereal crops. With the draft genome sequences in hand, there is a great scope to speed up the progress of genetic improvement of buckwheat. This article outlines the state of the art in buckwheat research and provides concrete perspectives how modern breeding approaches can be implemented to accelerate the genetic gain. Our suggestions are transferable to many minor and underutilized crops to address the issue of limited genetic gain and low productivity.

De novo transcriptome assembly and identification of genes related to seed size in common buckwheat (Fagopyrum esculentum M.)

Common buckwheat (Fagopyrum esculentum M.) belongs to the eudicot family Polygonaceae, Fagopyrum Mill, and its seeds have high nutritional value. The mechanism of seed development of common buckwheat remains unclear at the molecular level and no genes related to seed size have been identified. In this study, we performed genome-wide transcriptome sequencing and analysis using common buckwheat seeds at 5 days post anthesis (DPA) and 10 DPA from two cultivars (large-seeded and small-seeded). A total of 259,895 transcripts were assembled, resulting in 187,034 unigenes with average length of 1097 bp and N50 of 1538 bp. Based on gene expression profiles, 9127 differentially expressed genes (DEGs) were identified and analyzed in GO enrichment and KEGG analysis. In addition, genes related to seed size in the IKU pathway, ubiquitin-proteasome pathway, MAPK signaling pathway, TFs and phytohormones were identified and analyzed. AP2 and bZIP transcription factors, BR-signal and ABA were considered to be important regulators of seed size. This study provides a valuable genetic resource for future identification and functional analysis of candidate genes regulating seed size in common buckwheat and will be useful for improving seed yield in common buckwheat through molecular breeding in the future.

Transcriptome analysis reveals the different compatibility between LAAA × AA and LAAA × LL in Lilium

To unveil the mechanism of the compatibility of odd-allotetraploid lily (LAAA) as female with diploid male lily, the differences of expressed unigenes in the ovaries and leaves between LAAA × AA and LAAA × LL were investigated using transcriptome analysis. The results showed the fruits of LAAA × AA well developed, while those of LAAA × LL aborted. The number of differentially expressed genes was less in the ovaries of LAAA × AA than those of LAAA × LL, but it showed opposite trend in those of leaves. The unigenes related with auxins, cytokinins, gibberellins, antioxidants, expansins, chlorophylls, carbohydrates, transport proteins were usually up-expressed in the ovaries and leaves of LAAA × AA but not in LAAA × LL; while those of abscisic acid, ethylene, jasmonic acid, and salicylic acid were increased in the ovaries or leaves of LAAA × LL but not in LAAA × AA. The up-expressed unigenes in the ovaries and leaves of LAAA × AA played positive roles in its fruit development because the products of the genes, like phytohormones and antioxidants, had functions protecting leaves from senescence or scavenging ROS, and thus LAAA was compatible with AA, while those of LAAA × LL played negative roles and caused its fruits aborted, and hence LAAA was incompatible with LL.

Charting oat (Avena sativa) embryo and endosperm transcription factor expression reveals differential expression of potential importance for seed development

Uniquely, oat, among cereals, accumulates an appreciable amount of oil in the endosperm together with starch. Oat is also recognized for its soluble fibers in the form of β-glucans. Despite high and increasing interest in oat yield and quality, the genetic and molecular understanding of oat grain development is still very limited. Transcription factors (TFs) are important regulatory components for plant development, product quality and yield. This study aimed to develop a workflow to determine seed tissue specificity of transcripts encoding transcription factors to reveal differential expression of potential importance for storage compound deposition and quality characters in oat. We created a workflow through the de novo assembly of sequenced seed endosperm and embryo, and publicly available oat seed RNAseq dataset, later followed by TF identification. RNAseq data were assembled into 33,878 transcripts with approximately 90% completeness. A total of 3875 putative TF encoding transcripts were identified from the oat hybrid assemblies. Members of the B3, bHLH, bZIP, C3H, ERF, NAC, MYB and WRKY families were the most abundant TF transcripts. A total of 514 transcripts which were differentially expressed between embryo and endosperm were identified with a threshold of 16-fold expression difference. Among those, 36 TF transcript homologs, belonging to 7 TF families, could be identified through similarity search in wheat embryo and endosperm EST libraries of NCBI Unigene database, and almost all the closest homologs were specifically expressed in seed when explored in WheatExp database. We verified our findings by cloning, sequencing and finally confirming differential expression of two TF encoding transcripts in oat seed embryo and endosperm. The developed workflow for identifying tissue-specific transcription factors allows further functional characterization of specific genes to increase our understanding of grain filling and quality.

Central role of the LEAFY COTYLEDON1 transcription factor in seed development

Seed development is a complex period of the flowering plant life cycle. After fertilization, the three main regions of the seed, embryo, endosperm and seed coat, undergo a series of developmental processes that result in the production of a mature seed that is developmentally arrested, desiccated, and metabolically quiescent. These processes are highly coordinated, both temporally and spatially, to ensure the proper growth and development of the seed. The transcription factor, LEAFY COTYLEDON1 (LEC1), is a central regulator that controls several aspects of embryo and endosperm development, including embryo morphogenesis, photosynthesis, and storage reserve accumulation. Thus, LEC1 regulates distinct sets of genes at different stages of seed development. Despite its critical importance for seed development, an understanding of the mechanisms underlying LEC1's multifunctionality is only beginning to be obtained. Recent studies describe the roles of specific transcription factors and the hormones, gibberellic acid and abscisic acid, in controlling the activity and transcriptional specificity of LEC1 across seed development. Moreover, studies indicate that LEC1 acts as a pioneer transcription factor to promote epigenetic reprogramming during embryogenesis. In this review, we discuss the mechanisms that enable LEC1 to serve as a central regulator of seed development.