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

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.

Exploring the complexity of genome size reduction in angiosperms

The genome sizes of angiosperms decreased significantly more than the genome sizes of their ancestors (pteridophytes and gymnosperms). Decreases in genome size involve a highly complex process, with remnants of the genome size reduction scattered across the genome and not directly linked to specific genomic structures. This is because the associated mechanisms operate on a much smaller scale than the mechanisms mediating increases in genome size. This review thoroughly summarizes the available literature regarding the molecular mechanisms underlying genome size reductions and introduces Utricularia gibba and Arabidopsis thaliana as model species for the examination of the effects of these molecular mechanisms. Additionally, we propose that phosphorus deficiency and drought stress are the major external factors contributing to decreases in genome size. Considering these factors affect almost all land plants, angiosperms likely gained the mechanisms for genome size reductions. These environmental factors may affect the retention rates of deletions, while also influencing the mutation rates of deletions via the functional diversification of the proteins facilitating double-strand break repair. The biased retention and mutation rates of deletions may have synergistic effects that enhance deletions in intergenic regions, introns, transposable elements, duplicates, and repeats, leading to a rapid decrease in genome size. We suggest that these selection pressures and associated molecular mechanisms may drive key changes in angiosperms during recurrent cycles of genome size decreases and increases.

A low-cost dpMIG-seq method for elucidating complex inheritance in polysomic crops: a case study in tetraploid blueberry

Next-generation sequencing (NGS) library construction often requires high-quality DNA extraction, precise adjustment of DNA concentration, and restriction enzyme digestion to reduce genome complexity, which results in increased time and cost in sample preparation and processing. To address these challenges, a PCR-based method for rapid NGS library preparation, named dpMIG-seq, has been developed and proven effective for high-throughput genotyping. However, the application of dpMIG-seq has been limited to diploid and polyploid species with disomic inheritance. In this study, we obtained genome-wide single nucleotide polymorphism (SNP) markers for tetraploid blueberry to evaluate genotyping and downstream analysis outcomes. Comparison of genotyping qualities inferred across samples with different DNA concentrations and multiple bioinformatics approaches revealed high accuracy and reproducibility of dpMIG-seq-based genotyping, with Pearson's correlation coefficients between replicates in the range of 0.91 to 0.98. Furthermore, we demonstrated that dpMIG-seq enables accurate genotyping of samples with low DNA concentrations. Subsequently, we applied dpMIG-seq to a tetraploid F1 population to examine the inheritance probability of parental alleles. Pairing configuration analysis supported the random meiotic pairing of homologous chromosomes on a genome-wide level. On the other hand, preferential pairing was observed on chr-11, suggesting that there may be an exception to the random pairing. Genotypic data suggested quadrivalent formation within the population, although the frequency of quadrivalent formation varied by chromosome and cultivar. Collectively, the results confirmed applicability of dpMIG-seq for allele dosage genotyping and are expected to catalyze the adoption of this cost-effective and rapid genotyping technology in polyploid studies.

Chromosome-level genome assemblies for two quinoa inbred lines from northern and southern highlands of Altiplano where quinoa originated

