A haplotype map of genomic variations and genome-wide association studies of agronomic traits in foxtail millet (Setaria italica)

Assessment and phenotypic identification of millet germplasm (Setaria italica L.) in Liaoning, China

Aims

This study evaluated millet germplasms in Liaoning Province to support the collection, preservation and innovation of millet germplasm resources.

Methods

The study was conducted from 2018 to 2020, involved the selection of 105 millet germplasm resources from the Germplasm Bank of the Liaoning Academy of Agricultural Sciences (LAAS), the observation and recording of 31 traits, and the application of multivariate analysis methods to assess phenotypic diversity.

Results

From the diversity analysis and correlation analysis, it was found that the tested traits had abundant diversity and complex correlations among them. Principal component analysis (PCA) comprehensively analyzed all quantitative traits and extracted seven principal components. Grey relational analysis (GRA) highlighted the varied contributions of different traits to yield. Through systematic cluster analysis (SCA), the resources were categorized into six groups at Euclidean distance of 17.09. K-mean cluster analysis determined the distribution interval and central value of each trait, then identified resources with desirable traits.

Conclusion

The results revealed resources that possess characteristics such as upthrow seedling leaves, more tillers and branches, larger and well-formed ears, and lodging resistance prefer to higher grain yield. It was also discovered that the subear internode length (SIL) could be an indicator for maturity selection. Four specific resources, namely, Dungu No. 1, Xiao-li-xiang, Basen Shengu, and Yuhuanggu No. 1, were identified for further breeding and practical applications.

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

Genome-Wide Association and RNA-Seq Analyses Reveal a Potential Candidate Gene Related to Oil Content in Soybean Seeds

Soybean is a crucial crop globally, serving as a significant source of unsaturated fatty acids and protein in the human diet. However, further enhancements are required for the related genes that regulate soybean oil synthesis. In this study, 155 soybean germplasms were cultivated under three different environmental conditions, followed by phenotypic identification and genome-wide association analysis using simplified sequencing data. Genome-wide association analysis was performed using SLAF-seq data. A total of 36 QTLs were significantly associated with oil content (-log10(p) > 3). Out of the 36 QTLs associated with oil content, 27 exhibited genetic overlap with previously reported QTLs related to oil traits. Further transcriptome sequencing was performed on extreme high-low oil soybean varieties. Combined with transcriptome expression data, 22 candidate genes were identified (|log2FC| ≥ 3). Further haplotype analysis of the potential candidate genes showed that three potential candidate genes had excellent haplotypes, including Glyma.03G186200, Glyma.09G099500, and Glyma.18G248900. The identified loci harboring beneficial alleles and candidate genes likely contribute significantly to the molecular network's underlying marker-assisted selection (MAS) and oil content.

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

KEY MESSAGE

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

Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study providing fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

Omics-driven utilization of wild relatives for empowering pre-breeding in pearl millet

MAIN CONCLUSION

Pearl millet wild relatives harbour novel alleles which could be utilized to broaden genetic base of cultivated species. Genomics-informed pre-breeding is needed to speed up introgression from wild to cultivated gene pool in pearl millet. Rising episodes of intense biotic and abiotic stresses challenge pearl millet production globally. Wild relatives provide a wide spectrum of novel alleles which could address challenges posed by climate change. Pre-breeding holds potential to introgress novel diversity in genetically narrow cultivated Pennisetum glaucum from diverse gene pool. Practical utilization of gene pool diversity remained elusive due to genetic intricacies. Harnessing promising traits from wild pennisetum is limited by lack of information on underlying candidate genes/QTLs. Next-Generation Omics provide vast scope to speed up pre-breeding in pearl millet. Genomic resources generated out of draft genome sequence and improved genome assemblies can be employed to utilize gene bank accessions effectively. The article highlights genetic richness in pearl millet and its utilization with a focus on harnessing next-generation Omics to empower pre-breeding.

Genomic variation in weedy and cultivated broomcorn millet accessions uncovers the genetic architecture of agronomic traits

Large-scale genomic variations are fundamental resources for crop genetics and breeding. Here we sequenced 1,904 genomes of broomcorn millet to an average of 40× sequencing depth and constructed a comprehensive variation map of weedy and cultivated accessions. Being one of the oldest cultivated crops, broomcorn millet has extremely low nucleotide diversity and remarkably rapid decay of linkage disequilibrium. Genome-wide association studies identified 186 loci for 12 agronomic traits. Many causative candidate genes, such as PmGW8 for grain size and PmLG1 for panicle shape, showed strong selection signatures during domestication. Weedy accessions contained many beneficial variations for the grain traits that are largely lost in cultivated accessions. Weedy and cultivated broomcorn millet have adopted different loci controlling flowering time for regional adaptation in parallel. Our study uncovers the unique population genomic features of broomcorn millet and provides an agronomically important resource for cereal crops.

Leaf Senescence Regulation Mechanism Based on Comparative Transcriptome Analysis in Foxtail Millet

Leaf senescence, a pivotal process in plants, directly influences both crop yield and nutritional quality. Foxtail millet (Setaria italica) is a C4 model crop renowned for its exceptional nutritional value and stress tolerance characteristics. However, there is a lack of research on the identification of senescence-associated genes (SAGs) and the underlying molecular regulatory mechanisms governing this process. In this study, a dark-induced senescence (DIS) experimental system was applied to investigate the extensive physiological and transcriptomic changes in two foxtail millet varieties with different degrees of leaf senescence. The physiological and biochemical indices revealed that the light senescence (LS) variety exhibited a delayed senescence phenotype, whereas the severe senescence (SS) variety exhibited an accelerated senescence phenotype. The most evident differences in gene expression profiles between these two varieties during DIS included photosynthesis, chlorophyll, and lipid metabolism. Comparative transcriptome analysis further revealed a significant up-regulation of genes related to polysaccharide and calcium ion binding, nitrogen utilization, defense response, and malate metabolism in LS. In contrast, the expression of genes associated with redox homeostasis, carbohydrate metabolism, lipid homeostasis, and hormone signaling was significantly altered in SS. Through WGCNA and RT-qPCR analyses, we identified three SAGs that exhibit potential negative regulation towards dark-induced leaf senescence in foxtail millet. This study establishes the foundation for a further comprehensive examination of the regulatory network governing leaf senescence and provides potential genetic resources for manipulating senescence in foxtail millet.

