A high density genetic map and QTL for agronomic and yield traits in 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.

Crop Safety and Weed Control of Foliar Application of Penoxsulam in Foxtail Millet

Grass damage has become an important factor restricting foxtail millet production; chemical weeding can help resolve this issue. However, special herbicides in foxtail millet fields are lacking. Penoxsulam has a broad weed control spectrum and a good control effect. In this project, Jingu 21 was used as the test material, and five different concentrations of penoxsulam were used for spraying test in the three-five leaf stage. In this experiment, the effects on the growth of foxtail millet were discussed by measuring the agronomic characters and antioxidant capacity of foxtail millet after spraying penoxsulam. The results showed that: (1) penoxsulam is particularly effective in controlling Amaranthus retroflexus L. (A. retroflexus) and Echinochloa crus-galli (L.) Beauv. (E. crus-galli), but is ineffective in controlling Chenopodium album L. (C. album) and Digitaria sanguinalis (L.) Scop. (D. sanguinalis); (2) the stem diameter, fresh weight, and dry weight of the above-ground parts decreased with the increase in spraying amount; (3) as the spraying dosage increased, the superoxide (SOD), peroxidase (POD), and catalase (CAT) activities in the foxtail millet initially increased and subsequently decreased; the malonaldehyde (MDA) content increased. Our experiment found that 1/2X and 1X spraying dosages had certain application value in controlling gramineous weeds in foxtail millet field. Other spraying dosages are not recommended as they may harm the crops. Our findings provide reference for identifying new herbicides in the foxtail millet field.

Construction of a high-density genetic map and QTL localization of body weight and wool production related traits in Alpine Merino sheep based on WGR

BACKGROUND

The Alpine Merino is a new breed of fine-wool sheep adapted to the cold and arid climate of the plateau in the world. It has been popularized in Northwest China due to its superior adaptability as well as excellent production performance. Those traits related to body weight, wool yield, and wool fiber characteristics, which are economically essential traits in Alpine Merino sheep, are controlled by QTL (Quantitative Trait Loci). Therefore, the identification of QTL and genetic markers for these key economic traits is a critical step in establishing a MAS (Marker-Assisted Selection) breeding program.

RESULTS

In this study, we constructed the high-density genetic linkage map of Alpine Merino sheep by sequencing 110 F1 generation individuals using WGR (Whole Genome Resequencing) technology. 14,942 SNPs (Single Nucleotide Polymorphism) were identified and genotyped. The map spanned 2,697.86 cM, with an average genetic marker interval of 1.44 cM. A total of 1,871 high-quality SNP markers were distributed across 27 linkage groups, with an average of 69 markers per LG (Linkage Group). Among them, the smallest genetic distance is 19.62 cM for LG2, while the largest is 237.19 cM for LG19. The average genetic distance between markers in LGs ranged from 0.24 cM (LG2) to 3.57 cM (LG17). The marker density in the LGs ranged from LG14 (39 markers) to LG1 (150 markers).

CONCLUSIONS

The first genetic map of Alpine Merino sheep we constructed included 14,942 SNPs, while 46 QTLs associated with body weight, wool yield and wool fiber traits were identified, laying the foundation for genetic studies and molecular marker-assisted breeding. Notably, there were QTL intervals for overlapping traits on LG4 and LG8, providing potential opportunities for multi-trait co-breeding and further theoretical support for selection and breeding of ultra-fine and meaty Alpine Merino sheep.

Genetic linkage map construction and QTL analysis for plant height in proso millet (Panicum miliaceum L.)

KEY MESSAGE

A genetic linkage map representing proso millet genome was constructed with SSR markers, and a major QTL corresponding to plant height was mapped on chromosome 14 of this map. Proso millet (Panicum miliaceum L.) has the lowest water requirements of all cultivated cereal crops. However, the lack of a genetic map and the paucity of genomic resources for this species have limited the utility of proso millet for detailed genetic studies and hampered genetic improvement programs. In this study, 97,317 simple sequence repeat (SSR) markers were developed based on the genome sequence of the proso millet landrace Longmi 4. Using some of these markers in conjunction with previously identified SSRs, an SSR-based linkage map for proso millet was successfully constructed using a large mapping population (316 F2 offspring). In total, 186 SSR markers were assigned to 18 linkage groups corresponding to the haploid chromosomes. The constructed map had a total length of 3033.42 centimorgan (cM) covering 78.17% of the assembled reference genome. The length of the 18 linkage groups ranged from 88.89 cM (Chr. 15) to 274.82 cM (Chr. 16), with an average size of 168.17 cM. To our knowledge, this is the first genetic linkage map for proso millet based on SSR markers. Plant height is one of the most important traits in crop improvement. A major QTL was repeatedly detected in different environments, explaining 8.70-24.50% of the plant height variations. A candidate gene affecting auxin biosynthesis and transport, and ROS homeostasis regulation was predicted. Thus, the linkage map and QTL analysis provided herein will promote the development of gene mining and molecular breeding in proso millet.