Quinoa is emerging as a key seed crop for global food security due to its ability to grow in marginal environments and its excellent nutritional properties. Because quinoa is partially allogamous, we have developed quinoa inbred lines necessary for molecular genetic analysis. Our comprehensive genomic analysis showed that the quinoa inbred lines fall into three genetic subpopulations: northern highland, southern highland, and lowland. Lowland and highland quinoa are the same species, but have very different genotypes and phenotypes. Lowland quinoa has relatively small grains and a darker grain color, and is widely tested and grown around the world. In contrast, the white, large-grained highland quinoa is grown in the Andean highlands, including the region where quinoa originated, and is exported worldwide as high-quality quinoa. Recently, we have shown that viral vectors can be used to regulate endogenous genes in quinoa, paving the way for functional genomics to reveal the diversity of quinoa. However, although a high-quality assembly has recently been reported for a lowland quinoa line, genomic resources of the quality required for functional genomics are not available for highland quinoa lines. Here we present high-quality chromosome-level genome assemblies for two highland inbred quinoa lines, J075 representing the northern highland line and J100 representing the southern highland line, using PacBio HiFi sequencing and dpMIG-seq. In addition, we demonstrate the importance of verifying and correcting reference-based scaffold assembly with other approaches such as linkage maps. The assembled genome sizes of J075 and J100 are 1.29 and 1.32 Gb, with contigs N50 of 66.3 and 12.6 Mb, and scaffold N50 of 71.2 and 70.6 Mb, respectively, comprising 18 pseudochromosomes. The repetitive sequences of J075 and J100 represent 72.6% and 71.5% of the genome, the majority of which are long terminal repeats, representing 44.0% and 42.7% of the genome, respectively. The de novo assembled genomes of J075 and J100 were predicted to contain 65,303 and 64,945 protein-coding genes, respectively. The high quality genomes of these highland quinoa lines will facilitate quinoa functional genomics research on quinoa and contribute to the identification of key genes involved in environmental adaptation and quinoa domestication.

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.

Association study of morpho-phenological traits in quinoa (Chenopodium quinoa Willd.) using SSR markers

In this study, the genetic and molecular diversity of 60 quinoa accessions was assessed using agronomically important traits related to grain yield as well as microsatellite (SSR) markers, and informative markers linked to the studied traits were identified using association study. The results showed that most of the studied traits had a relatively high diversity, but grain saponin and protein content showed the highest diversity. High diversity was also observed in all SSR markers, but KAAT023, KAAT027, KAAT036, and KCAA014 showed the highest values for most of the diversity indices and can be introduced as the informative markers to assess genetic diversity in quinoa. Population structure analysis showed that the studied population probably includes two subclusters, so that out of 60 quinoa accessions, 29 (48%) and 23 (38%) accessions were assigned to the first and second subclusters, respectively, and eight (13%) accessions were considered as the mixed genotypes. The study of the population structure using Structure software showed two possible subgroups (K = 2) in the studied population and the results of the bar plot confirmed it. Association study using the general linear model (GLM) and mixed linear model (MLM) identified the number of 35 and 32 significant marker-trait associations (MTAs) for the first year (2019) and 37 and 35 significant MTAs for the second year (2020), respectively. Among the significant MTAs identified for different traits, the highest number of significant MTAs were obtained for grain yield and 1000-grain weight with six and five MTAs, respectively.

Upregulation of tandem duplicated BoFLC1 genes is associated with the non-flowering trait in Brassica oleracea var. capitata

Tandem duplicated BoFLC1 genes (BoFLC1a and BoFLC1b), which were identified as the candidate causal genes for the non-flowering trait in the cabbage mutant 'nfc', were upregulated during winter in 'nfc'. The non-flowering natural cabbage mutant 'nfc' was discovered from the breeding line 'T15' with normal flowering characteristics. In this study, we investigated the molecular basis underlying the non-flowering trait of 'nfc'. First, 'nfc' was induced to flower using the grafting floral induction method, and three F₂ populations were generated. The flowering phenotype of each F₂ population was widely distributed with non-flowering individuals appearing in two populations. QTL-seq analysis detected a genomic region associated with flowering date at approximately 51 Mb on chromosome 9 in two of the three F₂ populations. Subsequent validation and fine mapping of the candidate genomic region using QTL analysis identified the quantitative trait loci (QTL) at 50,177,696-51,474,818 bp on chromosome 9 covering 241 genes. Additionally, RNA-seq analysis in leaves and shoot apices of 'nfc' and 'T15' plants identified 19 and 15 differentially expressed genes related to flowering time, respectively. Based on these results, we identified tandem duplicated BoFLC1 genes, which are homologs of the floral repressor FLOWERING LOCUS C, as the candidate genes responsible for the non-flowering trait of 'nfc'. We designated the tandem duplicated BoFLC1 genes as BoFLC1a and BoFLC1b. Expression analysis revealed that the expression levels of BoFLC1a and BoFLC1b were downregulated during winter in 'T15' but were upregulated and maintained during winter in 'nfc'. Additionally, the expression level of the floral integrator BoFT was upregulated in the spring in 'T15' but hardly upregulated in 'nfc'. These results suggest that the upregulated levels of BoFLC1a and BoFLC1b contributed to the non-flowering trait of 'nfc'.

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.

