QTL-seq analysis identifies two genomic regions determining the heading date of foxtail millet, Setaria italica (L.) P.Beauv

Realizing visionary goals for the International Year of Millet (IYoM): accelerating interventions through advances in molecular breeding and multiomics resources

MAIN CONCLUSION

Leveraging advanced breeding and multi-omics resources is vital to position millet as an essential "nutricereal resource," aligning with IYoM goals, alleviating strain on global cereal production, boosting resilience to climate change, and advancing sustainable crop improvement and biodiversity. The global challenges of food security, nutrition, climate change, and agrarian sustainability demand the adoption of climate-resilient, nutrient-rich crops to support a growing population amidst shifting environmental conditions. Millets, also referred to as "Shree Anna," emerge as a promising solution to address these issues by bolstering food production, improving nutrient security, and fostering biodiversity conservation. Their resilience to harsh environments, nutritional density, cultural significance, and potential to enhance dietary quality index made them valuable assets in global agriculture. Recognizing their pivotal role, the United Nations designated 2023 as the "International Year of Millets (IYoM 2023)," emphasizing their contribution to climate-resilient agriculture and nutritional enhancement. Scientific progress has invigorated efforts to enhance millet production through genetic and genomic interventions, yielding a wealth of advanced molecular breeding technologies and multi-omics resources. These advancements offer opportunities to tackle prevailing challenges in millet, such as anti-nutritional factors, sensory acceptability issues, toxin contamination, and ancillary crop improvements. This review provides a comprehensive overview of molecular breeding and multi-omics resources for nine major millet species, focusing on their potential impact within the framework of IYoM. These resources include whole and pan-genome, elucidating adaptive responses to abiotic stressors, organelle-based studies revealing evolutionary resilience, markers linked to desirable traits for efficient breeding, QTL analysis facilitating trait selection, functional gene discovery for biotechnological interventions, regulatory ncRNAs for trait modulation, web-based platforms for stakeholder communication, tissue culture techniques for genetic modification, and integrated omics approaches enabled by precise application of CRISPR/Cas9 technology. Aligning these resources with the seven thematic areas outlined by IYoM catalyzes transformative changes in millet production and utilization, thereby contributing to global food security, sustainable agriculture, and enhanced nutritional consequences.

Latitudinal adaptation and dispersal pathway of foxtail millet suggested by geographical distribution of transposable elements inserted in the SiPRR37 gene

We investigated the variation and geographical distribution of the Pseudo-regulator response 37 (Setaria italica PRR37; SiPRR37) gene, which is involved in heading time (photoperiodism) in foxtail millet. An allele of the SiPRR37 gene, in which an approximately 4.9-kb transposable element (designated TE1) is inserted (a loss-of-function or reduction-of-function type), is distributed sporadically in East Asia and broadly in Southeast Asia and South Asia, implying that this gene is important in latitudinal adaptation. In addition, we found a new allele of SiPRR37 with an insertion of a 360-bp TE (TE2) at this locus and investigated the geographical distribution of this new type. This SiPRR37 allele with TE2 is distributed in Japan, Korea, Nepal, Iran and Turkey. Both TE1 and TE2 are useful markers for tracing foxtail millet dispersal pathways in Asia.

Crop Evolution of Foxtail Millet

Studies on the domestication, genetic differentiation, and crop evolution of foxtail millet are reviewed in this paper. Several genetic studies were carried out to elucidate the genetic relationships among foxtail millet accessions originating mainly from Eurasia based on intraspecific hybrid pollen semi-sterility, isozymes, DNA markers, and single-nucleotide polymorphisms. Most studies suggest that China is the center of diversity of foxtail millet, and landraces were categorized into geographical groups. These results indicate that this millet was domesticated in China and spread over Eurasia, but independent origin in other regions cannot be ruled out. Furthermore, the evolution of genes was reviewed (i.e., the Waxy gene conferring amylose content in the endosperm, the Si7PPO gene controlling polyphenol oxidase, the HD1 and SiPRR37 genes controlling heading time, the Sh1 and SvLes1 genes involved in grain shattering, and the C gene controlling leaf sheath pigmentation), and the variation and distribution of these genes suggested complex patterns of evolution under human and/or natural selection.

