Quinoa Genome Research Articles

Genome-Wide Analysis of GATA Transcription Factor Family in Quinoa (Chenopodium quinoa): Identification, Characterization, and Expression Profiles

GATA-binding factor (GATA) proteins are the transcription factor (TF) family that are commonly involved in plant growth and development. The GATA TF family has been successfully identified and characterized for various higher plant species, but there is little research on the GATA TF family in quinoa (Chenopodium quinoa). In this present study, a total of 32 CqGATA genes were identified and analyzed in the quinoa genome. While the general features of the CqGATA TFs in quinoa were slightly variable, the majority of genes encoding the CqGATA TFs contained two and three exons. Our phylogenetic analysis demonstrated that the CqGATA TFs could be classified into four different groups. Gene expression analysis indicated that the expression profiles of the CqGATA genes varied in different tissues. Overall, our study could provide a reference for further functional characterization of the CqGATA genes in quinoa.

Genomic variation induced by a low concentration of ethyl methanesulfonate (EMS) in quinoa ‘Longli-4’ variety

Quinoa (Chenopodium quinoa, 2n = 4x = 36), a super pseudocereal crop, has been introduced into China nearly 60 years. Many excellent varieties have been developed through massive selection; however, few are developed through mutagenesis breeding. In this study, the ‘Longli-4’ variety, locally cultivated in Gansu province, Northwest China, was selected for experimentation. The grains of ‘Longli-4’ were treated with ethyl methanesulfonate (EMS) at a concentration of 0.8% for 8 h. Nine plants from independent M2 families were randomly selected to investigate the mutagenesis effect of EMS on the quinoa genome. The results indicated that the single nucleotide polymorphisms (SNPs) induced by EMS were unevenly distributed across all 18 chromosomes, with an average mutation frequency of 91.2 SNPs/Mb, ranging from 4.5 to 203.5 SNPs/Mb. A significant positive correlation between the number of SNPs and chromosome length was identified through linear model analysis. Transitions from G/C to A/T were the most predominated in all variant categories, accounting for 34.4–67.2% of the mutations, and SNPs were significantly enriched in intergenic regions, representing 69.2–75.1% of the total mutations. This study provides empirical support for the application of low concentration EMS treatment in quinoa breeding.

Genome-wide identification and salt stress-expression analysis of the dehydrin gene family in Chenopodium quinoa

Dehydrins (DHNs) are essential proteins in the embryonic development and abiotic stress responses of plants. Due to their remarkable ability to confer tolerance to plants in conditions of drought, salinity and extreme temperatures, DHNs have garnered considerable interest. Quinoa (Chenopodium quinoa Willd.), a facultative halophyte plant, can thrive in a wide range of agroecosystems, making it a promising candidate for stress tolerance studies. In this study, we identified eleven DHN genes in the quinoa genome belonging to Y-, F- and H-orthologous groups found in angiosperms. Notably, the H-DHNs lack the K-segment, a feature observed in all Amaranthaceae species, but not in other angiosperms. We identified four DHN structural subgroups: FSKn, YnSKn, SKn-DHNs and the atypical HS-DHN. Phylogenetic analysis indicated that each structural subgroup, except for SK2-DHN, presents two paralogous genes, in accordance with the allotetraploid character of C. quinoa. Quantitative real-time PCR expression analysis revealed that DHN1s (FSK2) and DHN3s (Y2SK2) were expressed in all tissues, while DHN2s (FSK3) were predominant in roots and DHN4s (Y4SK2 and SK2) were predominant in flowers. Salt-response gene expression analysis in seedlings showed that CqDHN4s increase their expression in response to salt stress in all varieties studied, while CqDHN1s reduce their expression in a more salt stress-tolerant variety, suggesting a possible adaptive advantage. In silico analysis of the promoters of CqDHN1s and CqDHN4s supports the involvement of these DHNs in responding to abiotic stress.

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

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

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.

Genome-Wide Identification and Expression Analysis of 1-Aminocyclopropane-1-Carboxylate Synthase (ACS) Gene Family in Chenopodium quinoa