Genetic diversity of grain yield traits and identification of a grain weight gene SiTGW6 in foxtail millet

KEY MESSAGE

Agronomic traits were evaluated in 1250 foxtail millet accessions, and a crucial gene SiTGW6 governing grain yield was identified. Elite haplotypes and dCAPS markers developed for SiTGW6 facilitate molecular breeding. A comprehensive evaluation of phenotypic characteristics and genetic diversity in germplasm resources are important for gene discovery and breeding improvements. In this study, we conducted a comprehensive evaluation of 1250 foxtail millet varieties, assessing seven grain yield-related traits and fourteen common agronomic traits over two years. Principal component analysis, correlation analysis, and cluster analysis revealed a strong positive correlation between 1000-grain weight and grain width with grain yield, emphasizing their importance in foxtail millet breeding. Additionally, we found that panicle weight positively correlated with 1000-grain weight but negatively correlated with branch and tiller numbers, indicating selection factors during domestication and breeding. Using this information, we identified 27 germplasm resources suitable for high-yield foxtail millet breeding. Furthermore, through an integration of haplotype variations and phenotype association analysis, we pinpointed a crucial gene, SiTGW6, responsible for governing grain yield in foxtail millet. SiTGW6 encodes an IAA-glucose hydrolase, primarily localized in the cytoplasm and predominantly expressed in flowering panicles. Employing RNAseq analysis, we identified 1439 differentially expressed genes across various SiTGW6 haplotypes. Functional enrichment analysis indicating that SiTGW6 regulates grain yield through the orchestration of auxin and glucan metabolism, as well as plant hormone signaling pathways. Additionally, we have identified elite haplotypes and developed dCAPS markers for SiTGW6, providing valuable technical tools to facilitate molecular breeding efforts in foxtail millet.

A genomic variation map provides insights into peanut diversity in China and associations with 28 agronomic traits

Peanut (Arachis hypogaea L.) is an important allotetraploid oil and food legume crop. China is one of the world's largest peanut producers and consumers. However, genomic variations underlying the migration and divergence of peanuts in China remain unclear. Here we reported a genome-wide variation map based on the resequencing of 390 peanut accessions, suggesting that peanuts might have been introduced into southern and northern China separately, forming two cultivation centers. Selective sweep analysis highlights asymmetric selection between the two subgenomes during peanut improvement. A classical pedigree from South China offers a context for the examination of the impact of artificial selection on peanut genome. Genome-wide association studies identified 22,309 significant associations with 28 agronomic traits, including candidate genes for plant architecture and oil biosynthesis. Our findings shed light on peanut migration and diversity in China and provide valuable genomic resources for peanut improvement.

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.

Crop Evolution of Foxtail Millet

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

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

KEY MESSAGE

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

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.

The m6 A reader SiYTH1 enhances drought tolerance by affecting the messenger RNA stability of genes related to stomatal closure and reactive oxygen species scavenging in Setaria italica

Foxtail millet (Setaria italica), a vital drought-resistant crop, plays a significant role in ensuring food and nutritional security. However, its drought resistance mechanism is not fully understood. N6 -methyladenosine (m6 A) modification of RNA, a prevalent epi-transcriptomic modification in eukaryotes, provides a binding site for m6 A readers and affects plant growth and stress responses by regulating RNA metabolism. In this study, we unveiled that the YT521-B homology (YTH) family gene SiYTH1 positively regulated the drought tolerance of foxtail millet. Notably, the siyth1 mutant exhibited reduced stomatal closure and augmented accumulation of excessive H2 O2 under drought stress. Further investigations demonstrated that SiYTH1 positively regulated the transcripts harboring m6 A modification related to stomatal closure and reactive oxygen species (ROS) scavenging under drought stress. SiYTH1 was uniformly distributed in the cytoplasm of SiYTH1-GFP transgenic foxtail millet. It formed dynamic liquid-like SiYTH1 cytosol condensates in response to drought stress. Moreover, the cytoplasmic protein SiYTH1 was identified as a distinct m6 A reader, facilitating the stabilization of its directly bound SiARDP and ROS scavenging-related transcripts under drought stress. Furthermore, natural variation analysis revealed SiYTH1AGTG as the dominant allele responsible for drought tolerance in foxtail millet. Collectively, this study provides novel insights into the intricate mechanism of m6 A reader-mediated drought tolerance and presents a valuable genetic resource for improving drought tolerance in foxtail millet breeding.

Identification of Multiple Genetic Loci Related to Low-Temperature Tolerance during Germination in Maize (Zea maize L.) through a Genome-Wide Association Study