Comparative transcriptome analysis and genetic dissection of vegetative branching traits in foxtail millet (Setaria italica)

KEY MESSAGE

Two major genetic loci, qTN5.1 and qAB9.1, were identified and finely mapped to the 255 Kb region with one potential candidate gene for tiller number and the 521 Kb region with eight candidate genes for axillary branch number, respectively. Vegetative branching including tillering and axillary branching are vital traits affecting both the plant architecture and the biomass in cereal crops. However, the mechanism underlying the formation of vegetative branching in foxtail millet is largely unknown. Here, a foxtail millet cultivar and its bushy wild relative Setaria viridis accession were used to construct segregating populations to identify candidate genes regulating tiller number and axillary branch number. Transcriptome analysis using vegetative branching bud samples of parental accessions was performed, and key differentially expressed genes and pathways regulating vegetative branching were pointed out. Bulk segregant analysis on their F2:3 segregating population was carried out, and a major QTL for tiller number (qTN5.1) and two major QTLs for axillary branch number (qAB2.1 and qAB9.1) were detected. Fine-mapping strategy was further performed on F2:4 segregate population, and Seita.5G356600 encoding β-glucosidase 11 was identified as the promising candidate gene for qTN5.1, and eight genes, especially Seita.9G125300 and Seita.9G125400 annotated as B-S glucosidase 44, were finally identified as candidate genes for regulating axillary branching. Findings in this study will help to elucidate the genetic basis of the vegetative branching formation of foxtail millet and lay a foundation for breeding foxtail millet varieties with ideal vegetative branching numbers.

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.

Construction of a High-Density Paulownia Genetic Map and QTL Mapping of Important Phenotypic Traits Based on Genome Assembly and Whole-Genome Resequencing

Quantitative trait locus (QTL) mapping based on a genetic map is a very effective method of marker-assisted selection in breeding, and whole-genome resequencing is one of the useful methods to obtain high-density genetic maps. In this study, the hybrid assembly of Illumina, PacBio, and chromatin interaction mapping data was used to construct high-quality chromosomal genome sequences of Paulownia fortunei, with a size of 476.82 Mb, a heterozygosity of 0.52%, and a contig and scaffold N50s of 7.81 Mb and 21.81 Mb, respectively. Twenty scaffolds with a total length of 437.72 Mb were assembled into 20 pseudochromosomes. Repeat sequences with a total length of 243.96 Mb accounted for 51.16% of the entire genome. In all, 26,903 protein-coding gene loci were identified, and 26,008 (96.67%) genes had conserved functional motifs. Further comparative genomics analysis preliminarily showed that the split of P. fortunei with Tectona grandis likely occurred 38.8 (33.3-45.1) million years ago. Whole-genome resequencing was used to construct a merged genetic map of 20 linkage groups, with 2993 bin markers (3,312,780 SNPs), a total length of 1675.14 cm, and an average marker interval of 0.56 cm. In total, 73 QTLs for important phenotypic traits were identified (19 major QTLs with phenotypic variation explained ≥ 10%), including 10 for the diameter at breast height, 7 for the main trunk height, and 56 for branch-related traits. These results not only enrich P. fortunei genomic data but also form a solid foundation for fine QTL mapping and key marker/gene mining of Paulownia, which is of great significance for the directed genetic improvement of these species.

The Era of Plant Breeding: Conventional Breeding to Genomics-assisted Breeding for Crop Improvement

Plant breeding has made a significant contribution to increasing agricultural production. Conventional breeding based on phenotypic selection is not effective for crop improvement. Because phenotype is considerably influenced by environmental factors, which will affect the selection of breeding materials for crop improvement. The past two decades have seen tremendous progress in plant breeding research. Especially the availability of high-throughput molecular markers followed by genomic-assisted approaches significantly contributed to advancing plant breeding. Integration of speed breeding with genomic and phenomic facilities allowed rapid quantitative trait loci (QTL)/gene identifications and ultimately accelerated crop improvement programs. The advances in sequencing technology helps to understand the genome organization of many crops and helped with genomic selection in crop breeding. Plant breeding has gradually changed from phenotype-to-genotype-based to genotype-to-phenotype-based selection. High-throughput phenomic platforms have played a significant role in the modern breeding program and are considered an essential part of precision breeding. In this review, we discuss the rapid advance in plant breeding technology for efficient crop improvements and provide details on various approaches/platforms that are helpful for crop improvement. This review will help researchers understand the recent developments in crop breeding and improvements.

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.

