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

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.

Pipeline to Design Inbred Lines and F1 Hybrids of Leaf Chicory (Radicchio) Using Male Sterility and Genotyping-by-Sequencing

Chicory, a horticultural crop cultivated worldwide, presents many botanical varieties and local biotypes. Among these, cultivars of the Italian radicchio group of the pure species Cichorium intybus L. and its interspecific hybrids with Cichorium endivia L.-as the "Red of Chioggia" biotype-includes several phenotypes. This study uses a pipeline to address the marker-assisted breeding of F1 hybrids: it presents the genotyping-by-sequencing results of four elite inbred lines using a RADseq approach and an original molecular assay based on CAPS markers for screening mutants with nuclear male sterility in the radicchio of Chioggia. A total of 2953 SNP-carrying RADtags were identified and used to compute the actual estimates of homozygosity and overall genetic similarity and uniformity of the populations, as well as to determine their genetic distinctiveness and differentiation. Molecular data were further used to investigate the genomic distribution of the RADtags among the two Cichorium species, allowing their mapping in 1131 and 1071 coding sequences in chicory and endive, respectively. Paralleling this, an assay to screen the genotype at the male sterility locus Cims-1 was developed to discriminate wild-type and mutant alleles of the causative gene myb80-like. Moreover, a RADtag mapped close to this genomic region proved the potential application of this method for future marker-assisted selection tools. Finally, after combining the genotype information of the core collection, the best 10 individuals from each inbred line were selected to compute the observed genetic similarity as a measure of uniformity as well as the expected homozygosity and heterozygosity estimates scorable by the putative progenies derived from selfing (pollen parent) and full-sibling (seed parent) or pair-wise crossing (F1 hybrids). This predictive approach was conducted as a pilot study to understand the potential application of RADseq in the fine tuning of molecular marker-assisted breeding strategies aimed at the development of inbred lines and F1 hybrids in leaf chicory.

Soluble Protein Content, Bioactive Compounds and the Antioxidant Activity in Seeds of Ten Rheum tanguticum Lines from Qinghai-Tibet Plateau

Rheum tanguticum (Rh. tanguticum) is a Chinese medicinal plant traditionally used in the treatment of constipation. As a byproduct, the seeds of this plant are rich in nutrients and phytochemicals. This study aimed to determine and assess seed germination ability, seed physical characteristics, soluble protein content, chemical constituents and antioxidant capacity from different breeding lines, to promote the development and utilization of seed resources. Significant differences were observed for the soluble protein content and antioxidant assays among the ten lines. The contents of aloe-emodin, rhein and catechins accumulated in seeds were extremely low and significantly different from those in roots. In contrast, emodin and chrysophanol were abundant in seeds, and significant differences were observed between seeds and roots. It was found that associations between gallic acid and catechins were not significant for either soluble protein or antioxidant capacity. There was a significantly positive correlation between the contents of four anthraquinones (aloe-emodin, rhein, emodin and chrysophanol) and soluble protein. Seeds have potent antioxidative capacity and relatively high levels of soluble protein content. The rich chemical composition of seeds can be widely used in the medical industry for further development.

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.

Conuping BSA-Seq and RNA-Seq Reveal the Molecular Pathway and Genes Associated with the Plant Height of Foxtail Millet (Setaria italica)

Foxtail millet (Setaria italica) plays an important role in C4 crop research and agricultural development in arid areas due to its short growth period, drought tolerance, and barren tolerance. Exploration of the dwarfing mechanism and the dwarf genes of foxtail millet can provide a reference for dwarf breeding and dwarf research of other C4 crops. In this study, genetic analysis was performed using phenotypic data, candidate genes were screened by bulk segregant analysis sequencing (BSA-Seq); differentially expressed genes and metabolic pathways in different strains of high samples were analyzed by RNA sequencing (RNA-Seq). The association analysis of BSA-Seq and RNA-Seq further narrowed the candidate range. As a result, a total of three quantitative trait loci (QTLs) and nine candidate genes related to plant height were obtained on chromosomes I and IX. Based on the functional prediction of the candidate genes, we propose a hypothetical mechanism for the formation of millet dwarfing, in which, metabolism and MAPK signaling play important roles in the formation of foxtail millet plant height.

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.

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.

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.

