Reference genome sequence of the model plant Setaria

Genetic insights in pearl millet breeding in the genomic era: challenges and prospects

Pearl millet, a vital staple food and an important cereal, is emerging as crop having various end-uses as feed, food as well as fodder. Advancement in high-throughput sequencing technology has boosted up pearl millet genomic research in past few years. The available draft genome of pearl millet providing an insight into the advancement of several breeding lines. Comparative and functional genomics have untangled several loci and genes regulating adaptive and agronomic traits in pearl millet. Additionally, the knowledge achieved has far away from being applicable in real breeding practices. We believe that the best path ahead is to adopt genome-based approaches for tailored designing of pearl millet as multi-functional crop with outstanding agronomic traits for various end uses. Presently review highlight several novel concepts and techniques in crop breeding, and summarize the recent advances in pearl millet genomic research, peculiarly genome-wide association dissections of several novel alleles and genes for agronomically important traits.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

Heterologous Expression of SiFBP, a Folate-Binding Protein from Foxtail Millet, Confers Increased Folate Content and Altered Amino Acid Profiles with Nutritional Potential to Arabidopsis

The mechanism underlying folate degradation in foxtail millet grains remains unclear. Here, we identified SiFBP (Setaria italica folate-binding protein) from foxtail millet. A phylogenetic tree revealed that FBPs have close genetic relationships among cereal crop species. Docking analysis and heterologous expression of SiFBP in yeast showed that it could bind folic acid (FA). The SiFBP localized to the plasma membrane in tobacco mesophyll cells by transient expression. In Arabidopsis, it was expressed specifically in the roots and germinating seeds. Overexpressing SiFBP in yeast and Arabidopsis significantly increased folate contents. Untargeted metabolome analysis revealed differentially accumulated metabolites between the transgenic lines (TLs) and wild type (WT); these metabolites were mainly enriched in the amino acid metabolism pathway. The relative contents of lysine and leucine, threonine, and l-methionine were significantly higher in the TLs than in WT. Genes related to the folate and lysine synthesis pathways were upregulated in the TLs. Thus, SiFBP can be used for biofortification of folate and important amino acids in crops via genetic engineering.

Tissue Culture and Somatic Embryogenesis in Warm-Season Grasses—Current Status and Its Applications: A Review

Warm-season grasses are C4 plants and have a high capacity for biomass productivity. These grasses are utilized in many agricultural production systems with their greatest value as feeds for livestock, bioethanol, and turf. However, many important warm-season perennial grasses multiply either by vegetative propagation or form their seeds by an asexual mode of reproduction called apomixis. Therefore, the improvement of these grasses by conventional breeding is difficult and is dependent on the availability of natural genetic variation and its manipulation through breeding and selection. Recent studies have indicated that plant tissue culture system through somatic embryogenesis complements and could further develop conventional breeding programs by micropropagation, somaclonal variation, somatic hybridization, genetic transformation, and genome editing. This review summarizes the tissue culture and somatic embryogenesis in warm-season grasses and focus on current status and above applications including the author’s progress.

The Construction and Exploration of a Comprehensive MicroRNA Centered Regulatory Network in Foxtail Millet (Setaria italica L.)

MicroRNA (miRNA) is an essential endogenous post-transcriptional regulatory factor, and foxtail millet (Setaria italica L.) is an ideal C4 model cereal that is a highly valuable crop in semiarid and arid areas. The Research on comprehensive and high confidence identification and annotation of foxtail millet miRNAs needs to be strengthened, and to our knowledge, there is no information on the regulatory network of foxtail millet miRNA. In this study, 136 high confidence miRNAs were identified through high-throughput sequencing of the small RNAs in seven tissues at the shooting and grain filling stages of foxtail millet. A total of 2,417 target genes were obtained by combining computational biology software and degradome sequencing methods. Furthermore, an analysis using transcriptome sequencing revealed the relationships between miRNAs and their target genes and simultaneously explored key regulatory modules in panicles during the grain filling stage. An miRNA regulatory network was constructed to explore the functions of miRNA in more detail. This network, centered on miRNAs and combining upstream transcriptional factors and downstream target genes, is primarily composed of feed forward loop motifs, which greatly enhances our knowledge of the potential functions of miRNAs and uncovers numerous previously unknown regulatory links. This study provides a solid foundation for research on the function and regulatory network of miRNAs in foxtail millet.

Time-series transcriptomics reveals a drought-responsive temporal network and crosstalk between drought stress and the circadian clock in foxtail millet

Drought stress is a serious factor affecting crop growth and production worldwide. The circadian clock has been identified as key to improving regional adaptability of plants. However, our understanding of the contribution of the circadian clock to drought response and the impacts of drought stress on the circadian clock in plants is still limited. To explore the interactions between the circadian clock and drought stress, foxtail millet seedlings were treated with simulated drought (20% polyethylene glycol-6000) treatment starting at the day (DD) onset zeitgeber time 0 (ZT0, lights on) and at the night (DN) onset zeitgeber time 16 (ZT16, lights off). A high temporal-resolution transcriptomic investigation was performed using DD and DN samples collected at intervals of 2 or 4 h within a 24-h drought-treatment period. Overall, we identified 13 294 drought-responsive genes (DRGs). Among these DRGs, 7931 were common between DD and DN samples, 2638 were specific to DD, and 2725 were specific to DN. Additionally, we identified 1257 circadian genes, of which 67% were DRGs. Interestingly, with drought treatment starting at the day for 8, 12 or 16 h, the circadian phase shifted to 12 h. We also found that the circadian clock led to different day and night drought-responsive pathways. The identification of DRG_Clock (DRG and circadian clock) and DRG_NonClock (DRG and not circadian clock) genes provides a reference for selecting candidate drought resistance genes. Our work reveals the temporal drought-response process and crosstalk between drought stress and the circadian clock in foxtail millet.

Genome-Wide Characterization of DGATs and Their Expression Diversity Analysis in Response to Abiotic Stresses in Brassica napus

Triacylglycerol (TAG) is the most important storage lipid for oil plant seeds. Diacylglycerol acyltransferases (DGATs) are a key group of rate-limiting enzymes in the pathway of TAG biosynthesis. In plants, there are three types of DGATs, namely, DGAT1, DGAT2 and DGAT3. Brassica napus, an allotetraploid plant, is one of the most important oil plants in the world. Previous studies of Brassica napus DGATs (BnaDGATs) have mainly focused on BnaDGAT1s. In this study, four DGAT1s, four DGAT2s and two DGAT3s were identified and cloned from B. napus ZS11. The analyses of sequence identity, chromosomal location and collinearity, phylogenetic tree, exon/intron gene structures, conserved domains and motifs, and transmembrane domain (TMD) revealed that BnaDGAT1, BnaDGAT2 and BnaDGAT3 were derived from three different ancestors and shared little similarity in gene and protein structures. Overexpressing BnaDGATs showed that only four BnaDGAT1s can restore TAG synthesis in yeast H1246 and promote the accumulation of fatty acids in yeast H1246 and INVSc1, suggesting that the three BnaDGAT subfamilies had greater differentiation in function. Transcriptional analysis showed that the expression levels of BnaDGAT1s, BnaDGAT2s and BnaDGAT3s were different during plant development and under different stresses. In addition, analysis of fatty acid contents in roots, stems and leaves under abiotic stresses revealed that P starvation can promote the accumulation of fatty acids, but no obvious relationship was shown between the accumulation of fatty acids with the expression of BnaDGATs under P starvation. This study provides an extensive evaluation of BnaDGATs and a useful foundation for dissecting the functions of BnaDGATs in biochemical and physiological processes.

Evolutionary History and Functional Diversification of the JmjC Domain-Containing Histone Demethylase Gene Family in Plants

Histone demethylases containing JumonjiC (JmjC) domains regulate gene transcription and chromatin structure by changing the methylation status of lysine residues and play an important role in plant growth and development. In this study, a total of 332 JmjC family genes were identified from 21 different plant species. The evolutionary analysis results showed that the JmjC gene was detected in each species, that is, the gene has already appeared in algae. The phylogenetic analysis showed that the KDM3/JHDM2 subfamily genes may have appeared when plants transitioned from water to land, but were lost in lycophytes (Selaginella moellendorffii). During the evolutionary process, some subfamily genes may have been lost in individual species. According to the analysis of the conserved domains, all of the plant JmjC genes contained a typical JmjC domain, which was highly conserved during plant evolution. The analysis of cis-acting elements showed that the promoter region of the JmjC gene was rich in phytohormones and biotic and abiotic stress-related elements. The transcriptome data analysis and protein interaction analyses showed that JmjC genes play an important role in plant growth and development. The results clarified the evolutionary history of JmjC family genes in plants and lay the foundation for the analysis of the biological functions of JmjC family genes.

