Multi-omics analyses of 398 foxtail millet accessions reveal genomic regions associated with domestication, metabolite traits, and anti-inflammatory effects

Genomic and metabolomic insights into the selection and differentiation of bioactive compounds in citrus

Bioactive compounds play an increasingly prominent role in breeding functional and nutritive fruit crops such as citrus. However, the genomic and metabolic bases for the selection and differentiation underlying bioactive compound variations in citrus remain poorly understood. In this study, we constructed a species-level variation atlas of genomes and metabolomes using 299 citrus accessions. A total of 19 829 significant SNPs were targeted to 653 annotated metabolites, among which multiple significant signals were identified for secondary metabolites, especially flavonoids. Significant differential accumulation of bioactive compounds in the phenylpropane pathway, mainly flavonoids and coumarins, was unveiled across ancestral citrus species during differentiation, which is likely associated with the divergent haplotype distribution and/or expression profiles of relevant genes, including p-coumaroyl coenzyme A 2'-hydroxylases, flavone synthases, cytochrome P450 enzymes, prenyltransferases, and uridine diphosphate glycosyltransferases. Moreover, we systematically evaluated the beneficial bioactivities such as the antioxidant and anticancer capacities of 219 citrus varieties, and identified robust associations between distinct bioactivities and specific metabolites. Collectively, these findings provide citrus breeding options for enrichment of beneficial flavonoids and avoidance of potential risk of coumarins. Our study will accelerate the application of genomic and metabolic engineering strategies in developing modern healthy citrus cultivars.

Combined transcriptome and physiological analysis reveals exogenous sucrose enhances photosynthesis and source capacity in foxtail millet

Foxtail millet (Setaria italica (L.) P. Beauv.) is an environmentally friendly crop that meets the current requirements of international food security and is widely accepted as a photosynthesis research model. However, whether exogenous sucrose treatment has a positive effect on foxtail millet growth remains unknown. Here, we employed physiological and molecular approaches to identify photosynthesis and source capacity associated with exogenous sucrose during the growth of Jingu 21 seedlings. RNA-seq analysis showed that some differentially expressed genes (DEGs) related to photosynthesis and carotenoid biosynthesis were induced by exogenous sucrose and that most of these genes were up-regulated. An increase in gas exchange parameters, chlorophyll content, and chlorophyll fluorescence of Jingu 21 was noted after exogenous sucrose addition. Furthermore, exogenous sucrose up-regulated genes encoding sucrose and hexose transporters and enhanced starch and sucrose metabolism. More DEGs were up-regulated by sucrose, the nonstructural carbohydrate (NSC) content in the leaves increased and energy metabolism and sucrose loading subsequently improved, ultimately enhancing photosynthesis under normal and dark conditions. Further analysis revealed that WRKYs, ERFs, HY5, RAP2, and ABI5 could be key transcription factors involved in growth regulation. These results indicate that exogenous sucrose affects the normal photosynthetic performance of foxtail millet by increasing NSC transport and loading. They improve our understanding of the molecular mechanisms of the effects of exogenous sucrose on photosynthesis in foxtail millet, providing an effective measure to enhance source-sink relationships and improve yield.

Cell wall remodeling confers plant architecture with distinct wall structure in Nelumbo nucifera

Lotus (Nelumbo nucifera G.) is a perennial aquatic horticultural plant with diverse architectures. Distinct plant architecture (PA) has certain attractive and practical qualities, but its genetic morphogenesis in lotus remains elusive. In this study, we employ genome-wide association analysis (GWAS) for the seven traits of petiole length (PLL), leaf length (LL), leaf width (LW), peduncle length (PLF), flower diameter (FD), petal length (PeL), and petal width (PeW) in 301 lotus accessions. A total of 90 loci are identified to associate with these traits across 4 years of trials. Meanwhile, we perform RNA sequencing (RNA-seq) to analyze the differential expression of the gene (DEG) transcripts between large and small PA (LPA and SPA) of lotus stems (peduncles and petioles). As a result, eight key candidate genes are identified that are all primarily involved in plant cell wall remodeling significantly associated with PA traits by integrating the results of DEGs and GWAS. To verify this result, we compare the cell wall compositions and structures of LPA versus SPA in representative lotus germplasms. Intriguingly, compared with the SPA lotus, the LPA varieties have higher content of cellulose and hemicellulose, but less filling substrates of pectin and lignin. Additionally, we verified longer cellulose chains and higher cellulose crystallinity with less interference in LPA varieties. Taken together, our study illustrates how plant cell wall remodeling affects PA in lotus, shedding light on the genetic architecture of this significant ornamental trait and offering a priceless genetic resource for future genomic-enabled breeding.

Examining Carotenoid Metabolism Regulation and Its Role in Flower Color Variation in Brassica rapa L