Multi-environment GWAS identifies genomic regions underlying grain nutrient traits in foxtail millet (Setaria italica)

KEY MESSAGE

A total of 104 foxtail millet accessions were evaluated for 11 nutrients in three environments and 67 high-confidence marker-trait associations (MTAs) were identified. Six SNPs showed pleiotropic effect and associated with two or more nutrients, whereas 24 candidate genes were identified for 28 MTAs involving seven traits. Millets are known for their better nutritional profiles compared to major cereals. Foxtail millet (Setaria italica) is rich in nutrients essential to circumvent malnutrition and hidden hunger. However, the genetic determinants underlying this trait remain elusive. In this context, we evaluated 104 diverse foxtail millet accessions in three different environments (E1, E2, and E3) for 11 nutrients and genotyped with 30K SNPs. The genome-wide association study showed 67 high-confidence (Bonferroni-corrected) marker-trait associations (MTAs) for the nutrients except for phosphorus. Six pleiotropic SNPs were also identified, which were associated with two or more nutrients. Around 24 candidate genes (CGs) were identified for 28 MTAs involving seven nutrients. A total of 17 associated SNPs were present within the gene region, and five (5) were mapped in the exon of the CGs. Significant SNPs, desirable alleles and CGs identified in the present study will be useful in breeding programmes for trait improvement.

Revitalization of small millets for nutritional and food security by advanced genetics and genomics approaches

Small millets, also known as nutri-cereals, are smart foods that are expected to dominate food industries and diets to achieve nutritional security. Nutri-cereals are climate resilient and nutritious. Small millet-based foods are becoming popular in markets and are preferred for patients with celiac and diabetes. These crops once ruled as food and fodder but were pushed out of mainstream cultivation with shifts in dietary habits to staple crops during the green revolution. Nevertheless, small millets are rich in micronutrients and essential amino acids for regulatory activities. Hence, international and national organizations have recently aimed to restore these lost crops for their desirable traits. The major goal in reviving these crops is to boost the immune system of the upcoming generations to tackle emerging pandemics and disease infestations in crops. Earlier periods of civilization consumed these crops, which had a greater significance in ethnobotanical values. Along with nutrition, these crops also possess therapeutic traits and have shown vast medicinal use in tribal communities for the treatment of diseases like cancer, cardiovascular disease, and gastrointestinal issues. This review highlights the significance of small millets, their values in cultural heritage, and their prospects. Furthermore, this review dissects the nutritional and therapeutic traits of small millets for developing sustainable diets in near future.

Multi-omics intervention in Setaria to dissect climate-resilient traits: Progress and prospects

Millets constitute a significant proportion of underutilized grasses and are well known for their climate resilience as well as excellent nutritional profiles. Among millets, foxtail millet (Setaria italica) and its wild relative green foxtail (S. viridis) are collectively regarded as models for studying broad-spectrum traits, including abiotic stress tolerance, C₄ photosynthesis, biofuel, and nutritional traits. Since the genome sequence release, the crop has seen an exponential increase in omics studies to dissect agronomic, nutritional, biofuel, and climate-resilience traits. These studies have provided first-hand information on the structure, organization, evolution, and expression of several genes; however, knowledge of the precise roles of such genes and their products remains elusive. Several open-access databases have also been instituted to enable advanced scientific research on these important crops. In this context, the current review enumerates the contemporary trend of research on understanding the climate resilience and other essential traits in Setaria, the knowledge gap, and how the information could be translated for the crop improvement of related millets, biofuel crops, and cereals. Also, the review provides a roadmap for studying other underutilized crop species using Setaria as a model.

Exploring the SiCCT Gene Family and Its Role in Heading Date in Foxtail Millet

