Analysis of the quinoa genome reveals conservation and divergence of the flowering pathways

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

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

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.

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.

An overview of floral regulatory genes in annual and perennial plants

The floral initiation in angiosperms is a complex process influenced by endogenous and exogenous signals. With this approach, we aim to provide a comprehensive review to integrate this complex floral regulatory process and summarize the regulatory genes and their functions in annuals and perennials. Seven primary paths leading to flowering have been discovered in Arabidopsis under several growth condition that include; photoperiod, ambient temperature, vernalization, gibberellins, autonomous, aging and carbohydrates. These pathways involve a series of interlinked signaling pathways that respond to both internal and external signals, such as light, temperature, hormones, and developmental cues, to coordinate the expression of genes that are involved in flower development. Among them, the photoperiodic pathway was the most important and conserved as some of the fundamental loci and mechanisms are shared even by closely related plant species. The activation of floral regulatory genes such as FLC, FT, LFY, and SOC1 that determine floral meristem identity and the transition to the flowering stage result from the merging of these pathways. Recent studies confirmed that alternative splicing, antisense RNA and epigenetic modification play crucial roles by regulating the expression of genes related to blooming. In this review, we documented recent progress in the floral transition time in annuals and perennials, with emphasis on the specific regulatory mechanisms along with the application of various molecular approaches including overexpression studies, RNA interference and Virus-induced flowering. Furthermore, the similarities and differences between annual and perennial flowering will aid significant contributions to the field by elucidating the mechanisms of perennial plant development and floral initiation regulation.

Transcriptome analysis of critical genes related to flowering in Mikania micrantha at different altitudes provides insights for a potential control

BACKGROUND

Mikania micrantha is a vine with strong invasion ability, and its strong sexual reproduction ability is not only the main factor of harm, but also a serious obstacle to control. M. micrantha spreads mainly through seed production. Therefore, inhibiting the flowering and seed production of M. micrantha is an effective strategy to prevent from continuing to spread.

RESULT

The flowering number of M. micrantha is different at different altitudes. A total of 67.01 Gb of clean data were obtained from nine cDNA libraries, and more than 83.47% of the clean reads were mapped to the reference genome. In total, 5878 and 7686 significantly differentially expressed genes (DEGs) were found in E2 vs. E9 and E13 vs. E9, respectively. Based on the background annotation and gene expression, some candidate genes related to the flowering pathway were initially screened, and their expression levels in the three different altitudes in flower bud differentiation showed the same trend. That is, at an altitude of 1300 m, the flower integration gene and flower meristem gene were downregulated (such as SOC1 and AP1), and the flowering inhibition gene was upregulated (such as FRI and SVP). Additionally, the results showed that there were many DEGs involved in the hormone signal transduction pathway in the flower bud differentiation of M. micrantha at different altitudes.

CONCLUSIONS

Our results provide abundant sequence resources for clarifying the underlying mechanisms of flower bud differentiation and mining the key factors inhibiting the flowering and seed production of M. micrantha to provide technical support for the discovery of an efficient control method.

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.

Transcriptomic Analysis of Early Flowering Signals in ‘Royal’ Flax

Canada is one of the world’s leading producers and exporters of flax seed, with most production occurring in the Prairie Provinces. However, reduced season length and risk of frost restricts production in the northern grain belt of the Canadian Prairies. To expand the growing region of flax and increase production in Canada, flax breeders need to develop earlier-flowering varieties capable of avoiding the risk of abiotic stress. A thorough understanding of flowering control of flax is essential for the efficient breeding of such lines. We identified 722 putative flax flowering genes that span all major flowering-time pathways. Frequently, we found multiple flax homologues for a single Arabidopsis flowering gene. We used RNA sequencing to quantify the expression of genes in the shoot apical meristem (SAM) at 10, 15, 19, and 29 days after planting (dap) using the ‘Royal’ cultivar. We observed the expression of 80% of putative flax flowering genes and the differential expression of only 30%; these included homologues of major flowering regulators, such as SOC1, FUL, and AP1. We also found enrichment of differentially expressed genes (DEGs) in transcription factor (TF) families involved in flowering. Finally, we identified the candidates’ novel flowering genes amongst the uncharacterized flax genes. Our transcriptomic dataset provides a useful resource for investigating the regulatory control of the transition to flowering in flax and for the breeding of northern-adapted varieties.

Haplotype variations of major flowering time genes in quinoa unveil their role in the adaptation to different environmental conditions

Response to photoperiod is of major importance in crop production. It defines the adaptation of plants to local environments. Quinoa is a short-day plant which had been domesticated in the Andeans regions. We wanted to understand the adaptation to long-day conditions by studying orthologues of two major flowering time regulators of Arabidopsis, FLOWERING LOCUS T (FT) and CONSTANS (CO) in quinoa accessions with contrasting photoperiod response. By searching the quinoa reference genome sequence, we identified 24 FT and six CO homologs. CqFT genes displayed remarkably different expression patterns between long- and short-day conditions, whereas the influence of the photoperiod on CqCOL expressions was moderate. Cultivation of 276 quinoa accessions under short- and long-day conditions revealed great differences in photoperiod sensitivity. After sequencing their genomes, we identified large sequence variations in 12 flowering time genes. We found non-random distribution of haplotypes across accessions from different geographical origins, highlighting the role of CqFT and CqCOL genes in the adaptation to different day-length conditions. We identified five haplotypes causing early flowering under long days. This study provides assets for quinoa breeding because superior haplotypes can be assembled in a predictive breeding approach to produce well-adapted early flowering lines under long-day photoperiods.

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.