Carotenoids are vital organic pigments that determine the color of flowers, roots, and fruits in plants, imparting them yellow, orange, and red hues. This study comprehensively analyzes carotenoid accumulation in different tissues of the Brassica rapa mutant "YB1", which exhibits altered flower and root colors. Integrating physiological and biochemical assessments, transcriptome profiling, and quantitative metabolomics, we examined carotenoid accumulation in the flowers, roots, stems, and seeds of YB1 throughout its growth and development. The results indicated that carotenoids continued to accumulate in the roots and stems of YBI, especially in its cortex, throughout plant growth and development; however, the carotenoid levels in the petals decreased with progression of the flowering stage. In total, 54 carotenoid compounds were identified across tissues, with 30 being unique metabolites. Their levels correlated with the expression pattern of 22 differentially expressed genes related to carotenoid biosynthesis and degradation. Tissue-specific genes, including CCD8 and NCED in flowers and ZEP in the roots and stems, were identified as key regulators of color variations in different plant parts. Additionally, we identified genes in the seeds that regulated the conversion of carotenoids to abscisic acid. In conclusion, this study offers valuable insights into the regulation of carotenoid metabolism in B. rapa, which can guide the selection and breeding of carotenoid-rich varieties.

Comparative Transcriptomic Analysis Reveals Domestication and Improvement Patterns of Broomcorn Millet (Panicum miliaceum L.)

Broomcorn millet (Panicum miliaceum L.) is one of the earliest crops, domesticated nearly 8000 years ago in northern China. It gradually spread across the entire Eurasian continent, as well as to America and Africa, with recent improvement in various reproductive and vegetative traits. To identify the genes that were selected during the domestication and improvement processes, we performed a comparative transcriptome analysis based on wild types, landraces, and improved cultivars of broomcorn millet at both seeding and filling stages. The variations in gene expression patterns between wild types and landraces and between landraces and improved cultivars were further evaluated to explore the molecular mechanisms underlying the domestication and improvement of broomcorn millet. A total of 2155 and 3033 candidate genes involved in domestication and a total of 84 and 180 candidate genes related to improvement were identified at seedling and filling stages of broomcorn millet, respectively. The annotation results suggested that the genes related to metabolites, stress resistance, and plant hormones were widely selected during both domestication and improvement processes, while some genes were exclusively selected in either domestication or improvement stages, with higher selection pressure detected in the domestication process. Furthermore, some domestication- and improvement-related genes involved in stress resistance either lost their functions or reduced their expression levels due to the trade-offs between stress resistance and productivity. This study provided novel genetic materials for further molecular breeding of broomcorn millet varieties with improved agronomic traits.

Genome-Wide Identification and Characterization of Alternative Oxidase (AOX) Genes in Foxtail Millet (Setaria italica): Insights into Their Abiotic Stress Response

Alternative oxidase (AOX) serves as a critical terminal oxidase within the plant respiratory pathway, playing a significant role in cellular responses to various stresses. Foxtail millet (Setaria italica), a crop extensively cultivated across Asia, is renowned for its remarkable tolerance to abiotic stresses and minimal requirement for fertilizer. In this study, we conducted a comprehensive genome-wide identification of AOX genes in foxtail millet genome, discovering a total of five SiAOX genes. Phylogenetic analysis categorized these SiAOX members into two subgroups. Prediction of cis-elements within the promoter regions, coupled with co-expression network analysis, intimated that SiAOX proteins are likely involved in the plant's adaptive response to abiotic stresses. Employing RNA sequencing (RNA-seq) and real-time quantitative PCR (RT-qPCR), we scrutinized the expression patterns of the SiAOX genes across a variety of tissues and under multiple abiotic stress conditions. Specifically, our analysis uncovered that SiAOX1, SiAOX2, SiAOX4, and SiAOX5 display distinct tissue-specific expression profiles. Furthermore, SiAOX2, SiAOX3, SiAOX4, and SiAOX5 exhibit responsive expression patterns under abiotic stress conditions, with significant differences in expression levels observed between the shoot and root tissues of foxtail millet seedlings. Haplotype analysis of SiAOX4 and SiAOX5 revealed that these genes are in linkage disequilibrium, with Hap_2 being the superior haplotype for both, potentially conferring enhanced cold stress tolerance in the cultivar group. These findings suggest that both SiAOX4 and SiAOX5 may be targeted for selection in future breeding programs aimed at improving foxtail millet's resilience to cold stress.

An integration of genome-wide survey, homologous comparison and gene expression analysis provides a basic framework for the ZRT, IRT-like protein (ZIP) in foxtail millet

Essential mineral elements such as zinc and iron play a crucial role in maintaining crop growth and development, as well as ensuring human health. Foxtail millet is an ancient food crop rich in mineral elements and constitutes an important dietary supplement for nutrient-deficient populations. The ZIP (ZRT, IRT-like protein) transporters are primarily responsible for the absorption, transportation and accumulation of Zn, Fe and other metal ions in plants. Here, we identified 14 ZIP transporters in foxtail millet (SiZIP) and systematically characterized their phylogenetic relationships, expression characteristics, sequence variations, and responses to various abiotic stresses. As a result, SiZIPs display rich spatiotemporal expression characteristics in foxtail millet. Multiple SiZIPs demonstrated significant responses to Fe, Cd, Na, and K metal ions, as well as drought and cold stresses. Based on homologous comparisons, expression characteristics and previous studies, the functions of SiZIPs were predicted as being classified into several categories: absorption/efflux, transport/distribution and accumulation of metal ions. Simultaneously, a schematic diagram of SiZIP was drawn. In general, SiZIPs have diverse functions and extensively involve in the transport of metal ions and osmotic regulation under abiotic stresses. This work provides a fundamental framework for the transport and accumulation of mineral elements and will facilitate the quality improvement of foxtail millet.