Low-temperature stress during the germination stage is an important abiotic stress that affects the growth and development of northern spring maize and seriously restricts maize yield and quality. Although some quantitative trait locis (QTLs) related to low-temperature tolerance in maize have been detected, only a few can be commonly detected, and the QTL intervals are large, indicating that low-temperature tolerance is a complex trait that requires more in-depth research. In this study, 296 excellent inbred lines from domestic and foreign origins (America and Europe) were used as the study materials, and a low-coverage resequencing method was employed for genome sequencing. Five phenotypic traits related to low-temperature tolerance were used to assess the genetic diversity of maize through a genome-wide association study (GWAS). A total of 14 SNPs significantly associated with low-temperature tolerance were detected (-log10(P) > 4), and an SNP consistently linked to low-temperature tolerance in the field and indoors during germination was utilized as a marker. This SNP, 14,070, was located on chromosome 5 at position 2,205,723, which explained 4.84-9.68% of the phenotypic variation. The aim of this study was to enrich the genetic theory of low-temperature tolerance in maize and provide support for the innovation of low-temperature tolerance resources and the breeding of new varieties.

Genome-Wide Characterization and Haplotypic Variation Analysis of the YUC Gene Family in Foxtail Millet (Setaria italica)

Panicle development and grain production in crop species are essential breeding characteristics affected by the synthesis of auxin, which is influenced by flavin monooxygenase-encoding genes such as YUC (YUCCA) family members. In this trial, fourteen YUCs were identified and named uniformly in foxtail millet, an ancient crop species cultivated across the world. The phylogenetic analysis revealed that the SiYUCs were clustered into four subgroups; protein motif and gene structure analyses suggested that the closely clustered SiYUC genes were relatively conserved within each subgroup; while genome mapping analysis indicated that the SiYUC genes were unevenly distributed on foxtail millet chromosomes and colinear with other grass species. Transcription analysis revealed that the SiYUC genes differed greatly in expression pattern in different tissues and contained hormonal/light/stress-responding cis-elements. The haplotype characterization of SiYUC genes indicated many superior haplotypes of SiYUCs correlated with higher panicle and grain weight could be favorably selected by breeding. These results will be useful for the further study of the functional characteristics of SiYUC genes, particularly with regard to the marker-assisted pyramiding of beneficial haplotypes in foxtail millet breeding programs.

Association mapping in bambara groundnut [Vigna subterranea (L.) Verdc.] reveals loci associated with agro-morphological traits

BACKGROUND

Genome-wide association studies (GWAS) are important for the acceleration of crop improvement through knowledge of marker-trait association (MTA). This report used DArT SNP markers to successfully perform GWAS on agro-morphological traits using 270 bambara groundnut [Vigna subterranea (L.) Verdc.] landraces sourced from diverse origins. The study aimed to identify marker traits association for nine agronomic traits using GWAS and their candidate genes. The experiment was conducted at two different locations laid out in alpha lattice design. The cowpea [Vigna unguiculata (L.) Walp.] reference genome (i.e. legume genome most closely related to bambara groundnut) assisted in the identification of candidate genes.

RESULTS

The analyses showed that linkage disequilibrium was found to decay rapidly with an average genetic distance of 148 kb. The broadsense heritability was relatively high and ranged from 48.39% (terminal leaf length) to 79.39% (number of pods per plant). The GWAS identified a total of 27 significant marker-trait associations (MTAs) for the nine studied traits explaining 5.27% to 24.86% of phenotypic variations. Among studied traits, the highest number of MTAs was obtained from seed coat colour (6) followed by days to flowering (5), while the least is days to maturity (1), explaining 5.76% to 11.03%, 14.5% to 19.49%, and 11.66% phenotypic variations, respectively. Also, a total of 17 candidate genes were identified, varying in number for different traits; seed coat colour (6), days to flowering (3), terminal leaf length (2), terminal leaf width (2), number of seed per pod (2), pod width (1) and days to maturity (1).

CONCLUSION

These results revealed the prospect of GWAS in identification of SNP variations associated with agronomic traits in bambara groundnut. Also, its present new opportunity to explore GWAS and marker assisted strategies in breeding of bambara groundnut for acceleration of the crop improvement.

The osmotic stress-activated receptor-like kinase DPY1 mediates SnRK2 kinase activation and drought tolerance in Setaria

Plant genomes encode many receptor-like kinases (RLKs) that localize to the cell surface and perceive a wide variety of environmental cues to initiate downstream signaling cascades. Whether these RLKs participate in dehydration stress signaling in plants is largely unknown. DROOPY LEAF1 (DPY1), a leucine-rich repeat (LRR)-RLK, was recently shown to regulate plant architecture by orchestrating early brassinosteroid signaling in foxtail millet (Setaria italica). Here, we show that DPY1 is essential for the acclimation of foxtail millet to drought stress. DPY1 can be phosphorylated and activated in response to osmotic stress and is required for more than half of osmotic stress-induced global phosphorylation events, including the phosphorylation of sucrose nonfermenting kinase 2s (SnRK2s), the central kinases involved in osmotic stress. DPY1 acts upstream of STRESS-ACTIVATED PROTEIN KINASE 6 (SAPK6, a subclass I SnRK2) and is required for full SAPK6 activation, thereby allowing regulation of downstream genes to mount a response against drought stress. These signaling events are largely independent of DPY1-mediated brassinosteroid signaling. The DPY1-SAPK6 module is specific to seed plants and is absent in ancestral nonseed plants. Our findings reveal a dehydration stress-activated RLK that plays an indispensable role in osmotic stress signaling and mediates SnRK2 activation at the cell surface.

Association analysis of agronomic traits and construction of genetic networks by resequencing of 306 sugar beet (Beta vulgaris L.) lines

Due to the relatively brief domestication history of sugar beet (Beta vulgaris ssp. vulgaris), our understanding of the genomic diversity and functional genes in its cultivars is limited, resulting in slow breeding progress. To address this issue, a total of 306 germplasm materials of major cultivars and breeding lines from China, the USA, and Europe were selected for genome resequencing. We investigated population structure and genetic diversity and performed selective scanning of genomic regions, identifying six novel genes associated with important agronomic traits: the candidate genes DFAX2 and P5CS for skin roughness; the candidate genes FRO5, GL24, and PPR91 for root yield and sugar yield, and the pleiotropic candidate gene POLX for flourishing growth vigour, plant height, crown size, flesh coarseness, and sugar yield. In addition, we constructed a protein-protein interaction network map and a phenotype-gene network map, which provide valuable information for identifying and characterizing functional genes affecting agronomic traits in sugar beet. Overall, our study sheds light on the future improvement of sugar beet agronomic traits at the molecular level.

