High-Density Mapping of Quantitative Trait Loci Controlling Agronomically Important Traits in Quinoa (Chenopodium quinoa Willd.)

New crops on the block: effective strategies to broaden our food, fibre, and fuel repertoire in the face of increasingly volatile agricultural systems

Climate change poses significant challenges to our ability to keep a growing global population fed, clothed, and fuelled. This review sets the scene by summarizing the impacts of climate change on production of the major grain crop species rice, wheat, and maize, with a focus on yield reductions due to abiotic stresses and altered disease pressures. We discuss efforts to improve resilience, emphasizing traits such as water use efficiency, heat tolerance, and disease resistance. We move on to exploring production trends of established, re-emerging, and new crops, highlighting the challenges of developing and maintaining new arrivals in the global market. We analyse the potential of wild relatives for improving domesticated crops, or as candidates for de novo domestication. The importance of pangenomes for uncovering genetic variation for crop improvement is also discussed. We examine the impact of climate change on non-cereals, including fruit, nut, and fibre crops, and the potential of alternative multiuse crops to increase global sustainability and address climate change-related challenges. Agave is used as an exemplar to demonstrate the strategic pathway for developing a robust new crop option. There is a need for sustained investment in research and development across the entire value chain to facilitate the exploration of diverse species and genetic resources to enhance crop resilience and adaptability to future environmental conditions.

Root restriction accelerates genomic target identification in quinoa under controlled conditions

Quinoa (Chenopodium quinoa) is a nutritious and resilient crop that displays a high genetic and phenotypic variation. As the popularity of this crop increases, there is a growing need to integrate classic and modern breeding tools to favor its improvement. We tested root restriction as a method to reduce plant size and enable high-throughput phenotypic screening of large sets of quinoa plants under controlled conditions. We verified how increasing root restriction does not affect the prediction of field behavior with respect to other standard greenhouse cultivation procedures. We then combined the phenotypic information obtained with our root restriction system with whole-genome re-sequencing data to characterize a quinoa diversity panel of 100 accessions and showed that phenotypic data obtained from root-restricted plants provide real insights into quinoa genetics. Finally, we carried out a genome-wide association study (GWAS) and identified a previously described locus for betalain biosynthesis, as well as other candidate loci linked to betalain biosynthesis and seed size. Overall, we showed that a phenotyping system based on root restriction can aid the identification of genomic targets in quinoa, which can complement and inform field trials for certain traits. This work supports further breeding and faster improvement of quinoa.

Genome assembly of a diversity panel of Chenopodium quinoa

Quinoa (Chenopodium quinoa) is an important crop for the future challenges of food and nutrient security. Deep characterization of quinoa diversity is needed to support the agronomic improvement and adaptation of quinoa as its worldwide cultivation expands. In this study, we report the construction of chromosome-scale genome assemblies of eight quinoa accessions covering the range of phenotypic and genetic diversity of both lowland and highland quinoas. The assemblies were produced from a combination of PacBio HiFi reads and Bionano Saphyr optical maps, with total assembly sizes averaging 1.28 Gb with a mean N50 of 71.1 Mb. Between 43,733 and 48,564 gene models were predicted for the eight new quinoa genomes, and on average, 66% of each quinoa genome was classified as repetitive sequences. Alignment between the eight genome assemblies allowed the identification of structural rearrangements including inversions, translocations, and duplications. These eight novel quinoa genome assemblies provide a resource for association genetics, comparative genomics, and pan-genome analyses for the discovery of genetic components and variations underlying agriculturally important traits.

Early Flowering and Maturity Promote the Successful Adaptation and High Yield of Quinoa (Chenopodium quinoa Willd.) in Temperate Regions

Quinoa (Chenopodium quinoa Willd.) can offer an alternative for staple food considering its tolerance to abiotic stresses and high seed quality. However, its cultivation in temperate regions has not been successful due to its photoperiod sensitivity and low seed yield. This study investigated the agronomical performance and quality traits of 48 accessions for cultivation in northern Europe. We conducted two-year field trials and phenotyped traits related to phenological development, plant architecture, yield components, seed quality, and disease resistance. The major determinants of seed yield in this study were days to flowering, days to maturity, thousand-kernel weight, and panicle density, while downy mildew susceptibility and stem lodging showed a negative correlation with seed yield. We developed a selection index to enable simultaneous selection based on different important agronomical traits. We evaluated the stability of different accessions over the two years of the experiment. Finally, we provided a list of 10 selected accessions that can be directly integrated and serve as new crossing parents in quinoa breeding programs for temperate regions.

Epistasis and pleiotropy-induced variation for plant breeding

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability.

Site-directed genotype screening for elimination of antinutritional saponins in quinoa seeds identifies TSARL1 as a master controller of saponin biosynthesis selectively in seeds