A transposon insertion drove the loss of natural seed shattering during foxtail millet domestication

Loss of seed shattering was a key step during cereal domestication, and it greatly facilitated seed harvest of the staple cereal foxtail millet (Setaria italica) because the cereal has very small seeds. However, the genetic basis for this loss has been largely unknown. Here, we combined comparative and association mapping to identify an 855-bp Harbinger transposable element insertion in the second exon of the foxtail millet gene shattering1 (sh1) that was responsible for the loss of seed shattering. The sh1 gene encodes zinc finger and YABBY domains. The insert prevents transcription of the second exon, causing partial loss of the zinc finger domain and then loss of natural seed shattering. Specifically, sh1 functions as a transcription repressor and represses the transcription of genes associated with lignin synthesis in the abscission zone, including CAD2. The diversity of sh1 is highly reduced in foxtail millet, consistent with either a severe domestication bottleneck or a selective sweep. Phylogenetic analysis of sh1 further revealed a single origin of foxtail millet in China. Our results support the theories that transposons were the most active factors in genome evolution driving loss of natural seed shattering during foxtail millet domestication and that sh1 underwent parallel selection during domestication across different cereal species.

The role of AUX1 during lateral root development in the domestication of the model C4 grass Setaria italica

C4 photosynthesis increases the efficiency of carbon fixation by spatially separating high concentrations of molecular oxygen from Rubisco. The specialized leaf anatomy required for this separation evolved independently many times. The morphology of C4 root systems is also distinctive and adapted to support high rates of photosynthesis; however, little is known about the molecular mechanisms that have driven the evolution of C4 root system architecture. Using a mutant screen in the C4 model plant Setaria italica, we identify Siaux1-1 and Siaux1-2 as root system architecture mutants. Unlike in S. viridis, AUX1 promotes lateral root development in S. italica. A cell by cell analysis of the Siaux1-1 root apical meristem revealed changes in the distribution of cell volumes in all cell layers and a dependence of the frequency of protophloem and protoxylem strands on SiAUX1. We explore the molecular basis of the role of SiAUX1 in seedling development using an RNAseq analysis of wild-type and Siaux1-1 plants and present novel targets for SiAUX1-dependent gene regulation. Using a selection sweep and haplotype analysis of SiAUX1, we show that Hap-2412TT in the promoter region of SiAUX1 is an allele which is associated with lateral root number and has been strongly selected for during Setaria domestication.

Fine mapping of Rf5 region for a sorghum fertility restorer gene and microsynteny analysis across grass species

Cytoplasmic male sterility (CMS) is widely used to control pollination in the production of commercial F₁ hybrid seed in sorghum. So far, 6 major fertility restorer genes, Rf1 to Rf6, have been reported in sorghum. Here, we fine-mapped the Rf5 locus on sorghum chromosome 5 using descendant populations of a 'Nakei MS-3A' × 'JN43' cross. The Rf5 locus was narrowed to a 140-kb region in BTx623 genome (161-kb in JN43) with 16 predicted genes, including 6 homologous to the rice fertility restorer Rf1 (PPR.1 to PPR.6). These 6 homologs have tandem pentatricopeptide repeat (PPR) motifs. Many Rf genes encode PPR proteins, which bind RNA transcripts and modulate gene expression at the RNA level. No PPR genes were detected at the Rf5 locus on the corresponding homologous chromosome of rice, foxtail millet, or maize, so this gene cluster may have originated by chromosome translocation and duplication after the divergence of sorghum from these species. Comparison of the sequences of these genes between fertile and CMS lines identified PPR.4 as the most plausible candidate gene for Rf5.

Genome-wide identification, phylogenetic and expression pattern analysis of MADS-box family genes in foxtail millet (Setaria italica)

Foxtail millet (Setaria italica) is rich in nutrients and extremely beneficial to human health. We identified and comprehensively analyzed 89 MADS-box genes in the foxtail millet genome. According to the classification of MADS-box genes in Arabidopsis thaliana and rice, the SiMADS-box genes were divided into M-type (37) and MIKC-type (52). During evolution, the differentiation of MIKC-type MADS-box genes occurred before that of monocotyledons and dicotyledons. The SiMADS-box gene structure has undergone much differentiation, and the number of introns in the MIKC-type subfamily is much greater than that in the M-type subfamily. Analysis of gene duplication events revealed that MIKC-type MADS-box gene segmental duplication accounted for the vast majority of gene duplication events, and MIKC-type MADS-box genes played a major role in the amplification of SiMADS-box genes. Collinearity analysis showed highest collinearity between foxtail millet and maize MADS-box genes. Analysis of tissue-specific expression showed that SiMADS-box genes are highly expressed throughout the grain-filling process. Expression analysis of SiMADS-box genes under eight different abiotic stresses revealed many stress-tolerant genes, with induced expression of SiMADS33 and SiMADS78 under various stresses warranting further attention. Further, some SiMADS-box proteins may interact under external stress. This study provides insights for MADS-box gene mining and molecular breeding of foxtail millet in the future.

Elbaiomy R, Chen L, Wang Z, Jiao J, Zhu J, Zhou W, Chen B, Soaud SA, Abbas M, Lin N, El-Sappah AH. Evolution analysis of FRIZZY PANICLE (FZP) orthologs explored the mutations in DNA coding sequences in the grass family (Poaceae)

FRIZZY PANICLE (FZP), an essential gene that controls spikelet differentiation and development in the grass family (Poaceae), prevents the formation of axillary bud meristems and is closely associated with crop yields. It is unclear whether the FZP gene or its orthologs were selected during the evolutionary process of grass species, which possess diverse spike morphologies. In the present study, we adopted bioinformatics methods for the evolutionary analysis of FZP orthologs in species of the grass family. Thirty-five orthologs with protein sequences identical to that of the FZP gene were identified from 29 grass species. Analysis of conserved domains revealed that the AP2/ERF domains were highly conserved with almost no amino acid mutations. However, species of the tribe Triticeae, genus Oryza, and C4 plants exhibited more significant amino acid mutations in the acidic C-terminus region. Results of the phylogenetic analysis showed that the 29 grass species could be classified into three groups, namely, Triticeae, Oryza, and C4 plants. Within the Triticeae group, the FZP genes originating from the same genome were classified into the same sub-group. When selection pressure analysis was performed, significant positive selection sites were detected in species of the Triticeae and Oryza groups. Our results show that the FZP gene was selected during the grass family's evolutionary process, and functional divergence may have already occurred among the various species. Therefore, researchers investigating the FZP gene's functions should take note of the possible presence of various roles in other grass species.

Comparative metabolomic and transcriptomic analysis reveals a coexpression network of the carotenoid metabolism pathway in the panicle of Setaria italica

BACKGROUND

The grains of foxtail millet are enriched in carotenoids, which endow this plant with a yellow color and extremely high nutritional value. However, the underlying molecular regulation mechanism and gene coexpression network remain unclear.

METHODS

The carotenoid species and content were detected by HPLC for two foxtail millet varieties at three panicle development stages. Based on a homologous sequence BLAST analysis, these genes related to carotenoid metabolism were identified from the foxtail millet genome database. The conserved protein domains, chromosome locations, gene structures and phylogenetic trees were analyzed using bioinformatics tools. RNA-seq was performed for these samples to identify differentially expressed genes (DEGs). A Pearson correlation analysis was performed between the expression of genes related to carotenoid metabolism and the content of carotenoid metabolites. Furthermore, the expression levels of the key DEGs were verified by qRT-PCR. The gene coexpression network was constructed by a weighted gene coexpression network analysis (WGCNA).

RESULT