Mapping of major QTL and candidate gene analysis for hull colour in foxtail millet (Setaria italica (L.) P. Beauv.)

BACKGROUND

Hull colour is an important morphological marker for selection in seed production of foxtail millet. However, the molecular mechanisms underlying hull colour variation remain unknown.

RESULTS

An F7 recombinant inbred line (RIL) population containing 215 lines derived from Hongjiugu × Yugu18 was used to analyze inheritance and detect the quantitative trait loci (QTL) for four hull colour traits using major gene plus polygene mixed inheritance analysis and composite interval mapping (CIM) in four environments. Genetic analysis revealed that the hull colour L* value (HCL*) was controlled by two major genes plus additive polygenes, the hull colour a* value (HCa*) was controlled by three major genes, the hull colour b* value (HCb*) was controlled by two major genes plus polygenes, and the hull colour C* value (HCC*) was controlled by four major genes. A high-density genetic linkage map covering 1227.383 cM of the foxtail millet genome, with an average interval of 0.879 cM between adjacent bin markers, was constructed using 1420 bin markers. Based on the genetic linkage map and the phenotypic data, a total of 39 QTL were detected for these four hull colour traits across four environments, each explaining 1.50%-49.20% of the phenotypic variation. Of these, six environmentally stable major QTL were co-localized to regions on chromosomes 1 and 9, playing a major role in hull colour. There were 556 annotated genes within the two QTL regions. Based on the functions of homologous genes in Arabidopsis and the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) gene annotations, five genes were predicted as candidate genes for further studies.

CONCLUSIONS

This is the first study to use an inheritance model and QTL mapping to determine the genetic mechanisms of hull colour trait in foxtail millet. We identified six major environmentally stable QTL and predicted five potential candidate genes to be associated with hull colour. These results advance the current understanding of the genetic mechanisms underlying hull colour traits in foxtail millet and provide additional resources for application in genomics-assisted breeding and potential isolation and functional characterization of the candidate genes.

Role of SiPHR1 in the Response to Low Phosphate in Foxtail Millet via Comparative Transcriptomic and Co-Expression Network Analyses

Enhancing the absorption and utilization of phosphorus by crops is an important aim for ensuring food security worldwide. However, the gene regulatory network underlying phosphorus use in foxtail millet remains unclear. In this study, the molecular mechanism underlying low-phosphorus (LP) responsiveness in foxtail millet was evaluated using a comparative transcriptome analysis. LP reduced the chlorophyll content in shoots, increased the anthocyanin content in roots, and up-regulated purple acid phosphatase and phytase activities as well as antioxidant systems (CAT, POD, and SOD). Finally, 13 differentially expressed genes related to LP response were identified and verified using transcriptomic data and qRT-PCR. Two gene co-expression network modules related to phosphorus responsiveness were positively correlated with POD, CAT, and PAPs. Of these, SiPHR1, functionally annotated as PHOSPHATE STARVATION RESPONSE 1, was identified as an MYB transcription factor related to phosphate responsiveness. SiPHR1 overexpression in Arabidopsis significantly modified the root architecture. LP stress caused cellular, physiological, and phenotypic changes in seedlings. SiPHR1 functioned as a positive regulator by activating downstream genes related to LP tolerance. These results improve our understanding of the molecular mechanism underlying responsiveness to LP stress, thereby laying a theoretical foundation for the genetic modification and breeding of new LP-tolerant foxtail millet varieties.

Plant genome resequencing and population genomics: Current status and future prospects

Advances in DNA sequencing technology have sparked a genomics revolution, driving breakthroughs in plant genetics and crop breeding. Recently, the focus has shifted from cataloging genetic diversity in plants to exploring their functional significance and delivering beneficial alleles for crop improvement. This transformation has been facilitated by the increasing adoption of whole-genome resequencing. In this review, we summarize the current progress of population-based genome resequencing studies and how these studies affect crop breeding. A total of 187 land plants from 163 countries have been resequenced, comprising 54 413 accessions. As part of resequencing efforts 367 traits have been surveyed and 86 genome-wide association studies have been conducted. Economically important crops, particularly cereals, vegetables, and legumes, have dominated the resequencing efforts, leaving a gap in 49 orders, including Lycopodiales, Liliales, Acorales, Austrobaileyales, and Commelinales. The resequenced germplasm is distributed across diverse geographic locations, providing a global perspective on plant genomics. We highlight genes that have been selected during domestication, or associated with agronomic traits, and form a repository of candidate genes for future research and application. Despite the opportunities for cross-species comparative genomics, many population genomic datasets are not accessible, impeding secondary analyses. We call for a more open and collaborative approach to population genomics that promotes data sharing and encourages contribution-based credit policy. The number of plant genome resequencing studies will continue to rise with the decreasing DNA sequencing costs, coupled with advances in analysis and computational technologies. This expansion, in terms of both scale and quality, holds promise for deeper insights into plant trait genetics and breeding design.

Association mapping in multiple yam species (Dioscorea spp.) of quantitative trait loci for yield-related traits

BACKGROUND

Yam (Dioscorea spp.) is multiple species with various ploidy level and considered as cash crop in many producing areas. Selection based phenotyping for yield and its related traits such as mosaic virus and anthracnose diseases resistance and plant vigor in multiple species of yam is lengthy however, marker information has proven to enhance selection efficiency.

