Draft genome sequence of an inbred line of Chenopodium quinoa, an allotetraploid crop with great environmental adaptability and outstanding nutritional properties

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.

Current trends and prospects in quinoa research: An approach for strategic knowledge areas

Currently, the demand for healthy consumption and the use of alternatives to dairy proteins for the development of foods with good nutritional value are growing. Quinoa has received much attention because it contains a high content of proteins, essential amino acids, essential fatty acids, minerals, vitamins, dietary fibers, and bioactive compounds. Nevertheless, this content and the bioavailability of specific compounds of interest are related to the genotype, the agri-environmental conditions, and management practices where quinoa is grown and postharvest management. This article aimed to analyze the research trends for three knowledge areas: quinoa plant breeding for nutraceutical properties, plant-soil relations focused on abiotic stresses, and postharvest and value-added transformation activities. To this end, a specific methodological design based on bibliometrics and scientometrics methods was used. Through these analyses based on publications' keywords, titles, abstracts, and conclusions sections, for each knowledge area, the key research trends (scope and main topics), the classification of trends based on their development and relevance degree, and the core of knowledge were established. The trends comprise the current state of research. Finally, analyzing the conclusions, recommendations, and future research sections of key publications, a strong correlation among plant breeding research to obtain varieties with tolerance to biotic and abiotic stresses, nutritional and functional compounds of interest for food safety, and the development of products with higher added value established interest in further research on the potential bioactivity of quinoa and the verification of health benefits to humans.

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 resources and breeding approaches for improvement of amaranth (Amaranthus spp.) and quinoa (Chenopodium quinoa)

Nowadays, the human population is more concerned about their diet and very specific in choosing their food sources to ensure a healthy lifestyle and avoid diseases. So people are shifting to more smart nutritious food choices other than regular cereals and staple foods they have been eating for a long time. Pseudocereals, especially, amaranth and quinoa, are important alternatives to traditional cereals due to comparatively higher nutrition, essential minerals, amino acids, and zero gluten. Both Amaranchaceae crops are low-input demanding and hardy plants tolerant to stress, drought, and salinity conditions. Thus, these crops may benefit developing countries that follow subsistence agriculture and have limited farming resources. However, these are underutilized orphan crops, and the efforts to improve them by reducing their saponin content remain ignored for a long time. Furthermore, these crops have very rich variability, but the progress of their genetic gain for getting high-yielding genotypes is slow. Realizing problems in traditional cereals and opting for crop diversification to tackle climate change, research should be focused on the genetic improvement for low saponin, nutritionally rich, tolerant to biotic and abiotic stresses, location-specific photoperiod, and high yielding varietal development of amaranth and quinoa to expand their commercial cultivation. The latest technologies that can accelerate the breeding to improve yield and quality in these crops are much behind and slower than the already established major crops of the world. We could learn from past mistakes and utilize the latest trends such as CRISPR/Cas, TILLING, and RNA interference (RNAi) technology to improve these pseudocereals genetically. Hence, the study reviewed important nutrition quality traits, morphological descriptors, their breeding behavior, available genetic resources, and breeding approaches for these crops to shed light on future breeding strategies to develop superior genotypes.

Genome-wide identification and characterization of glutathione S-transferase gene family in quinoa (Chenopodium quinoa Willd.)

The present investigation was envisaged for large scale in-silico genome wide identification and characterization of glutathione S-transferases (GSTs) in Chenopodium quinoa. In this study, a total of 120 GST genes (CqGSTs) were identified and divided into 11 classes of which tau and phi were highest in numbers. The average protein length of protein was found to be 279.06 with their corresponding average molecular weight of 31,819.4 kDa. The subcellular localization analysis results showed that proteins were centrally localized in the cytoplasm followed by chloroplast, mitochondria and plastids. Structural analysis revealed the presence of 2 -14 exons in CqGST genes. Most of the proteins possessed two exon one intron organization. MEME analysis identified 15 significantly conserved motifs with a width of 6-50 amino acids. Motifs 1, 3, 2, 5, 6, 8, 9 and 13 were found specifically in tau class family; motifs 3, 4, 5, 6, 7 and 9 were found in phi class gene family, while motifs 3, 4, 13 and 14 were found in metaxin class. Multiple sequence alignment revealed highly conserved N-terminus with active site serine (Ser; S) or cysteine (Cys; C) residue for the activation of GSH binding and GST catalytic activity. The gene loci were found to be unevenly distributed across 18 different chromosomes with a maximum of 17 genes located on chromosome number 7. Dominance of alpha helix was followed by coil, extended strand and beta turns. Gene duplication analysis revealed that segmental duplication and purifying type selection were highest in number and found to be main source of expansion of GST gene family. Cis acting regulatory elements analysis showed the presence of 21 different elements involved in stress, hormone and light response and cellular development. The evolutionary relationship of CqGST proteins carried out using maximum likelihood method revealed that all the tau and phi class GSTs were closely associated with those of G. max, O. sativa and A. thaliana. Molecular docking of GST molecules with the fungicide metalaxyl showed that the CqGSTF1 had the lowest binding energy. The comprehensive study of CqGST gene family in quinoa provides groundwork for further functional analysis of CqGST genes in the species at molecular level and has potential applications in plant breeding.

Quinoa (Chenopodium quinoa Willd) Bran Saponins Alleviate Hyperuricemia and Inhibit Renal Injury by Regulating the PI3K/AKT/NFκB Signaling Pathway and Uric Acid Transport

Triterpenoids derived from natural products can exert antihyperuricemic effects. Here, we investigated the antihyperuricemic activity and mechanism of quinoa bran saponins (QBSs) in hyperuricemic mouse and cell models. The QBS4 fraction, with the highest saponin content, was used. Fourier-transform infrared, high-performance liquid chromatography, and ultrahigh-performance liquid chromatography-mass spectrometry identified 11 individual saponins in QBS4, of which the main components were hederagenin and oleanolic acid. The QBS4 effects on hyperuricemic mice (induced by adenine and potassium oxonate) were then studied. QBS4 reduced the levels of uric acid (UA), serum urea nitrogen, creatinine, and lipids in mice with hyperuricemia (HUA) and decreased renal inflammation and renal damage. Molecular analysis revealed that QBS4 may alleviate HUA by regulating the expression of key genes involved in the transport of UA and by inhibiting the activation of the PI3K/AKT/NFκB inflammatory signaling pathway. In conclusion, QBS4 has promise for using as a natural dietary supplement to treat and prevent HUA.

Shotgun proteomics of quinoa seeds reveals chitinases enrichment under rainfed conditions

Quinoa is an Andean crop whose cultivation has been extended to many different parts of the world in the last decade. It shows a great capacity for adaptation to diverse climate conditions, including environmental stressors, and, moreover, the seeds are very nutritious in part due to their high protein content, which is rich in essential amino acids. They are gluten-free seeds and contain good amounts of other nutrients such as unsaturated fatty acids, vitamins, or minerals. Also, the use of quinoa hydrolysates and peptides has been linked to numerous health benefits. Altogether, these aspects have situated quinoa as a crop able to contribute to food security worldwide. Aiming to deepen our understanding of the protein quality and function of quinoa seeds and how they can vary when this crop is subjected to water-limiting conditions, a shotgun proteomics analysis was performed to obtain the proteomes of quinoa seeds harvested from two different water regimes in the field: rainfed and irrigated conditions. Differentially increased levels of proteins determined in seeds from each field condition were analysed, and the enrichment of chitinase-related proteins in seeds harvested from rainfed conditions was found. These proteins are described as pathogen-related proteins and can be accumulated under abiotic stress. Thus, our findings suggest that chitinase-like proteins in quinoa seeds can be potential biomarkers of drought. Also, this study points to the need for further research to unveil their role in conferring tolerance when coping with water-deficient conditions.

Snapshot of four mature quinoa (Chenopodium quinoa) seeds: a shotgun proteomics analysis with emphasis on seed maturation, reserves and early germination