Milletomics: a metabolomics centered integrated omics approach toward genetic progression

Producing alternative staple foods like millet will be essential to feeding ten billion people by 2050. The increased demand for millet is driving researchers to improve its genetic variation. Millets include protein, dietary fiber, phenolic substances, and flavonoid components. Its climate resilience makes millet an appealing crop for agronomic sustainability. Integrative omics technologies could potentially identify and develop millets with desirable phenotypes that may have high agronomic value. Millets' salinity and drought tolerance have been enhanced using transcriptomics. In foxtail, finger, and pearl millet, proteomics has discovered salt-tolerant protein, phytohormone-focused protein, and drought tolerance. Metabolomics studies have revealed that certain metabolic pathways including those involving lignin, flavonoids, phenylpropanoid, and lysophospholipids are critical for many processes, including seed germination, photosynthesis, energy metabolism, and the synthesis of bioactive chemicals necessary for drought tolerance. Metabolomics integration with other omics revealed metabolome engineering and trait-specific metabolite creation. Integrated metabolomics and ionomics are still in the development stage, but they could potentially assist in comprehending the pathway of ionomers to control nutrient levels and biofortify millet. Epigenomic analysis has shown alterations in DNA methylation patterns and chromatin structure in foxtail and pearl millets in response to abiotic stress. Whole-genome sequencing utilizing next-generation sequencing is the most proficient method for finding stress-induced phytoconstituent genes. New genome sequencing enables novel biotechnological interventions including genome-wide association, mutation-based research, and other omics approaches. Millets can breed more effectively by employing next-generation sequencing and genotyping by sequencing, which may mitigate climate change. Millet marker-assisted breeding has advanced with high-throughput markers and combined genotyping technologies.

Dissecting the genetic basis of UV-B responsive metabolites in rice

BACKGROUND

UV-B, an important environmental factor, has been shown to affect the yield and quality of rice (Oryza sativa) worldwide. However, the molecular mechanisms underlying the response to UV-B stress remain elusive in rice.

RESULTS

We perform comprehensive metabolic profiling of leaves from 160 diverse rice accessions under UV-B and normal light conditions using a widely targeted metabolomics approach. Our results reveal substantial differences in metabolite accumulation between the two major rice subspecies indica and japonica, especially after UV-B treatment, implying the possible role and mechanism of metabolome changes in subspecies differentiation and the stress response. We next conduct a transcriptome analysis from four representative rice varieties under UV-B stress, revealing genes from amino acid and flavonoid pathways involved in the UV-B response. We further perform a metabolite-based genome-wide association study (mGWAS), which reveals 3307 distinct loci under UV-B stress. Identification and functional validation of candidate genes show that OsMYB44 regulates tryptamine accumulation to mediate UV-B tolerance, while OsUVR8 interacts with OsMYB110 to promote flavonoid accumulation and UV-B tolerance in a coordinated manner. Additionally, haplotype analysis suggests that natural variation of OsUVR8groupA contributes to UV-B resistance in rice.

CONCLUSIONS

Our study reveals the complex biochemical and genetic foundations that govern the metabolite dynamics underlying the response, tolerance, and adaptive strategies of rice to UV-B stress. These findings provide new insights into the biochemical and genetic basis of the metabolome underlying the crop response, tolerance, and adaptation to UV-B stress.

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.

Genome-wide association study reveals genetic loci for ten trace elements in foxtail millet (Setaria italica)

KEY MESSAGE

One hundred and fifty-five QTL for trace element concentrations in foxtail millet were identified using a genome-wide association study, and a candidate gene associated with Ni-Co-Cr concentrations was detected. Foxtail millet (Setaria italica) is an important regional crop known for its rich mineral nutrient content, which has beneficial effects on human health. We assessed the concentrations of ten trace elements (Ba, Co, Cr, Cu, Fe, Mn, Ni, Pb, Sr, and Zn) in the grain of 408 foxtail millet accessions. Significant differences in the concentrations of five elements (Ba, Co, Ni, Sr, and Zn) were observed between two subpopulations of spring- and summer-sown foxtail millet varieties. Moreover, 84.4% of the element pairs exhibited significant correlations. To identify the genetic factors influencing trace element accumulation, a comprehensive genome-wide association study was conducted, identifying 155 quantitative trait locus (QTL) for the ten trace elements across three different environments. Among them, ten QTL were consistently detected in multiple environments, including qZn2.1, qZn4.4, qCr4.1, qFe6.3, qFe6.5, qCo6.1, qPb7.3, qPb7.5, qBa9.1, and qNi9.1. Thirteen QTL clusters were detected for multiple elements, which partially explained the correlations between elements. Additionally, the different concentrations of five elements between foxtail millet subpopulations were caused by the different frequencies of high-concentration alleles associated with important marker-trait associations. Haplotype analysis identified a candidate gene SETIT_036676mg associated with Ni accumulation, with the GG haplotype significantly increasing Ni-Co-Cr concentrations in foxtail millet. A cleaved amplified polymorphic sequence marker (cNi6676) based on the two haplotypes of SETIT_036676mg was developed and validated. Results of this study provide valuable reference information for the genetic research and improvement of trace element content in foxtail millet.