METHODOLOGY

In this study, a panel of 182 yam accessions distributed across six yam species were assessed for diversity and marker-traits association study using SNP markers generated from Diversity Array Technology platform. For the traits association analysis, the relation matrix alongside the population structure were used as co-factor to avoid false discovery using Multiple random Mixed Linear Model (MrMLM) followed by gene annotation.

RESULTS

Accessions performance were significantly different (p < 0.001) across all the traits with high broad-sense heritability (H²). Phenotypic and genotypic correlations showed positive relationships between yield and vigor but negative for yield and yam mosaic disease severity. Population structure revealed k = 6 as optimal clusters-based species. A total of 22 SNP markers were identified to be associated with yield, vigor, mosaic and anthracnose diseases resistance. Gene annotation for the significant SNP loci identified some putative genes associated with primary metabolism, pest and resistance to anthracnose disease, maintenance of NADPH in biosynthetic reaction especially those involving nitro-oxidative stress for resistance to mosaic virus, and seed development, photosynthesis, nutrition use efficiency, stress tolerance, vegetative and reproductive development for tuber yield.

CONCLUSION

This study provides valuable insights into the genetic control of plant vigor, anthracnose, mosaic virus resistance, and tuber yield in yam and thus, opens an avenue for developing additional genomic resources for markers-assisted selection focusing on multiple yam species.

The potentialities of omics resources for millet improvement

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

A graph-based genome and pan-genome variation of the model plant Setaria

Setaria italica (foxtail millet), a founder crop of East Asian agriculture, is a model plant for C4 photosynthesis and developing approaches to adaptive breeding across multiple climates. Here we established the Setaria pan-genome by assembling 110 representative genomes from a worldwide collection. The pan-genome is composed of 73,528 gene families, of which 23.8%, 42.9%, 29.4% and 3.9% are core, soft core, dispensable and private genes, respectively; 202,884 nonredundant structural variants were also detected. The characterization of pan-genomic variants suggests their importance during foxtail millet domestication and improvement, as exemplified by the identification of the yield gene SiGW3, where a 366-bp presence/absence promoter variant accompanies gene expression variation. We developed a graph-based genome and performed large-scale genetic studies for 68 traits across 13 environments, identifying potential genes for millet improvement at different geographic sites. These can be used in marker-assisted breeding, genomic selection and genome editing to accelerate crop improvement under different climatic conditions.

An E2-E3 pair contributes to seed size control in grain crops

Understanding the molecular mechanisms that regulate grain yield is important for improving agricultural productivity. Protein ubiquitination controls various aspects of plant growth but lacks understanding on how E2-E3 enzyme pairs impact grain yield in major crops. Here, we identified a RING-type E3 ligase SGD1 and its E2 partner SiUBC32 responsible for grain yield control in Setaria italica. The conserved role of SGD1 was observed in wheat, maize, and rice. Furthermore, SGD1 ubiquitinates the brassinosteroid receptor BRI1, stabilizing it and promoting plant growth. Overexpression of an elite SGD1 haplotype improved grain yield by about 12.8% per plant, and promote complex biological processes such as protein processing in endoplasmic reticulum, stress responses, photosystem stabilization, and nitrogen metabolism. Our research not only identifies the SiUBC32-SGD1-BRI1 genetic module that contributes to grain yield improvement but also provides a strategy for exploring key genes controlling important traits in Poaceae crops using the Setaria model system.

Insights into the genetic determination of tuber shape and eye depth in potato natural population based on autotetraploid potato genome

Potato is one of the world's most important food crops, with a time-consuming breeding process. In this study, we performed a genome-wide association (GWAS) analysis of the two important traits of potato tuber shape and eye depth, using the tetraploid potato genome (2n=4x=48) as a reference. A total of 370 potatoes were divided into three subgroups based on the principal component analysis and evolutionary tree analysis. The genetic diversity within subgroups is low (5.18×10^-5, 4.36×10^-5 and 4.24×10^-5). Genome-wide linkage disequilibrium (LD) analysis showed that their LD is about 60 Kb. GWAS analysis identified that 146 significant single nucleotide polymorphism (SNP) loci at Chr01A1:34.44-35.25 Mb and Chr02A1:28.35-28.54 Mb regions are significantly associated with potato tuber shape, and that three candidate genes that might be related to potato tuber traits, PLATZ transcription factor, UTP-glucose-1-phosphate uridylyltransferase and FAR1 DNA-binding domain, are in the association region of Chr02A1. GWAS analysis identified 53 significant SNP loci at Chr05A2: 49.644-50.146 Mb and Chr06A2: 25.866-26.384 Mb regions with robust associations with potato tuber eye depth. Hydrolase and methyltransferases are present in the association region of Chr05A2, and three CYPs are present in the association region of Chr06A2. Our findings suggested that these genes are closely associated with potato tuber shape and eye depth. Our study identified molecular markers and candidate genes for improving tetraploid potato tuber shape and eye depth and provided ideas and insights for tetraploid potato breeding.

How ambient temperature affects the heading date of foxtail millet (Setaria italica)