Effects of exogenous auxin on yield in foxtail millet (Setaria italica L.) when applied at the grain-filling stage

Foxtail millet (Setaria italica L.) is of high nutritious value, which is an important crop in arid and semi-arid regions. The objective of this experiment was to explore the effects of the synthetic auxin naphthalene acetic acid (NAA) on the physiological processes of foxtail millet, and to provide a theoretical basis and technical approaches for its efficient use in millet cultivation. Two foxtail millet varieties ('Jingu 21' and 'Zhangzagu 5') were treated with six concentrations of NAA from 0-144 mg L^-1 at the grain-filling stage in field experiments. The photosynthetic pigment contents, gas exchange parameters, chlorophyll fluorescence parameters, and grain yield were measured in foxtail millet. The results showed that low concentrations of NAA (18-36 mg L^-1) increased the contents of photosynthetic pigments, and increased the activities of antioxidant enzymes, the photosynthetic rate, and the activity of photosystem system II (PS II). At higher NAA concentrations, the facilitation effect of the treatments diminished, showing a clear concentration effect. In this study, yield was significantly and positively correlated with PS II effective quantum yield (Y(II)) and the PSII electron transport rate (ETR), and the net photosynthetic rate (Pn) was significantly and positively correlated with chlorophyll content, stomatal conductance (Gs), Y(II), and ETR. These results also indicated that exogenous NAA application promotes the production of ATP and NADPH by increasing the efficiency of electron transfer within the photosystems and also improved photochemical utilization, which facilitates the fixation and reduction of carbon, ultimately leading to an increase in Pn and increasing grain yield in foxtail millet.

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

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

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.

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.

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.

Fine mapping of qPH9, a major quantitative trait locus, responsible for plant height in foxtail millet [Setaria italica (L.) P. Beauv.]

Plant height is vital for crop yield by influencing plant architecture and resistance to lodging. Although lots of quantitative trait loci (QTLs) controlling plant height had been mapped in foxtail millet, their contributions to phenotypic variation were generally small and mapping regions were relatively large, indicating the difficult application in molecular breeding using marker-assisted selection (MAS). In the present paper, a total of 23 QTLs involving in 15 traits were identified via a high-density Bin map containing 3024 Bin markers with an average distance of 0.48 cM through an F2 population. Among them, qPH9, with a large phenotypic variation explained (51.6%) related to plant height, was one of the major QTLs. Furthermore, qPH9 was repeatedly detected in multi-environments under field conditions using two new developed F2 populations from the same F1 plant, and was narrowed down to a smaller interval of 281 kb using 1024 recessive F2 individuals from the same F1 plant. Finally, we found that there was an extremely significant correlation between marker MRI1016 and plant height, and further speculated that Seita.9G088900 and Seita.9G089700 could be key candidates of qPH9. This study laid an important foundation for the cloning of qPH9 and molecular breeding of dwarf varieties via MAS.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01261-w.

Genetic control and phenotypic characterization of panicle architecture and grain yield-related traits in foxtail millet (Setaria italica)

Multi-environment QTL mapping identified 23 stable loci and 34 co-located QTL clusters for panicle architecture and grain yield-related traits, which provide a genetic basis for foxtail millet yield improvement. Panicle architecture and grain weight, both of which are influenced by genetic and environmental factors, have significant effects on grain yield potential. Here, we used a recombinant inbred line (RIL) population of 333 lines of foxtail millet, which were grown in 13 trials with varying environmental conditions, to identify quantitative trait loci (QTL) controlling nine agronomic traits related to panicle architecture and grain yield. We found that panicle weight, grain weight per panicle, panicle length, panicle diameter, and panicle exsertion length varied across different geographical locations. QTL mapping revealed 159 QTL for nine traits. Of the 159 QTL, 34 were identified in 2 to 12 environments, suggesting that the genetic control of panicle architecture in foxtail millet is sensitive to photoperiod and/or other environmental factors. Eighty-eight QTL controlling different traits formed 34 co-located QTL clusters, including the triple QTL cluster qPD9.2/qPL9.5/qPEL9.3, which was detected 23 times in 13 environments. Several candidate genes, including Seita.2G388700, Seita.3G136000, Seita.4G185300, Seita.5G241500, Seita.5G243100, Seita.9G281300, and Seita.9G342700, were identified in the genomic intervals of multi-environmental QTL or co-located QTL clusters. Using available phenotypic and genotype data, we conducted haplotype analysis for Seita.2G002300 and Seita.9G064000,which showed high correlations with panicle weight and panicle exsertion length, respectively. These results not only provided a basis for further fine mapping, functional studies and marker-assisted selection of traits related to panicle architecture in foxtail millet, but also provide information for comparative genomics analyses of cereal crops.

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.