The major carotenoid metabolites in the panicles of DHD and JG21 were lutein and β-carotene. These carotenoid metabolite contents sharply decreased during the panicle development stage. The lutein and β-carotene contents were highest at the S1 stage of DHD, with values of 11.474 μg /100 mg and 12.524 μg /100 mg, respectively. Fifty-four genes related to carotenoid metabolism were identified in the foxtail millet genome. Cis-acting element analysis showed that these gene promoters mainly contain 'plant hormone', 'drought stress resistance', 'MYB binding site', 'endosperm specific' and 'seed specific' cis-acting elements and especially the 'light-responsive' and 'ABA-responsive' elements. In the carotenoid metabolic pathways, SiHDS, SiHMGS3, SiPDS and SiNCED1 were more highly expressed in the panicle of foxtail millet. The expression of SiCMT, SiAACT3, SiPSY1, SiZEP1/2, and SiCCD8c/8d was significantly correlated with the lutein content. The expression of SiCMT, SiHDR, SiIDI2, SiAACT3, SiPSY1, and SiZEP1/2 was significantly correlated with the content of β-carotene. WGCNA showed that the coral module was highly correlated with lutein and β-carotene, and 13 structural genes from the carotenoid biosynthetic pathway were identified. Network visualization revealed 25 intramodular hub genes that putatively control carotenoid metabolism.

CONCLUSION

Based on the integrative analysis of the transcriptomics and carotenoid metabonomics, we found that DEGs related to carotenoid metabolism had a stronger correlation with the key carotenoid metabolite content. The correlation analysis and WGCNA identified and predicted the gene regulation network related to carotenoid metabolism. These results lay the foundation for exploring the key target genes regulating carotenoid metabolism flux in the panicle of foxtail millet. We hope that these target genes could be used to genetically modify millet to enhance the carotenoid content in the future.

A high-density linkage map of finger millet provides QTL for blast resistance and other agronomic traits

Finger millet [Eleusine coracana (L.) Gaertn.] is a critical subsistence crop in eastern Africa and southern Asia but has few genomic resources and modern breeding programs. To aid in the understanding of finger millet genomic organization and genes underlying disease resistance and agronomically important traits, we generated a F_2:3 population from a cross between E. coracana (L.) Gaertn. subsp. coracana accession ACC 100007 and E. coracana (L.) Gaertn. subsp. africana , accession GBK 030647. Phenotypic data on morphology, yield, and blast (Magnaporthe oryzae) resistance traits were taken on a subset of the F_2:3 population in a Kenyan field trial. The F_2:3 population was genotyped via genotyping-by-sequencing (GBS) and the UGbS-Flex pipeline was used for sequence alignment, nucleotide polymorphism calling, and genetic map construction. An 18-linkage-group genetic map consisting of 5,422 markers was generated that enabled comparative genomic analyses with rice (Oryza sativa L.), foxtail millet [Setaria italica (L.) P. Beauv.], and sorghum [Sorghum bicolor (L.) Moench]. Notably, we identified conserved acrocentric homoeologous chromosomes (4A and 4B in finger millet) across all species. Significant quantitative trait loci (QTL) were discovered for flowering date, plant height, panicle number, and blast incidence and severity. Sixteen putative candidate genes that may underlie trait variation were identified. Seven LEUCINE-RICH REPEAT-CONTAINING PROTEIN genes, with homology to nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance proteins, were found on three chromosomes under blast resistance QTL. This high-marker-density genetic map provides an important tool for plant breeding programs and identifies genomic regions and genes of critical interest for agronomic traits and blast resistance.

Improvement of millets in the post-genomic era

Millets are food and nutrient security crops in the semi-arid tropics of developing countries. Crop improvement using modern tools is one of the priority areas of research in millets. The whole-genome sequence (WGS) of millets provides new insight into understanding and studying the genes, genome organization and genomic-assisted improvement of millets. The WGS of millets helps to carry out genome-wide comparison and co-linearity studies among millets and other cereal crops. This approach might lead to the identification of genes underlying biotic and abiotic stress tolerance in millets. The available genome sequence of millets can be used for SNP identification, allele discovery, association and linkage mapping, identification of valuable candidate genes, and marker-assisted breeding (MAB) programs. Next generation sequencing (NGS) technology provides opportunities for genome-assisted breeding (GAB) through genomic selection (GS) and genome-wide association studies (GAWS) for crop improvement. Clustered, regularly interspaced, short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) genome editing (GE) system provides new opportunities for millet improvement. In this review, we discuss the details on the WGS available for millets and highlight the importance of utilizing such resources in the post-genomic era for millet improvement. We also draw inroads on the utilization of various approaches such as GS, GWAS, functional genomics, gene validation and GE for millet improvement. This review might be helpful for understanding the developments in the post-genomic era of millet improvement.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Natural variation in Glume Coverage 1 causes naked grains in sorghum

One of the most critical steps in cereal threshing is the ease with which seeds are detached from sticky glumes. Naked grains with low glume coverage have dramatically increased threshing efficiency and seed quality. Here, we demonstrate that GC1 (Glume Coverage 1), encoding an atypical G protein γ subunit, negatively regulates sorghum glume coverage. Naturally truncated variations of GC1 C-terminus accumulate at higher protein levels and affect the stability of a patatin-related phospholipase SbpPLAII-1. A strong positive selection signature around the GC1 genic region is found in the naked sorghum cultivars. Our findings reveal a crucial event during sorghum domestication through a subtle regulation of glume development by GC1 C-terminus variation, and establish a strategy for future breeding of naked grains.

Low glume coverage is the preferred for easy threshing in grain production, but the genetic basis remains unclear. Here, the authors report the gene GC1, which encodes an atypical G protein γ subunit, negatively regulates sorghum glume coverage and the naturally truncated alleles can be useful in the naked grain breeding.

Orphan Crops: A Best Fit for Dietary Enrichment and Diversification in Highly Deteriorated Marginal Environments

Orphan crops are indigenous and invariably grown by small and marginal farmers under subsistence farming systems. These crops, which are common and widely accepted by local farmers, are highly rich in nutritional profile, good for medicinal purposes, and well adapted to suboptimal growing conditions. However, these crops have suffered neglect and abandonment from the scientific community because of very low or no investments in research and genetic improvement. A plausible reason for this is that these crops are not traded internationally at a rate comparable to that of the major food crops such as wheat, rice, and maize. Furthermore, marginal environments have poor soils and are characterized by extreme weather conditions such as heat, erratic rainfall, water deficit, and soil and water salinity, among others. With more frequent extreme climatic events and continued land degradation, orphan crops are beginning to receive renewed attention as alternative crops for dietary diversification in marginal environments and, by extension, across the globe. Increased awareness of good health is also a major contributor to the revived attention accorded to orphan crops. Thus, the introduction, evaluation, and adaptation of outstanding varieties of orphan crops for dietary diversification will contribute not only to sustained food production but also to improved nutrition in marginal environments. In this review article, the concept of orphan crops vis-à-vis marginality and food and nutritional security is defined for a few orphan crops. We also examined recent advances in research involving orphan crops and the potential of these crops for dietary diversification within the context of harsh marginal environments. Recent advances in genomics coupled with molecular breeding will play a pivotal role in improving the genetic potential of orphan crops and help in developing sustainable food systems. We concluded by presenting a potential roadmap to future research engagement and a policy framework with recommendations aimed at facilitating and enhancing the adoption and sustainable production of orphan crops under agriculturally marginal conditions.

Comparative transcriptome analysis provides insights into grain filling commonalities and differences between foxtail millet [Setaria italica (L.) P. Beauv.] varieties with different panicle types

Grain filling affects grain weight and quality and is among the most critical factors in determining the yield and quality of cereal crops. Though hybrids have larger panicles and numerous spikelets with a larger sink capacity than conventional varieties, data on the grain filling commonalities and differences between foxtail millet varieties with different panicle types remain sparse. In this study, we found that "Zhang Gu 13" (ZG, large panicle) exhibits a significantly higher panicle weight than "Yu Gu 18" (YG, conventional panicle) at the early stage of grain filling, but the weight of YG increased rapidly and gradually overtook ZG during the middle stages. A temporal expression pattern analysis demonstrated that the genes involved in photosynthesis, metabolic pathways, and phenylpropanoid biosynthesis were downregulated, while those related to peroxisome function, purine metabolism, and zeatin biosynthesis were upregulated during grain filling in both varieties. A total of 6,832 differentially expressed genes (DEGs) were identified in both varieties, with the majority identified at the early and late stages. A Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis further revealed that the upregulated DEGs in YG were associated with gibberellin (GA) biosynthesis, ATP-binding cassette (ABC) transporters, and plant hormone signal transduction. Photosynthesis-related DEGs, such as photosystem and antenna proteins, were significantly upregulated in ZG. This study provides preliminary insights into the differences in gene expression and molecular mechanisms of grain filling between ZG and YG in the North China summer-sowing region.