Foxtail millet (Setaria italica), a short-day plant, is one of the important crops for food security encountering climate change, particularly in regions where it is a staple food. Under the short-day condition in Taiwan, the heading dates (HDs) of foxtail millet accessions varied by genotypes and ambient temperature (AT). The allelic polymorphisms in flowering time (FT)-related genes were associated with HD variations. AT, in the range of 13°C-30°C that was based on field studies at three different latitudes in Taiwan and observations in the phytotron at four different AT regimes, was positively correlated with growth rate, and high AT promoted HD. To elucidate the molecular mechanism of foxtail millet HD, the expression of 14 key FT-related genes in four accessions at different ATs was assessed. We found that the expression levels of SiPRR95, SiPRR1, SiPRR59, SiGhd7-2, SiPHYB, and SiGhd7 were negatively correlated with AT, whereas the expression levels of SiEhd1, SiFT11, and SiCO4 were positively correlated with AT. Furthermore, the expression levels of SiGhd7-2, SiEhd1, SiFT, and SiFT11 were significantly associated with HD. A coexpression regulatory network was identified that shown genes involved in the circadian clock, light and temperature signaling, and regulation of flowering, but not those involved in photoperiod pathway, interacted and were influenced by AT. The results reveal how gene × temperature and gene × gene interactions affect the HD in foxtail millet and could serve as a foundation for breeding foxtail millet cultivars for shift production to increase yield in response to global warming.

Genome-wide identification and expression analysis of the SAUR gene family in foxtail millet (Setaria italica L.)

BACKGROUND

Auxin performs important functions in plant growth and development processes, as well as abiotic stress. Small auxin-up RNA (SAUR) is the largest gene family of auxin-responsive factors. However, the knowledge of the SAUR gene family in foxtail millet is largely obscure.

RESULTS

In the current study, 72 SiSAUR genes were identified and renamed according to their chromosomal distribution in the foxtail millet genome. These SiSAUR genes were unevenly distributed on nine chromosomes and were classified into three groups through phylogenetic tree analysis. Most of the SiSAUR members from the same group showed similar gene structure and motif composition characteristics. Analysis of cis-acting elements showed that many hormone and stress response elements were identified in the promoter region of SiSAURs. Gene replication analysis revealed that many SiSAUR genes were derived from gene duplication events. We also found that the expression of 10 SiSAURs was induced by abiotic stress and exogenous hormones, which indicated that SiSAUR genes may participated in complex physiological processes.

CONCLUSIONS

Overall, these results will be valuable for further studies on the biological role of SAUR genes in foxtail development and response to stress conditions and may shed light on the improvement of the genetic breeding of foxtail millet.

Domestication-associated PHYTOCHROME C is a flowering time repressor and a key factor determining Setaria as a short-day plant

Phytochromes play vital roles in the regulation of flowering time, but little is known in Panicoideae species, especially the C₄ model Setaria. Here, genomic variations of PHYTOCHROME C (PHYC) between wild and cultivated Setaria gene pools were analysed and three SiphyC mutants were identified. The function of SiPHYC was verified by CRISPR-Cas9 approach and transcriptome sequencing. Furthermore, efficiency of indoor cultivation of SiphyC mutants were systematically evaluated. An extreme purified selection of PHYC was detected in wild to cultivated domestication process of Setaria. SiphyC mutants and knockout transgenic plants showed an early heading date and a loss of response to short-day photoperiod. Furthermore, variable expression of SiFTa, SiMADS14 and SiMADS15 might be responsible for promoting flowering of SiphyC mutants. Moreover, SiphyC mutant was four times that of the indoor plot ratio of wild-type and produced over 200 seeds within 45 d per individual. Our results suggest that domestication-associated SiPHYC repressed flowering and determined Setaria as a short-day plant, and SiphyC mutants possess the potential for creating efficient indoor cultivation system suitable for research on Setaria as a model, and either for maize or sorghum as well.

Convergence in domesticated fungi used for cheese and dry-cured meat maturation: beneficial traits, genomic mechanisms, and degeneration

Humans have domesticated genetically distant fungi for similar uses, the fermentation of lipid-rich and sugar-rich food to generate attractive aspects, odor and aroma, and to improve shelf life and product safety. Multiple independent domestication events also occurred within species. We review recent evidence of phenotypic convergence during the domestication of fungi for making cheese (Saccharomyces cerevisiae, Penicillium roqueforti, P. camemberti, and Geotrichum candidum) and for dry-cured meat making (P. nalgiovense and P. salamii). Convergence following adaptation to similar ecological niches involved colony aspect (fluffiness and color), lipolysis, proteolysis, volatile compound production, and competitive ability against food spoilers. We review evidence for convergence in genetic diversity loss in domesticated populations and in the degeneration of unused traits, such as toxin production and sexual reproduction. Phenotypic convergence sometimes occurred by similar mechanisms of genomic adaptation, in particular horizontal gene transfers and loss of genes.

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.

Genome-wide association study for resistance in bread wheat (Triticum aestivum L.) to stripe rust (Puccinia striiformis f. sp. tritici) races in Argentina

BACKGROUND

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most devastating diseases of the wheat crop. It causes significant reductions in both grain yield and grain quality. In recent years, new and more virulent races have overcome many of the known resistance genes in Argentinian germplasm. In order to identify loci conferring resistance to the local races of Pst for effective utilization in future breeding programs, a genome-wide association study (GWAS) was performed using a collection of 245 bread wheat lines genotyped with 90 K SNPs.

RESULTS

To search for adult plant resistance (APR) the panel was evaluated for disease severity (DS) and area under disease progress curve (AUDPC) in field trials during two years under natural infection conditions. To look for seedling or all-stage resistance (ASR) the panel was evaluated to determine infection type (IT) under greenhouse conditions against two prevalent races in Argentina. The phenotypic data showed that the panel possessed enough genetic variability for searching for sources of resistance to Pst. Significant correlations between years were observed for Pst response in the field and high heritability values were found for DS (H² = 0.89) and AUDPC (H² = 0.93). Based on GWAS, eight markers associated with Pst resistance (FDR < 0.01) were identified, of these, five were associated with ASR (on chromosomes 1B, 2A, 3A and 5B) and three with APR (on chromosomes 3B and 7A). These markers explained between 2% and 32.62% of the phenotypic variation. Five of the markers corresponded with previously reported Yr genes/QTL, while the other three (QYr.Bce.1B.sd.1, QYr.Bce.3A.sd and QYr.Bce.3B.APR.2) might be novel resistance loci.