Proximal and remote sensing in plant phenomics: 20 years of progress, challenges, and perspectives

Plant phenomics (PP) has been recognized as a bottleneck in studying the interactions of genomics and environment on plants, limiting the progress of smart breeding and precise cultivation. High-throughput plant phenotyping is challenging owing to the spatio-temporal dynamics of traits. Proximal and remote sensing (PRS) techniques are increasingly used for plant phenotyping because of their advantages in multi-dimensional data acquisition and analysis. Substantial progress of PRS applications in PP has been observed over the last two decades and is analyzed here from an interdisciplinary perspective based on 2972 publications. This progress covers most aspects of PRS application in PP, including patterns of global spatial distribution and temporal dynamics, specific PRS technologies, phenotypic research fields, working environments, species, and traits. Subsequently, we demonstrate how to link PRS to multi-omics studies, including how to achieve multi-dimensional PRS data acquisition and processing, how to systematically integrate all kinds of phenotypic information and derive phenotypic knowledge with biological significance, and how to link PP to multi-omics association analysis. Finally, we identify three future perspectives for PRS-based PP: (1) strengthening the spatial and temporal consistency of PRS data, (2) exploring novel phenotypic traits, and (3) facilitating multi-omics communication.

Photosynthesis is not the unique useful trait for discriminating salt tolerance capacity between sensitive and tolerant quinoa varieties

Growth was not strictly linked to photosynthesis performance under salinity conditions in quinoa. Other key traits, which were varieties-specific, rather than photosynthesis explained better growth performance. Phenotyping for salinity stress tolerance in quinoa is of great interest to select traits contributing to overall salinity tolerance and to understand the response mechanisms to salinity at a whole plant level. The objective of this work was to dissect the responses of specific traits and analyse relations between these traits to better understand growth response under salinity conditions in quinoa. Growth response to salinity was mostly related to differences in basal values of biomass, being reduced the most in plants with higher basal biomass. Regarding the relationship between growth and specific traits, in Puno variety, better photosynthetic performance was related to a better maintenance of growth. Nevertheless, in the rest of the varieties other traits rather than photosynthesis could better explain growth response. In this way, the development of succulence in F-16 and Collana varieties, also the osmotic adjustment but in smaller dimensions in Pasankalla, Marisma and S-15-15 helped to maintain better growth. Besides, smaller increases of Cl^- could have caused a limited nitrate uptake reducing more growth in Vikinga. Ascorbate was considered a key trait as a noticeable fall of it was also related to higher reductions in growth in Titicaca. These results suggest that, due to the genetic variability of quinoa and the complexity of salinity tolerance, no unique and specific traits should be taken into consideration when using phenotyping for analysing salinity tolerance in quinoa.

Assessment of Phenotypic Diversity in the USDA Collection of Quinoa Links Genotypic Adaptation to Germplasm Origin