Chenopodium quinoa Willd. is a crop species domesticated over 5000 years ago. This species is highly diverse, with a geographical distribution that covers more than 5000 km from Colombia to Chile, going through a variety of edaphoclimatic conditions. Quinoa grains have great nutritional quality, raising interest at a worldwide level. In this work, by using shotgun proteomics and in silico analysis, we present an overview of mature quinoa seed proteins from a physiological context and considering the process of seed maturation and future seed germination. For this purpose, we selected grains from four contrasting quinoa cultivars (Amarilla de Maranganí, Chadmo, Sajama and Nariño) with different edaphoclimatic and geographical origins. The results give insight on the most important metabolic pathways for mature quinoa seeds including: starch synthesis, protein bodies and lipid bodies composition, reserves and their mobilization, redox homeostasis, and stress related proteins like heat-shock proteins (HSPs) and late embryogenesis abundant proteins (LEAs), as well as evidence for capped and uncapped mRNA translation. LEAs present in our analysis show a specific pattern of expression matching that of other species. Overall, this work presents a complete snapshot of quinoa seeds physiological context, providing a reference point for further studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01295-8.

Functional analysis of CqPORB in the regulation of chlorophyll biosynthesis in Chenopodium quinoa

Protochlorophyllide oxidoreductase (POR) plays a key role in catalyzing the light-dependent reduction of protochlorophyllide (Pchlide) to chlorophyllide (Chlide), and thus promotes the transit from etiolated seedlings to green plants. In this study, by exploring ethyl methanesulfonate (EMS)-mediated mutagenesis in Chenopodium quinoa NL-6 variety, we identified a mutant nl6-35 that displays faded green leaf and reduced chlorophyll (Chl) and carotenoid contents. Bulk segregant analysis (BSA) revealed that a mutation in CqPORB gene is genetically associated with the faded green leaf of the nl6-35 mutant. Further study indicates that the nl6-35 mutant exhibits abnormal grana stacks and compromised conversion of Pchlide to Chlide upon illumination, suggesting the important role of CqPORB in producing photoactive Pchlide. Totally three CqPOR isoforms, including CqPORA, CqPORA-like, and CqPORB are identified in NL-6 variety. Transcriptional analysis shows that the expression of all these three CqPOR isoforms is regulated in light- and development-dependent manners, and in mature quinoa plants only CqPORB isoform is predominantly expressed. Subcellular localization analysis indicates that CqPORB is exclusively localized in chloroplast. Together, our study elucidates the important role of CqPORB in the regulation of Chl biosynthesis and chloroplast development in quinoa.

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Quinoa (Chenopodium quinoa)

BACKGROUND

Squamous promoter binding protein-like (SPL) proteins are a class of transcription factors that play essential roles in plant growth and development, signal transduction, and responses to biotic and abiotic stresses. The rapid development of whole genome sequencing has enabled the identification and characterization of SPL gene families in many plant species, but to date this has not been performed in quinoa (Chenopodium quinoa).

RESULTS

This study identified 23 SPL genes in quinoa, which were unevenly distributed on 18 quinoa chromosomes. Quinoa SPL genes were then classified into eight subfamilies based on homology to Arabidopsis thaliana SPL genes. We selected three dicotyledonous and monocotyledonous representative species, each associated with C. quinoa, for comparative sympatric mapping to better understand the evolution of the developmental mechanisms of the CqSPL family. Furthermore, we also used 15 representative genes from eight subfamilies to characterize CqSPLs gene expression in different tissues and at different fruit developmental stages under six different abiotic stress conditions.

CONCLUSIONS

This study, the first to identify and characterize SPL genes in quinoa, reported that CqSPL genes, especially CqSPL1, play a critical role in quinoa development and in its response to various abiotic stresses.

Genomic Variation Underlying the Breeding Selection of Quinoa Varieties Longli-4 and CA3-1 in China

Quinoa (Chenopodium quinoa) is a well-known climate-resilient crop and has been introduced into multiple marginal lands across the world, including China, to improve food security and/or balanced nutrient supplies. Conventional breeding has been widely applied in the selection and breeding of quinoa varieties in China since 1980s; however, few studies have been implemented on the genetic variances among different varieties developed by diversity breeding objectives. In this study, the phenotypic and genetic differences between two varieties (Longli-4 and CA3-1) from China were systematically analyzed. A total of 407,651 and 2,731,411 single nucleotide polymorphisms (SNPs) and 212,724 and 587,935 small insertion and deletion (INDELs) were detected for Longli-4 and CA3-1, respectively, when compared with the reference genome of PI614886. The SNPs/INDELs were unevenly distributed across each chromosome for both varieties. There were 143,996 SNPs and 83,410 INDELs shared between Longli-4 and CA3-1, accounting for 4% of the total variances. The variation was then screened based on the SNP effects. There were 818 and 73 genes with the variety-specific non-synonymous and stop-gain variation in Longli-4, whereas there were 13,701 and 733 genes in CA3-1. Specifically, 3501 genes with the non-synonymous variation and 74 genes with the stop-gain variation were found in both Longli-4 and CA3-1. These results suggest that convergent selection occurred during the different breeding processes. A set of candidate genes related to agronomic traits and domestication were further selected to detect the genetic divergence in detail in the two varieties. Only one domestication gene was identified having Longli-4-specific stop-gain variation. Twelve candidate genes related to betalain (1), flowering (4), seed size (2), domestication (1), and saponin (4) were identified having CA3-1-specific stop-gain variation. Interestingly, one seed size gene homologous of CKX1 (cytokinin oxidase/dehydrogenase 1) had the stop-gain variation in both varieties. This research will therefore provide guidance for the molecular-assisted breeding in quinoa.

The Disease Progression and Molecular Defense Response in Chenopodium Quinoa Infected with Peronospora Variabilis, the Causal Agent of Quinoa Downy Mildew

Downy mildew disease, caused by the biotrophic oomycete Peronospora variabilis, is the largest threat to the cultivation of quinoa (Chenopodium quinoa Willd.) in the Andean highlands, and occurs worldwide. However, so far, no molecular study of the quinoa-Peronospora interaction has been reported. Here, we developed tools to study downy mildew disease in quinoa at the gene expression level. P. variabilis was isolated and maintained, allowing the study of downy mildew disease progression in two quinoa cultivars under controlled conditions. Quinoa gene expression changes induced by P. variabilis were analyzed by qRT-PCR, for quinoa homologues of A. thaliana pathogen-associated genes. Overall, we observed a slower disease progression and higher tolerance in the quinoa cultivar Kurmi than in the cultivar Maniqueña Real. The quinoa orthologs of putative defense genes such as the catalase CqCAT2 and the endochitinase CqEP3 showed no changes in gene expression. In contrast, quinoa orthologs of other defense response genes such as the transcription factor CqWRKY33 and the chaperone CqHSP90 were significantly induced in plants infected with P. variabilis. These genes could be used as defense response markers to select quinoa cultivars that are more tolerant to P. variabilis infection.

Genome-Wide Analysis and Expression Profiles of the VOZ Gene Family in Quinoa (Chenopodium quinoa)

Vascular plant one zinc-finger (VOZ) proteins are a plant-specific transcription factor family and play important roles in plant development and stress responses. However, little is known about the VOZ genes in quinoa. In the present study, a genome-wide investigation of the VOZ gene family in quinoa was performed, including gene structures, conserved motifs, phylogeny, and expression profiles. A total of four quinoa VOZ genes distributed on three chromosomes were identified. Based on phylogenetic analysis, CqVOZ1 and CqVOZ3 belong to subfamily II, and CqVOZ2 and CqVOZ4 belong to subfamily III. Furthermore, the VOZ transcription factors of quinoa and sugarbeet were more closely related than other species. Except for CqVOZ3, all the other three CqVOZs have four exons and four introns. Analysis of conserved motifs indicated that each CqVOZ member contained seven common motifs. Multiple sequence alignment showed that the CqVOZ genes were highly conserved with consensus sequences, which might be plausibly significant for the preservation of structural integrity of the family proteins. Tissue expression analysis revealed that four CqVOZ genes were highly expressed in inflorescence and relatively low in leaves and stems, suggesting that these genes had obvious tissue expression specificity. The expression profiles of the quinoa CqVOZs under various abiotic stresses demonstrated that these genes were differentially induced by cold stress, salt stress, and drought stress. The transcript level of CqVOZ1 and CqVOZ4 were down-regulated by salt stress and drought stress, while CqVOZ2 and CqVOZ3 were up-regulated by cold, salt, and drought stress, which could be used as abiotic stress resistance candidate genes. This study systematically identifies the CqVOZ genes at the genome-wide level, contributing to a better understanding of the quinoa VOZ transcription factor family and laying a foundation for further exploring the molecular mechanism of development and stress resistance of quinoa.