Ethylene plays an important role in plant development and stress resistance. The rate-limiting enzyme in ethylene biosynthesis is 1-aminocyclopropane-1-carboxylic acid synthase (ACS). C. quinoa (Chenopodium quinoa) is an important food crop known for its strong tolerance to abiotic stresses. However, knowledge regarding the ACS gene family in C. quinoa remains restricted. In this study, we successfully identified 12 ACS genes (CqACSs) from the C. quinoa genome. Through thorough analysis of their sequences and phylogenetic relationships, it was verified that 8 out of these 12 CqACS isozymes exhibited substantial resemblance to ACS isozymes possessing ACS activity. Furthermore, these eight isozymes could be categorized into three distinct groups. The four remaining CqACS genes grouped under category IV displayed notable similarities with AtACS10 and AtACS12, known as amido transferases lacking ACS activity. The CqACS proteins bore resemblance to the AtACS proteins and had the characteristic structural features typically observed in plant ACS enzymes. Twelve CqACS genes were distributed across 8 out of the 18 chromosomes of C. quinoa. The CqACS genes were expanded from segment duplication. Many cis-regulatory elements related with various abiotic stresses, phytohormones, and light were found. The expression patterns of ACS genes varied across different tissues of C. quinoa. Furthermore, the analysis of gene expression patterns under abiotic stress showed that CqACS genes can be responsive to various stresses, implying their potential functions in adapting to various abiotic stresses. The findings from this research serve as a foundation for delving deeper into the functional roles of CqACS genes.

Molecular Characterization and Expression Analysis of YABBY Genes in Chenopodium quinoa.

Plant-specific YABBY transcription factors play an important role in lateral organ development and abiotic stress responses. However, the functions of the YABBY genes in quinoa remain elusive. In this study, twelve YABBY (CqYAB) genes were identified in the quinoa genome, and they were distributed on nine chromosomes. They were classified into FIL/YAB3, YAB2, YAB5, INO, and CRC clades. All CqYAB genes consist of six or seven exons, and their proteins contain both N-terminal C2C2 zinc finger motifs and C-terminal YABBY domains. Ninety-three cis-regulatory elements were revealed in CqYAB gene promoters, and they were divided into six groups, such as cis-elements involved in light response, hormone response, development, and stress response. Six CqYAB genes were significantly upregulated by salt stress, while one was downregulated. Nine CqYAB genes were upregulated under drought stress, whereas six CqYAB genes were downregulated under cadmium treatment. Tissue expression profiles showed that nine CqYAB genes were expressed in seedlings, leaves, and flowers, seven in seeds, and two specifically in flowers, but no CqYAB expression was detected in roots. Furthermore, CqYAB4 could rescue the ino mutant phenotype in Arabidopsis but not CqYAB10, a paralog of CqYAB4, indicative of functional conservation and divergence among these YABBY genes. Taken together, these results lay a foundation for further functional analysis of CqYAB genes in quinoa growth, development, and abiotic stress responses.

Identification of CDK gene family and functional analysis of CqCDK15 under drought and salt stress in quinoa

as one of the oldest cultivated crops in the world, quinoa has been widely valued for its rich nutritional value and green health. In this study, 22 CDK genes (CqCDK01-CqCDK22) were identified from quinoa genome using bioinformatics method. The number of amino acids was 173–811, the molecular weight was 19,554.89 Da-91,375.70 Da, and the isoelectric point was 4.57–9.77. The phylogenetic tree divided 21 CqCDK genes into six subfamilies, the gene structure showed that 12 (54.5%) CqCDK genes (CqCDK03, CqCDK04, CqCDK05, CqCDK06, CqCDK07, CqCDK11, CqCDK14, CqCDK16, CqCDK18, CqCDK19, CqCDK20 and CqCDK21) had UTR regions at 5’ and 3’ ends. Each CDK protein had different motifs (3–9 motifs), but the genes with the same motifs were located in the same branch. Promoter analysis revealed 41 cis-regulatory elements related to plant hormones, abiotic stresses, tissue-specific expression and photoresponse. The results of real-time fluorescence quantitative analysis showed that the expression level of some CDK genes was higher under drought and salt stress, which indicated that CDK genes could help plants to resist adverse environmental effects. Subcellular localization showed that CqCDK15 gene was localized to the nucleus and cytoplasm, and transgenic plants overexpressing CqCDK15 gene showed higher drought and salt tolerance compared to the controls. Therefore, CDK genes are closely related to quinoa stress resistance. In this study, the main functions of quinoa CDK gene family and its expression level in different tissues and organs were analyzed in detail, which provided some theoretical support for quinoa stress-resistant breeding. Meanwhile, this study has important implications for further understanding the function of the CDK gene family in quinoa and our understanding of the CDK family in vascular plant.