Genome wide association study of agronomic and seed traits in a world collection of proso millet (Panicum miliaceum L.)

BACKGROUND

The climate crisis threatens sustainability of crop production worldwide. Crop diversification may enhance food security while reducing the negative impacts of climate change. Proso millet (Panicum milaceum L.) is a minor cereal crop which holds potential for diversification and adaptation to different environmental conditions. In this study, we assembled a world collection of proso millet consisting of 88 varieties and landraces to investigate its genomic and phenotypic diversity for seed traits, and to identify marker-trait associations (MTA).

RESULTS

Sequencing of restriction-site associated DNA fragments yielded 494 million reads and 2,412 high quality single nucleotide polymorphisms (SNPs). SNPs were used to study the diversity in the collection and perform a genome wide association study (GWAS). A genotypic diversity analysis separated accessions originating in Western Europe, Eastern Asia and Americas from accessions sampled in Southern Asia, Western Asia, and Africa. A Bayesian structure analysis reported four cryptic genetic groups, showing that landraces accessions had a significant level of admixture and that most of the improved proso millet materials clustered separately from landraces. The collection was highly diverse for seed traits, with color varying from white to dark brown and width spanning from 1.8 to 2.6 mm. A GWAS study for seed morphology traits identified 10 MTAs. In addition, we identified three MTAs for agronomic traits that were previously measured on the collection.

CONCLUSION

Using genomics and automated seed phenotyping, we elucidated phylogenetic relationships and seed diversity in a global millet collection. Overall, we identified 13 MTAs for key agronomic and seed traits indicating the presence of alleles with potential for application in proso breeding programs.

Genome-wide identification of PTI1 family in Setaria italica and salinity-responsive functional analysis of SiPTI1-5

BACKGROUND

PTI1 (Pto-interacting 1) protein kinase belongs to the receptor-like cytoplasmic kinase (RLCK) group of receptor-like protein kinases (RLK), but lack extracellular and transmembrane domains. PTI1 was first identified in tomato (Solanum lycopersicum) and named SlPTI1, which has been reported to interact with bacterial effector Pto, a serine/threonine protein kinase involved in plant resistance to bacterial disease. Briefly, the host PTI1 specifically recognizes and interacts with the bacterial effector AvrPto, which triggers hypersensitive cell death to inhibit the pathogen growth in the local infection site. Previous studies have demonstrated that PTI1 is associated with oxidative stress and hypersensitivity.

RESULTS

We identified 12 putative PTI1 genes from the genome of foxtail millet (Setaria italica) in this study. Gene replication analysis indicated that both segmental replication events played an important role in the expansion of PTI1 gene family in foxtail millet. The PTI1 family members of model plants, i.e. S. italica, Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), maize (Zea mays), S. lycopersicum, and soybean (Glycine max), were classified into six major categories according to the phylogenetic analysis, among which the PTI1 family members in foxtail millet showed higher degree of homology with those of rice and maize. The analysis of a complete set of SiPTI1 genes/proteins including classification, chromosomal location, orthologous relationships and duplication. The tissue expression characteristics revealed that SiPTI1 genes are mainly expressed in stems and leaves. Experimental qRT-PCR results demonstrated that 12 SiPTI1 genes were induced by multiple stresses. Subcellular localization visualized that all of foxtail millet SiPTI1s were localized to the plasma membrane. Additionally, heterologous expression of SiPTI1-5 in yeast and E. coli enhanced their tolerance to salt stress.

CONCLUSIONS

Our results contribute to a more comprehensive understanding of the roles of PTI1 protein kinases and will be useful in prioritizing particular PTI1 for future functional validation studies in foxtail millet.

QTL mapping for foxtail millet plant height in multi-environment using an ultra-high density bin map

Using a fixed RIL population derived from a widely used foxtail millet backbone breeding line and an elite cultivar, we constructed a high-density bin map and identified six novel multi-environment effect QTLs and seven candidate genes for dwarf phenotype. Plant height is an important trait that determines tradeoffs between competition and resource allocation, which is crucial for yield potential. To improve the C4 model plant foxtail millet (Setaria italica) productivity, it is necessary to isolate plant height-related genes that contribute to ideal plant architecture in breeding. In the present study, we generated a foxtail millet population of 333 recombinant inbred lines (RILs) derived from a cross between a backbone line Ai 88 and an elite cultivar Liaogu 1. We evaluated plant height in 13 environmental conditions across 4 years, the mean plant height of the RIL population ranged from 89.5 to 149.9 cm. Using deep re-sequencing data, we constructed a high-density bin map with 3744 marker bins. Quantitative trait locus (QTL) mapping identified 26 QTLs significantly associated with plant height. Of these, 13 QTLs were repeatedly detected under multiple environments, including six novel QTLs that have not been reported before. Seita.1G242300, a gene encodes gibberellin 2-oxidase-8, which was detected in nine environments in a 1.54-Mb interval of qPH1.3, was considered as an important candidate gene. Moreover, other six genes involved in GA biosynthesis or signaling pathways, and fifteen genes encode F-box domain proteins which might function as E3 ligases, were also considered as candidate genes in different QTLs. These QTLs and candidate genes identified in this study will help to elucidate the genetic basis of foxtail millet plant height, and the linked markers will be useful for marker-assistant selection of varieties with ideal plant architecture and high yield potential.