Stalk cell polar ion transport provide for bladder-based salinity tolerance in Chenopodium quinoa

Chenopodium quinoa uses epidermal bladder cells (EBCs) to sequester excess salt. Each EBC complex consists of a leaf epidermal cell, a stalk cell, and the bladder. Under salt stress, sodium (Na⁺ ), chloride (Cl^- ), potassium (K⁺ ) and various metabolites are shuttled from the leaf lamina to the bladders. Stalk cells operate as both a selectivity filter and a flux controller. In line with the nature of a transfer cell, advanced transmission electron tomography, electrophysiology, and fluorescent tracer flux studies revealed the stalk cell's polar organization and bladder-directed solute flow. RNA sequencing and cluster analysis revealed the gene expression profiles of the stalk cells. Among the stalk cell enriched genes, ion channels and carriers as well as sugar transporters were most pronounced. Based on their electrophysiological fingerprint and thermodynamic considerations, a model for stalk cell transcellular transport was derived.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

Chromosome-Level Genome Assembly Provides New Insights into Genome Evolution and Tuberous Root Formation of Potentilla anserina

Potentilla anserina is a perennial stoloniferous plant with edible tuberous roots in Rosaceae, served as important food and medicine sources for Tibetans in the Qinghai-Tibetan Plateau (QTP), China, over thousands of years. However, a lack of genome information hindered the genetic study. Here, we presented a chromosome-level genome assembly using single-molecule long-read sequencing, and the Hi-C technique. The assembled genome was 454.28 Mb, containing 14 chromosomes, with contig N50 of 2.14 Mb. A total of 46,495 protein-coding genes, 169.74 Mb repeat regions, and 31.76 Kb non-coding RNA were predicted. P. anserina diverged from Potentilla micrantha ∼28.52 million years ago (Mya). Furthermore, P. anserina underwent a recent tetraploidization ∼6.4 Mya. The species-specific genes were enriched in Starch and sucrose metabolism and Galactose metabolism pathways. We identified the sub-genome structures of P. anserina, with A sub-genome was larger than B sub-genome and closer to P. micrantha phylogenetically. Despite lacking significant genome-wide expression dominance, the A sub-genome had higher homoeologous gene expression in shoot apical meristem, flower and tuberous root. The resistance genes was contracted in P. anserina genome. Key genes involved in starch biosynthesis were expanded and highly expressed in tuberous roots, which probably drives the tuber formation. The genomics and transcriptomics data generated in this study advance our understanding of the genomic landscape of P. anserina, and will accelerate genetic studies and breeding programs.

Transcriptomic Analysis of Quinoa Reveals a Group of Germin-Like Proteins Induced by Trichoderma

Symbiotic strains of fungi in the genus Trichoderma affect growth and pathogen resistance of many plant species, but the interaction is not known in molecular detail. Here we describe the transcriptomic response of two cultivars of the crop Chenopodium quinoa to axenic co-cultivation with Trichoderma harzianum BOL-12 and Trichoderma afroharzianum T22. The response of C. quinoa roots to BOL-12 and T22 in the early phases of interaction was studied by RNA sequencing and RT-qPCR verification. Interaction with the two fungal strains induced partially overlapping gene expression responses. Comparing the two plant genotypes, a broad spectrum of putative quinoa defense genes were found activated in the cultivar Kurmi but not in the Real cultivar. In cultivar Kurmi, relatively small effects were observed for classical pathogen response pathways but instead a C. quinoa-specific clade of germin-like genes were activated. Germin-like genes were found to be more rapidly induced in cultivar Kurmi as compared to Real. The same germin-like genes were found to also be upregulated systemically in the leaves. No strong correlation was observed between any of the known hormone-mediated defense response pathways and any of the quinoa-Trichoderma interactions. The differences in responses are relevant for the capabilities of applying Trichoderma agents for crop protection of different cultivars of C. quinoa.

Quinoa genome assembly employing genomic variation for guided scaffolding

We propose to use the natural variation between individuals of a population for genome assembly scaffolding. In today's genome projects, multiple accessions get sequenced, leading to variant catalogs. Using such information to improve genome assemblies is attractive both cost-wise as well as scientifically, because the value of an assembly increases with its contiguity. We conclude that haplotype information is a valuable resource to group and order contigs toward the generation of pseudomolecules. Quinoa (Chenopodium quinoa) has been under cultivation in Latin America for more than 7500 years. Recently, quinoa has gained increasing attention due to its stress resistance and its nutritional value. We generated a novel quinoa genome assembly for the Bolivian accession CHEN125 using PacBio long-read sequencing data (assembly size 1.32 Gbp, initial N50 size 608 kbp). Next, we re-sequenced 50 quinoa accessions from Peru and Bolivia. This set of accessions differed at 4.4 million single-nucleotide variant (SNV) positions compared to CHEN125 (1.4 million SNV positions on average per accession). We show how to exploit variation in accessions that are distantly related to establish a genome-wide ordered set of contigs for guided scaffolding of a reference assembly. The method is based on detecting shared haplotypes and their expected continuity throughout the genome (i.e., the effect of linkage disequilibrium), as an extension of what is expected in mapping populations where only a few haplotypes are present. We test the approach using Arabidopsis thaliana data from different populations. After applying the method on our CHEN125 quinoa assembly we validated the results with mate-pairs, genetic markers, and another quinoa assembly originating from a Chilean cultivar. We show consistency between these information sources and the haplotype-based relations as determined by us and obtain an improved assembly with an N50 size of 1079 kbp and ordered contig groups of up to 39.7 Mbp. We conclude that haplotype information in distantly related individuals of the same species is a valuable resource to group and order contigs according to their adjacency in the genome toward the generation of pseudomolecules.

Developing Chenopodium ficifolium as a potential B genome diploid model system for genetic characterization and improvement of allotetraploid quinoa (Chenopodium quinoa)

BACKGROUND

Quinoa (Chenopodium quinoa) is a high-value grain known for its excellent nutritional balance. It is an allotetraploid species (AABB, 2n = 4x = 36) formed by the hybridization between AA and BB genome diploid (2n = 2x = 18) species. This study reports genetic studies in Chenopodium ficifolium as a potential B genome diploid model system to simplify the genetic studies of quinoa including gene identification and marker-assisted breeding.

RESULTS

Portsmouth, New Hampshire and Quebec City, Quebec accessions of C. ficifolium were used to develop an F2 population segregating for agronomically relevant traits including flowering time, plant height, the number of branches, branch angle, and internode length. Marker-trait associations were identified for the FLOWERING LOCUS T-LIKE 1 (FTL1) marker gene, where the alternate alleles (A1/A2) were segregating among the F2 generation plants in association with flowering time, plant height, and the number of branches. There was a strong correlation of the flowering time trait with both plant height and the number of branches. Thus, a possible multifaceted functional role for FTL1 may be considered. The parental Portsmouth and Quebec City accessions were homozygous for the alternate FTL1 alleles, which were found to be substantially diverged. SNPs were identified in the FTL1 coding sequence that could have some functional significance in relation to the observed trait variation.

CONCLUSION

These results draw further attention to the possible functional roles of the FTL1 locus in Chenopodium and justify continued exploration of C. ficifolium as a potential diploid model system for the genetic study of quinoa. We expect our findings to aid in quinoa breeding as well as to any studies related to the Chenopodium genus.

Genome-Wide Identification and Analysis of the Phosphoenolpyruvate Carboxylase Gene Family in Suaeda aralocaspica, an Annual Halophyte With Single-Cellular C4 Anatomy