Genotyping Analysis by RAD-Seq Reads Is Useful to Assess the Genetic Identity and Relationships of Breeding Lines in Lavender Species Aimed at Managing Plant Variety Protection

Lavender species are widely distributed in their wild forms around the Mediterranean Basin and they are also cultivated worldwide as improved and registered clonal varieties. The economic interest of the species belonging to the Lavandula genus is determined by their use as ornamental plants and important source of essential oils that are destinated to the production of cosmetics, pharmaceuticals and foodstuffs. Because of the increasing number of cases of illegal commercialization of selected varieties, the protection of plant breeders’ rights has become of main relevance for the recognition of breeding companies’ royalties. With this aim, genomic tools based on molecular markers have been demonstrated to be very reliable and transferable among laboratories, and also much more informative than morphological descriptors. With the rising of the next-generation sequencing (NGS) technologies, several genotyping-by-sequencing approaches are now available. This study deals with a deep characterization of 15 varietal clones, belonging to two distinct Lavandula species, by means of restriction-site associated DNA sequencing (RAD-Seq). We demonstrated that this technology screens single nucleotide variants that enable to assess the genetic identity of individual accessions, to reconstruct genetic relationships among related breeding lines, to group them into genetically distinguishable main subclusters, and to assign their molecular lineages to distinct ancestors. Moreover, a number of polymorphic sites were identified within genes putatively involved in biosynthetic pathways related to both tissue pigmentation and terpene production, useful for breeding and/or protecting newly registered varieties. Overall, the results highlighted the presence of pure ancestries and interspecific hybrids for the analyzed Lavandula species, and demonstrated that RAD-Seq analysis is very informative and highly reliable for characterizing Lavandula clones and managing plant variety protection.

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.

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.

Genetic Mapping by Sequencing More Precisely Detects Loci Responsible for Anaerobic Germination Tolerance in Rice

Direct seeded rice (DSR) is a mainstay for planting rice in the Americas, and it is rapidly becoming more popular in Asia. It is essential to develop rice varieties that are suitable for this type of production system. ASD1, a landrace from India, possesses several traits desirable for direct-seeded fields, including tolerance to anaerobic germination (AG). To map the genetic basis of its tolerance, we examined a population of 200 F_2:3 families derived from a cross between IR64 and ASD1 using the restriction site-associated DNA sequencing (RAD-seq) technology. This genotyping platform enabled the identification of 1921 single nucleotide polymorphism (SNP) markers to construct a high-resolution genetic linkage map with an average interval of 0.9 cM. Two significant quantitative trait loci (QTLs) were detected on chromosomes 7 and 9, qAG7 and qAG9, with LOD scores of 7.1 and 15.0 and R² values of 15.1 and 29.4, respectively. Here, we obtained more precise locations of the QTLs than traditional simple sequence repeat and low-density SNP genotyping methods and may help further dissect the genetic factors of these QTLs.

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.

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 assessment of population structure and genetic diversity of Chinese Lou onion using specific length amplified fragment (SLAF) sequencing

Lou onion (Allium fistulosum L. var. viviparum) is an abundant source of flavonols which provides additional health benefits to diseases. Genome-wide specific length amplified fragment (SLAF) sequencing method is a rapidly developed deep sequencing technologies used for selection and identification of genetic loci or markers. This study aimed to elucidate the genetic diversity of 122 onion accessions in China using the SLAF-seq method. A set of 122 onion accessions including 107 A.fistulosum L. var. viviparum Makino, 3 A.fistulosum L. var. gigantum Makino, 3 A.mongolicum Regel and 9 A.cepa L. accessions (3 whites, 3 reds and 3 yellows) from different regions in China were enrolled. Genomic DNA was isolated from young leaves and prepared for the SLAF-seq, which generated a total of 1,387.55 M reads and 162,321 high quality SNPs (integrity >0.5 and MAF >0.05). These SNPs were used for the construction of neighbor-joining phylogenetic tree, in which 10 A.fistulosum L. var. viviparum Makino accessions from Yinchuan (Ningxia province) and Datong (Qinghai province) had close genetic relationship. The 3 A.cepa L. clusters (red, white and yellow) had close genetic relationship especially with the 97 A.fistulosum L. var. viviparum Makino accessions. Population structure analysis suggested entire population could be clustered into 3 groups, while principal component analysis (PCA) showed there were 4 genetic groups. We confirmed the SLAF-seq approach was effective in genetic diversity analysis in red onion accessions. The key findings would provide a reference to the Lou onion germplasm in China.