Small Millets for Enduring Food Security Amidst Pandemics

Food security is an ongoing problem, and current staple foods are not sufficient to overcome challenges such as the present COVID-19 pandemic. We propose here that small millets have the potential to become new staple crops, especially in hunger hotspots. Currently, the absence of intensification of millet farming, lack of deployment of genetic tools for trait improvement, and the need for optimization of storage and supply chains limit crop production. We highlight a roadmap to strengthen small millet cultivation, such as identifying varieties suitable for particular environments and targeting trait improvement using genetic and genomic approaches. These approaches will help to combat hunger and malnutrition and also economically benefit the farmers and stakeholders involved in small millet cultivation amidst the ongoing pandemic.

SSRgenotyper: A simple sequence repeat genotyping application for whole-genome resequencing and reduced representational sequencing projects

PREMISE

Many programs can identify simple sequence repeat (SSR) motifs in genomic data. SSRgenotyper extends SSR identification to en masse genotyping from resequencing data for diversity panels and linkage mapping populations.

METHODS AND RESULTS

SSRgenotyper will find and genotype SSRs from SAM files and an SSR reference FASTA. Several outputs are possible, including a simple table with the SSR marker name, position, and SSR alleles, defined by the repeat number of the repeat motif. Specific output files include a GENEPOP-formatted file for downstream genetic diversity analyses and a traditional A, H, B mapping file output that is phased to the parents of the population for biparental linkage map construction. Linkage maps produced using SSRgenotyper genotypes were highly collinear with physical maps and correctly inferred known phylogenies.

CONCLUSIONS

SSRgenotyper provides an easy-to-use, accurate, and scalable SSR genotyping platform for whole-genome resequencing data. SSRgenotyper is freely available at https://github.com/dlewis27/SSRgenotyper.

High-density genetic map construction and identification of loci controlling flower-type traits in Chrysanthemum (Chrysanthemum × morifolium Ramat.)

Flower type is an important and extremely complicated trait of chrysanthemum. The corolla tube merged degree (CTMD) and the relative number of ray florets (RNRF) are the two key factors affecting chrysanthemum flower type. However, few reports have clarified the inheritance of these two complex traits, which limits directed breeding for flower-type improvement. In this study, 305 F1 hybrids were obtained from two parents with obvious differences in CTMD and RNRF performance. Using specific-locus amplified fragment sequencing (SLAF-seq) technology, we constructed a high-density genetic linkage map with an average map distance of 0.76 cM. Three major QTLs controlling CTMD and four major QTLs underlying RNRF were repeatedly detected in the 2 years. Moreover, the synteny between the genetic map and other Compositae species was investigated, and weak collinearity was observed. In QTL regions with a high degree of genomic collinearity, eight annotated genes were probed in the Helianthus annuus L. and Lactuca sativa L. var. ramosa Hort. genomes. Furthermore, 20 and 11 unigenes were identified via BLAST searches between the SNP markers of the QTL regions and the C. vestitum and C. lavandulifolium transcriptomes, respectively. These results lay a foundation for molecular marker-assisted breeding and candidate gene exploration in chrysanthemum without a reference assembly.

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.

Genomics-assisted breeding in minor and pseudo-cereals

Minor and pseudo-cereals, which can grow with lower input and often produce specific nutrients compared to major cereal crops, are attracting worldwide attention. Since these crops generally have a large genetic diversity in a breeding population, rapid genetic improvement can be possible by the application of genomics-assisted breeding methods. In this review, we discuss studies related to biparental quantitative trait locus (QTL) mapping, genome-wide association study, and genomic selection for minor and pseudo-cereals. Especially, we focus on the current progress in a pseudo-cereal, buckwheat. Prospects for the practical utilization of genomics-assisted breeding in minor and pseudo-cereals are discussed including the issues to overcome especially for these crops.

High-Resolution Linkage Map With Allele Dosage Allows the Identification of Regions Governing Complex Traits and Apospory in Guinea Grass (Megathyrsus maximus)