Phosphoenolpyruvate carboxylase (PEPC) plays pivotal roles in the carbon fixation of photosynthesis and a variety of metabolic and stress pathways. Suaeda aralocaspica belongs to a single-cellular C₄ species and carries out a photosynthetic pathway in an unusually elongated chlorenchyma cell, which is expected to have PEPCs with different characteristics. To identify the different isoforms of PEPC genes in S. aralocaspica and comparatively analyze their expression and regulation patterns as well as the biochemical and enzymatic properties in this study, we characterized a bacterial-type PEPC (BTPC; SaPEPC-4) in addition to the two plant-type PEPCs (PTPCs; SaPEPC-1 and SaPEPC-2) using a genome-wide identification. SaPEPC-4 presented a lower expression level in all test combinations with an unknown function; two SaPTPCs showed distinct subcellular localizations and different spatiotemporal expression patterns but positively responded to abiotic stresses. Compared to SaPEPC-2, the expression of SaPEPC-1 specifically in chlorenchyma cell tissues was much more active with the progression of development and under various stresses, particularly sensitive to light, implying the involvement of SaPEPC-1 in a C₄ photosynthetic pathway. In contrast, SaPEPC-2 was more like a non-photosynthetic PEPC. The expression trends of two SaPTPCs in response to light, development, and abiotic stresses were also matched with the changes in PEPC activity in vivo (native) or in vitro (recombinant), and the biochemical properties of the two recombinant SaPTPCs were similar in response to various effectors while the catalytic efficiency, substrate affinity, and enzyme activity of SaPEPC-2 were higher than that of SaPEPC-1 in vitro. All the different properties between these two SaPTPCs might be involved in transcriptional (e.g., specific cis-elements), posttranscriptional [e.g., 5'-untranslated region (5'-UTR) secondary structure], or translational (e.g., PEPC phosphorylation/dephosphorylation) regulatory events. The comparative studies on the different isoforms of the PEPC gene family in S. aralocaspica may help to decipher their exact role in C₄ photosynthesis, plant growth/development, and stress resistance.

Integrating Omics and Gene Editing Tools for Rapid Improvement of Traditional Food Plants for Diversified and Sustainable Food Security

Indigenous communities across the globe, especially in rural areas, consume locally available plants known as Traditional Food Plants (TFPs) for their nutritional and health-related needs. Recent research shows that many TFPs are highly nutritious as they contain health beneficial metabolites, vitamins, mineral elements and other nutrients. Excessive reliance on the mainstream staple crops has its own disadvantages. Traditional food plants are nowadays considered important crops of the future and can act as supplementary foods for the burgeoning global population. They can also act as emergency foods in situations such as COVID-19 and in times of other pandemics. The current situation necessitates locally available alternative nutritious TFPs for sustainable food production. To increase the cultivation or improve the traits in TFPs, it is essential to understand the molecular basis of the genes that regulate some important traits such as nutritional components and resilience to biotic and abiotic stresses. The integrated use of modern omics and gene editing technologies provide great opportunities to better understand the genetic and molecular basis of superior nutrient content, climate-resilient traits and adaptation to local agroclimatic zones. Recently, realizing the importance and benefits of TFPs, scientists have shown interest in the prospection and sequencing of TFPs for their improvements, cultivation and mainstreaming. Integrated omics such as genomics, transcriptomics, proteomics, metabolomics and ionomics are successfully used in plants and have provided a comprehensive understanding of gene-protein-metabolite networks. Combined use of omics and editing tools has led to successful editing of beneficial traits in several TFPs. This suggests that there is ample scope for improvement of TFPs for sustainable food production. In this article, we highlight the importance, scope and progress towards improvement of TFPs for valuable traits by integrated use of omics and gene editing techniques.

The genotype-dependent phenotypic landscape of quinoa in salt tolerance and key growth traits

Cultivation of quinoa (Chenopodium quinoa), an annual pseudocereal crop that originated in the Andes, is spreading globally. Because quinoa is highly nutritious and resistant to multiple abiotic stresses, it is emerging as a valuable crop to provide food and nutrition security worldwide. However, molecular analyses have been hindered by the genetic heterogeneity resulting from partial outcrossing. In this study, we generated 136 inbred quinoa lines as a basis for the molecular identification and characterization of gene functions in quinoa through genotyping and phenotyping. Following genotyping-by-sequencing analysis of the inbred lines, we selected 5,753 single-nucleotide polymorphisms (SNPs) in the quinoa genome. Based on these SNPs, we show that our quinoa inbred lines fall into three genetic sub-populations. Moreover, we measured phenotypes, such as salt tolerance and key growth traits in the inbred quinoa lines and generated a heatmap that provides a succinct overview of the genotype–phenotype relationship between inbred quinoa lines. We also demonstrate that, in contrast to northern highland lines, most lowland and southern highland lines can germinate even under high salinity conditions. These findings provide a basis for the molecular elucidation and genetic improvement of quinoa and improve our understanding of the evolutionary process underlying quinoa domestication.

Virus-Mediated Transient Expression Techniques Enable Functional Genomics Studies and Modulations of Betalain Biosynthesis and Plant Height in Quinoa

Quinoa (Chenopodium quinoa), native to the Andean region of South America, has been recognized as a potentially important crop in terms of global food and nutrition security since it can thrive in harsh environments and has an excellent nutritional profile. Even though challenges of analyzing the complex and heterogeneous allotetraploid genome of quinoa have recently been overcome, with the whole genome-sequencing of quinoa and the creation of genotyped inbred lines, the lack of technology to analyze gene function in planta is a major limiting factor in quinoa research. Here, we demonstrate that two virus-mediated transient expression techniques, virus-induced gene silencing (VIGS) and virus-mediated overexpression (VOX), can be used in quinoa. We show that apple latent spherical virus (ALSV) can induce gene silencing of quinoa phytoene desaturase (CqPDS1) in a broad range of quinoa inbred lines derived from the northern and southern highland and lowland sub-populations. In addition, we show that ALSV can be used as a VOX vector in roots. Our data also indicate that silencing a quinoa 3,4-dihydroxyphenylalanine 4,5-dioxygenase gene (CqDODA1) or a cytochrome P450 enzyme gene (CqCYP76AD1) inhibits betalain production and that knockdown of a reduced-height gene homolog (CqRHT1) causes an overgrowth phenotype in quinoa. Moreover, we show that ALSV can be transmitted to the progeny of quinoa plants. Thus, our findings enable functional genomics in quinoa, ushering in a new era of quinoa research.

The genome of Cleistogenes songorica provides a blueprint for functional dissection of dimorphic flower differentiation and drought adaptability

Cleistogenes songorica (2n = 4x = 40) is a desert grass with a unique dimorphic flowering mechanism and an ability to survive extreme drought. Little is known about the genetics underlying drought tolerance and its reproductive adaptability. Here, we sequenced and assembled a high-quality chromosome-level C. songorica genome (contig N50 = 21.28 Mb). Complete assemblies of all telomeres, and of ten chromosomes were derived. C. songorica underwent a recent tetraploidization (~19 million years ago) and four major chromosomal rearrangements. Expanded genes were significantly enriched in fatty acid elongation, phenylpropanoid biosynthesis, starch and sucrose metabolism, and circadian rhythm pathways. By comparative transcriptomic analysis we found that conserved drought tolerance related genes were expanded. Transcription of CsMYB genes was associated with differential development of chasmogamous and cleistogamous flowers, as well as drought tolerance. Furthermore, we found that regulation modules encompassing miRNA, transcription factors and target genes are involved in dimorphic flower development, validated by overexpression of CsAP2_9 and its targeted miR172 in rice. Our findings enable further understanding of the mechanisms of drought tolerance and flowering in C. songorica, and provide new insights into the adaptability of native grass species in evolution, along with potential resources for trait improvement in agronomically important species.

A multilayered cross-species analysis of GRAS transcription factors uncovered their functional networks in plant adaptation to the environment

Introduction

Environmental stress is both a major force of natural selection and a prime factor affecting crop qualities and yields. The impact of the GRAS [gibberellic acid-insensitive (GAI), repressor of GA1-3 mutant (RGA), and scarecrow (SCR)] family on plant development and the potential to resist environmental stress needs much emphasis.