CCT transcription factors are involved in the regulation of photoperiod and abiotic stress in Arabidopsis and rice. It is not clear that how CCT gene family expand and regulate heading date in foxtail millet. In this study, we conducted a systematic analysis of the CCT gene family in foxtail millet. Thirty-nine CCT genes were identified and divided into four subfamilies based on functional motifs. Analysis showed that dispersed duplication played a predominant role in the expansion of CCT genes during evolution. Nucleotide diversity analysis suggested that genes in CONSTANS (COL)-like, CCT MOTIF FAMILY (CMF)-like, and pseudoresponse response regulator (PRR)-like subfamilies were subjected to selection. Fifteen CCT genes were colocalized with previous heading date quantitative trait loci (QTL) and genome-wide association analysis (GWAS) signals. Transgenic plants were then employed to confirm that overexpression of the CCT gene SiPRR37 delayed the heading date and increased plant height. Our study first investigated the characterization and expansion of the CCT family in foxtail millet and demonstrated the role of SiPRR37. These results lay a significant foundation for further research on the function of CCT genes and provide a cue for the regulation of heading date.

Genome-wide identification of quantitative trait loci for morpho-agronomic and yield-related traits in foxtail millet (Setaria italica) across multi-environments

Foxtail millet (Setaria italica) is an ideal model of genetic system for functional genomics of the Panicoideae crop. Identification of QTL responsible for morpho-agronomic and yield-related traits facilitates dissection of genetic control and breeding in cereal crops. Here, based on a Yugu1 × Longgu7 RIL population and genome-wide resequencing data, an updated linkage map harboring 2297 bin and 74 SSR markers was constructed, spanning 1315.1 cM with an average distance of 0.56 cM between adjacent markers. A total of 221 QTL for 17 morpho-agronomic and yield-related traits explaining 5.5 ~ 36% of phenotypic variation were identified across multi-environments. Of these, 109 QTL were detected in two to nine environments, including the most stable qLMS6.1 harboring a promising candidate gene Seita.6G250500, of which 70 were repeatedly identified in different trials in the same geographic location, suggesting that foxtail millet has more identical genetic modules under the similar ecological environment. One hundred-thirty QTL with overlapping intervals formed 22 QTL clusters. Furthermore, six superior recombinant inbred lines, RIL35, RIL48, RIL77, RIL80, RIL115 and RIL125 with transgressive inheritance and enrichment of favorable alleles in plant height, tiller, panicle morphology and yield related-traits were screened by hierarchical cluster. These identified QTL, QTL clusters and superior lines lay ground for further gene-trait association studies and breeding practice in foxtail millet.

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.

Rapid identification of candidate genes controlling male-sterility in Foxtail millet (Setaria italica)

Photo-(thermo-) sensitive genic male-sterile line is the key component of two-line hybridization system in foxtail millet (Setaria italica), but the genetic basis of male sterility in most male-sterile lines is still unclear. In the present study, a large F2 population was developed derived from a cross between the photo-(thermo-) sensitive male-sterile line A2 and the fertile-line 1484-5. Thirty plants with extreme high and extreme low fertility were selected from the population to construct a sterile DNA pool and a fertile DNA pool, respectively. Sequencing both DNA pools and data analysis revealed that two QTLs conferred male-sterility, qSiMS6.1 with a major effect and qSiMS6.2 with a minor effect, on chromosome 6. Both QTLs exhibited complete dominance. The major QTL, qSiMS6.1, was delimited to a 186-kb interval between the markers SiM20 and SiM9 by the joint analysis of QTL-seq and QTL mapping with SSR and structure variation markers. Millet_GLEAN_10020454 in this region is the most likely candidate gene for qSiMS6.1 since it is predicted to encode a male-sterile 5 like protein. These results lay a solid foundation for qSiMS6.1 cloning and provided gene resources for breeding new male-sterile lines.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01269-2.

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.

High-depth resequencing of 312 accessions reveals the local adaptation of foxtail millet