Forage grasses are mainly used in animal feed to fatten cattle and dairy herds, and guinea grass (Megathyrsus maximus) is considered one of the most productive of the tropical forage crops that reproduce by seeds. Due to the recent process of domestication, this species has several genomic complexities, such as autotetraploidy and aposporous apomixis. Consequently, approaches that relate phenotypic and genotypic data are incipient. In this context, we built a linkage map with allele dosage and generated novel information of the genetic architecture of traits that are important for the breeding of M. maximus. From a full-sib progeny, a linkage map containing 858 single nucleotide polymorphism (SNP) markers with allele dosage information expected for an autotetraploid was obtained. The high genetic variability of the progeny allowed us to map 10 quantitative trait loci (QTLs) related to agronomic traits, such as regrowth capacity and total dry matter, and 36 QTLs related to nutritional quality, which were distributed among all homology groups (HGs). Various overlapping regions associated with the quantitative traits suggested QTL hotspots. In addition, we were able to map one locus that controls apospory (apo-locus) in HG II. A total of 55 different gene families involved in cellular metabolism and plant growth were identified from markers adjacent to the QTLs and APOSPORY locus using the Panicum virgatum genome as a reference in comparisons with the genomes of Arabidopsis thaliana and Oryza sativa. Our results provide a better understanding of the genetic basis of reproduction by apomixis and traits important for breeding programs that considerably influence animal productivity as well as the quality of meat and milk.

QTL mapping of yield component traits on bin map generated from resequencing a RIL population of foxtail millet (Setaria italica)

Background

Foxtail millet (Setaria italica) has been developed into a model genetical system for deciphering architectural evolution, C₄ photosynthesis, nutritional properties, abiotic tolerance and bioenergy in cereal grasses because of its advantageous characters with the small genome size, self-fertilization, short growing cycle, small growth stature, efficient genetic transformation and abundant diverse germplasm resources. Therefore, excavating QTLs of yield component traits, which are closely related to aspects mentioned above, will further facilitate genetic research in foxtail millet and close cereal species.

Results

Here, 164 Recombinant inbreed lines from a cross between Longgu7 and Yugu1 were created and 1,047,978 SNPs were identified between both parents via resequencing. A total of 3413 bin markers developed from SNPs were used to construct a binary map, containing 3963 recombinant breakpoints and totaling 1222.26 cM with an average distance of 0.36 cM between adjacent markers. Forty-seven QTLs were identified for four traits of straw weight, panicle weight, grain weight per plant and 1000-grain weight. These QTLs explained 5.5–14.7% of phenotypic variance. Thirty-nine favorable QTL alleles were found to inherit from Yugu1. Three stable QTLs were detected in multi-environments, and nine QTL clusters were identified on Chromosome 3, 6, 7 and 9.

Conclusions

A high-density genetic map with 3413 bin markers was constructed and three stable QTLs and 9 QTL clusters for yield component traits were identified. The results laid a powerful foundation for fine mapping, identifying candidate genes, elaborating molecular mechanisms and application in foxtail millet breeding programs by marker-assisted selection.

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.

Generation and Characterization of a Foxtail Millet (Setaria italica) Mutant Library

Foxtail millet (Setaria italica) is attractive to plant scientists as a model plant because of several distinct characteristics, such as its short stature, rapid life cycle, sufficient seed production per plant, self-compatibility, true diploid nature, high photosynthetic efficiency, small genome size, and tolerance to abiotic and biotic stress. However, the study on the genetic resources of foxtail millet largely lag behind those of the other model plants such as Arabidopsis, rice and maize. Mutagenized populations cannot only create new germplasm resources, but also provide materials for gene function research. In this manuscript, an ethyl methanesulfonate (EMS)-induced foxtail millet population comprising ∼15,000 individual M₁ lines was established. Total 1353 independent lines with diverse abnormal phenotypes of leaf color, plant morphologies and panicle shapes were identified in M₂. Resequencing of sixteen randomly selected M₂ plants showed an average estimated mutation density of 1 loci/213 kb. Moreover, we provided an example for rapid cloning of the WP1 gene by a map-based cloning method. A white panicle mutant, named as wp1.a, exhibited significantly reduced chlorophyll (Chl) and carotenoid contents in leaf and panicle. Map-based cloning results showed an eight-base pair deletion located at the sixth exon of wp1.a in LOC101786849, which caused the premature termination. WP1 encoded phytoene synthase. Moreover, the sequencing analysis and cross test verified that a white panicle mutant wp1.b was an allelic mutant of wp1.a. The filed phenotypic observation and gene cloning example showed that our foxtail millet EMS-induced mutant population would be used as an important resource for functional genomics studies of foxtail millet.

QTL Mapping Combined With Comparative Analyses Identified Candidate Genes for Reduced Shattering in Setaria italica