Objectives

This study aims to investigate the evolution, expansion, and adaptive mechanisms of GRASs of important representative plants during polyploidization.

Methods

We explored the evolutionary characteristics of GRASs in 15 representative plant species by systematic biological analysis of the genome, transcriptome, metabolite, protein complex map and phenotype.

Results

The GRAS family was systematically identified from 15 representative plant species of scientific and agricultural importance. The detection of gene duplication types of GRASs in all species showed that the widespread expansion of GRASs in these species was mainly contributed by polyploidization events. Evolutionary analysis reveals that most species experience independent genome-wide duplication (WGD) events and that interspecies GRAS functions may be broadly conserved. Polyploidy-related Chenopodium quinoa GRASs (CqGRASs) and Arabidopsis thaliana GRASs (AtGRASs) formed robust networks with flavonoid pathways by crosstalk with auxin and photosynthetic pathways. Furthermore, Arabidopsis thaliana population transcriptomes and the 1000 Plants (OneKP) project confirmed that GRASs are components of flavonoid biosynthesis, which enables plants to adapt to the environment by promoting flavonoid accumulation. More importantly, the GRASs of important species that may potentially improve important agronomic traits were mapped through TAIR and RARGE-II publicly available phenotypic data. Determining protein interactions and target genes contributes to determining GRAS functions.

Conclusion

The results of this study suggest that polyploidy-related GRASs in multiple species may be a target for improving plant growth, development, and environmental adaptation.

Mainstreaming Barahnaja cultivation for food and nutritional security in the Himalayan region

Selective production of input intensive crops in the present scenario have resulted in productivity stagnation or even decline due to excessive usage of chemicals, affecting the farmers economically. Sustainable agriculture is the way to increase agricultural productivity and economic prosperity by protecting all natural resources. It maintains a balance of soil fertility with crop productivity and nutritional quality. The mixed cropping systems followed earlier in different regions according to their tradition, climatic zone, soil and water conditions were climate-smart approaches to sustainable food production based on practical experiences over the years of old generations. The life style changes, imbalance in farming system in last 70 years and demand for more food as well as declining land resources resulted in intensive agriculture. Besides, least returns and less demand of ethnic crops gave more preference to major staple food crops. Barahnaja is a traditional orphan crops based mixed cropping system practiced in Himalayan region due to its sustainability and assured crop harvest during erratic weather conditions. This traditional farming method is an exemplary scientific approach to derive innovations with respect to productivity, quality, plant soil interactions and organic agriculture. The main focus of the review is to substantiate the characteristics of the traditional mixed cropping system by describing the advantages of the system and opportunities for scientific innovation towards new knowledge and sustainability.

Orphan Crops and their Wild Relatives in the Genomic Era

More than half of the calories consumed by humans are provided by three major cereal crops (rice, maize, and wheat). Orphan crops are usually well adapted to low-input agricultural conditions, and they not only play vital roles in local areas but can also contribute to food and nutritional needs worldwide. Interestingly, many wild relatives of orphan crops are important weeds of major crops. Although orphan crops and their wild relatives have received little attentions from researchers for many years, genomic studies have recently been performed on these plants. Here, we provide an overview of genomic studies on orphan crops, with a focus on orphan cereals and their wild relatives. The genomes of at least 12 orphan cereals and/or their wild relatives have been sequenced. In addition to genomic benefits for orphan crop breeding, we discuss the potential ways for mutual utilization of genomic data from major crops, orphan crops, and their wild relatives (including weeds) and provide perspectives on genetic improvement of both orphan and major crops (including de novo domestication of orphan crops) in the coming genomic era.

The Cassandra retrotransposon landscape in sugar beet (Beta vulgaris) and related Amaranthaceae: recombination and re-shuffling lead to a high structural variability

BACKGROUND AND AIMS

Plant genomes contain many retrotransposons and their derivatives, which are subject to rapid sequence turnover. As non-autonomous retrotransposons do not encode any proteins, they experience reduced selective constraints leading to their diversification into multiple families, usually limited to a few closely related species. In contrast, the non-coding Cassandra terminal repeat retrotransposons in miniature (TRIMs) are widespread in many plants. Their hallmark is a conserved 5S rDNA-derived promoter in their long terminal repeats (LTRs). As sugar beet (Beta vulgaris) has a well-described LTR retrotransposon landscape, we aim to characterize TRIMs in beet and related genomes.

METHODS

We identified Cassandra retrotransposons in the sugar beet reference genome and characterized their structural relationships. Genomic organization, chromosomal localization, and distribution of Cassandra-TRIMs across the Amaranthaceae were verified by Southern and fluorescent in situ hybridization.

KEY RESULTS

All 638 Cassandra sequences in the sugar beet genome contain conserved LTRs and thus constitute a single family. Nevertheless, variable internal regions required a subdivision into two Cassandra subfamilies within B. vulgaris. The related Chenopodium quinoa harbours a third subfamily. These subfamilies vary in their distribution within Amaranthaceae genomes, their insertion times and the degree of silencing by small RNAs. Cassandra retrotransposons gave rise to many structural variants, such as solo LTRs or tandemly arranged Cassandra retrotransposons. These Cassandra derivatives point to an interplay of template switch and recombination processes - mechanisms that likely caused Cassandra's subfamily formation and diversification.

CONCLUSIONS

We traced the evolution of Cassandra in the Amaranthaceae and detected a considerable variability within the short internal regions, whereas the LTRs are strongly conserved in sequence and length. Presumably these hallmarks make Cassandra a prime target for unequal recombination, resulting in the observed structural diversity, an example of the impact of LTR-mediated evolutionary mechanisms on the host genome.

Transcriptome analysis and differential gene expression profiling of two contrasting quinoa genotypes in response to salt stress

BACKGROUND

Soil salinity is one of the major abiotic stress factors that affect crop growth and yield, which seriously restricts the sustainable development of agriculture. Quinoa is considered as one of the most promising crops in the future for its high nutrition value and strong adaptability to extreme weather and soil conditions. However, the molecular mechanisms underlying the adaptive response to salinity stress of quinoa remain poorly understood. To identify candidate genes related to salt tolerance, we performed reference-guided assembly and compared the gene expression in roots treated with 300 mM NaCl for 0, 0.5, 2, and 24 h of two contrasting quinoa genotypes differing in salt tolerance.

RESULTS

The salt-tolerant (ST) genotype displayed higher seed germination rate and plant survival rate, and stronger seedling growth potential as well than the salt-sensitive (SS) genotype under salt stress. An average of 38,510,203 high-quality clean reads were generated. Significant Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified to deeper understand the differential response. Transcriptome analysis indicated that salt-responsive genes in quinoa were mainly related to biosynthesis of secondary metabolites, alpha-Linolenic acid metabolism, plant hormone signal transduction, and metabolic pathways. Moreover, several pathways were significantly enriched amongst the differentially expressed genes (DEGs) in ST genotypes, such as phenylpropanoid biosynthesis, plant-pathogen interaction, isoquinoline alkaloid biosynthesis, and tyrosine metabolism. One hundred seventeen DEGs were common to various stages of both genotypes, identified as core salt-responsive genes, including some transcription factor members, like MYB, WRKY and NAC, and some plant hormone signal transduction related genes, like PYL, PP2C and TIFY10A, which play an important role in the adaptation to salt conditions of this species. The expression patterns of 21 DEGs were detected by quantitative real-time PCR (qRT-PCR) and confirmed the reliability of the RNA-Seq results.

CONCLUSIONS

We identified candidate genes involved in salt tolerance in quinoa, as well as some DEGs exclusively expressed in ST genotype. The DEGs common to both genotypes under salt stress may be the key genes for quinoa to adapt to salinity environment. These candidate genes regulate salt tolerance primarily by participating in reactive oxygen species (ROS) scavenging system, protein kinases biosynthesis, plant hormone signal transduction and other important biological processes. These findings provide theoretical basis for further understanding the regulation mechanism underlying salt tolerance network of quinoa, as well establish foundation for improving its tolerance to salinity in future breeding programs.

A novel WD40-repeat protein involved in formation of epidermal bladder cells in the halophyte quinoa