Identification, Molecular Characteristics, and Evolution of GRF Gene Family in Foxtail Millet (Setaria italica L.)

Growth-regulating factor (GRF) is a multigene family that plays a vital role in the growth and development of plants. In the past, the GRF family of many plants has been studied. However, there is not a report about identification and evolution of GRF in foxtail millet (Setaria italia). Here, we identified 10 GRF genes in foxtail millet. Seven (70.00%) were regulated by Sit-miR396, and there were 19 optimal codons in GRFs of foxtail millet. Additionally, we found that WGD or segmental duplication have affected GRFs in foxtail millet between 15.07 and 45.97 million years ago. Regarding the GRF gene family of land plants, we identified a total of 157 GRF genes in 15 representative land plants. We found that GRF gene family originated from Group E, and the GRF gene family in monocots was gradually shrinking. Also, more loss resulted from the small number of GRF genes in lower plants. Exploring the evolution of GRF and functional analysis in the foxtail millet help us to understand GRF better and make a further study about the mechanism of GRF. These results provide a basis for the genetic improvement of foxtail millet and indicate an improvement of the yield.

Plant Growth-promoting Effects of Viable and Dead Spores of Bacillus pumilus TUAT1 on Setaria viridis

Spores are a stress-resistant form of Bacillus spp., which include species that are plant growth-promoting rhizobacteria (PGPR). Previous studies showed that the inoculation of plants with vegetative cells or spores exerted different plant growth-promoting effects. To elucidate the spore-specific mechanism, we compared the effects of viable vegetative cells, autoclaved dead spores, and viable spores of Bacillus pumilus TUAT1 inoculated at 10⁷ CFU plant^–1 on the growth of the C4 model plant, Setaria viridis A10.1. B. pumilus TUAT1 spores exerted stronger growth-promoting effects on Setaria than on control plants 14 days after the inoculation. Viable spores increased shoot weight, root weight, shoot length, root length, and nitrogen uptake efficiency 21 days after the inoculation. These increases involved primary and crown root formation. Additionally, autoclaved dead spores inoculated at 10⁸ or 10⁹ CFU plant^–1 had a positive impact on crown root differentiation, which increased total lateral root length, resulting in a greater biomass and more efficient nitrogen uptake. The present results indicate that an inoculation with viable spores of B. pumilus TUAT1 is more effective at enhancing the growth of Setaria than that with vegetative cells. The plant response to dead spores suggests that the spore-specific plant growth-promoting mechanism is at least partly independent of symbiotic colonization.

Response of Multiple Tissues to Drought Revealed by a Weighted Gene Co-Expression Network Analysis in Foxtail Millet [Setaria italica (L.) P. Beauv.]

Characterization of drought-tolerance mechanisms during the jointing stage in foxtail millet under water-limited conditions is essential for improving the grain yield of this C₄ crop species. In this trial, two drought-tolerant and two drought-sensitive cultivars were examined using transcriptomic dissections of three tissues (root, stem, and leaf) under naturally occurring water-limited conditions. We detected a total of 32,170 expressed genes and characterized 13,552 differentially expressed genes (DEGs) correlated with drought treatment. The majority of DEGs were identified in the root tissue, followed by leaf and stem tissues, and the number of DEGs identified in the stems of drought-sensitive cultivars was about two times higher than the drought-tolerant ones. A total of 127 differentially expressed transcription factors (DETFs) with different drought-responsive patterns were identified between drought-tolerant and drought-sensitive genotypes (including MYB, b-ZIP, ERF, and WRKY). Furthermore, a total of 34 modules were constructed for all expressed genes using a weighted gene co-expression network analysis (WGCNA), and seven modules were closely related to the drought treatment. A total of 1,343 hub genes (including RAB18, LEA14, and RD22) were detected in the drought-related module, and cell cycle and DNA replication-related transcriptional pathways were identified as vital regulators of drought tolerance in foxtail millet. The results of this study provide a comprehensive overview of how Setaria italica copes with drought-inflicted environments during the jointing stage through transcriptional regulating strategies in different organs and lays a foundation for the improvement of drought-tolerant cereal cultivars through genomic editing approaches in the future.

Insights Into the Genetic Architecture of Complex Traits in Napier Grass (Cenchrus purpureus) and QTL Regions Governing Forage Biomass Yield, Water Use Efficiency and Feed Quality Traits

Napier grass is the most important perennial tropical grass native to Sub-Saharan Africa and widely grown in tropical and subtropical regions around the world, primarily as a forage crop for animal feed, but with potential as an energy crop and in a wide range of other areas. Genomic resources have recently been developed for Napier grass that need to be deployed for genetic improvement and molecular dissection of important agro-morphological and feed quality traits. From a diverse set of Napier grass genotypes assembled from two independent collections, a subset of 84 genotypes (although a small population size, the genotypes were selected to best represent the genetic diversity of the collections) were selected and evaluated for 2 years in dry (DS) and wet (WS) seasons under three soil moisture conditions: moderate water stress in DS (DS-MWS); severe water stress in DS (DS-SWS) and, under rainfed (RF) conditions in WS (WS-RF). Data for agro-morphological and feed quality traits, adjusted for the spatial heterogeneity in the experimental blocks, were collected over a 2-year period from 2018 to 2020. A total of 135,706 molecular markers were filtered, after removing markers with missing values >10% and a minor allele frequency (MAF) <5%, from the high-density genome-wide markers generated previously using the genotyping by sequencing (GBS) method of the DArTseq platform. A genome-wide association study (GWAS), using two different mixed linear model algorithms implemented in the GAPIT R package, identified more than 35 QTL regions and markers associated with agronomic, morphological, and water-use efficiency traits. QTL regions governing purple pigmentation and feed quality traits were also identified. The identified markers will be useful in the genetic improvement of Napier grass through the application of marker-assisted selection and for further characterization and map-based cloning of the QTLs.

Recombinant inbred lines and next-generation sequencing enable rapid identification of candidate genes involved in morphological and agronomic traits in foxtail millet

We constructed recombinant inbred lines (RILs) between a Japanese and a Taiwanese landrace of foxtail millet and employed next-generation sequencing, such as flexible ddRAD-seq and Nanopore sequencing to identify the candidate genes involved in the crop evolution of foxtail millet. We successfully constructed a linkage map using flexible ddRAD-seq with parents and RILs and detected major QTLs for each of three traits: leaf sheath colors, spikelet-tipped bristles (stb), and days to heading (DTH). (1) For leaf sheath colors, we identified the C gene on chromosome IV. (2) We identified a homeobox (HOX14) gene for stb on chromosome II, which shows homology with HvVrs1 in barley. (3) Finally, we identified a QTL with a large effect on DTH on chromosome II. A parent of the RILs from Taiwan and Yugu1 had a Harbinger-like TE in intron 3 of this gene. We also investigated the geographical distribution of the TE insertion type of this gene and found that the insertion type is distributed in the northern part of East Asia and intensively in South and Southeast Asia, suggesting that loss/reduction of function of this gene plays an important role in spreading into the northern part of East Asia and subtropical and tropical zones.

Gramene: A Resource for Comparative Analysis of Plants Genomes and Pathways

Gramene is an integrated bioinformatics resource for accessing, visualizing, and comparing plant genomes and biological pathways. Originally targeting grasses, Gramene has grown to host annotations for over 90 plant genomes including agronomically important cereals (e.g., maize, sorghum, wheat, teff), fruits and vegetables (e.g., apple, watermelon, clementine, tomato, cassava), specialty crops (e.g., coffee, olive tree, pistachio, almond), and plants of special or emerging interest (e.g., cotton, tobacco, cannabis, or hemp). For some species, the resource includes multiple varieties of the same species, which has paved the road for the creation of species-specific pan-genome browsers. The resource also features plant research models, including Arabidopsis and C4 warm-season grasses and brassicas, as well as other species that fill phylogenetic gaps for plant evolution studies. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. This chapter outlines system requirements for end-users and database hosting, data types and basic navigation within Gramene, and provides examples of how to (1) explore Gramene's search results, (2) explore gene-centric comparative genomics data visualizations in Gramene, and (3) explore genetic variation associated with a gene locus. This is the first publication describing in detail Gramene's integrated search interface-intended to provide a simplified entry portal for the resource's main data categories (genomic location, phylogeny, gene expression, pathways, and external references) to the most complete and up-to-date set of plant genome and pathway annotations.