Characterization of the CqCAMTA gene family reveals the role of CqCAMTA03 in drought tolerance

BackgroundCalmodulin-binding transcription activators (CAMTAs) are relatively conserved calmodulin-binding transcription factors widely found in eukaryotes and play important roles in plant growth and stress response. CAMTA transcription factors have been identified in several plant species, but the family members and functions have not yet been identified and analyzed in quinoa.ResultsIn this study, we identified seven CAMTA genes across the whole quinoa genome and analyzed the expression patterns of CqCAMTAs in root and leaf tissues. Gene structure, protein domain, and phylogenetic analyses showed that the quinoa CAMTAs were structurally similar and clustered into the same three major groups as other plant CAMTAs. A large number of stress response-related cis-elements existed in the 2 kb promoter region upstream of the transcription start site of the CqCAMTA genes. qRT-PCR indicated that CqCAMTA genes were expressed differentially under PEG treatments in leaves, and responded to drought stress in leaves and roots. In particular, the CqCAMTA03 gene strongly responded to drought. The transient expression of CqCAMTA03-GFP fusion protein in the tobacco leaf showed that CqCAMTA03 was localized in the nucleus. In addition, transgenic Arabidopsis lines exhibited higher concentration levels of the antioxidant enzymes measured, including POD, SOD, and CAT, under drought conditions with very low levels of H2O2 and MDA. Moreover, relative water content and the degree of stomatal opening showed that the transgenic Arabidopsis lines were more tolerant of both stress factors as compared to their wild types.ConclusionIn this study, the structures and functions of the CAMTA family in quinoa were systematically explored. Many CAMTAs may play vital roles in the regulation of organ development, growth, and responses to drought stress. The results of the present study serve as a basis for future functional studies on the quinoa CAMTA family.

Genome-Wide Identification and Characterization of SPL Family Genes in Chenopodium quinoa.

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes encode a large family of plant-specific transcription factors that play important roles in plant growth, development, and stress responses. However, there is little information available on SPL genes in Chenopodiaceae. Here, 23 SPL genes were identified and characterized in the highly nutritious crop Chenopodium quinoa. Chromosome localization analysis indicated that the 23 CqSPL genes were unevenly distributed on 12 of 18 chromosomes. Two zinc finger-like structures and a nuclear location signal were present in the SBP domains of all CqSPLs, with the exception of CqSPL21/22. Phylogenetic analysis revealed that these genes were classified into eight groups (group I–VIII). The exon–intron structure and motif composition of the genes in each group were similar. Of the 23 CqSPLs, 13 were potential targets of miR156/7. In addition, 5 putative miR156-encoding loci and 13 putative miR157-encoding loci were predicted in the quinoa genome, and they were unevenly distributed on chromosome 1–4. The expression of several Cqu-MIR156/7 loci was confirmed by reverse transcription polymerase chain reaction in seedlings. Many putative cis-elements associated with light, stress, and phytohormone responses were identified in the promoter regions of CqSPLs, suggesting that CqSPL genes are likely involved in the regulation of key developmental processes and stress responses. Expression analysis revealed highly diverse expression patterns of CqSPLs among tissues. Many CqSPLs were highly expressed in leaves, flowers, and seeds, and their expression levels were low in the roots, suggesting that CqSPLs play distinct roles in the development and growth of quinoa. The expression of 13 of 23 CqSPL genes responded to salt treatment (11 up-regulated and 2 down-regulated). A total of 22 of 23 CqSPL genes responded to drought stress (21 up-regulated and 1 down-regulated). Moreover, the expression of 14 CqSPL genes was significantly altered following cadmium treatment (3 up-regulated and 11 down-regulated). CqSPL genes are thus involved in quinoa responses to salt/drought and cadmium stresses. These findings provide new insights that will aid future studies of the biological functions of CqSPLs in C. quinoa.