Halophytes are plants that grow in high-salt environments and form characteristic epidermal bladder cells (EBCs) that are important for saline tolerance. To date, however, little has been revealed about the formation of these structures. To determine the genetic basis for their formation, we applied ethylmethanesulfonate mutagenesis and obtained two mutants with reduced levels of EBCs (rebc) and abnormal chloroplasts. In silico subtraction experiments revealed that the rebc phenotype was caused by mutation of REBC, which encodes a WD40 protein that localizes to the nucleus and chloroplasts. Phylogenetic and transformant analyses revealed that the REBC protein differs from TTG1, a WD40 protein involved in trichome formation. Furthermore, rebc mutants displayed damage to their shoot apices under abiotic stress, suggesting that EBCs may protect the shoot apex from such stress. These findings will help clarify the mechanisms underlying EBC formation and function.

Genome-wide survey, characterization, and expression analysis of bZIP transcription factors in Chenopodium quinoa

BACKGROUND

Chenopodium quinoa Willd. (quinoa) is a pseudocereal crop of the Amaranthaceae family and represents a promising species with the nutritional content and high tolerance to stressful environments, such as soils affected by high salinity. The basic leucine zipper (bZIP) transcription factor represents exclusively in eukaryotes and can be related to many biological processes. So far, the genomes of quinoa and 3 other Amaranthaceae crops (Spinacia oleracea, Beta vulgaris, and Amaranthus hypochondriacus) have been fully sequenced. However, information about the bZIPs in these Amaranthaceae species is limited, and genome-wide analysis of the bZIP family is lacking in quinoa.

RESULTS

We identified 94 bZIPs in quinoa (named as CqbZIP1-CqbZIP94). All the CqbZIPs were phylogenetically splitted into 12 distinct subfamilies. The proportion of CqbZIPs was different in each subfamily, and members within the same subgroup shared conserved exon-intron structures and protein motifs. Besides, 32 duplicated CqbZIP gene pairs were investigated, and the duplicated CqbZIPs had mainly undergone purifying selection pressure, which suggested that the functions of the duplicated CqbZIPs might not diverge much. Moreover, we identified the bZIP members in 3 other Amaranthaceae species, and 41, 32, and 16 orthologous gene pairs were identified between quinoa and S. oleracea, B. vulgaris, and A. hypochondriacus, respectively. Among them, most were a single copy being present in S. oleracea, B. vulgaris, and A. hypochondriacus, and two copies being present in allotetraploid quinoa. The function divergence within the bZIP orthologous genes might be limited. Additionally, 11 selected CqbZIPs had specific spatial expression patterns, and 6 of 11 CqbZIPs were up-regulated in response to salt stress. Among the selected CqbZIPs, 3 of 4 duplicated gene pairs shared similar expression patterns, suggesting that these duplicated genes might retain some essential functions during subsequent evolution.

CONCLUSIONS

The present study provided the first systematic analysis for the phylogenetic classification, motif and gene structure, expansion pattern, and expression profile of the bZIP family in quinoa. Our results would lay an important foundation for functional and evolutionary analysis of CqbZIPs, and provide promising candidate genes for further investigation in tissue specificity and their functional involvement in quinoa's resistance to salt stress.

Evolution of l-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales

The evolution of l-DOPA 4,5-dioxygenase activity, encoded by the gene DODA, was a key step in the origin of betalain biosynthesis in Caryophyllales. We previously proposed that l-DOPA 4,5-dioxygenase activity evolved via a single Caryophyllales-specific neofunctionalisation event within the DODA gene lineage. However, this neofunctionalisation event has not been confirmed and the DODA gene lineage exhibits numerous gene duplication events, whose evolutionary significance is unclear. To address this, we functionally characterised 23 distinct DODA proteins for l-DOPA 4,5-dioxygenase activity, from four betalain-pigmented and five anthocyanin-pigmented species, representing key evolutionary transitions across Caryophyllales. By mapping these functional data to an updated DODA phylogeny, we then explored the evolution of l-DOPA 4,5-dioxygenase activity. We find that low l-DOPA 4,5-dioxygenase activity is distributed across the DODA gene lineage. In this context, repeated gene duplication events within the DODA gene lineage give rise to polyphyletic occurrences of elevated l-DOPA 4,5-dioxygenase activity, accompanied by convergent shifts in key functional residues and distinct genomic patterns of micro-synteny. In the context of an updated organismal phylogeny and newly inferred pigment reconstructions, we argue that repeated convergent acquisition of elevated l-DOPA 4,5-dioxygenase activity is consistent with recurrent specialisation to betalain synthesis in Caryophyllales.

Satellite DNA landscapes after allotetraploidization of quinoa (Chenopodium quinoa) reveal unique A and B subgenomes

If two related plant species hybridize, their genomes may be combined and duplicated within a single nucleus, thereby forming an allotetraploid. How the emerging plant balances two co-evolved genomes is still a matter of ongoing research. Here, we focus on satellite DNA (satDNA), the fastest turn-over sequence class in eukaryotes, aiming to trace its emergence, amplification, and loss during plant speciation and allopolyploidization. As a model, we used Chenopodium quinoa Willd. (quinoa), an allopolyploid crop with 2n = 4x = 36 chromosomes. Quinoa originated by hybridization of an unknown female American Chenopodium diploid (AA genome) with an unknown male Old World diploid species (BB genome), dating back 3.3-6.3 million years. Applying short read clustering to quinoa (AABB), C. pallidicaule (AA), and C. suecicum (BB) whole genome shotgun sequences, we classified their repetitive fractions, and identified and characterized seven satDNA families, together with the 5S rDNA model repeat. We show unequal satDNA amplification (two families) and exclusive occurrence (four families) in the AA and BB diploids by read mapping as well as Southern, genomic, and fluorescent in situ hybridization. Whereas the satDNA distributions support C. suecicum as possible parental species, we were able to exclude C. pallidicaule as progenitor due to unique repeat profiles. Using quinoa long reads and scaffolds, we detected only limited evidence of intergenomic homogenization of satDNA after allopolyploidization, but were able to exclude dispersal of 5S rRNA genes between subgenomes. Our results exemplify the complex route of tandem repeat evolution through Chenopodium speciation and allopolyploidization, and may provide sequence targets for the identification of quinoa's progenitors.

Genomics-assisted breeding in minor and pseudo-cereals

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

Healthy and Resilient Cereals and Pseudo-Cereals for Marginal Agriculture: Molecular Advances for Improving Nutrient Bioavailability

With the ever-increasing world population, an extra 1.5 billion mouths need to be fed by 2050 with continuously dwindling arable land. Hence, it is imperative that extra food come from the marginal lands that are expected to be unsuitable for growing major staple crops under the adverse climate change scenario. Crop diversity provides right alternatives for marginal environments to improve food, feed, and nutritional security. Well-adapted and climate-resilient crops will be the best fit for such a scenario to produce seed and biomass. The minor millets are known for their high nutritional profile and better resilience for several abiotic stresses that make them the suitable crops for arid and salt-affected soils and poor-quality waters. Finger millet (Eleucine coracana) and foxtail millet (Setaria italica), also considered as orphan crops, are highly tolerant grass crop species that grow well in marginal and degraded lands of Africa and Asia with better nutritional profile. Another category of grains, called pseudo-cereals, is considered as rich foods because of their protein quality and content, high mineral content, and healthy and balance food quality. Quinoa (Chenopodium quinoa), amaranth (Amaranthus sp.), and buckwheat (Fagopyrum esculentum) fall under this category. Nevertheless, both minor millets and pseudo-cereals are morphologically different, although similar for micronutrient bioavailability, and their grains are gluten-free. The cultivation of these millets can make dry lands productive and ensure future food as well as nutritional security. Although the natural nutrient profile of these crop plant species is remarkably good, little development has occurred in advances in molecular genetics and breeding efforts to improve the bioavailability of nutrients. Recent advances in NGS have enabled the genome and transcriptome sequencing of these millets and pseudo-cereals for the faster development of molecular markers and application in molecular breeding. Genomic information on finger millet (1,196 Mb with 85,243 genes); S. italica, a model small millet (well-annotated draft genome of 420 Mb with 38,801 protein-coding genes); amaranth (466 Mb genome and 23,059 protein-coding genes); buckwheat (genome size of 1.12 Gb with 35,816 annotated genes); and quinoa (genome size of 1.5 Gb containing 54,438 protein-coding genes) could pave the way for the genetic improvement of these grains. These genomic resources are an important first step toward genetic improvement of these crops. This review highlights the current advances and available resources on genomics to improve nutrient bioavailability in these five suitable crops for the sustained healthy livelihood.