Combined analysis of the metabolome and transcriptome reveals the metabolic characteristics and candidate genes involved in alkaloid metabolism in Heuchera micrantha Douglas ex Lindl

BACKGROUND

Alkaloids, important secondary metabolites produced by plants, play a crucial role in responding to environmental stress. Heuchera micrantha, a well-known plant used in landscaping, has the ability to purify air, and absorb toxic and radioactive substances, showing strong environmental adaptability. However, there is still limited understanding of the accumulation characteristics and metabolic mechanism of alkaloids in H. micrantha.

RESULTS

In this study, four distinct varieties of H. micrantha were used to investigate the accumulation and metabolic traits of alkaloids in its leaves. We conducted a combined analysis of the plant's metabolome and transcriptome. Our analysis identified 44 alkaloids metabolites in the leaves of the four H. micrantha varieties, with 26 showing different levels of accumulation among the groups. The HT and JQ varieties exhibited higher accumulation of differential alkaloid metabolites compared to YH and HY. We annotated the differential alkaloid metabolites to 22 metabolic pathways, including several alkaloid metabolism. Transcriptome data revealed 5064 differentially expressed genes involved in these metabolic pathways. Multivariate analysis showed that four key metabolites (N-hydroxytryptamine, L-tyramine, tryptamine, and 2-phenylethylamine) and three candidate genes (Cluster-15488.116815, Cluster-15488.146268, and Cluster-15488.173297) that merit further investigation.

CONCLUSIONS

This study provided preliminarily insight into the molecular mechanism of the biosynthesis of alkaloids in H. micrantha. However, further analysis is required to elucidate the specific regulatory mechanisms of the candidate gene involved in the synthesis of key alkaloid metabolites. In summary, our findings provide important information about how alkaloid metabolites build up and the metabolic pathways involved in H. micrantha varieties. This gives us a good starting point for future research on the regulation mechanism, and development, and utilization of alkaloids in H. micrantha.

Assembly and comparative analysis of the first complete mitochondrial genome of Setaria italica

MAIN CONCLUSION

In this study, we assembled the first complete mitochondrial genome of Setaria italica and confirmed the multi-branched architecture. The foxtail millet (Setaria italica) holds significant agricultural importance, particularly in arid and semi-arid regions. It plays a pivotal role in diversifying dietary patterns and shaping planting strategies. Although the chloroplast genome of S. italica has been elucidated in recent studies, the complete mitochondrial genome remains largely unexplored. In this study, we employed PacBio HiFi sequencing platforms to sequence and assemble the complete mitochondrial genome. The mitochondrial genome spans a total length of 446,614 base pairs and harbors a comprehensive set of genetic elements, including 33 unique protein-coding genes (PCGs), encompassing 24 unique mitochondrial core genes and 9 variable genes, along with 20 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. Our analysis of mitochondrial PCGs revealed a pronounced codon usage preference. For instance, the termination codon exhibits a marked preference for UAA, while alanine (Ala) exhibits a preference for GCU, and glutamine (Gln) favors CAA. Notably, the maximum Relative Synonymous Codon Usage (RSCU) values for cysteine (Cys) and phenylalanine (Phe) are both below 1.2, indicating a lack of strong codon usage preference for these amino acids. Phylogenetic analyses consistently place S. italica in close evolutionary proximity to Chrysopogon zizanioides, relative to other Panicoideae plants. Collinearity analysis showed that a total of 39 fragments were identified to display homology with both the mitochondrial and chloroplast genomes. A total of 417 potential RNA-editing sites were discovered across the 33 mitochondrial PCGs. Notably, all these editing events involved the conversion of cytosine (C) to uracil (U). Through the employment of PCR validation coupled with Sanger sequencing for the anticipated editing sites of these codons, RNA-editing events were conclusively identified at two specific loci: nad4L-2 and atp6-1030. The results of this study provide a pivotal foundation for advanced genomic breeding research in foxtail millet. Furthermore, they impart essential insights that will be instrumental for forthcoming investigations into the evolutionary and molecular dynamics of Panicoideae species.

Identification, molecular evolution, codon bias, and expansion analysis of NLP transcription factor family in foxtail millet (Setaria italica L.) and closely related crops

The NODULE-INCEPTION-like protein (NLP) family is a plant-specific transcription factor (TF) family involved in nitrate transport and assimilation in plants, which are essential for improving plant nitrogen use efficiency. Currently, the molecular nature and evolutionary trajectory of NLP genes in the C4 model crop foxtail millet are unknown. Therefore, we performed a comprehensive analysis of NLP and molecular evolution in foxtail millet by scanning the genomes of foxtail millet and representative species of the plant kingdom. We identified seven NLP genes in the foxtail millet genome, all of which are individually and separately distributed on different chromosomes. They were not structurally identical to each other and were mainly expressed on root tissues. We unearthed two key genes (Si5G004100.1 and Si6G248300.1) with a variety of excellent characteristics. Regarding its molecular evolution, we found that NLP genes in Gramineae mainly underwent dispersed duplication, but maize NLP genes were mainly generated via WGD events. Other factors such as base mutations and natural selection have combined to promote the evolution of NLP genes. Intriguingly, the family in plants showed a gradual expansion during evolution with more duplications than losses, contrary to most gene families. In conclusion, this study advances the use of NLP genetic resources and the understanding of molecular evolution in cereals.