Involvement of Auxin-Mediated CqEXPA50 Contributes to Salt Tolerance in Quinoa (Chenopodium quinoa) by Interaction with Auxin Pathway Genes.

Soil salinization is a global problem that limits crop yields and threatens agricultural development. Auxin-induced expansins contribute to plant salt tolerance through cell wall loosening. However, how auxins and expansins contribute to the adaptation of the halophyte quinoa (Chenopodium quinoa) to salt stress has not yet been reported. Here, auxin was found to contribute to the salt tolerance of quinoa by promoting the accumulation of photosynthetic pigments under salt stress, maintaining enzymatic and nonenzymatic antioxidant systems and scavenging excess reactive oxygen species (ROS). The Chenopodium quinoa expansin (Cqexpansin) family and the auxin pathway gene family (Chenopodium quinoa auxin response factor (CqARF), Chenopodium quinoa auxin/indoleacetic acid (CqAux/IAA), Chenopodium quinoa Gretchen Hagen 3 (CqGH3) and Chenopodium quinoa small auxin upregulated RNA (CqSAUR)) were identified from the quinoa genome. Combined expression profiling identified Chenopodium quinoa α-expansin 50 (CqEXPA50) as being involved in auxin-mediated salt tolerance. CqEXPA50 enhanced salt tolerance in quinoa seedlings was revealed by transient overexpression and physiological and biochemical analyses. Furthermore, the auxin pathway and salt stress-related genes regulated by CqEXPA50 were identified. The interaction of CqEXPA50 with these proteins was demonstrated by bimolecular fluorescence complementation (BIFC). The proteins that interact with CqEXPA50 were also found to improve salt tolerance. In conclusion, this study identified some genes potentially involved in the salt tolerance regulatory network of quinoa, providing new insights into salt tolerance.

Malate-mediated CqMADS68 enhances aluminum tolerance in quinoa seedlings through interaction with CqSTOP6, CqALMT6 and CqWRKY88

Aluminum (Al) stress in acidic soils has severe negative effects on crop productivity. In this study, the alleviating effect and related mechanism of malate on Al stress in quinoa (Chenopodium quinoa) seedlings were investigated. The findings indicated that malate alleviated the growth inhibition of quinoa seedlings under Al stress, maintained the enzymatic and nonenzymatic antioxidant systems, and aided resistance to the damage caused by excessive reactive oxygen species (ROS). Under Al stress, malate significantly increased the contents of chlorophyll and carotenoids in quinoa shoots by 103.8% and 240.7%, and significantly increased the ratios of glutathione (GSH)/oxidized glutathione (GSSG), and ascorbate (AsA)/dehydroascorbate (DHA) in roots by 59.9% and 699.2%, respectively. However, malate significantly decreased the superoxide radical (O2•-), hydrogen peroxide (H2O2), malondialdehyde (MDA) and Al contents in quinoa roots under Al stress by 32.7%, 60.9%, 63.1% and 49%, respectively. Moreover, the CqMADS family and the Al stress-responsive gene families (CqSTOP, CqALMT, and CqWRKY) were identified from the quinoa genome. Comprehensive expression profiling identified CqMADS68 as being involved in malate-mediated Al resistance. Transient overexpression of CqMADS68 increased Al tolerance in quinoa seedlings. More importantly, we found that CqMADS68 regulated the expression of CqSTOP6, CqALMT6 and CqWRKY88 and further demonstrated the interaction of CqMADS68 with CqSTOP6, CqALMT6 and CqWRKY88 by bimolecular fluorescence complementation (BIFC) experiments. Moreover, transient overexpression and physiological and biochemical analyses demonstrated that CqSTOP6, CqALMT6 and CqWRKY88 could also improve Al tolerance by maintaining the antioxidant capacity of quinoa seedlings. Taken together, these findings reveal that CqMADS68, CqSTOP6, CqALMT6 and CqWRKY88 may be important contributors to the Al tolerance regulatory network in quinoa, providing new insights into Al stress resistance.

Genome-wide identification, phylogenetic analysis, and expression profiles of trihelix transcription factor family genes in quinoa (Chenopodium quinoa Willd.) under abiotic stress conditions