Setaria (L.) P. Beauv is a genus of grasses that belongs to the Poaceae (grass) family, subfamily Panicoideae. Two members of the Setaria genus, Setaria italica (foxtail millet) and S. viridis (green foxtail), have been studied extensively over the past few years as model species for C4-photosynthesis and to facilitate genome studies in complex Panicoid bioenergy grasses. We exploited the available genetic and genomic resources for S. italica and its wild progenitor, S. viridis, to study the genetic basis of seed shattering. Reduced shattering is a key trait that underwent positive selection during domestication. Phenotyping of F2:3 and recombinant inbred line (RIL) populations generated from a cross between S. italica accession B100 and S. viridis accession A10 identified the presence of additive main effect quantitative trait loci (QTL) on chromosomes V and IX. As expected, enhanced seed shattering was contributed by the wild S. viridis. Comparative analyses pinpointed Sh1 and qSH1, two shattering genes previously identified in sorghum and rice, as potentially underlying the QTL on Setaria chromosomes IX and V, respectively. The Sh1 allele in S. italica was shown to carry a PIF/Harbinger MITE in exon 2, which gave rise to an alternatively spliced transcript that lacked exon 2. This MITE was universally present in S. italica accessions around the world and absent from the S. viridis germplasm tested, strongly suggesting a single origin of foxtail millet domestication. The qSH1 gene carried two MITEs in the 5'UTR. Presence of one or both MITEs was strongly associated with cultivated germplasm. If the MITE insertion(s) in qSH1 played a role in reducing shattering in S. italica accessions, selection for the variants likely occurred after the domestication of foxtail millet.

Domestication and Improvement in the Model C4 Grass, Setaria

Setaria viridis (green foxtail) and its domesticated relative S. italica (foxtail millet) are diploid C4 panicoid grasses that are being developed as model systems for studying grass genomics, genetics, development, and evolution. According to archeological evidence, foxtail millet was domesticated from green foxtail approximately 9,000 to 6,000 YBP in China. Under long-term human selection, domesticated foxtail millet developed many traits adapted to human cultivation and agricultural production. In comparison with its wild ancestor, foxtail millet has fewer vegetative branches, reduced grain shattering, delayed flowering time and less photoperiod sensitivity. Foxtail millet is the only present-day crop in the genus Setaria, although archeological records suggest that other species were domesticated and later abandoned in the last 10,000 years. We present an overview of domestication in foxtail millet, by reviewing recent studies on the genetic regulation of several domesticated traits in foxtail millet and discuss how the foxtail millet and green foxtail system could be further developed to both better understand its domestication history, and to provide more tools for future breeding efforts.

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

Heading date is an important event to ensure successful seed production. Although foxtail millet (Setaria italica (L.) P.Beauv.) is an important foodstuff in semiarid regions around the world, the genetic basis determining heading date is unclear. To identify genomic regions regulating days to heading (DTH), we conducted a QTL-seq analysis based on combining whole-genome re-sequencing and bulked-segregant analysis of an F₂ population derived from crosses between the middle-heading cultivar Shinanotsubuhime and the early-heading cultivar Yuikogane. Under field conditions, transgressive segregation of DTH toward late heading was observed in the F₂ population. We made three types of bulk samples: Y-bulk (early-heading), S-bulk (late-heading) and L-bulk (extremely late-heading). By genome-wide comparison of SNPs in the Y-bulk vs. the S-bulk and the Y-bulk vs. the L-bulk, we identified two QTLs associated with DTH. The first QTL, qDTH2, was detected on chromosome 2 from the Y-bulk and S-bulk comparison. The second QTL, qDTH7, was detected on chromosome 7 from the Y-bulk and L-bulk comparison. The Shinanotsubuhime allele for qDTH2 caused late heading in the F₂ population, whereas the Yuikogane allele for qDTH7 led to extremely late heading. These results suggest that allelic differences in both qDTH2 and qDTH7 determine regional adaptability in S. italica.

Detection of genomic loci associated with chromosomal recombination using high-density linkage mapping in Setaria

Meiotic recombination is essential to sexual reproduction and the generation of genetic diversity. Variation in recombination rates is presently of particular interest due to efforts being made to increase the rate of genetic gain in agricultural crops by breaking up large linkage disequilibrium blocks containing both beneficial and detrimental alleles. Here, a high-density genetic linkage map of Setaria was constructed using tunable genotyping by sequencing (tGBS) analysis of a population of recombinant inbred lines (RILs). Several regions of the Setaria genome exhibited significant levels of segregation distortion (SD), and recombination crossovers (COs) were also detected. The regions with high SD generally tended to have fewer COs, particularly for pericentromeric chromosomal areas. Recombination crossovers detected in Setaria were unevenly distributed across the genome and occurred more often in intergenic regions. Quantitative trait loci (QTLs) contributing towards the recombination frequency (Type I) and occurrence of COs in designated loci (Type II) were identified, and Type II QTLs garnered higher statistical power. The result of this study suggest that QTLs analysis of Type II traits using RILs might provide an opportunity to further understand meiotic recombination using high throughput genome sequencing and genotyping technologies.