Quinoa’s germplasm evaluation is the first step towards determining its suitability under new environmental conditions. The aim of this study was to introduce suitable germplasm to the lowland areas of the Faisalabad Plain that could then be used to introduce quinoa more effectively to that region. A set of 117 quinoa genotypes belonging to the USDA quinoa collection was evaluated for 11 phenotypic quantitative traits (grain yield (Y), its biological and numerical components plus phenological variables) in a RCBD during two consecutive growing seasons at the University of Agriculture, Faisalabad, Pakistan under mid-autumn sowings. Genotypic performance changed across the years, however most phenotypic traits showed high heritability, from 0.75 for Harvest Index (HI) to 0.97 for aerial biomass (B) and Y. Ordination and cluster analyses differentiated four groups dominated by genotypes from: Peru and the Bolivian Highlands (G1); the Bolivian Highlands (G2); the Ballón collection (regarded as a cross between Bolivian and Sea Level (Chilean) genotypes) plus Bolivian Highlands (G3); and Ballón plus Sea Level (G4), this latter group being the most differentiated one. This genetic structure shared similarities with previous groups identified using SSR markers and G×E data from an international quinoa test. G4 genotypes showed the highest Y associated with higher B and seed numbers (SN), while HI made a significant contribution to yield determination in G2 and seed weight (SW) in G3. G1 and G2 showed the lowest Y associated with a lower B and SN. Moreover, SW showed a strongly negative association with SN in G2. Accordingly, G4 followed by G3 are better suited to the lowland areas of Faisalabad plain and the physiological traits underlying yield determination among genotypic groups should be considered in future breeding programs.

MIG-seq is an effective method for high-throughput genotyping in wheat (Triticum spp.)

MIG-seq (Multiplexed inter-simple sequence repeats genotyping by sequencing) has been developed as a low cost genotyping technology, although the number of polymorphisms obtained is assumed to be minimal, resulting in the low application of this technique to analyses of agricultural plants. We applied MIG-seq to 12 plant species that include various crops and investigated the relationship between genome size and the number of bases that can be stably sequenced. The genome size and the number of loci, which can be sequenced by MIG-seq, are positively correlated. This is due to the linkage between genome size and the number of simple sequence repeats (SSRs) through the genome. The applicability of MIG-seq to population structure analysis, linkage mapping, and quantitative trait loci (QTL) analysis in wheat, which has a relatively large genome, was further evaluated. The results of population structure analysis for tetraploid wheat showed the differences among collection sites and subspecies, which agreed with previous findings. Additionally, in wheat biparental mapping populations, over 3,000 SNPs/indels with low deficiency were detected using MIG-seq, and the QTL analysis was able to detect recognized flowering-related genes. These results revealed the effectiveness of MIG-seq for genomic analysis of agricultural plants with large genomes, including wheat.

R. Patiranage D, Johansen K, Schmöckel SM, Bertero D, Oakey H, Colque-Little C, Afzal I, Raubach S, Miller N, Streich J, Amby DB, Emrani N, Warmington M, Mousa MAA, Wu D, Jacobson D, Andreasen C, Jung C, Murphy K, Bazile D, Tester M. Quinoa Phenotyping Methodologies: An International Consensus

Quinoa is a crop originating in the Andes but grown more widely and with the genetic potential for significant further expansion. Due to the phenotypic plasticity of quinoa, varieties need to be assessed across years and multiple locations. To improve comparability among field trials across the globe and to facilitate collaborations, components of the trials need to be kept consistent, including the type and methods of data collected. Here, an internationally open-access framework for phenotyping a wide range of quinoa features is proposed to facilitate the systematic agronomic, physiological and genetic characterization of quinoa for crop adaptation and improvement. Mature plant phenotyping is a central aspect of this paper, including detailed descriptions and the provision of phenotyping cards to facilitate consistency in data collection. High-throughput methods for multi-temporal phenotyping based on remote sensing technologies are described. Tools for higher-throughput post-harvest phenotyping of seeds are presented. A guideline for approaching quinoa field trials including the collection of environmental data and designing layouts with statistical robustness is suggested. To move towards developing resources for quinoa in line with major cereal crops, a database was created. The Quinoa Germinate Platform will serve as a central repository of data for quinoa researchers globally.

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.

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.