Genotyping by RAD Sequencing Analysis Assessed the Genetic Distinctiveness of Experimental Lines and Narrowed Down the Genomic Region Responsible for Leaf Shape in Endive (Cichorium endivia L.)

The characterization of genetic diversity in elite breeding stocks is crucial for the registration and protection of new varieties. Moreover, experimental population structure analysis and information about the genetic distinctiveness of commercial materials are essential for crop breeding programs. The purpose of our research was to assess the genetic relationships of 32 endive (Cichorium endivia L.) breeding lines, 18 from var. latifolium (escarole) and 14 from var. crispum (curly), using heterologous Cichorium intybus-derived simple sequence repeats (SSR) markers and single-nucleotide polymorphisms (SNP) markers. We found that 14 out of 29 SSR markers were successfully amplified, but only 8 of them were related to polymorphic loci. To overcome the limitation of the low number of informative SSR marker loci, an alternative SNP-based approach was employed. The 4621 SNPs produced by a restriction site-associated DNA marker sequencing approach were able to fully discriminate the 32 endive accessions; most importantly, as many as 50 marker loci were found to distinguish the curly group from the escarole group. Interestingly, 24 of the marker loci mapped within a peripheral segment of chromosome 8 of lettuce (Lactuca sativa L.), spanning a chromosomal region of 49.6 Mb. Following Sanger sequencing-based validation, three genes were determined to carry nonsynonymous SNPs, and one of them matched a putative ortholog of AtELP1, subunit 1 of the Elongator complex. Considering that several previously characterized Elongator complex subunit mutants exhibited elongated and/or curly leaf phenotypes, this gene should be taken into consideration for a better understanding of the underlying mechanism controlling leaf shape in endive.

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.

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.

Artificially designed hybrids facilitate efficient generation of high-resolution linkage maps

When sequencing eukaryotic genomes, linkage maps are indispensable for building scaffolds to assemble and/or to validate chromosomes. However, current approaches to constructing linkage maps are limited by marker density and cost-effectiveness, especially for wild organisms. We have now devised a new strategy based on artificially generated hybrid organisms to acquire ultrahigh-density genomic markers at reduced cost and build highly accurate linkage maps. We have also developed the novel analysis pipeline Scaffold Extender with Low Depth Linkage Analysis (SELDLA) for data processing to generate linkage maps and draft genomes. Using SELDLA, linkage maps and improved genomes for two species of pufferfish, Takifugu rubripes and Takifugu stictonotus, were obtained simultaneously. The strategy is applicable to a wide range of sexually reproducing organisms, and could, therefore, accelerate the whole genome analysis of various organisms including fish, mollusks, amphibians, insects, plants, and even mammals.

Acid phosphatase gene GmHAD1 linked to low phosphorus tolerance in soybean, through fine mapping

KEY MESSAGE

Map-based cloning identified GmHAD1, a gene which encodes a HAD-like acid phosphatase, associated with soybean tolerance to low phosphorus stress. Phosphorus (P) deficiency in soils is a major limiting factor for crop growth worldwide. Plants may adapt to low phosphorus (LP) conditions via changes to root morphology, including the number, length, orientation, and branching of the principal root classes. To elucidate the genetic mechanisms for LP tolerance in soybean, quantitative trait loci (QTL) related to root morphology responses to LP were identified via hydroponic experiments. In total, we identified 14 major loci associated with these traits in a RIL population. The log-likelihood scores ranged from 2.81 to 7.43, explaining 4.23-13.98% of phenotypic variance. A major locus on chromosome 08, named qP8-2, was co-localized with an important P efficiency QTL (qPE8), containing phosphatase genes GmACP1 and GmACP2. Another major locus on chromosome 10 named qP10-2 explained 4.80-13.98% of the total phenotypic variance in root morphology. The qP10-2 contains GmHAD1, a gene which encodes an acid phosphatase. In the transgenic soybean hairy roots, GmHAD1 overexpression increased P efficiency by 8.4-16.5% relative to the control. Transgenic Arabidopsis plants had higher biomass than wild-type plants across both short- and long-term P reduction. These results suggest that GmHAD1, an acid phosphatase gene, improved the utilization of organic phosphate by soybean and Arabidopsis plants.

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.