Climate change may result in a drier climate and increased salinization, threatening agricultural productivity worldwide. Quinoa (Chenopodium quinoa) produces highly nutritious seeds and tolerates abiotic stresses such as drought and high salinity, making it a promising future food source. However, the presence of antinutritional saponins in their seeds is an undesirable trait. We mapped genes controlling seed saponin content to a genomic region that includes TSARL1. We isolated desired genetic variation in this gene by producing a large mutant library of a commercial quinoa cultivar and screening the library for specific nucleotide substitutions using droplet digital PCR. We were able to rapidly isolate two independent tsarl1 mutants, which retained saponins in the leaves and roots for defence, but saponins were undetectable in the seed coat. We further could show that TSARL1 specifically controls seed saponin biosynthesis in the committed step after 2,3-oxidosqualene. Our work provides new important knowledge on the function of TSARL1 and represents a breakthrough for quinoa breeding.

Leaf and shoot apical meristem transcriptomes of quinoa (Chenopodium quinoa Willd.) in response to photoperiod and plant development

Understanding the regulation of flowering time is crucial for adaptation of crops to new environment. In this study, we examined the timing of floral transition and analysed transcriptomes in leaf and shoot apical meristems of photoperiod-sensitive and -insensitive quinoa accessions. Histological analysis showed that floral transition in quinoa initiates 2-3 weeks after sowing. We found four groups of differentially expressed genes in quinoa genome that responded to plant development and floral transition: (i) 222 genes responsive to photoperiod in leaves, (ii) 1812 genes differentially expressed between accessions under long-day conditions in leaves, (iii) 57 genes responding to developmental changes under short-day conditions in leaves and (iv) 911 genes responding to floral transition within the shoot apical meristem. Interestingly, among numerous candidate genes, two putative FT orthologs together with other genes (e.g. SOC1, COL, AP1) were previously reported as key regulators of flowering time in other species. Additionally, we used coexpression networks to associate novel transcripts to a putative biological process based on the annotated genes within the same coexpression cluster. The candidate genes in this study would benefit quinoa breeding by identifying and integrating their beneficial haplotypes in crossing programs to develop adapted cultivars to diverse environmental conditions.

Mining genomic regions associated with agronomic and biochemical traits in quinoa through GWAS

Quinoa (Chenopodium quinoa Willd.), an Andean crop, is a facultative halophyte food crop recognized globally for its high nutritional value and plasticity to adapt to harsh conditions. We conducted a genome-wide association study on a diverse set of quinoa germplasm accessions. These accessions were evaluated for the following agronomic and biochemical traits: days to 50% flowering (DTF), plant height (PH), panicle length (PL), stem diameter (SD), seed yield (SY), grain diameter (GD), and thousand-grain weight (TGW). These accessions underwent genotyping-by-sequencing using the DNBSeq-G400R platform. Among all evaluated traits, TGW represented maximum broad-sense heritability. Our study revealed average SNP density of ≈ 3.11 SNPs/10 kb for the whole genome, with the lowest and highest on chromosomes Cq1B and Cq9A, respectively. Principal component analysis clustered the quinoa population in three main clusters, one clearly representing lowland Chilean accessions, whereas the other two groups corresponded to germplasm from the highlands of Peru and Bolivia. In our germplasm set, we estimated linkage disequilibrium decay to be ≈ 118.5 kb. Marker-trait analyses revealed major and consistent effect associations for DTF on chromosomes 3A, 4B, 5B, 6A, 7A, 7B and 8B, with phenotypic variance explained (PVE) as high as 19.15%. Nine associations across eight chromosomes were also found for saponin content with 20% PVE by qSPN5A.1. More QTLs were identified for PL and TGW on multiple chromosomal locations. We identified putative candidate genes in the genomic regions associated with DTF and saponin content. The consistent and major-effect genomic associations can be used in fast-tracking quinoa breeding for wider adaptation across marginal environments.

Empirical phenotyping and genome-wide association study reveal the association of panicle architecture with yield in Chenopodium quinoa

Chenopodium quinoa manifests adaptability to grow under varying agro-climatic scenarios. Assessing quinoa germplasm's phenotypic and genetic variability is a prerequisite for introducing it as a potential candidate in cropping systems. Adaptability is the basic outcome of ecological genomics of crop plants. Adaptive variation predicted with a genome-wide association study provides a valuable basis for marker-assisted breeding. Hence, a panel of 72 quinoa plants was phenotyped for agro morphological attributes and association-mapping for distinct imperative agronomic traits. Inter simple sequence repeat (ISSR) markers were employed to assess genetic relatedness and population structure. Heatmap analysis showed three genotypes were early maturing, and six genotypes were attributed for highest yield. The SD-121-07 exhibited highest yield per plant possessing green, glomerulate shaped, compact density panicle with less leaves. However, SJrecm-03 yielded less exhibiting pink, intermediate shape, intermediate density panicles with less leaves. The phenotyping revealed strong correlation of panicle architecture with yield in quinoa. A genome-wide association study unraveled the associations between ISSR makers and agro-morphological traits. Mixed linear modes analysis yielded nine markers associated with eight traits at p ≤ 0.01. Moreover, ISSR markers significantly associated with panicle shape and leafiness were also associated with yield per plant. These findings contribute to the provision of authenticity for marker-assisted selection that ultimately would support quinoa breeding programs.