The role of omics in improving the orphan crop tef

Tef or teff [Eragrostis tef (Zucc.) Trotter] is a cereal crop indigenous to the Horn of Africa, where it is a staple food for a large population. The popularity of tef arises from its resilience to environmental stresses and its nutritional value. For many years, tef has been considered an orphan crop, but recent research initiatives from across the globe are helping to unravel its undisclosed potential. Advanced omics tools and techniques have been directed toward the exploration of tef's diversity with the aim of increasing its productivity. In this review, we report on the most recent advances in tef omics that brought the crop into the spotlight of international research.

Efficient identification of QTL for agronomic traits in foxtail millet (Setaria italica) using RTM- and MLM-GWAS

KEY MESSAGE

Eleven QTLs for agronomic traits were identified by RTM- and MLM-GWAS, putative candidate genes were predicted and two markers for grain weight were developed and validated. Foxtail millet (Setaria italica), the second most cultivated millet crop after pearl millet, is an important grain crop in arid regions. Seven agronomic traits of 408 diverse foxtail millet accessions from 15 provinces in China were evaluated in three environments. They were clustered into two divergent groups based on genotypic data using ADMIXTURE, which was highly consistent with their geographical distribution. Two models for genome-wide association studies (GWAS), namely restricted two-stage multi-locus multi-allele (RTM)-GWAS and mixed linear model (MLM)-GWAS, were used to dissect the genetic architecture of the agronomic traits based on 13,723 SNPs. Eleven quantitative trait loci (QTLs) for seven traits were identified using two models (RTM- and MLM-GWAS). Among them, five were considered stable QTLs that were identified in at least two environments using MLM-GWAS. One putative candidate gene (SETIT_006045mg, Chr4: 744,701-746,852) that can enhance grain weight per panicle was identified based on homologous gene comparison and gene expression analysis and was validated by haplotype analysis of 330 accessions with high-depth (10×) resequencing data (unpublished). In addition, homologous gene comparison and haplotype analysis identified one putative foxtail millet ortholog (SETIT_032906mg, Chr2: 5,020,600-5,029,771) with rice affecting the target traits. Two markers (cGWP6045 and kTGW2906) were developed and validated and can be used for marker-assisted selection of foxtail millet with high grain weight. The results provide a fundamental resource for foxtail millet genetic research and breeding and demonstrate the power of integrating RTM- and MLM-GWAS approaches as a complementary strategy for investigating complex traits in foxtail millet.

Time-Course Transcriptome Analysis of Aquilegia vulgaris Root Reveals the Cell Wall's Roles in Salinity Tolerance

Salt stress has a considerable impact on the development and growth of plants. The soil is currently affected by salinisation, a problem that is becoming worse every year. This means that a significant amount of salt-tolerant plant material needs to be added. Aquilegia vulgaris has aesthetically pleasing leaves, unique flowers, and a remarkable tolerance to salt. In this study, RNA-seq technology was used to sequence and analyse the transcriptome of the root of Aquilegia vulgaris seedlings subjected to 200 mM NaCl treatment for 12, 24, and 48 h. In total, 12 Aquilegia vulgaris seedling root transcriptome libraries were constructed. At the three time points of salt treatment compared with the control, 3888, 1907, and 1479 differentially expressed genes (DEGs) were identified, respectively. Various families of transcription factors (TFs), mainly AP2, MYB, and bHLH, were identified and might be linked to salt tolerance. Gene Ontology (GO) analysis of DEGs revealed that the structure and composition of the cell wall and cytoskeleton may be crucial in the response to salt stress. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the DEGs showed a significant enrichment of the pentose and glucuronate interconversion pathway, which is associated with cell wall metabolism after 24 and 48 h of salt treatment. Based on GO and KEGG analyses of DEGs, the pentose and glucuronate interconversion pathway was selected for further investigation. AP2, MYB, and bHLH were found to be correlated with the functional genes in this pathway based on a correlation network. This study provides the groundwork for understanding the key pathways and gene networks in response to salt stress, thereby providing a theoretical basis for improving salt tolerance in Aquilegia vulgaris.

Multiomics Analysis Reveals the Chemical and Genetic Bases of Pigmented Potato Tuber