The draft nuclear genome assembly of Eucalyptus pauciflora: a pipeline for comparing de novo assemblies

BACKGROUND

Eucalyptus pauciflora (the snow gum) is a long-lived tree with high economic and ecological importance. Currently, little genomic information for E. pauciflora is available. Here, we sequentially assemble the genome of Eucalyptus pauciflora with different methods, and combine multiple existing and novel approaches to help to select the best genome assembly.

FINDINGS

We generated high coverage of long- (Nanopore, 174×) and short- (Illumina, 228×) read data from a single E. pauciflora individual and compared assemblies from 5 assemblers (Canu, SMARTdenovo, Flye, Marvel, and MaSuRCA) with different read lengths (1 and 35 kb minimum read length). A key component of our approach is to keep a randomly selected collection of ∼10% of both long and short reads separated from the assemblies to use as a validation set for assessing assemblies. Using this validation set along with a range of existing tools, we compared the assemblies in 8 ways: contig N50, BUSCO scores, LAI (long terminal repeat assembly index) scores, assembly ploidy, base-level error rate, CGAL (computing genome assembly likelihoods) scores, structural variation, and genome sequence similarity. Our result showed that MaSuRCA generated the best assembly, which is 594.87 Mb in size, with a contig N50 of 3.23 Mb, and an estimated error rate of ∼0.006 errors per base.

CONCLUSIONS

We report a draft genome of E. pauciflora, which will be a valuable resource for further genomic studies of eucalypts. The approaches for assessing and comparing genomes should help in assessing and choosing among many potential genome assemblies from a single dataset.

RNA-seq Analysis of Salt-Stressed Versus Non Salt-Stressed Transcriptomes of Chenopodium quinoa Landrace R49

Quinoa (Chenopodium quinoa Willd.), a model halophytic crop species, was used to shed light on salt tolerance mechanisms at the transcriptomic level. An RNA-sequencing analysis of genotype R49 at an early vegetative stage was performed by Illumina paired-ends method comparing high salinity and control conditions in a time-course pot experiment. Genome-wide transcriptional salt-induced changes and expression profiling of relevant salt-responsive genes in plants treated or not with 300 mM NaCl were analyzed after 1 h and 5 days. We obtained up to 49 million pairs of short reads with an average length of 101 bp, identifying a total of 2416 differentially expressed genes (DEGs) based on the treatment and time of sampling. In salt-treated vs. control plants, the total number of up-regulated and down-regulated genes was 945 and 1471, respectively. The number of DEGs was higher at 5 days than at 1 h after salt treatment, as reflected in the number of transcription factors, which increased with time. We report a strong transcriptional reprogramming of genes involved in biological processes like oxidation-reduction, response to stress and response to abscisic acid (ABA), and cell wall organization. Transcript analyses by real-time RT- qPCR supported the RNA-seq results and shed light on the contribution of roots and shoots to the overall transcriptional response. In addition, it revealed a time-dependent response in the expression of the analyzed DEGs, including a quick (within 1 h) response for some genes, suggesting a "stress-anticipatory preparedness" in this highly salt-tolerant genotype.

Shotgun proteomic analysis of quinoa seeds reveals novel lysine-rich seed storage globulins

Quinoa seeds have high protein content and an exceptional balance of amino acids, with higher contents of lysine, methionine and cysteine than common cereals. To date, only three globulins, all of which have a content of lysine mass that does not exceed 3.8%, have been identified in quinoa. To address the protein present in quinoa seeds, TCA/Acetone protein extraction was performed using four different quinoa seed genotypes with contrasting edaphoclimatic origins. Proteins were identified and analyzed using label-free shotgun proteomics followed by in silico analysis, using the three published quinoa genomes. This analysis allowed us to identify sixteen globulins, thirteen of which are novel: nine legumin-like proteins and seven vicilin-like proteins. Seven of the novel proteins contain 7.5% or more of lysine mass, justifying the high content of lysine repeatedly reported in quinoa seeds. No significant differences were found between the four genotypes here analyzed.

Genome-Wide Identification, Characterization, and Expression Analysis of the NAC Transcription Factor in Chenopodium quinoa

The NAC (NAM, ATAF, and CUC) family is one of the largest families of plant-specific transcription factors. It is involved in many plant growth and development processes, as well as abiotic/biotic stress responses. So far, little is known about the NAC family in Chenopodium quinoa. In the present study, a total of 90 NACs were identified in quinoa (named as CqNAC1-CqNAC90) and phylogenetically divided into 14 distinct subfamilies. Different subfamilies showed diversities in gene proportions, exon-intron structures, and motif compositions. In addition, 28 CqNAC duplication events were investigated, and a strong subfamily preference was found during the NAC expansion in quinoa, indicating that the duplication event was not random across NAC subfamilies during quinoa evolution. Moreover, the analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios suggested that the duplicated CqNACs might have mainly experienced purifying selection pressure with limited functional divergence. Additionally, 11 selected CqNACs showed significant tissue-specific expression patterns, and all the CqNACs were positively regulated in response to salt stress. The result provided evidence for selecting candidate genes for further characterization in tissue/organ specificity and their functional involvement in quinoa's strong salinity tolerance.

Isolation of amaranthin synthetase from Chenopodium quinoa and construction of an amaranthin production system using suspension-cultured tobacco BY-2 cells

Betalains are plant pigments primarily produced by plants of the order Caryophyllales. Because betalain possesses anti-inflammatory and anticancer activities, it may be useful as a pharmaceutical agent and dietary supplement. Recent studies have identified the genes involved in the betalain biosynthesis of betanin. Amaranthin and celosianin II are abundant in the quinoa (Chenopodium quinoa Willd.) hypocotyl, and amaranthin comprises glucuronic acid bound to betanin; therefore, this suggests the existence of a glucuronyltransferase involved in the synthesis of amaranthin in the quinoa hypocotyl. To identify the gene involved in amaranthin biosynthesis, we performed a BLAST analysis and phylogenetic tree analysis based on sequences homologous to flavonoid glycosyltransferase, followed by expression analysis on the quinoa hypocotyl to obtain three candidate proteins. Production of amaranthin in a transient Nicotiana benthamiana expression system was evaluated for these candidates and one was identified as having the ability to produce amaranthin. The gene encoding this protein was quinoa amaranthin synthetase 1 (CqAmaSy1). We also created a transgenic tobacco bright yellow-2 (BY-2) cell line wherein four betalain biosynthesis genes were introduced to facilitate amaranthin production. This transgenic cell line produced 13.67 ± 4.13 μm (mean ± SEM) amaranthin and 26.60 ± 1.53 μm betanin, whereas the production of isoamaranthin and isobetanin could not be detected. Tests confirmed the ability of amaranthin and betanin to slightly suppress cancer cell viability. Furthermore, amaranthin was shown to significantly inhibit HIV-1 protease activity, whereas betanin did not.

Primula vulgaris (primrose) genome assembly, annotation and gene expression, with comparative genomics on the heterostyly supergene

Primula vulgaris (primrose) exhibits heterostyly: plants produce self-incompatible pin- or thrum-form flowers, with anthers and stigma at reciprocal heights. Darwin concluded that this arrangement promotes insect-mediated cross-pollination; later studies revealed control by a cluster of genes, or supergene, known as the S (Style length) locus. The P. vulgaris S locus is absent from pin plants and hemizygous in thrum plants (thrum-specific); mutation of S locus genes produces self-fertile homostyle flowers with anthers and stigma at equal heights. Here, we present a 411 Mb P. vulgaris genome assembly of a homozygous inbred long homostyle, representing ~87% of the genome. We annotate over 24,000 P. vulgaris genes, and reveal more genes up-regulated in thrum than pin flowers. We show reduced genomic read coverage across the S locus in other Primula species, including P. veris, where we define the conserved structure and expression of the S locus genes in thrum. Further analysis reveals the S locus has elevated repeat content (64%) compared to the wider genome (37%). Our studies suggest conservation of S locus genetic architecture in Primula, and provide a platform for identification and evolutionary analysis of the S locus and downstream targets that regulate heterostyly in diverse heterostylous species.