A chromosome-scale assembly of the quinoa genome provides insights into the structure and dynamics of its subgenomes

Quinoa (Chenopodium quinoa Willd.) is an allotetraploid seed crop with the potential to help address global food security concerns. Genomes have been assembled for four accessions of quinoa; however, all assemblies are fragmented and do not reflect known chromosome biology. Here, we use in vitro and in vivo Hi-C data to produce a chromosome-scale assembly of the Chilean accession PI 614886 (QQ74). The final assembly spans 1.326 Gb, of which 90.5% is assembled into 18 chromosome-scale scaffolds. The genome is annotated with 54,499 protein-coding genes, 96.9% of which are located on the 18 largest scaffolds. We also report an updated genome assembly for the B-genome diploid C. suecicum and use it, together with the A-genome diploid C. pallidicaule, to identify genomic rearrangements within the quinoa genome, including a large pericentromeric inversion representing 71.7% of chromosome Cq3B. Repetitive sequences comprise 65.2%, 48.6%, and 57.9% of the quinoa, C. pallidicaule, and C. suecicum genomes, respectively. Evidence suggests that the B subgenome is more dynamic and has expanded more than the A subgenome. These genomic resources will enable more accurate assessments of genome evolution within the Amaranthaceae and will facilitate future efforts to identify variation in genes underlying important agronomic traits in quinoa.

Genetic mapping of QTLs controlling brown seed coat traits by genome resequencing in sesame (Sesamum indicum L.)

Introduction

Sesame seeds have become an irreplaceable source of edible oils and food products with rich nutrients and a unique flavor, and their metabolite contents and physiological functions vary widely across different seed coat colors. Although the quantitative trait loci (QTLs) for genetic variation in seed coat color have been extensively investigated, the identification of unique genetic loci for intermediate colors such as brown has not been reported due to their complexity.

Methods

Here, we crossed the white sesame 'Yuzhi No. 8' (YZ8) and the brown sesame 'Yanzhou Erhongpi' (YZEHP) to construct a recombinant inbred line (RIL) population with consecutive self-fertilization for ten generations.

Results

The selfed F1 seeds were brown which was controlled by a dominant gene. Based on the genotyping by whole-genome resequencing of the RILs, a major-effect QTL for brown coat color was identified through both bulk segregant analysis (BSA) and genetic linkage mapping in sesame, which was located within a 1.19 Mb interval on chromosome 6 (qBSCchr6). Moreover, we found that the YZEHP seed coat initially became pigmented at 20 days post-anthesis (DPA) and was substantially colored at 30 DPA. We screened 13 possible candidate genes based on the effects of genetic variants on protein coding and predicted gene functions. Furthermore, qRT‒PCR was used to verify the expression patterns of these genes in different post-anthesis developmental periods. We noted that in comparison to YZ8 seeds, YZEHP seeds had expression of SIN_1023239 that was significantly up-regulated 2.5-, 9.41-, 6.0-, and 5.9-fold at 15, 20, 25, and 30 DPA, respectively, which was consistent with the pattern of brown seed coat pigment accumulation.

Discussion

This study identified the first major-effect QTL for the control of the brown seed coat trait in sesame. This finding lays the foundation for further fine mapping and cloning as well as investigating the regulatory mechanism of seed coat color in sesame.

A comprehensive characterization of agronomic and end-use quality phenotypes across a quinoa world core collection

Quinoa (Chenopodium quinoa Willd.), a pseudocereal with high protein quality originating from the Andean region of South America, has broad genetic variation and adaptability to diverse agroecological conditions, contributing to the potential to serve as a global keystone protein crop in a changing climate. However, the germplasm resources currently available to facilitate quinoa expansion worldwide are restricted to a small portion of quinoa's total genetic diversity, in part because of day-length sensitivity and issues related to seed sovereignty. This study aimed to characterize phenotypic relationships and variation within a quinoa world core collection. The 360 accessions were planted in a randomized complete block design with four replicates in each of two greenhouses in Pullman, WA during the summer of 2018. Phenological stages, plant height, and inflorescence characteristics were recorded. Seed yield, composition, thousand seed weight, nutritional composition, shape, size, and color were measured using a high-throughput phenotyping pipeline. Considerable variation existed among the germplasm. Crude protein content ranged from 11.24% to 17.81% (fixed at 14% moisture). We found that protein content was negatively correlated with yield and positively correlated with total amino acid content and days to harvest. Mean essential amino acids values met adult daily requirements but not leucine and lysine infant requirements. Yield was positively correlated with thousand seed weight and seed area, and negatively correlated with ash content and days to harvest. The accessions clustered into four groups, with one-group representing useful accessions for long-day breeding programs. The results of this study establish a practical resource for plant breeders to leverage as they strategically develop germplasm in support of the global expansion of quinoa.