Anthocyanins and carotenoids determine the diversity of potato tuber flesh pigmentation; here, the underlying chemical and genetic bases were elucidated by multiomics analyses. A total of 31 anthocyanins and 30 carotenoids were quantified in five differently pigmented tubers. Cyanidin and pelargonidin derivatives determined the redness, while malvidin, petunidin, and delphinidin derivatives contributed to purpleness. Violaxanthin derivatives determined the light-yellow color, while zeaxanthin and antheraxanthin derivatives further enhanced the deep-yellow deposition. Integrated transcriptome and proteome analyses identified that F3'5'H highly enhanced anthocyanin biosynthesis in purple flesh and was responsible for metabolic divergence between red and purple samples. BCH2 significantly enhanced carotenoid biosynthesis in yellow samples and along with ZEP, NCED1, and CCD1 genes determined metabolic divergence between light and deep-yellow samples. The weighted correlation network analysis constructed a regulatory network revealing the central role of AN1 in regulating anthocyanin biosynthesis, and 10 new transcription factors related to anthocyanin and carotenoid metabolism regulation were identified. Our findings provide targeted genes controlling tuber pigmentation, which will be meaningful for the genetic manipulation of tuber quality improvement.

The potentialities of omics resources for millet improvement

Millets are nutrient-rich (nutri-rich) cereals with climate resilience attributes. However, its full productive potential is not realized due to the lack of a focused yield improvement approach, as evidenced by the available literature. Also, the lack of well-characterized genomic resources significantly limits millet improvement. But the recent availability of genomic data and advancement in omics tools has shown its enormous potential to enhance the efficiency and precision faced by conventional breeding in millet improvement. The development of high throughput genotyping platforms based on next-generation sequencing (NGS) has provided a low-cost method for genomic information, specifically for neglected nutri-rich cereals with the availability of a limited number of reference genome sequences. NGS has created new avenues for millet biotechnological interventions such as mutation-based study, GWAS, GS, and other omics technologies. The simultaneous discovery of high-throughput markers and multiplexed genotyping platform has aggressively aided marker-assisted breeding for millet improvement. Therefore, omics technology offers excellent opportunities to explore and combine useful variations for targeted traits that could impart high nutritional value to high-yielding cultivars under changing climatic conditions. In millet improvement, an in-depth account of NGS, integrating genomics data with different biotechnology tools, is reviewed in this context.

MDSi: Multi-omics Database for Setaria italica

BACKGROUND

Foxtail millet (Setaria italica) harbors the small diploid genome (~ 450 Mb) and shows the high inbreeding rate and close relationship to several major foods, feed, fuel and bioenergy grasses. Previously, we created a mini foxtail millet, xiaomi, with an Arabidopsis-like life cycle. The de novo assembled genome data with high-quality and an efficient Agrobacterium-mediated genetic transformation system made xiaomi an ideal C₄ model system. The mini foxtail millet has been widely shared in the research community and as a result there is a growing need for a user-friendly portal and intuitive interface to perform exploratory analysis of the data.

RESULTS