BackgroundThe trihelix family of transcription factors plays essential roles in the growth, development, and abiotic stress response of plants. Although several studies have been performed on the trihelix gene family in several dicots and monocots, this gene family is yet to be studied in Chenopodium quinoa (quinoa).ResultsIn this study, 47 C. quinoa trihelix (CqTH) genes were in the quinoa genome. Phylogenetic analysis of the CqTH and trihelix genes from Arabidopsis thaliana and Beta vulgaris revealed that the genes were clustered into five subfamilies: SIP1, GTγ, GT1, GT2, and SH4. Additionally, synteny analysis revealed that the CqTH genes were located on 17 chromosomes, with the exception of chromosomes 8 and 11, and 23 pairs of segmental duplication genes were detected. Furthermore, expression patterns of 10 CqTH genes in different plant tissues and at different developmental stages under abiotic stress and phytohormone treatment were examined. Among the 10 genes, CqTH02, CqTH25, CqTH18, CqTH19, CqTH25, CqTH31, and CqTH36, were highly expressed in unripe achenes 21 d after flowering and in mature achenes compared with other plant tissues. Notably, the 10 CqTH genes were upregulated in UV-treated leaves, whereas CqTH36 was consistently upregulated in the leaves under all abiotic stress conditions.ConclusionsThe findings of this study suggest that gene duplication could be a major driver of trihelix gene evolution in quinoa. These findings could serve as a basis for future studies on the roles of CqTH transcription factors and present potential genetic markers for breeding stress-resistant and high-yielding quinoa varieties.

Genome-Wide Analysis of CqCrRLK1L and CqRALF Gene Families in Chenopodium quinoa and Their Roles in Salt Stress Response.

Chenopodium quinoa is a halophyte with exceptional nutritional qualities, and therefore it is potentially an ideal crop to grow in saline soils, not only addressing the problem of land salinization, but also providing nutrient food for the health of humans. Currently, the molecular mechanisms underlying salt tolerance in quinoa are still largely unknown. In Arabidopsis thaliana, Catharanthus roseus receptor-like kinase (CrRLK1Ls) FERONIA (FER) and its ligands rapid alkalinization factors (RALFs) have been reported that participate in the regulation of salt tolerance. Here, we performed a genome-wide analysis and identified 26 CqCrRLK1L and 18 CqRALF family genes in quinoa genome. Transcriptomic profiling of the leaf, root, stamen, and pistil tissues of quinoa reveals that different CqCrRLK1L and CqRALF genes exhibit tissue-specific expression patterns, which is consistent with that observed in other plant species. RNA-seq data show that three CqCrRLK1L genes are highly up-regulated after salt treatment, suggesting that some CqCrRLK1L family genes are transcriptionally responsive to salt stress in quinoa. Biochemical study indicates that CqRALF15, a paralog of Arabidopsis RALF22, is physically associated with CrRLK1L proteins CqFER and AtFER. CqRALF15 and AtRALF22 are functionally conserved in inducing the internalization of AtFER and in triggering root growth inhibition in both quinoa and Arabidopsis. Moreover, overexpression of CqRALF15 in Arabidopsis results in enhanced leaf bleaching under salt stress, indicating that CqRALF15 is involved in salt stress response. Together, our study characterizes CqCrRLK1L and CqRALF family genes in quinoa at genomic, transcriptional, and protein levels, and provides evidence to support their roles in salt stress response.

Identification and expression analysis of the CqSnRK2 gene family and a functional study of the CqSnRK2.12 gene in quinoa (Chenopodium quinoa Willd.)