CONCLUSION

Our results revealed high genetic variation for resistance to Argentinian stripe rust races in the germplasm used here. It constitutes a very promising step towards the improvement of Pst resistance of bread wheat in Argentina. Also, the identification of new resistance loci would represent a substantial advance for diversifying the current set of resistance genes and to advance in the improvement of the durable resistance to the disease.

Modulating physiological and transcriptional regulatory mechanisms for enhanced climate resilience in cereal crops

Climate change adversely affects the yield and productivity of cereal crops, which consequently impacts food security. Therefore, studying stress acclimation, particularly transcriptional patterns and morpho-physiological responses of cereal crops to different stresses, will provide insights into the molecular determinants underlying climate resilience. The availability of advanced tools and approaches has enabled the characterization of plants at morphological, physiological, biochemical, and molecular levels, which will lead to the identification of genomic regions regulating the stress responses at these levels. This will further facilitate using transgenic, breeding, or genome editing approaches to manipulate the identified regions (genes, alleles, or QTLs) to enhance stress resilience. Next-generation sequencing approaches have advanced the identification of causal genes and markers in the genomes through forward or reverse genetics. In this context, the review enumerates the progress of dissecting the molecular mechanisms underlying transcriptional and physiological responses of major cereals to climate-induced stresses. The review systematically discusses different tools and approaches available to study the response of plants to various stresses and identify the molecular determinants regulating stress-resilience. Further, the application of genomics-assisted breeding, transgene-, and targeted editing-based approaches for modulating the genetic determinants for enhanced climate resilience has been elaborated.

Conventional and Omics Approaches for Understanding the Abiotic Stress Response in Cereal Crops-An Updated Overview

Cereals have evolved various tolerance mechanisms to cope with abiotic stress. Understanding the abiotic stress response mechanism of cereal crops at the molecular level offers a path to high-yielding and stress-tolerant cultivars to sustain food and nutritional security. In this regard, enormous progress has been made in the omics field in the areas of genomics, transcriptomics, and proteomics. Omics approaches generate a massive amount of data, and adequate advancements in computational tools have been achieved for effective analysis. The combination of integrated omics and bioinformatics approaches has been recognized as vital to generating insights into genome-wide stress-regulation mechanisms. In this review, we have described the self-driven drought, heat, and salt stress-responsive mechanisms that are highlighted by the integration of stress-manipulating components, including transcription factors, co-expressed genes, proteins, etc. This review also provides a comprehensive catalog of available online omics resources for cereal crops and their effective utilization. Thus, the details provided in the review will enable us to choose the appropriate tools and techniques to reduce the negative impacts and limit the failures in the intensive crop improvement study.

Genome-Wide Identification and Characterization of the CCT Gene Family in Foxtail Millet (Setaria italica) Response to Diurnal Rhythm and Abiotic Stress

The CCT gene family plays important roles in diurnal rhythm and abiotic stress response, affecting crop growth and development, and thus yield. However, little information is available on the CCT family in foxtail millet (Setaria italica). In the present study, we identified 37 putative SiCCT genes from the foxtail millet genome. A phylogenetic tree was constructed from the predicted full-length SiCCT amino acid sequences, together with CCT proteins from rice and Arabidopsis as representatives of monocotyledonous and dicotyledonous plants, respectively. Based on the conserved structure and phylogenetic relationships, 13, 5, and 19 SiCCT proteins were classified in the COL, PRR, and CMF subfamilies, respectively. The gene structure and protein conserved motifs analysis exhibited highly similar compositions within the same subfamily. Whole-genome duplication analysis indicated that segmental duplication events played an important role in the expansion of the CCT gene family in foxtail millet. Analysis of transcriptome data showed that 16 SiCCT genes had significant diurnal rhythm oscillations. Under abiotic stress and exogenous hormonal treatment, the expression of many CMF subfamily genes was significantly changed. Especially after drought treatment, the expression of CMF subfamily genes except SiCCT32 was significantly up-regulated. This work provides valuable information for further study of the molecular mechanism of diurnal rhythm regulation, abiotic stress responses, and the identification of candidate genes for foxtail millet molecular breeding.

GWAS, MWAS and mGWAS provide insights into precision agriculture based on genotype-dependent microbial effects in foxtail millet

Genetic and environmental factors collectively determine plant growth and yield. In the past 20 years, genome-wide association studies (GWAS) have been conducted on crops to decipher genetic loci that contribute to growth and yield, however, plant genotype appears to be insufficient to explain the trait variations. Here, we unravel the associations between genotypic, phenotypic, and rhizoplane microbiota variables of 827 foxtail millet cultivars by an integrated GWAS, microbiome-wide association studies (MWAS) and microbiome genome-wide association studies (mGWAS) method. We identify 257 rhizoplane microbial biomarkers associated with six key agronomic traits and validated the microbial-mediated growth effects on foxtail millet using marker strains isolated from the field. The rhizoplane microbiota composition is mainly driven by variations in plant genes related to immunity, metabolites, hormone signaling and nutrient uptake. Among these, the host immune gene FLS2 and transcription factor bHLH35 are widely associated with the microbial taxa of the rhizoplane. We further uncover a plant genotype-microbiota interaction network that contributes to phenotype plasticity. The microbial-mediated growth effects on foxtail millet are dependent on the host genotype, suggesting that precision microbiome management could be used to engineer high-yielding cultivars in agriculture systems.

CRISPR for accelerating genetic gains in under-utilized crops of the drylands: Progress and prospects