Developing naturally stress-resistant crops for a sustainable agriculture

A major problem facing humanity is that our numbers are growing but the availability of land and fresh water for agriculture is not. This problem is being exacerbated by climate change-induced increases in drought, and other abiotic stresses. Stress-resistant crops are needed to ensure yield stability under stress conditions and to minimize the environmental impacts of crop production. Evolution has created thousands of species of naturally stress-resistant plants (NSRPs), some of which have already been subjected to human domestication and are considered minor crops. Broader cultivation of these minor crops will diversify plant agriculture and the human diet, and will therefore help improve global food security and human health. More research should be directed toward understanding and utilizing NSRPs. Technologies are now available that will enable researchers to rapidly improve the genetics of NSRPs, with the goal of increasing NSRP productivity while retaining NSRP stress resistance and nutritional value.

Quinoa Abiotic Stress Responses: A Review

Quinoa (Chenopodium quinoa Willd.) is a genetically diverse Andean crop that has earned special attention worldwide due to its nutritional and health benefits and its ability to adapt to contrasting environments, including nutrient-poor and saline soils and drought stressed marginal agroecosystems. Drought and salinity are the abiotic stresses most studied in quinoa; however, studies of other important stress factors, such as heat, cold, heavy metals, and UV-B light irradiance, are severely limited. In the last few decades, the incidence of abiotic stress has been accentuated by the increase in unpredictable weather patterns. Furthermore, stresses habitually occur as combinations of two or more. The goals of this review are to: (1) provide an in-depth description of the existing knowledge of quinoa's tolerance to different abiotic stressors; (2) summarize quinoa's physiological responses to these stressors; and (3) describe novel advances in molecular tools that can aid our understanding of the mechanisms underlying quinoa's abiotic stress tolerance.

Understanding the Molecular Basis of Salt Sequestration in Epidermal Bladder Cells of Chenopodium quinoa

Soil salinity is destroying arable land and is considered to be one of the major threats to global food security in the 21st century. Therefore, the ability of naturally salt-tolerant halophyte plants to sequester large quantities of salt in external structures, such as epidermal bladder cells (EBCs), is of great interest. Using Chenopodium quinoa, a pseudo-cereal halophyte of great economic potential, we have shown previously that, upon removal of salt bladders, quinoa becomes salt sensitive. In this work, we analyzed the molecular mechanism underlying the unique salt dumping capabilities of bladder cells in quinoa. The transporters differentially expressed in the EBC transcriptome and functional electrophysiological testing of key EBC transporters in Xenopus oocytes revealed that loading of Na⁺ and Cl^- into EBCs is mediated by a set of tailored plasma and vacuole membrane-based sodium-selective channel and chloride-permeable transporter.

Isolation and characterization of the betalain biosynthesis gene involved in hypocotyl pigmentation of the allotetraploid Chenopodium quinoa

In quinoa seedlings, the pigment betalain accumulates in the hypocotyl. To isolate the genes involved in betalain biosynthesis in the hypocotyl, we performed ethyl methanesulfonate (EMS) mutagenesis on the CQ127 variety of quinoa seedlings. While putative amaranthin and celosianin II primarily accumulate in the hypocotyls, this process produced a green hypocotyl mutant (ghy). This MutMap+ method using the quinoa draft genome revealed that the causative gene of the mutant is CqCYP76AD1-1. Our results indicated that the expression of CqCYP76AD1-1 was light-dependent. In addition, the transient expression of CqCYP76AD1-1 in Nicotiana benthamiana leaves resulted in the accumulation of betanin but not isobetanin, and the presence of a polymorphism in CqCYP76A1-2 in the CQ127 variety was shown to have resulted in its loss of function. These findings suggested that CqCYP76AD1-1 is involved in betalain biosynthesis during the hypocotyl pigmentation process in quinoa. To our knowledge, CqCYP76AD1-1 is the first quinoa gene identified by EMS mutagenesis using a draft gene sequence.

Influence of natural and synthetic vitamin C (ascorbic acid) on primary and secondary metabolites and associated metabolism in quinoa (Chenopodium quinoa Willd.) plants under water deficit regimes

Phytoextracts are being widely used these days as a source of bioactive compounds for mitigating the harmful effects of abiotic stresses including drought stress. In this study, it was assessed how far foliar applied pure synthetic ascorbic acid (AsA) or natural sweet orange juice (OJ) enriched with AsA could mitigate the drought stress induced adverse effects on growth and some key metabolic processes in quinoa (Chenopodium quinoa Willd.; cultivar V₉) plants. Two weeks old quinoa seedlings were subjected to varying irrigation regimes as control [100% field capacity (FC)] and drought stress (60% FC, 40% FC and 20% FC). After one month of water deficit treatments, various levels of ascorbic acid (150 mg L^-1 AsA or 25% OJ) besides control [distilled water (DW) and no spray (NS)] were applied as a foliar spray. After 15 days of AsA application, different physio-biochemical attributes were measured. The results showed that water deficit markedly decreased plant growth, relative water content (RWC), photosynthetic rate, total carotenoids (CAR) and total flavonoids, while it increased relative membrane permeability (RMP), intrinsic AsA content, hydrogen peroxide (H₂O₂), malondialdehyde (MDA), glycinebetaine (GB), total phenolics, total soluble proteins (TSP), total free amino acids, activities of key antioxidant enzymes [superoxide dismutase (SOD), peroxidase (POD)], total soluble sugars (TSS), reducing (RS) and non-reducing sugars (NRS). Most obvious results of most of these parameters were observed at 40% and 20% FC. Foliar-applied pure 150 mg L^-1 AsA and 25% OJ were found to be very effective in improving plant growth, RMP, photosynthetic rate, CAR, proline, AsA, MDA, GB, TSP, free amino acids, SOD, POD, TSS, RS, NRS and total flavonoids. It was noticed that 25% OJ enriched with AsA and other essential nutrients and biomolecules was as efficient as 150 mg L^-1 AsA in reducing the adverse effects of drought stress on quinoa plants. So, it was concluded that OJ, a cheaper source of vitamin C, can be used as a mitigating agent for improving drought tolerance in plants under drought-prone environments.

Molecular dynamics in germinating, endophyte-colonized quinoa seeds

AIMS

The pseudo-cereal quinoa has an outstanding nutritional value. Seed germination is unusually fast, and plant tolerance to salt stress exceptionally high. Seemingly all seeds harbor bacterial endophytes. This work examines mitogen-activated protein kinase (MAPK) activities during early development. It evaluates possible contribution of endophytes to rapid germination and plant robustness.

METHODS

MAPK activities were monitored in water- and NaCl-imbibed seeds over a 4-h-period using an immunoblot-based approach. Cellulolytic and pectinolytic abilities of bacteria were assessed biochemically, and cellular movement, biofilm, elicitor and antimicrobial compound synthesis genes sequenced. GyrA-based, cultivation-independent studies provided first insight into endophyte diversity.

RESULTS

Quinoa seeds and seedlings exhibit remarkably complex and dynamic MAPK activity profiles. Depending on seed origin, variances exist in MAPK patterns and probably also in endophyte assemblages. Mucilage-degrading activities enable endophytes to colonize seed surfaces of a non-host species, chia, without apparent adverse effects.

CONCLUSIONS

Owing to their motility, cell wall-loosening and elicitor-generating abilities, quinoa endophytes have the potential to drive cell expansion, move across cell walls, generate damage-associated molecular patterns and activate MAPKs in their host. Bacteria may thus facilitate rapid germination and confer a primed state directly upon seed rehydration. Transfer into non-native crops appears both desirable and feasible.