A high-density genetic map and QTL analysis of agronomic traits in foxtail millet [Setaria italica (L.) P. Beauv.] using RAD-seq

Foxtail millet (Setaria italica), a very important grain crop in China, has become a new model plant for cereal crops and biofuel grasses. Although its reference genome sequence was released recently, quantitative trait loci (QTLs) controlling complex agronomic traits remains limited. The development of massively parallel genotyping methods and next-generation sequencing technologies provides an excellent opportunity for developing single-nucleotide polymorphisms (SNPs) for linkage map construction and QTL analysis of complex quantitative traits. In this study, a high-throughput and cost-effective RAD-seq approach was employed to generate a high-density genetic map for foxtail millet. A total of 2,668,587 SNP loci were detected according to the reference genome sequence; meanwhile, 9,968 SNP markers were used to genotype 124 F₂ progenies derived from the cross between Hongmiaozhangu and Changnong35; a high-density genetic map spanning 1648.8 cM, with an average distance of 0.17 cM between adjacent markers was constructed; 11 major QTLs for eight agronomic traits were identified; five co-dominant DNA markers were developed. These findings will be of value for the identification of candidate genes and marker-assisted selection in foxtail millet.

Utilization of UPLC/Q-TOF-MS-Based Metabolomics and AFLP-Based Marker-Assisted Selection to Facilitate/Assist Conventional Breeding of Polygala tenuifolia

As one of the most important traditional Chinese medicine, the quality of Polygala tenuifolia is difficult to control and a new method must be established to facilitate/assist the breeding of P. tenuifolia. In this study, UPLC/Q-TOF-MS-based metabolomics analysis was performed to determine the chemical composition and screen metabolite biomarkers according to agronomic traits. A total of 29 compounds and 18 metabolite biomarkers were found. AFLP-based marker-assisted selection (MAS) was used to identify molecular marker bands and screen characteristic bands associated with specific agronomic traits. 184 bands and 76 characteristic AFLP bands were found. The correlation network between compounds and characteristic AFLP bands was built, so we may directly breed certain P. tenuifolia herbs with special agronomic traits (or characteristic AFLP bands), which exhibit specific pharmacological functions depending on the content of the active compounds. The proposed method of metabolomics coupled with MAS could facilitate/assist the breeding of P. tenuifolia.

Identification of QTLs for 14 Agronomically Important Traits in Setaria italica Based on SNPs Generated from High-Throughput Sequencing

Foxtail millet (Setaria italica) is an important crop possessing C4 photosynthesis capability. The S. italica genome was de novo sequenced in 2012, but the sequence lacked high-density genetic maps with agronomic and yield trait linkages. In the present study, we resequenced a foxtail millet population of 439 recombinant inbred lines (RILs) and developed high-resolution bin map and high-density SNP markers, which could provide an effective approach for gene identification. A total of 59 QTL for 14 agronomic traits in plants grown under long- and short-day photoperiods were identified. The phenotypic variation explained ranged from 4.9 to 43.94%. In addition, we suggested that there may be segregation distortion on chromosome 6 that is significantly distorted toward Zhang gu. The newly identified QTL will provide a platform for sequence-based research on the S. italica genome, and for molecular marker-assisted breeding.

Setaria viridis as a Model System to Advance Millet Genetics and Genomics

Millet is a common name for a group of polyphyletic, small-seeded cereal crops that include pearl, finger and foxtail millet. Millet species are an important source of calories for many societies, often in developing countries. Compared to major cereal crops such as rice and maize, millets are generally better adapted to dry and hot environments. Despite their food security value, the genetic architecture of agronomically important traits in millets, including both morphological traits and climate resilience remains poorly studied. These complex traits have been challenging to dissect in large part because of the lack of sufficient genetic tools and resources. In this article, we review the phylogenetic relationship among various millet species and discuss the value of a genetic model system for millet research. We propose that a broader adoption of green foxtail (Setaria viridis) as a model system for millets could greatly accelerate the pace of gene discovery in the millets, and summarize available and emerging resources in S. viridis and its domesticated relative S. italica. These resources have value in forward genetics, reverse genetics and high throughput phenotyping. We describe methods and strategies to best utilize these resources to facilitate the genetic dissection of complex traits. We envision that coupling cutting-edge technologies and the use of S. viridis for gene discovery will accelerate genetic research in millets in general. This will enable strategies and provide opportunities to increase productivity, especially in the semi-arid tropics of Asia and Africa where millets are staple food crops.