Technologies and innovations are critical for addressing the future food system needs where genetic resources are an essential component of the change process. Advanced breeding tools like "genome editing" are vital for modernizing crop breeding to provide game-changing solutions to some of the "must needed" traits in agriculture. CRISPR/Cas-based tools have been rapidly repurposed for editing applications based on their improved efficiency, specificity and reduced off-target effects. Additionally, precise gene-editing tools such as base editing, prime editing, and multiplexing provide precision in stacking of multiple traits in an elite variety, and facilitating specific and targeted crop improvement. This has helped in advancing research and delivery of products in a short time span, thereby enhancing the rate of genetic gains. A special focus has been on food security in the drylands through crops including millets, teff, fonio, quinoa, Bambara groundnut, pigeonpea and cassava. While these crops contribute significantly to the agricultural economy and resilience of the dryland, improvement of several traits including increased stress tolerance, nutritional value, and yields are urgently required. Although CRISPR has potential to deliver disruptive innovations, prioritization of traits should consider breeding product profiles and market segments for designing and accelerating delivery of locally adapted and preferred crop varieties for the drylands. In this context, the scope of regulatory environment has been stated, implying the dire impacts of unreasonable scrutiny of genome-edited plants on the evolution and progress of much-needed technological advances.

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

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

Genomic analyses of rice bean landraces reveal adaptation and yield related loci to accelerate breeding

Rice bean (Vigna umbellata) is an underexploited domesticated legume crop consumed for dietary protein in Asia, yet little is known about the genetic diversity of this species. Here, we present a high-quality reference genome for a rice bean landrace (FF25) built using PacBio long-read data and a Hi-C chromatin interaction map, and assess the phylogenetic position and speciation time of rice bean within the Vigna genus. We sequence 440 landraces (two core collections), and GWAS based on data for growth sites at three widely divergent latitudes reveal loci associated with flowering and yield. Loci harboring orthologs of FUL (FRUITFULL), FT (FLOWERING LOCUS T), and PRR3 (PSEUDO-RESPONSE REGULATOR 3) contribute to the adaptation of rice bean from its low latitude center of origin towards higher latitudes, and the landraces which pyramid early-flowering alleles for these loci display maximally short flowering times. We also demonstrate that copy-number-variation for VumCYP78A6 can regulate seed-yield traits. Intriguingly, 32 landraces collected from a mountainous region in South-Central China harbor a recently acquired InDel in TFL1 (TERMINAL FLOWER1) affecting stem determinacy; these materials also have exceptionally high values for multiple human-desired traits and could therefore substantially advance breeding efforts to improve rice bean.

Rice bean is an underexploited legume crop that has many desirable properties against bio and abiotic stresses. Here, the authors report the genome assembly of this species, conduct population genetics studies and reveal the genetic variations associated with adaptation and yield traits.

Integrative genomics analysis of the ever-shrinking pectin methylesterase (PME) gene family in foxtail millet (Setaria italica)

Pectin methylesterase (PME) plays a vital role in the growth and development of plants. Their genes can be classified into two types, with Type-1 having an extra domain, PMEI. PME genes in foxtail millet (Setaria italica L.) have not been identified, and their sequence features and evolution have not been explored. Here, we identified 41 foxtail millet PME genes. Decoding the pro-region, containing the PMEI domain, revealed its more active nature than the DNA encoding PME domain, easier to be lost to produce Type-2 PME genes. We inferred that the active nature of the pro-region could be related to its harbouring more repetitive DNA sequences. Further, we revealed that though whole-genome duplication and tandem duplication contributed to producing new copies of PME genes, phylogenetic analysis provided clear evidence of ever-shrinking gene family size in foxtail millet and the other grasses in the past 100 million years. Phylogenetic analysis also supports the existence of two gene groups, Group I and Group II, with genes in Group II being more conservative. Our research contributes to understanding how DNA sequence structure affects the functional innovation and evolution of PME genes.

Genome-wide association study in quinoa reveals selection pattern typical for crops with a short breeding history

Quinoa germplasm preserves useful and substantial genetic variation, yet it remains untapped due to a lack of implementation of modern breeding tools. We have integrated field and sequence data to characterize a large diversity panel of quinoa. Whole-genome sequencing of 310 accessions revealed 2.9 million polymorphic high confidence SNP loci. Highland and Lowland quinoa were clustered into two main groups, with F_ST divergence of 0.36 and LD decay of 6.5 and 49.8 Kb, respectively. A genome-wide association study using multi-year phenotyping trials uncovered 600 SNPs stably associated with 17 traits. Two candidate genes are associated with thousand seed weight, and a resistance gene analog is associated with downy mildew resistance. We also identified pleiotropically acting loci for four agronomic traits important for adaptation. This work demonstrates the use of re-sequencing data of an orphan crop, which is partially domesticated to rapidly identify marker-trait association and provides the underpinning elements for genomics-enabled quinoa breeding.

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

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

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

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

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.

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

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

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

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

Fine mapping of a novel QTL CsFSG1 for fruit skin gloss in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.) is an important vegetable crop that is popular with many people. Peel gloss is a highly valued external quality trait that affects the market value of cucumbers, and it directly influences the purchasing psychology of consumers. Nonetheless, the candidate genes and underlying genetic mechanism for this important cucumber trait are not well understood. In this study, we successfully mapped a fruit skin gloss QTL interval to chromosome 3 (26.04-26.14 Mb) using BSA and GWAS methods. Among the eleven candidate genes in the interval, the cytochrome P450 family gene CsCYP86B1 was identified as the candidate for control of fruit skin gloss in cucumber. The expression of CsCYP86B1 in 0-day fruit skin was significantly lower in the low-gloss isogenic line NIL-1334 than in the high-gloss isogenic line NIL-1325. Our findings provide new insights for improving fruit skin gloss in cucumber breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01291-y.