BackgroundSucrose non-fermenting 1 (SNF1)-associated protein kinase 2 (SnRK2) proteins belong to a relatively small family of plant-specific serine/threonine (Ser/Thr) protein kinases. SnRK2s participate in the abscisic acid (ABA) signaling pathway and play important roles in many biotic and abiotic stresses. At present, no SnRK2 gene has been reported in quinoa, and the recently published genome for this species provides an opportunity to identify and characterize the SnRK2 gene family.ResultsWe identified 13 SnRK2 genes in the C. quinoa genome by bioinformatics analysis. Based on their phylogenetic relationships, these genes were divided into three subfamilies, similar to the situation in other plant species. Gene duplication analysis showed that there were seven pairs of homologous genes in the CqSnRK2 family, and that purifying selection played an important role in the evolution of SnRK2 genes. Gene structure analysis showed that the first exon in the SnRK2 family genes has the same length as the last exon, and that CqSnRK2 genes in the same subfamily have similar gene structures. Sequence analysis showed that the N-terminal region contains three highly conserved motifs. In addition, many kinds of cis-elements were identified in the promoter region of CqSnRK2, including those for hormone responses, stress responses, and tissue-specific expression. Transcription data analysis and qRT-PCR results showed that CqSnRK2 has different expression patterns in roots, stems, and leaves, and responded to biotic and abiotic stresses such as low temperature, salt, drought, and abscisic acid (ABA). In addition, we found that the protein encoded by CqSnRK2.12 was localized to the cytoplasm and nucleus, and there was no self-activation. The results of CqSnRK2.12 overexpression showed that transgenic Arabidopsis thaliana lines had increased drought tolerance compared to the controls.ConclusionThe results of our study provide references for further studies on the evolution, function, and expression of the SnRK2 gene family in quinoa.

Full-Length Transcriptome Sequencing Reveals the Impact of Cold Stress on Alternative Splicing in Quinoa.

Quinoa is a cold-resistant and nutrient-rich crop. To decipher the cold stress response of quinoa, the full-length transcriptomes of the cold-resistant quinoa variety CRQ64 and the cold-sensitive quinoa variety CSQ5 were compared. We identified 55,389 novel isoforms and 6432 novel genes in these transcriptomes. Under cold stress, CRQ64 had more differentially expressed genes (DEGs) and differentially alternative splicing events compared to non-stress conditions than CSQ5. DEGs that were specifically present only in CRQ64 were significantly enriched in processes which contribute to osmoregulation and ROS homeostasis in plants, such as sucrose metabolism and phenylpropanoid biosynthesis. More genes with differential alternative splicing under cold stress were enriched in peroxidase functions in CRQ64. In total, 5988 transcription factors and 2956 long non-coding RNAs (LncRNAs) were detected in this dataset. Many of these had altered expression patterns under cold stress compared to non-stress conditions. Our transcriptome results demonstrate that CRQ64 undergoes a wider stress response than CSQ5 under cold stress. Our results improved the annotation of the quinoa genome and provide new insight into the mechanisms of cold resistance in quinoa.

Diversity Control in Early Transcribed ROS/RNS-balancing Genes: A Common Mechanism for Healthy Resilience?

Diversity Control in Early Transcribed ROS/RNS-balancing Genes: A Common Mechanism for Healthy Resilience?

Quinoa genome assembly employing genomic variation for guided scaffolding

Key messageWe 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.

Involvement of several putative transporters of different families in β-cyclocitral-induced alleviation of cadmium toxicity in quinoa (Chenopodium quinoa) seedlings

Cadmium (Cd) has a serious negative impact on crop growth and human food security. This study investigated the alleviating effect of β-cyclocitral, a potential heavy metal barrier, on Cd stress in quinoa seedlings and the associated mechanisms. Our results showed that β-cyclocitral alleviated Cd stress-induced growth inhibition in quinoa seedlings and promoted quinoa seedling root development under Cd stress. Moreover, it maintained the antioxidant system of quinoa seedlings, including the enzymatic, i.e., superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), ascorbate peroxidase (APX), and nonenzymatic, i.e., reduced glutathione (GSH) and ascorbic acid (ASA), antioxidants, which eliminate the damage from excessive reactive oxygen species (ROS). Our results showed that β-cyclocitral could reduce the amount of Cd absorbed by roots. Furthermore, we systematically identified five transporter families from the quinoa genome, and the RT-qPCR results showed that ZIP, Nramp and YSL gene families were downregulated by β-cyclocitral to reduce Cd uptake by roots. Thus, β-cyclocitral promoted the growth, photosynthetic capacity and antioxidant capacity of the aboveground parts of quinoa seedlings. Taken together, these results suggested that the β-cyclocitral-induced decrease in Cd uptake may be caused by the downregulation of several selected transporter genes.

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.