Isolation and identification of SiCOL5, which is involved in photoperiod response, based on the quantitative trait locus mapping of Setaria italica

Foxtail millet (Setaria italica) is a versatile grain and fodder crop grown in arid and semi-arid regions. It is an especially important crop for combating malnutrition in certain poverty-stricken areas of the world. Photoperiod sensitivity is a major constraint to the distribution and utilization of foxtail millet germplasm resources. Foxtail millet may be suitable as a model species for studying the photoperiod sensitivity of C₄ crops. However, the genetic basis of the photoperiod response of foxtail millet remains poorly studied. To detect the genetic basis of photoperiod sensitivity-related traits, a recombinant inbred line (RIL) population consisting of 313 lines derived from a cross between the spring-sown cultivar "Longgu 3" and the summer-sown cultivar "Canggu 3" was established. The RIL population was genotyped using whole-genome re-sequencing and was phenotyped in four environments. A high-density genetic linkage map was constructed with an average distance between adjacent markers of 0.69 cM. A total of 21 quantitative trait loci (QTLs) were identified by composite interval mapping, and 116 candidate genes were predicted according to gene annotations and variations between parents, among which three genes were considered important candidate genes by the integration and overall consideration of the results from gene annotation, SNP and indel analysis, cis-element analysis, and the expression pattern of different genes in different varieties, which have different photoperiod sensitivities. A putative candidate gene, SiCOL5, was isolated based on QTL mapping analysis. The expression of SiCOL5 was sensitive to photoperiod and was regulated by biological rhythm-related genes. Function analysis suggested that SiCOL5 positively regulated flowering time. Yeast two-hybrid and bimolecular fluorescence complementation assays showed that SiCOL5 was capable of interacting with SiNF-YA1 in the nucleus.

The C4 cycle and beyond: diverse metabolic adaptations accompany dual-cell photosynthetic functions in Setaria

C4 photosynthesis is typically characterized by the spatial compartmentalization of the photosynthetic reactions into mesophyll (M) and bundle sheath (BS) cells. Initial carbon fixation within M cells gives rise to C4 acids, which are transported to the BS cells. There, C4 acids are decarboxylated so that the resulting CO2 is incorporated into the Calvin cycle. This work is focused on the study of Setaria viridis, a C4 model plant, closely related to several major feed and bioenergy grasses. First, we performed the heterologous expression and biochemical characterization of Setaria isoforms for chloroplastic NADP-malic enzyme (NADP-ME) and mitochondrial NAD-malic enzyme (NAD-ME). The kinetic parameters obtained agree with a major role for NADP-ME in the decarboxylation of the C4 acid malate in the chloroplasts of BS cells. In addition, mitochondria-located NAD-ME showed regulatory properties that could be important in the context of the operation of the C4 carbon shuttle. Secondly, we compared the proteomes of M and BS compartments and found 825 differentially accumulated proteins that could support different metabolic scenarios. Most interestingly, we found evidence of metabolic strategies to insulate the C4 core avoiding the leakage of intermediates by either up-regulation or down-regulation of chloroplastic, mitochondrial, and peroxisomal proteins. Overall, the results presented in this work provide novel data concerning the complexity of C4 metabolism, uncovering future lines of research that will undoubtedly contribute to the expansion of knowledge on this topic.

Evaluation of the molecular and physiological response to dehydration of two accessions of the model plant Setaria viridis

Water deficits are responsible for countless agricultural losses. Among the affected crops, C4 plants are of special interest due to their high water and nitrogen use efficiency. Two accessions of Setaria viridis (Ast-1 and A10.1) with contrasting responses to water deficit were used in the current work to better understand the mechanisms behind drought tolerance in C4 species. Our results showed that although the A10.1 accession exhibited a reduced size and lower Rfd values in comparison to Ast-1, it had overall higher Fv/Fm and lower NPQ values in well-watered conditions. The water deficit induction was performed with PEG-8000 at the grain-filling stage using dehydration cycles. Analysis of physiological measurements showed the A10.1 accession as being more tolerant to multiple water deficit exposures. In addition, PCA identified a clear difference in the pattern of drought response of the accessions. Four drought marker genes previously described in the literature were chosen to evaluate the response at the molecular level: SvP5CS2, SvDHN1, SvNAC6, and SvWRKY1. Besides confirming that Ast-1 is a more sensitive accession, the expression analysis revealed that SvNAC1 might better monitor drought stress, while SvWRKY1 was able to differentiate the two accessions. Distinct evolutionary histories of each accession may be behind their differences in response to water deficits.

Genome-editing in millets: current knowledge and future perspectives

Millets are small seeded cereal crops predominantly cultivated and consumed by resource-poor farmers in the semi-arid tropics of Asia and Africa. Millets possess rich nutrients and a climate resilience property when compared to the other cereals such as rice and wheat. Millet improvement using modern genetic and genomic tools is falling behind other cereal crops due to their cultivation being restricted to less developed countries. Genome editing tools have been successfully applied to major cereal crops and, as a result, many key traits have been introduced into rice, wheat and maize. However, genome editing tools have not yet been used for most millets although they possess rich nutrients. The foxtail millet is the only millet utilised up to now for genome editing works. Limited genomic resources and lack of efficient transformation systems may slow down genome editing in millets. As millets possess many important traits of agricultural importance, high resolution studies with genome editing tools will help to understand the specific mechanism and transfer such traits to major cereals in the future. This review covers the current status of genome editing studies in millets and discusses the future prospects of genome editing in millets to understand key traits of nutrient fortification and develop climate resilient crops in the future.

Genome-Wide DNA Polymorphism Analysis and Molecular Marker Development for the Setaria italica Variety "SSR41" and Positional Cloning of the Setaria White Leaf Sheath Gene SiWLS1

Genome-wide DNA polymorphism analysis and molecular marker development are important for forward genetics research and DNA marker-assisted breeding. As an ideal model system for Panicoideae grasses and an important minor crop in East Asia, foxtail millet (Setaria italica) has a high-quality reference genome as well as large mutant libraries based on the "Yugu1" variety. However, there is still a lack of genetic and mutation mapping tools available for forward genetics research on S. italica. Here, we screened another S. italica genotype, "SSR41", which is morphologically similar to, and readily cross-pollinates with, "Yugu1". High-throughput resequencing of "SSR41" identified 1,102,064 reliable single nucleotide polymorphisms (SNPs) and 196,782 insertions/deletions (InDels) between the two genotypes, indicating that these two genotypes have high genetic diversity. Of the 8,361 high-quality InDels longer than 20 bp that were developed as molecular markers, 180 were validated with 91.5% accuracy. We used "SSR41" and these developed molecular markers to map the white leaf sheath gene SiWLS1. Further analyses showed that SiWLS1 encodes a chloroplast-localized protein that is involved in the regulation of chloroplast development in bundle sheath cells in the leaf sheath in S. italica and is related to sensitivity to heavy metals. Our study provides the methodology and an important resource for forward genetics research on Setaria.

Genome-wide investigation of the GRAS transcription factor family in foxtail millet (Setaria italica L.)

BACKGROUND

GRAS transcription factors perform indispensable functions in various biological processes, such as plant growth, fruit development, and biotic and abiotic stress responses. The development of whole-genome sequencing has allowed the GRAS gene family to be identified and characterized in many species. However, thorough in-depth identification or systematic analysis of GRAS family genes in foxtail millet has not been conducted.

RESULTS