Here, we built a Multi-omics Database for Setaria italica (MDSi, http://sky.sxau.edu.cn/MDSi.htm ), that contains xiaomi genome of 161,844 annotations, 34,436 protein-coding genes and their expression information in 29 different tissues of xiaomi (6) and JG21 (23) samples that can be showed as an Electronic Fluorescent Pictograph (xEFP) in-situ. Moreover, the whole-genome resequencing (WGS) data of 398 germplasms, including 360 foxtail millets and 38 green foxtails and the corresponding metabolic data were available in MDSi. The SNPs and Indels of these germplasms were called in advance and can be searched and compared in an interactive manner. Common tools including BLAST, GBrowse, JBrowse, map viewer, and data downloads were implemented in MDSi.

CONCLUSION

The MDSi constructed in this study integrated and visualized data from three levels of genomics, transcriptomics and metabolomics, and also provides information on the variation of hundreds of germplasm resources that can satisfies the mainstream requirements and supports the corresponding research community.

CandiHap: a haplotype analysis toolkit for natural variation study

Haplotype blocks greatly assist association-based mapping of casual candidate genes by significantly reducing genotyping effort. The gene haplotype could be used to evaluate variants of affected traits captured from the gene region. While there is a rising interest in gene haplotypes, much of the corresponding analysis was carried out manually. CandiHap allows rapid and robust haplotype analysis and candidate identification preselection of candidate causal single-nucleotide polymorphisms and InDels from Sanger or next-generation sequencing data. Investigators can use CandiHap to specify a gene or linkage sites based on genome-wide association studies and explore favorable haplotypes of candidate genes for target traits. CandiHap can be run on computers with Windows, Mac, or UNIX platforms in a graphical user interface or command line, and applied to any species, such as plant, animal, and microbial. The CandiHap software, user manual, and example datasets are freely available at BioCode (https://ngdc.cncb.ac.cn/biocode/tools/BT007080) or GitHub (https://github.com/xukaili/CandiHap).

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-023-01366-4.

Large-scale metabolome analysis reveals dynamic changes of metabolites during foxtail millet grain filling

Compared with traditional staple crops, foxtail millet grain is rich in nutrition and beneficial to human health. Foxtail millet is also tolerance to various abiotic stresses, including drought, making it a good plant for growing in barren land. The study on the composition of metabolites and its dynamics changes during grain development is helpful to understand the process of foxtail millet grain formation. In our study, metabolic and transcriptional analysis were used to uncover the metabolic processes that could influence grain filling in foxtail millet. A total of 2104 known metabolites, belonging to 14 categories, were identified during grain filling. Functional analysis of DAMs and DEGs revealed a stage-specific metabolic properties in foxtail millet grain filling. Some important metabolic processes, such as flavonoid biosynthesis, glutathione metabolism, linoleic acid metabolism, starch and sucrose metabolism and valine, leucine and isoleucine biosynthesis were co-mapped for DEGs and DAMs. Thus, we constructed a gene-metabolite regulatory network of these metabolic pathways to explain their potential functions during grain filling. Our study showed the important metabolic processes during grain filling and focused on the dynamic changes of related metabolites and genes at different stages, which provided a reference for us to better understand and improve foxtail millet grain development and yield.

5M approach to decipher starch-lipid interaction in minor millets

The 5M approach can be applied to understand genetic complexity underlying nutritional traits of minor millets. It will help to systematically identify genomic regions/candidate genes imprinting metabolite profiles.

Protective Effect of Foxtail Millet Protein Hydrolysate on Ethanol and Pyloric Ligation-Induced Gastric Ulcers in Mice

Foxtail millet has been traditionally considered to possess gastroprotective effects, but studies evaluating its use as a treatment for gastric ulcers are lacking. Here, we assessed the antiulcer effects of foxtail millet protein hydrolysate (FPH) and explored its mechanism by using blocking agents. In a mouse model of ethanol-induced gastric ulcers, pretreatment with FPH reduced the ulcerative lesion index, downregulated the expression of inflammatory cytokines in the gastric tissue, increased the activity of antioxidant enzymes, and improved the oxidative status. FPH increased constitutive the activity of nitric oxide synthase (cNOS), NO levels, and mucin expression in gastric mucosa, and inhibited the activation of the ET-1/PI3K/Akt pathway. In a mouse model of pyloric ligation-induced gastric ulcers, FPH inhibited gastric acid secretion and decreased the activity of gastric protease. Pretreatment of mice with the sulfhydryl blocker NEM and the NO synthesis inhibitor L-NAME abolished the gastroprotective effect of FPH, but not the KATP channel blocker glibenclamide and the PGE2 synthesis blocker indomethacin. Among the peptides identified in FPH, 10 peptides were predicted to have regulatory effects on the gastric mucosa, and the key sequences were GP and PG. The results confirmed the gastroprotective effect of FPH and revealed that its mechanism was through the regulation of gastric mucosal mucus and NO synthesis. This study supports the health effects of a millet-enriched diet and provides a basis for millet protein as a functional food to improve gastric ulcers and its related oxidative stress.

Genome-wide identification and characterization of the NPF genes provide new insight into low nitrogen tolerance in Setaria

Introduction

Nitrogen (N) is essential for plant growth and yield production and can be taken up from soil in the form of nitrate or peptides. The NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER family (NPF) genes play important roles in the uptake and transportation of these two forms of N.

Methods

Bioinformatic analysis was used to identify and characterize the NPF genes in Setaria. RNA-seq was employed to analyze time-series low nitrate stress response of the SiNPF genes. Yeast and Arabidopsis mutant complementation were used to test the nitrate transport ability of SiNRT1.1B1 and SiNRT1.1B2.

Results

We identified 92 and 88 putative NPF genes from foxtail millet (Setaria italica L.) and its wild ancestor green foxtail (Setaria viridis L.), respectively. These NPF genes were divided into eight groups according to their sequence characteristics and phylogenetic relationship, with similar intron-exon structure and motifs in the same subfamily. Twenty-six tandem duplication and 13 segmental duplication events promoted the expansion of SiNPF gene family. Interestingly, we found that the tandem duplication of the SiNRT1.1B gene might contribute to low nitrogen tolerance of foxtail millet. The gene expression atlas showed that the SiNPFs were divided into two major clusters, which were mainly expressed in root and the above ground tissues, respectively. Time series transcriptomic analysis further revealed the response of these SiNPF genes to short- and long- time low nitrate stress. To provide natural variation of gene information, we carried out a haplotype analysis of these SiNPFs and identified 2,924 SNPs and 400 InDels based on the re-sequence data of 398 foxtail millet accessions. We also predicted the three-dimensional structure of the 92 SiNPFs and found that the conserved proline 492 residues were not in the substrate binding pocket. The interactions of SiNPF proteins withNO 3 - were analyzed using molecular docking and the pockets were then identified. We found that the SiNPFs-NO 3 - binding energy ranged from -3.8 to -2.7 kcal/mol.

Discussion

Taken together, our study provides a comprehensive understanding of the NPF gene family in Setaria and will contribute to function dissection of these genes for crop breeding aimed at improving high nitrogen use efficiency.

Multi-omics revolution to promote plant breeding efficiency

Crop production is the primary goal of agricultural activities, which is always taken into consideration. However, global agricultural systems are coming under increasing pressure from the rising food demand of the rapidly growing world population and changing climate. To address these issues, improving high-yield and climate-resilient related-traits in crop breeding is an effective strategy. In recent years, advances in omics techniques, including genomics, transcriptomics, proteomics, and metabolomics, paved the way for accelerating plant/crop breeding to cope with the changing climate and enhance food production. Optimized omics and phenotypic plasticity platform integration, exploited by evolving machine learning algorithms will aid in the development of biological interpretations for complex crop traits. The precise and progressive assembly of desire alleles using precise genome editing approaches and enhanced breeding strategies would enable future crops to excel in combating the changing climates. Furthermore, plant breeding and genetic engineering ensures an exclusive approach to developing nutrient sufficient and climate-resilient crops, the productivity of which can sustainably and adequately meet the world's food, nutrition, and energy needs. This review provides an overview of how the integration of omics approaches could be exploited to select crop varieties with desired traits.

The Integration of Genome-Wide Association Study and Homology Analysis to Explore the Genomic Regions and Candidate Genes for Panicle-Related Traits in Foxtail Millet

Panicle traits are important factors affecting yield, and their improvement has long been a critical goal in foxtail millet breeding. In order to understand the genetic basis of panicle formation, a large-scale genome-wide association study (GWAS) was performed in this study for six panicle-related traits based on 706,646 high-polymorphism SNP loci in 407 accessions. As a result, 87 quantitative trait loci (QTL) regions with a physical distance of less than 100 kb were detected to be associated with these traits in three environments. Among them, 27 core regions were stably detected in at least two environments. Based on rice-foxtail millet homologous comparison, expression, and haplotype analysis, 27 high-confidence candidate genes in the QTL regions, such as Si3g11200 (OsDER1), Si1g27910 (OsMADS6), Si7g27560 (GS5), etc., affected panicle-related traits by involving multiple plant growth regulator pathways, a photoperiod response, as well as panicle and grain development. Most of these genes showed multiple effects on different panicle-related traits, such as Si3g11200 affecting all six traits. In summary, this study clarified a strategy based on the integration of GWAS, a homologous comparison, and haplotype analysis to discover the genomic regions and candidate genes for important traits in foxtail millet. The detected QTL regions and candidate genes could be further used for gene clone and marker-assisted selection in foxtail millet breeding.

Genetic characterization of hull color using BSR-Seq and genome re-sequencing approaches in foxtail millet

Hull color of foxtail millet is an important indicator of certain nutritional quality parameters. An F_2:6 recombinant inbred line (RIL) population developed by crossing a yellow-hulled cultivar Yugu 5 and a brown-hulled cultivar Jigu 31 was used to determine the genetic control of the hull color trait. This population segregated for yellow and brown hull colors in a ratio of 2:1, indicating that hull color is regulated by multiple genetic loci. A bulk segregant analysis-RNA sequencing (BSR-Seq) approach performed using the RNA bulks from 30 lines with brown and yellow hull colors each identified three genomic regions on chromosomes 1 (4,570,517-10,698,955 bp), 2 (40,301,380-46,168,003 bp), and 3 (44,469,860-50,532,757 bp). A new QTL for brown hull color of Jigu 31, QHC.czas1, was detected between bin markers Block43 and Block697 on chromosome 1 with the genetic linkage map constructed by re-sequencing a subset of the 147 RILs. This QTL explained a high level of phenotypic variation ranging from 28.0% to 47.0%. The corresponding genomic region of this QTL in the foxtail millet reference genome overlapped with that detected on chromosome 1 by the BSR-Seq analysis. Nineteen genes associated with biosynthesis of anthocyanin were annotated in this genomic region. Gene Si1g06530 encoding a SANT/Myb domain protein was highly expressed in developing panicles and seeds, which warrants further verification as the candidate gene for the brown color hull of Jigu 31. Moreover, several annotated genes for biosynthesis of anthocyanin were identified in the genomic regions of chromosomes 2 and 3.

An optimistic future of C4 crop broomcorn millet (Panicum miliaceum L.) for food security under increasing atmospheric CO2 concentrations

Broomcorn millet, a C₄ cereal, has better tolerance to environmental stresses. Although elevated atmospheric CO₂ concentration has led to grain nutrition reduction in most staple crops, studies evaluating its effects on broomcorn millet are still scarce. The yield, nutritional quality and metabolites of broomcorn millet were investigated under ambient CO₂ (aCO₂, 400 µmol mol^-1) and elevated CO₂ (eCO₂, aCO₂+ 200 µmol mol^-1) for three years using open-top chambers (OTC). The results showed that the yield of broomcorn millet was markedly increased under eCO₂ compared with aCO₂. On average, eCO₂ significantly increased the concentration of Mg (27.3%), Mn (14.6%), and B (21.2%) over three years, whereas it did not affect the concentration of P, K, Fe, Ca, Cu or Zn. Protein content was significantly decreased, whereas starch and oil concentrations were not changed by eCO₂. With the greater increase in grain yield, eCO₂ induced increase in the grain accumulations of P (23.87%), K (29.5%), Mn (40.08%), Ca (22.58%), Mg (51.31%), Zn (40.95%), B (48.54%), starch (16.96%) and oil (28.37%) on average for three years. Flavonoids such as kaempferol, apigenin, eriodictyol, luteolin, and chrysoeriol were accumulated under eCO₂. The reduction in L-glutamine and L-lysine metabolites, which were the most representative amino acid in grain proteins, led to a reduction of protein concentration under eCO₂. Broomcorn millet has more desirable nutritional traits for combating hidden hunger. This may potentially be useful for breeding more nutritious plants in the era of climate change.