Based on the high-density variation map, we identified genome-level evidence for local adaptation and demonstrated that Siprr37 with transposon insertion contributes to the fitness of foxtail millet in the northeastern ecoregion. Adaptation is a robust way through which plants are able to overcome environmental constraints. The mechanisms of adaptation in heterogeneous natural environments are largely unknown. Deciphering the genomic basis of local adaptation will contribute to further improvement in domesticated plants. To this end, we describe a high-depth (19.4 ×) haplotype map of 3.02 million single nucleotide polymorphisms in foxtail millet (Setaria italica) from whole-genome resequencing of 312 accessions. In the genome-wide scan, we identified a set of improvement signals (including the homologous gene of OsIPA1, a key gene controlling ideal plant architecture) related to the geographical adaptation to four ecoregions in China. In particular, based on the genome-wide association analysis results, we identified the contribution of a pseudo-response regulator gene, SiPRR37, to heading date adaptation in foxtail millet. We observed the expression changes of SiPRR37 resulted from a key Tc1-Mariner transposon insertion in the first intron. Positive selection analyses revealed that SiPRR37 mainly contributed to the adaptation of northeastern ecoregions. Taken together, foxtail millet adapted to the northeastern region by regulating the function of SiPRR37, which sheds lights on genome-level evidence for adaptive geographical divergence. Besides, our data provide a nearly complete catalog of genomic variation aiding the identification of functionally important variants.

Barnyard Millet for Food and Nutritional Security: Current Status and Future Research Direction

Barnyard millet (Echinochloa species) has become one of the most important minor millet crops in Asia, showing a firm upsurge in world production. The genus Echinochloa comprises of two major species, Echinochloa esculenta and Echinochloa frumentacea, which are predominantly cultivated for human consumption and livestock feed. They are less susceptible to biotic and abiotic stresses. Barnyard millet grain is a good source of protein, carbohydrate, fiber, and, most notably, contains more micronutrients (iron and zinc) than other major cereals. Despite its nutritional and agronomic benefits, barnyard millet has remained an underutilized crop. Over the past decades, very limited attempts have been made to study the features of this crop. Hence, more concerted research efforts are required to characterize germplasm resources, identify trait-specific donors, develop mapping population, and discover QTL/gene (s). The recent release of genome and transcriptome sequences of wild and cultivated Echinochloa species, respectively has facilitated in understanding the genetic architecture and decoding the rapport between genotype and phenotype of micronutrients and agronomic traits in this crop. In this review, we highlight the importance of barnyard millet in the current scenario and discuss the up-to-date status of genetic and genomics research and the research gaps to be worked upon by suggesting directions for future research to make barnyard millet a potential crop in contributing to food and nutritional security.

Genome-Wide Association Study of Major Agronomic Traits in Foxtail Millet (Setaria italica L.) Using ddRAD Sequencing

Foxtail millet (Setaria italica), the second largest cultivated millet crop after pearl millet, is utilized for food and forage globally. Further, it is also considered as a model crop for studying agronomic, nutritional and biofuel traits. In the present study, a genome-wide association study (GWAS) was performed for ten important agronomic traits in 142 foxtail millet core eco-geographically diverse genotypes using 10 K SNPs developed through GBS-ddRAD approach. Number of SNPs on individual chromosome ranged from 844 (chromosome 5) to 2153 (chromosome 8) with an average SNP frequency of 25.9 per Mb. The pairwise linkage disequilibrium (LD) estimated using the squared-allele frequency correlations was found to decay rapidly with the genetic distance of 177 Kb. However, for individual chromosome, LD decay distance ranged from 76 Kb (chromosome 6) to 357 Kb (chromosome 4). GWAS identified 81 MTAs (marker-trait associations) for ten traits across the genome. High confidence MTAs for three important agronomic traits including FLW (flag leaf width), GY (grain yield) and TGW (thousand-grain weight) were identified. Significant pyramiding effect of identified MTAs further supplemented its importance in breeding programs. Desirable alleles and superior genotypes identified in the present study may prove valuable for foxtail millet improvement through marker-assisted selection.