In this study, 57 GRAS genes of foxtail millet (SiGRASs) were identified and renamed according to the chromosomal distribution of the SiGRAS genes. Based on the number of conserved domains and gene structure, the SiGRAS genes were divided into 13 subfamilies via phylogenetic tree analysis. The GRAS genes were unevenly distributed on nine chromosomes, and members of the same subfamily had similar gene structures and motif compositions. Genetic structure analysis showed that most SiGRAS genes lacked introns. Some SiGRAS genes were derived from gene duplication events, and segmental duplications may have contributed more to GRAS gene family expansion than tandem duplications. Quantitative polymerase chain reaction showed significant differences in the expression of SiGRAS genes in different tissues and stages of fruits development, which indicated the complexity of the physiological functions of SiGRAS. In addition, exogenous paclobutrazol treatment significantly altered the transcription levels of DELLA subfamily members, downregulated the gibberellin content, and decreased the plant height of foxtail millet, while it increased the fruit weight. In addition, SiGRAS13 and SiGRAS25 may have the potential for genetic improvement and functional gene research in foxtail millet.

CONCLUSIONS

Collectively, this study will be helpful for further analysing the biological function of SiGRAS. Our results may contribute to improving the genetic breeding of foxtail millet.

Cytochrome P450-catalyzed biosynthesis of furanoditerpenoids in the bioenergy crop switchgrass (Panicum virgatum L.)

Specialized diterpenoid metabolites are important mediators of plant-environment interactions in monocot crops. To understand metabolite functions in plant environmental adaptation that ultimately can enable crop improvement strategies, a deeper knowledge of the underlying species-specific biosynthetic pathways is required. Here, we report the genomics-enabled discovery of five cytochrome P450 monooxygenases (CYP71Z25-CYP71Z29) that form previously unknown furanoditerpenoids in the monocot bioenergy crop Panicum virgatum (switchgrass). Combinatorial pathway reconstruction showed that CYP71Z25-CYP71Z29 catalyze furan ring addition directly to primary diterpene alcohol intermediates derived from distinct class II diterpene synthase products. Transcriptional co-expression patterns and the presence of select diterpenoids in switchgrass roots support the occurrence of P450-derived furanoditerpenoids in planta. Integrating molecular dynamics, structural analysis and targeted mutagenesis identified active site determinants that contribute to the distinct catalytic specificities underlying the broad substrate promiscuity of CYP71Z25-CYP71Z29 for native and non-native diterpenoids.

Phylogenetic relationships in the Sorghum genus based on sequencing of the chloroplast and nuclear genes

Sorghum [Sorghum bicolor (L.) Moench] is an important food crop with a diverse gene pool residing in its wild relatives. A total of 15 sorghum accessions from the unexploited wild gene pool of the Sorghum genus, representing the five subgenera, were sequenced, and the complete chloroplast genomes and 99 common single-copy concatenated nuclear genes were assembled. Annotation of the chloroplast genomes identified a total of 81 protein-coding genes, 38 tRNA, and four rRNA genes. The gene content and gene order among the species was identical. A total of 153 nonsynonymous amino acid changes in 40 genes were identified across the species. Phylogenetic analysis of both the whole chloroplast genome and nuclear genes revealed a similar topology with two distinct clades within the genus. The species within the subgenera Eusorghum, Chaetosorghum, and Heterosorghum clustered in one clade, whereas the species within the subgenera Parasorghum and Stiposorghum clustered in a second clade. However, the subgenera Parasorghum and Stiposorghum were not monophyletic, suggesting the need for further research to resolve the relationships within this group. The close relationship between the two monotypic subgenera Chaetosorghum and Heterosorghum suggests that species within these subgenera could be considered as one group. This analysis provides an improved understanding of the genetic relationships within the Sorghum genus and defines diversity in wild sorghum species that may be useful for crop improvement.

Genome-wide identification and expression analysis of the bHLH transcription factor family and its response to abiotic stress in foxtail millet (Setaria italica L.)

Background

Members of the basic helix-loop-helix (bHLH) transcription factor family perform indispensable functions in various biological processes, such as plant growth, seed maturation, and abiotic stress responses. However, the bHLH family in foxtail millet (Setaria italica), an important food and feed crop, has not been thoroughly studied.

Results

In this study, 187 bHLH genes of foxtail millet (SibHLHs) were identified and renamed according to the chromosomal distribution of the SibHLH genes. Based on the number of conserved domains and gene structure, the SibHLH genes were divided into 21 subfamilies and two orphan genes via phylogenetic tree analysis. According to the phylogenetic tree, the subfamilies 15 and 18 may have experienced stronger expansion in the process of evolution. Then, the motif compositions, gene structures, chromosomal spread, and gene duplication events were discussed in detail. A total of sixteen tandem repeat events and thirty-eight pairs of segment duplications were identified in bHLH family of foxtail millet. To further investigate the evolutionary relationship in the SibHLH family, we constructed the comparative syntenic maps of foxtail millet associated with representative monocotyledons and dicotyledons species. Finally, the gene expression response characteristics of 15 typical SibHLH genes in different tissues and fruit development stages, and eight different abiotic stresses were analysed. The results showed that there were significant differences in the transcription levels of some SibHLH members in different tissues and fruit development stages, and different abiotic stresses, implying that SibHLH members might have different physiological functions.

Conclusions

In this study, we identified 187 SibHLH genes in foxtail millet and further analysed the evolution and expression patterns of the encoded proteins. The findings provide a comprehensive understanding of the bHLH family in foxtail millet, which will inform further studies on the functional characteristics of SibHLH genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08095-y.

Genetics and Genomics Interventions for Promoting Millets as Functional Foods

Several crops, including millets with immense nutritional and therapeutic values, were once a part of our regular diet. However, due to domestication and selection pressures, many of them have become marginally cultivated crops confined to a particular region, race, or locality. Millets are a perfect example of neglected species that have the potential to address both food and nutritional insecurities prevalent among the ever-growing global population. Starvation and malnutrition contribute to a large number of health-related issues, being the main reason behind the occurrence of most of the severe diseases worldwide. These constraints are repeatedly disturbing both the social and economic health of global society. Naturally, millets are rich in minerals, nutrients, and bioactive compounds, and these crops are less dependent on synthetic fertilizers, systemic irrigation, and pest/weed control. Given this, the review emphasizes the nutritional values, health benefits, processing techniques, and genomic advancements of millets. In addition, it proposes a roadmap for enhancing the utility and commercialization of millets.

Improvement of small seed for big nutritional feed

Exploding global population, rapid urbanization, salinization of soils, decreasing arable land availability, groundwater resources, and dynamic climatic conditions pose impending damage to our food security by reducing the grain quality and quantity. This issue is further compounded in arid and semi-arid regions due to the shortage of irrigation water and erratic rainfalls. Millets are gluten (a family of proteins)-free and cultivated all over the globe for human consumption, fuel, feed, and fodder. They provide nutritional security for the under- and malnourished. With the deployment of strategies like foliar spray, traditional/marker-assisted breeding, identification of candidate genes for the translocation of important minerals, and genome-editing technologies, it is now tenable to biofortify important millets. Since the bioavailability of iron and zinc has been proven in human trials, the challenge is to make such grains accessible. This review encompasses nutritional benefits, progress made, challenges being encountered, and prospects of enriching millet crops with essential minerals.

Genome-wide identification and characterization of the BES/BZR gene family in wheat and foxtail millet

Background

BES/BZR family genes have vital roles in plant growth, development, and adaptation to environmental stimuli. However, they have not yet been characterized and systematically analyzed in wheat and foxtail millet.

Results

In the current study, five common and two unique BES/BZR genes were identified by genome-wide analysis in wheat and foxtail millet, respectively. These genes were unevenly distributed on 14 and five chromosomes of wheat and foxtail millet, respectively, and clustered in two subgroups in a phylogenetic analysis. The BES/BZR gene family members in each subgroup contained similar conserved motifs. Investigation of cis-acting elements and expression profile analysis revealed that the BES/BZR genes were predominantly expressed in leaf tissues of wheat and foxtail millet seedlings and responded to brassinosteroid, abscisic acid, and NaCl treatments.

Conclusions

Our results provide a basis for future studies on the function and molecular mechanisms of the BES/BZR gene family in wheat, foxtail millet, and other plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08002-5.

DNA Methylation and RNA-Sequencing Analysis Show Epigenetic Function During Grain Filling in Foxtail Millet (Setaria italica L.)