Identification of the specific long-noncoding RNAs involved in night-break mediated flowering retardation in Chenopodium quinoa

Background

Night-break (NB) has been proven to repress flowering of short-day plants (SDPs). Long-noncoding RNAs (lncRNAs) play key roles in plant flowering. However, investigation of the relationship between lncRNAs and NB responses is still limited, especially in Chenopodium quinoa, an important short-day coarse cereal.

Results

In this study, we performed strand-specific RNA-seq of leaf samples collected from quinoa seedlings treated by SD and NB. A total of 4914 high-confidence lncRNAs were identified, out of which 91 lncRNAs showed specific responses to SD and NB. Based on the expression profiles, we identified 17 positive- and 7 negative-flowering lncRNAs. Co-expression network analysis indicated that 1653 mRNAs were the common targets of both types of flowering lncRNAs. By mapping these targets to the known flowering pathways in model plants, we found some pivotal flowering homologs, including 2 florigen encoding genes (FT (FLOWERING LOCUS T) and TSF (TWIN SISTER of FT) homologs), 3 circadian clock related genes (EARLY FLOWERING 3 (ELF3), LATE ELONGATED HYPOCOTYL (LHY) and ELONGATED HYPOCOTYL 5 (HY5) homologs), 2 photoreceptor genes (PHYTOCHROME A (PHYA) and CRYPTOCHROME1 (CRY1) homologs), 1 B-BOX type CONSTANS (CO) homolog and 1 RELATED TO ABI3/VP1 (RAV1) homolog, were specifically affected by NB and competed by the positive and negative-flowering lncRNAs. We speculated that these potential flowering lncRNAs may mediate quinoa NB responses by modifying the expression of the floral homologous genes.

Conclusions

Together, the findings in this study will deepen our understanding of the roles of lncRNAs in NB responses, and provide valuable information for functional characterization in future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07605-2.

Prospects for the accelerated improvement of the resilient crop quinoa

Crops tolerant to drought and salt stress may be developed by two approaches. First, major crops may be improved by introducing genes from tolerant plants. For example, many major crops have wild relatives that are more tolerant to drought and high salinity than the cultivated crops, and, once deciphered, the underlying resilience mechanisms could be genetically manipulated to produce crops with improved tolerance. Secondly, some minor (orphan) crops cultivated in marginal areas are already drought and salt tolerant. Improving the agronomic performance of these crops may be an effective way to increase crop and food diversity, and an alternative to engineering tolerance in major crops. Quinoa (Chenopodium quinoa Willd.), a nutritious minor crop that tolerates drought and salinity better than most other crops, is an ideal candidate for both of these approaches. Although quinoa has yet to reach its potential as a fully domesticated crop, breeding efforts to improve the plant have been limited. Molecular and genetic techniques combined with traditional breeding are likely to change this picture. Here we analyse protein-coding sequences in the quinoa genome that are orthologous to domestication genes in established crops. Mutating only a limited number of such genes by targeted mutagenesis appears to be a promising route for accelerating the improvement of quinoa and generating a nutritious high-yielding crop that can meet the future demand for food production in a changing climate.

Morphological Characterization of Flower Buds Development and Related Gene Expression Profiling at Bud Break Stage in Heterodichogamous Cyclocarya paliurus (Batal.) lljinskaja

Cyclocarya paliurus (Batal.) Iljinskaja, a unique species growing in southern China, is a multi-function tree species with medicinal, healthcare, material, and ornamental values. So far, sexual reproduction is the main method for extensive cultivation of C. paliurus plantations, but this is limited by low seed plumpness resulted from the character of heterodichogamy. Phenological observations have revealed the asynchronism of flower development in this species. However, its molecular mechanism remains largely unknown. To reveal molecular mechanism of heterodichogamy in C. paliurus, transcriptome of female (F) and male (M) buds from two mating types (protandry, PA; protogyny, PG) at bud break stage were sequenced using Illumina Hiseq 4000 platform. The expression patterns of both 32 genes related to flowering and 58 differentially expressed transcription factors (DETFs) selected from 6 families were divided four groups (PG-F, PG-M, PA-F, and PA-M) into two categories: first flowers (PG-F and PA-M) and later flowers (PA-F and PG-M). The results indicated that genes related to plant hormones (IAA, ABA, and GA) synthesis and response, glucose metabolism, and transcription factors (especially in MIKC family) played significant roles in regulating asynchronism of male and female flowers in the same mating type. The expression of DETFs showed two patterns. One contained DETFs up-regulated in first flowers in comparison to later flowers, and the other was the reverse. Nine genes related to flowering were selected for qRT-PCR to confirm the accuracy of RNA-seq, and generally, the RPKM values of these genes were consistent with the result of qRT-PCR. The results of this work could improve our understanding in asynchronism of floral development within one mating type in C. paliurus at transcriptional level, as well as lay a foundation for further study in heterodichogamous plants.