Grain filling is a crucial process for crop yield and quality. Certain studies already gained insight into the molecular mechanism of grain filling. However, it is unclear whether epigenetic modifications are associated with grain filling in foxtail millet. Global DNA methylation and transcriptome analysis were conducted in foxtail millet spikelets during different stages to interpret the epigenetic effects of the grain filling process. The study employed the whole-genome bisulfite deep sequencing and advanced bioinformatics to sequence and identify all DNA methylation during foxtail millet grain filling; the DNA methylation-mediated gene expression profiles and their involved gene network and biological pathway were systematically studied. One context of DNA methylation, namely, CHH methylation, was accounted for the largest percentage, and it was gradually increased during grain filling. Among all developmental stages, the methylation levels were lowest at T2, followed by T4, which mainly occurred in CHG. The distribution of differentially methylated regions (DMR) was varied in the different genetic regions for three contexts. In addition, gene expression was negatively associated with DNA methylation. Evaluation of the interconnection of the DNA methylome and transcriptome identified some stage-specific differentially expressed genes associated with the DMR at different stages compared with the T1 developmental stage, indicating the potential function of epigenetics on the expression regulation of genes related to the specific pathway at different stages of grain development. The results demonstrated that the dynamic change of DNA methylation plays a crucial function in gene regulation, revealing the potential function of epigenetics in grain development in foxtail millet.

Predicting transcriptional responses to cold stress across plant species

The same gene is often regulated differently in response to stress in even closely related plant species. Directly measuring stress-responsive gene expression can be financially and logistically challenging in nonmodel species. Here, we show that models trained using data on which genes respond to cold in one species can predict which genes will respond to cold in related species, even when the training and target species vary in their degree of tolerance to cold. The prediction models we used require only genomic sequence and gene models. As a result, data from well-studied model species may be used to predict which genes will respond to stress in less-studied species with sequenced genomes.

Although genome-sequence assemblies are available for a growing number of plant species, gene-expression responses to stimuli have been cataloged for only a subset of these species. Many genes show altered transcription patterns in response to abiotic stresses. However, orthologous genes in related species often exhibit different responses to a given stress. Accordingly, data on the regulation of gene expression in one species are not reliable predictors of orthologous gene responses in a related species. Here, we trained a supervised classification model to identify genes that transcriptionally respond to cold stress. A model trained with only features calculated directly from genome assemblies exhibited only modest decreases in performance relative to models trained by using genomic, chromatin, and evolution/diversity features. Models trained with data from one species successfully predicted which genes would respond to cold stress in other related species. Cross-species predictions remained accurate when training was performed in cold-sensitive species and predictions were performed in cold-tolerant species and vice versa. Models trained with data on gene expression in multiple species provided at least equivalent performance to models trained and tested in a single species and outperformed single-species models in cross-species prediction. These results suggest that classifiers trained on stress data from well-studied species may suffice for predicting gene-expression patterns in related, less-studied species with sequenced genomes.

Genetic Diversity and Population Structure of a Rhodes Grass (Chloris gayana) Collection

Rhodes grass (Chloris gayana Kunth) is one of the most important forage grasses used throughout the tropical and subtropical regions of the world. Enhancing the conservation and use of genetic resources requires the development of knowledge and understanding about the existing global diversity of the species. In this study, 104 Rhodes grass accessions, held in trust in the ILRI forage genebank, were characterized using DArTSeq markers to evaluate the genetic diversity and population structure, and to develop representative subsets, of the collection. The genotyping produced 193,988 SNP and 142,522 SilicoDArT markers with an average polymorphic information content of 0.18 and 0.26, respectively. Hierarchical clustering using selected informative markers showed the presence of two and three main clusters using SNP and SilicoDArT markers, respectively, with a cophenetic correction coefficient of 82%. Bayesian population structure analysis also showed the presence of two main subpopulations using both marker types indicating the existence of significant genetic variation in the collection. A representative subset, containing 21 accessions from diverse origins, was developed using the SNP markers. In general, the results revealed substantial genetic diversity in the Rhodes grass collection, and the generated molecular information, together with the developed subset, should help enhance the management, use and improvement of Rhodes grass germplasm in the future.

CHH Methylation Islands: A Nonconserved Feature of Grass Genomes That Is Positively Associated with Transposable Elements but Negatively Associated with Gene-Body Methylation

Methylated CHH (mCHH) islands are peaks of CHH methylation that occur primarily upstream to genes. These regions are actively targeted by the methylation machinery, occur at boundaries between heterochromatin and euchromatin, and tend to be near highly expressed genes. Here we took an evolutionary perspective by studying upstream mCHH islands across a sample of eight grass species. Using a statistical approach to define mCHH islands as regions that differ from genome-wide background CHH methylation levels, we demonstrated that mCHH islands are common and associate with 39% of genes, on average. We hypothesized that islands should be more frequent in genomes of large size, because they have more heterochromatin and hence more need for defined boundaries. We found, however, that smaller genomes tended to have a higher proportion of genes associated with 5' mCHH islands. Consistent with previous work suggesting that islands reflect the silencing of the edge of transposable elements (TEs), genes with nearby TEs were more likely to have mCHH islands. However, the presence of mCHH islands was not a function solely of TEs, both because the underlying sequences of islands were often not homologous to TEs and because genic properties also predicted the presence of 5' mCHH islands. These genic properties included length and gene-body methylation (gbM); in fact, in three of eight species, the absence of gbM was a stronger predictor of a 5' mCHH island than TE proximity. In contrast, gene expression level was a positive but weak predictor of the presence of an island. Finally, we assessed whether mCHH islands were evolutionarily conserved by focusing on a set of 2,720 orthologs across the eight species. They were generally not conserved across evolutionary time. Overall, our data establish additional genic properties that are associated with mCHH islands and suggest that they are not just a consequence of the TE silencing machinery.

The first draft genome of feather grasses using SMRT sequencing and its implications in molecular studies of Stipa

The Eurasian plant Stipa capillata is the most widespread species within feather grasses. Many taxa of the genus are dominants in steppe plant communities and can be used for their classification and in studies related to climate change. Moreover, some species are of economic importance mainly as fodder plants and can be used for soil remediation processes. Although large-scale molecular data has begun to appear, there is still no complete or draft genome for any Stipa species. Thus, here we present a single-molecule long-read sequencing dataset generated using the Pacific Biosciences Sequel System. A draft genome of about 1004 Mb was obtained with a contig N50 length of 351 kb. Importantly, here we report 81,224 annotated protein-coding genes, present 77,614 perfect and 58 unique imperfect SSRs, reveal the putative allopolyploid nature of S. capillata, investigate the evolutionary history of the genus, demonstrate structural heteroplasmy of the chloroplast genome and announce for the first time the mitochondrial genome in Stipa. The assembled nuclear, mitochondrial and chloroplast genomes provide a significant source of genetic data for further works on phylogeny, hybridisation and population studies within Stipa and the grass family Poaceae.

Metabolome-Based Genome-Wide Association Study Provides Genetic Insights Into the Natural Variation of Foxtail Millet

The plant metabolome is considered as a bridge between the genome and the phenome and is essential for the interaction between plant growth and the plant environment. Here, we used the liquid chromatography-tandem mass spectrometry method to perform a widely targeted metabolomics analysis of 150 millet germplasm and simultaneous identification and quantification of 330 annotated metabolites. Comparing the metabolic content of different millets revealed significant natural variation of both primary and secondary metabolites, including flavonoids, phenolamides, hydroxycinnamoyl derivatives, nucleotides, and lipids, in the millets from India and the north and south of China; among them, some of the flavonoids are the most prominent. A total of 2.2 TB sequence data were obtained by sequencing 150 accessions of foxtail millet using the Illumina platform. Further digging into the genetic basis of metabolites by mGWAS analysis found that cyanidin 3-O-glucoside and quercetin O-acetylhexside are concentratedly located at 43.55 Mb on chromosome 5 and 26.9 Mb on chromosome 7, and two Lc were mined as candidate genes, respectively. However, the signals of luteolin 7-O-glucoside and kaempferol 3-O-glucoside were also detected at 14.36 Mb on chromosome 3, and five glycosyltransferase genes on this loci were deemed to regulate their content. Our work is the first research to use mGWAS in millet, and it paves the way for future dissection of complex physiological traits in millet.

Identification of QTL related to anther color and hull color by RAD sequencing in a RIL population of Setaria italica

Background

Foxtail millet (Setaria italica) is one of the oldest domesticated crops and has been considered as an ideal model plant for C₄ grasses. It has abundant type of anther and hull colors which is not only a most intuitive morphological marker for color selection in seed production, but also has very important biological significance for the study of molecular mechanism of regulating the synthesis and metabolism of flavonoids and lignin. However, only a few genetic studies have been reported for anther color and hull color in foxtail millet.

Results

Quantitative trait loci (QTL) analysis for anther color and hull color was conducted using 400 F₆ and F₇ recombinant inbreed lines (RILs) derived from a cross between parents Yugu18 and Jigu19. Using restriction-site associated DNA sequencing, 43,001 single-nucleotide polymorphisms (SNPs) and 3,022 indels were identified between both the parents and the RILs. A total of 1,304 bin markers developed from the SNPs and indels were used to construct a genetic map that spanned 2196 cM of the foxtail millet genome with an average of 1.68 cM/bin. Combined with this genetic map and the phenotypic data observed in two locations for two years, two QTL located on chromosome 6 (Chr6) in a 1.215-Mb interval (33,627,819–34,877,940 bp) for anther color (yellow - white) and three QTL located on Chr1 in a 6.23-Mb interval (1–6,229,734 bp) for hull color (gold-reddish brown) were detected. To narrow the QTL regions identified from the genetic map and QTL analysis, we developed a new method named “inconsistent rate analysis” and efficiently narrowed the QTL regions of anther color into a 60-kb interval (34.13–34.19 Mb) in Chr6, and narrowed the QTL regions of hull color into 70-kb (5.43–5.50 Mb) and 30-kb (5.69–5.72 Mb) intervals in Chr1. Two genes (Seita.6G228600.v2.2 and Seita.6G228700.v2.2) and a cinnamyl alcohol dehydrogenase (CAD) gene (Seita.1G057300.v2.2) with amino acid changes between the parents detected by whole-genome resequencing were identified as candidate genes for anther and hull color, respectively.

Conclusions

This work presents the related QTL and candidate genes of anther and hull color in foxtail millet and developed a new method named inconsistent rate analysis to detect the chromosome fragments linked with the quality trait in RILs. This is the first study of the QTL related to hull color in foxtail millet and clarifying that the CAD gene (Seita.1G057300.v2.2) is the key gene responsible for this trait. It lays the foundation for further cloning of the functional genes and provides a powerful tool to detect the chromosome fragments linked with quality traits in RILs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07882-x.

Efficient genome editing and gene knockout in Setaria viridis with CRISPR/Cas9 directed gene editing by the non-homologous end-joining pathway

The CRISPR/Cas9 system has been used for genome editing in several organisms, including higher plants. This system induces site-specific mutations in the genome based on the nucleotide sequence of engineered guide RNAs. The complex genomes of C4 grasses makes genome editing a challenge in key grass crops like maize (Zea mays), sorghum (Sorghum bicolor), Brachiaria spp., switchgrass (Panicum virgatum), and sugarcane (Saccharum spp.). Setaria viridis is a diploid C4 grass widely used as a model for these C4 crop plants. Here, an optimized CRISPR/Cas9 binary vector that exploits the non-homologous end joining (NHEJ) system was used to knockout a green fluorescent protein (gfp) transgene in S. viridis accession A10.1. Transformation of embryogenic callus by A. tumefaciens generated ten glufosinate-ammonium resistant transgenic events. In the T0 generation, 60% of the events were biallelic mutants in the gfp transgene with no detectable accumulation of GFP protein and without insertions or deletions in predicted off-target sites. The gfp mutations generated by CRISPR/Cas9 were stable and displayed Mendelian segregation in the T1 generation. Altogether, the system described here is a highly efficient genome editing system for S. viridis, an important model plant for functional genomics studies in C4 grasses. Also, this system is a potential tool for improvement of agronomic traits in C4 crop plants with complex genomes.

Correlation and co-localization of QTL for stomatal density, canopy temperature, and productivity with and without drought stress in Setaria

Mechanistic modeling indicates that stomatal conductance could be reduced to improve water use efficiency (WUE) in C4 crops. Genetic variation in stomatal density and canopy temperature was evaluated in the model C4 genus, Setaria. Recombinant inbred lines (RILs) derived from a Setaria italica×Setaria viridis cross were grown with ample or limiting water supply under field conditions in Illinois. An optical profilometer was used to rapidly assess stomatal patterning, and canopy temperature was measured using infrared imaging. Stomatal density and canopy temperature were positively correlated but both were negatively correlated with total above-ground biomass. These trait relationships suggest a likely interaction between stomatal density and the other drivers of water use such as stomatal size and aperture. Multiple quantitative trait loci (QTL) were identified for stomatal density and canopy temperature, including co-located QTL on chromosomes 5 and 9. The direction of the additive effect of these QTL on chromosome 5 and 9 was in accordance with the positive phenotypic relationship between these two traits. This, along with prior experiments, suggests a common genetic architecture between stomatal patterning and WUE in controlled environments with canopy transpiration and productivity in the field, while highlighting the potential of Setaria as a model to understand the physiology and genetics of WUE in C4 species.

Comparative Analysis of Flavonoid Metabolites in Foxtail Millet (Setaria italica) with Different Eating Quality

Foxtail millet (Setaria italica) is an important minor cereal crop in China. The yellow color of the de-husked grain is the most direct aspect for evaluating the foxtail millet quality. The yellow pigment mainly includes carotenoids (lutein and zeaxanthin) and flavonoids. To reveal the diversity and specificity of flavonoids in foxtail millet, we chose three high eating quality and two poor eating quality varieties as research materials. A total of 116 flavonoid metabolites were identified based on Ultra Performance Liquid Chromatography-Electrospray Ionization-Tandem Mass Spectrometry (UPLC-ESI-MS/MS) system. The tested varieties contained similar levels of flavonoid metabolites, but with each variety accumulating its unique flavonoid metabolites. A total of 33 flavonoid metabolites were identified as significantly discrepant between high eating quality and poor eating quality varieties, which were mainly in the flavonoid biosynthesis pathway and one of its branches, the flavone and flavonol biosynthesis pathway. These results showed the diversified components of flavonoids accumulated in foxtail millets and laid the foundation for further research on flavonoids and the breeding for high-quality foxtail millet varieties.

De novo genome assembly of a foxtail millet cultivar Huagu11 uncovered the genetic difference to the cultivar Yugu1, and the genetic mechanism of imazethapyr tolerance

BACKGROUND

Setaria italica is the second-most widely planted species of millets in the world and an important model grain crop for the research of C4 photosynthesis and abiotic stress tolerance. Through three genomes assembly and annotation efforts, all genomes were based on next generation sequencing technology, which limited the genome continuity.

RESULTS

Here we report a high-quality whole-genome of new cultivar Huagu11, using single-molecule real-time sequencing and High-throughput chromosome conformation capture (Hi-C) mapping technologies. The total assembly size of the Huagu11 genome was 408.37 Mb with a scaffold N50 size of 45.89 Mb. Compared with the other three reported millet genomes based on the next generation sequencing technology, the Huagu11 genome had the highest genomic continuity. Intraspecies comparison showed about 94.97 and 94.66% of the Yugu1 and Huagu11 genomes, respectively, were able to be aligned as one-to-one blocks with four chromosome inversion. The Huagu11 genome contained approximately 19.43 Mb Presence/absence Variation (PAV) with 627 protein-coding transcripts, while Yugu1 genomes had 20.53 Mb PAV sequences encoding 737 proteins. Overall, 969,596 Single-nucleotide polymorphism (SNPs) and 156,282 insertion-deletion (InDels) were identified between these two genomes. The genome comparison between Huagu11 and Yugu1 should reflect the genetic identity and variation between the cultivars of foxtail millet to a certain extent. The Ser-626-Aln substitution in acetohydroxy acid synthase (AHAS) was found to be relative to the imazethapyr tolerance in Huagu11.

CONCLUSIONS

A new improved high-quality reference genome sequence of Setaria italica was assembled, and intraspecies genome comparison determined the genetic identity and variation between the cultivars of foxtail millet. Based on the genome sequence, it was inferred that the Ser-626-Aln substitution in AHAS was responsible for the imazethapyr tolerance in Huagu11. The new improved reference genome of Setaria italica will promote the genic and genomic studies of this species and be beneficial for cultivar improvement.