Transcriptional Responses of Chilean Quinoa (Chenopodium quinoa Willd.) Under Water Deficit Conditions Uncovers ABA-Independent Expression Patterns

Quinoa: A Promising Crop for Resolving the Bottleneck of Cultivation in Soils Affected by Multiple Environmental Abiotic Stresses

Quinoa (Chenopodium quinoa Willd.) has gained worldwide recognition for its nutritional values, adaptability to diverse environments, and genetic diversity. This review explores the current understanding of quinoa tolerance to environmental stress, focusing on drought, salinity, heat, heavy metals, and UV-B radiation. Although drought and salinity have been extensively studied, other stress factors remain underexplored. The ever-increasing incidence of abiotic stress, exacerbated by unpredictable weather patterns and climate change, underscores the importance of understanding quinoa's responses to these challenges. Global gene banks safeguard quinoa's genetic diversity, supporting breeding efforts to develop stress-tolerant varieties. Recent advances in genomics and molecular tools offer promising opportunities to improve stress tolerance and increase the yield potential of quinoa. Transcriptomic studies have shed light on the responses of quinoa to drought and salinity, yet further studies are needed to elucidate its resilience to other abiotic stresses. Quinoa's ability to thrive on poor soils and limited water resources makes it a sustainable option for land restoration and food security enterprises. In conclusion, quinoa is a versatile and robust crop with the potential to address food security challenges under environmental constraints.

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.

Evaluating relationships between seed morphological traits and seed dormancy in Chenopodium quinoa Willd

Introduction

Quinoa is a high-value, nutritious crop that performs well in variable environments, marginal soils, and in diverse crop rotations. Quinoa's many attributes make it an ideal crop for supporting human health in global communities and economies. To date, quinoa research has largely focused on traits in adult plants important for enhancing plant phenotypic plasticity, abiotic stress, disease resistance, and yield. Fewer studies have evaluated quinoa seed dormancy and suggest that most modern quinoa varieties have weak or no seed dormancy, and a narrow window of seed viability post-harvest. In other crops, diminished seed dormancy is a major risk factor for preharvest sprouting (PHS; germination on the panicle due to rain prior to harvest) and may also pose a similar risk for quinoa.

Methods

This study (1) developed a dormancy screening assay to characterize seed dormancy strength in a large collection of quinoa varieties, (2) investigated if morphological variables including seed coat color, seed coat thickness, seed shape including eccentricity which evaluates the roundness or flatness of a seed, and other agronomic traits like crude protein content and seed moisture, contribute to quinoa seed dormancy, and (3) evaluated the use of a phenetic modeling approach to explore relationships between seed morphology and seed dormancy.

Results

Dormancy screening indicated seed dormancy ranges in quinoa varieties from none to strong dormancy. Further, phenetic modeling approaches indicate that seed coat thickness and eccentricity are important morphological variables that impact quinoa seed dormancy strength.

Conclusions

While dormancy screening and phenetic modeling approaches do not provide a direct solution to preventing PHS in quinoa, they do provide new tools for identifying dormant varieties as well as morphological variables contributing to seed dormancy.

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.

An Evidence Theory and Fuzzy Logic Combined Approach for the Prediction of Potential ARF-Regulated Genes in Quinoa

Quinoa constitutes among the tolerant plants to the challenging and harmful abiotic environmental factors. Quinoa was selected as among the model crops destined for bio-saline agriculture that could contribute to the staple food security for an ever-growing worldwide population under various climate change scenarios. The auxin response factors (ARFs) constitute the main contributors in the plant adaptation to severe environmental conditions. Thus, the determination of the ARF-binding sites represents the major step that could provide promising insights helping in plant breeding programs and improving agronomic traits. Hence, determining the ARF-binding sites is a challenging task, particularly in species with large genome sizes. In this report, we present a data fusion approach based on Dempster-Shafer evidence theory and fuzzy set theory to predict the ARF-binding sites. We then performed an "In-silico" identification of the ARF-binding sites in Chenopodium quinoa. The characterization of some known pathways implicated in the auxin signaling in other higher plants confirms our prediction reliability. Furthermore, several pathways with no or little available information about their functions were identified to play important roles in the adaptation of quinoa to environmental conditions. The predictive auxin response genes associated with the detected ARF-binding sites may certainly help to explore the biological roles of some unknown genes newly identified 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.

Yield, phenology and triterpene saponins in Colombian quinoa

Exploring yield, phenology and their relationship with secondary metabolites in seeds provides a fundamental analysis that expands knowledge on the nutritional quality of seeds and the effect on productive potential. This knowledge is fundamental when improving or selecting nutritionally important crops, including Chenopodium quinoa Willd, which has excellent nutritional properties and contributes to global food security. This species contains saponins, a metabolite that imparts a bitter taste and can be highly toxic to consumers in large quantities. Therefore, the identification and selection of genotypes according to their saponin contents and outstanding agronomic characteristics are fundamental objectives for the genetic improvement programs of these species. Therefore, the objective of this research was to evaluate the characteristics of the grain, the phenology and the saponin content of 30 C. quinoa accessions with an aim to select or relate genotypes according to their yield and grain quality. The accessions were sown using randomized complete blocks (RCB) with nine repetitions for each material. Seven FAO-defined descriptors were evaluated to characterize the grain and physiological maturity. Saponin was extracted using microwave, and the quantification was done with high-performance liquid chromatography (HPLC) which a UV-VIS detector at 277 nm wavelength. The accessions were classified according to their phenology: semi-late (56.7%), late (36.7%), and semi-early (3.3%). The total triterpene saponin content varied from 0.018 to 0.537%. The multivariate and cluster analyses formed groups of accessions with good yields (>62.02 g of seeds per plant) and desirable grain morphological characteristics. The more suitable accessions for the production of saponins are Quinoa semiamarga (0.537%), Quinoa peruana (0.412%) and Amarilla de maranganí (0.305%). Quinoa real and Quinoa primavera are more suitable for food products, which can be used as parents in future quinoa genetic improvement programs in Colombia.

Transcriptome dynamics uncovers long non-coding RNAs response to salinity stress in Chenopodium quinoa

Chenopodium quinoa is a crop with outstanding tolerance to saline soil, but long non-coding RNAs (LncRNAs) expression profile driven by salt stress in quinoa has rarely been observed yet. Based on the high-quality quinoa reference genome and high-throughput RNA sequencing (RNA-seq), genome-wide identification of LncRNAs was performed, and their dynamic response under salt stress was then investigated. In total, 153,751 high-confidence LncRNAs were discovered and dispersed intensively in chromosomes. Expression profile analysis demonstrated significant differences between LncRNAs and coding RNAs. Under salt stress conditions, 4,460 differentially expressed LncRNAs were discovered, of which only 54 were differentially expressed at all the stress time points. Besides, strongly significantly correlation was observed between salt-responsive LncRNAs and their closest neighboring genes (r = 0.346, p-value < 2.2e-16). Furthermore, a weighted co-expression network was then constructed to infer the potential biological functions of LncRNAs. Seven modules were significantly correlated with salt treatments, resulting in 210 hub genes, including 22 transcription factors and 70 LncRNAs. These results indicated that LncRNAs might interact with transcription factors to respond to salinity stress. Gene ontology enrichment of the coding genes of these modules showed that they were highly related to regulating metabolic processes, biological regulation and response to stress. This study is the genome-wide analysis of the LncRNAs responding to salt stress in quinoa. The findings will provide a solid framework for further functional research of salt responsive LncRNAs, contributing to quinoa genetic improvement.

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

Background

Sucrose 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.

Results

We 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.

Conclusion

The results of our study provide references for further studies on the evolution, function, and expression of the SnRK2 gene family in quinoa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08626-1.

Morphological and Physiological Traits Associated with Yield under Reduced Irrigation in Chilean Coastal Lowland Quinoa

Quinoa (Chenopodium quinoa Willd.) is a genetically diverse crop that has gained popularity in recent years due to its high nutritional content and ability to tolerate abiotic stresses such as salinity and drought. Varieties from the coastal lowland ecotype are of particular interest due to their insensitivity to photoperiod and their potential to be cultivated in higher latitudes. We performed a field experiment in the southern Atacama Desert in Chile to investigate the responses to reduced irrigation of nine previously selected coastal lowland self-pollinated (CLS) lines and the commercial cultivar Regalona. We found that several lines exhibited a yield and seed size superior to Regalona, also under reduced irrigation. Plant productivity data were analyzed together with morphological and physiological traits measured at the visible inflorescence stage to estimate the contribution of these traits to differences between the CLS lines and Regalona under full and reduced irrigation. We applied proximal sensing methods and found that thermal imaging provided a promising means to estimate variation in plant water use relating to yield, whereas hyperspectral imaging separated lines in a different way, potentially related to photosynthesis as well as water use.

Transcriptional Regulation of Quinoa Seed Quality: Identification of Novel Candidate Genetic Markers for Increased Protein Content

Quinoa (Chenopodium quinoa Willd.) is a crop that has great potential for increased cultivation in diverse climate regions. The seed protein quality obtained from this crop is high concerning the requirements to meet human nutritional needs, but the seed protein content is relatively low if compared to crops such as grain legumes. Increased seed protein content is desirable for increasing the economic viability of this crop in order for it to be used as a protein crop. In this study, we characterized three genotypes of quinoa with different levels of seed protein content. By performing RNA sequencing of developing seeds, we determined the genotype differences in gene expression and identified genetic polymorphisms that could be associated with increased protein content. Storage nutrient analyses of seeds of three quinoa genotypes (Titicaca, Pasankalla, and Regalona) from different ecoregions grown under controlled climate conditions showed that Pasankalla had the highest protein content (20%) and the lowest starch content (46%). Our seed transcriptome analyses revealed highly differentially expressed transcripts (DETs) in Pasankalla as compared to the other genotypes. These DETs encoded functions in sugar transport, starch and protein synthesis, genes regulating embryo size, and seed transcription factors. We selected 60 genes that encode functions in the central carbon metabolism and transcription factors as potential targets for the development of high-precision markers. Genetic polymorphisms, such as single nucleotide polymorphisms (SNPs) and base insertions and deletions (InDels), were found in 19 of the 60 selected genes, which can be further evaluated for the development of genetic markers for high seed protein content in quinoa. Increased cultivation of quinoa can contribute to a more diversified agriculture and support the plant protein diet shift. The identification of quinoa genotypes with contrasting seed quality can help establish a model system that can be used for the identification of precise breeding targets to improve the seed quality of quinoa. The data presented in this study based on nutrient and transcriptome analyses contribute to an enhanced understanding of the genetic regulation of seed quality traits in quinoa and suggest high-precision candidate markers for such traits.

Systematic Identification and Validation of Housekeeping and Tissue-Specific Genes in Allotetraploid Chenopodium quinoa

Quinoa is a gluten-free food crop that contains all the essential amino acids and vitamins. The selection of proper housekeeping and tissue-specific genes is the crucial prerequisite for gene expression analysis using the common approach, real-time quantitative PCR (RT-qPCR). In this study, we identified 40 novel candidate housekeeping genes by the minimum transcript per million (TPM), coefficient of variation (CV) and maximum fold change (MFC) methods and 19 candidate tissue-specific genes by the co-expression network method based on an RNA-seq dataset that included 53 stem, leaf, flower and seed samples, as well as additional shoot and root samples under different stresses. The expression stability of 12 housekeeping and tissue-specific genes, as well as that of another two traditionally used housekeeping genes, was further evaluated using qPCR and ranked using NormFinder, BestKeeper and the comparative delta-Ct method. The results demonstrated that MIF, RGGA, VATE and UBA2B were ranked as the top four most stable candidate housekeeping genes. qPCR analysis also revealed three leaf-specific genes and five root-specific genes, but no stem-specific gene was identified. Gene Ontology (GO) enrichment analysis identified that housekeeping genes were mainly enriched in the small molecule metabolic process, organonitrogen compound metabolic process, NAD binding and ligase activity. In addition, tissue-specific genes are closely associated with the major functions of a specific tissue. Specifically, GO terms “photosynthesis” and “thylakoid” were most significantly overrepresented in candidate leaf-specific genes. The novel housekeeping and tissue-specific genes in our study will enable better normalization and quantification of transcript levels in quinoa.

Quinoa for the Brazilian Cerrado: Agronomic Characteristics of Elite Genotypes under Different Water Regimes

Quinoa stands out as an excellent crop in the Cerrado region for cultivation in the off-season or irrigated winter season. Here, we tested the effects of different water regimes on the agronomic characteristics, physiology, and grain quality of different elite quinoa genotypes under field conditions. The experiment was conducted under field conditions at Embrapa Cerrados (Planaltina, DF, Brazil). The experimental design was in randomized blocks, in a split-plot scheme, with four replications. The plots were composed of 18 quinoa genotypes and modified BRS Piabiru (the currently used genotype), and the split-plots were divided into 4 different water regimes. The following variables were evaluated: productivity and productivity per unit of applied water (PUAA), plant height, flavonoids, anthocyanins, gas exchange, chlorophyll, leaf proline, and relative water content. Our results showed that water regimes between 309 and 389 mm can be recommended for quinoa in the Cerrado region. CPAC6 and CPAC13 presented the highest yield and PUAA under high and intermediate WRs, and hence were the most suitable for winter growth under irrigation. CPAC17 is most suitable for off-season growth under rainfed conditions, as it presented the highest PUAA under the low WRs (247 and 150). CPAC9 stood out in terms of accumulation of flavonoids and anthocyanins in all WRs. Physiological analyses revealed different responses of the genotypes to water restriction, together with symptoms of stress under lower water regimes. Our study reinforces the importance of detailed analyses of the relationship between productivity, physiology, and water use when choosing genotypes for planting and harvest in different seasons.

The role of auxin during early berry development in grapevine as revealed by transcript profiling from pollination to fruit set

Auxin is a key phytohormone that modulates fruit formation in many fleshy fruits through the regulation of cell division and expansion. Auxin content rapidly increases after pollination and the manipulation in its levels may lead to the parthenocarpic development. ln Vitis vinifera L., little is known about the early fruit development that encompasses from pollination to fruit set. Pollination/fertilization events trigger fruit formation, and auxin treatment mimics their effect in grape berry set. However, the role of auxin in this process at the molecular level is not well understood. To elucidate the participation of auxin in grapevine fruit formation, morphological, reproductive, and molecular events from anthesis to fruit set were described in sequential days after pollination. Exploratory RNA-seq analysis at four time points from anthesis to fruit set revealed that the highest percentage of genes induced/repressed within the hormone-related gene category were auxin-related genes. Transcript profiling showed significant transcript variations in auxin signaling and homeostasis-related genes during the early fruit development. Indole acetic acid and several auxin metabolites were present during this period. Finally, application of an inhibitor of auxin action reduced cell number and the mesocarp diameter, similarly to unpollinated berries, further confirming the key role of auxin during early berry development. This work sheds light into the molecular features of the initial fruit development and highlights the auxin participation during this stage in grapevine.

Hormonal Regulation in Different Varieties of Chenopodium quinoa Willd. Exposed to Short Acute UV-B Irradiation

Increased ultraviolet-B (UV-B) due to global change can affect plant development and metabolism. Quinoa tolerates extreme conditions including high UV levels. However, the physiological mechanisms behind its abiotic stress tolerance are unclear, especially those related to UV-B. We previously demonstrated that 9.12 kJ m⁻² d⁻¹ may induce UV-B-specific signaling while 18.24 kJ m⁻² d⁻¹ promotes a UV-B-independent response. Here, we explored the effects of these UV-B doses on hormonal regulation linked to plant morphology and defense among diverse varieties. Changes in fluorescence parameters of photosystem II, flavonoids and hormones (indoleacetic acid (IAA), jasmonic acid (JA), abscisic acid (ABA) and salicylic acid (SA)) were surveyed under controlled conditions. Here, we showed that the sensitivity to short acute UV-B doses in varieties from different habitats is influenced by their parental lines and breeding time. UV-B sensitivity does not necessarily correlate with quinoa’s geographical distribution. The role of flavonoids in the UV-B response seems to be different depending on varieties. Moreover, we found that the extent of changes in JA and SA correlate with UV-B tolerance, while the increase of ABA was mainly related to UV-B stress.

Seed Dormancy and Preharvest Sprouting in Quinoa (Chenopodium quinoa Willd.)

Quinoa (Chenopodium quinoa Willd.) is a culturally significant staple food source that has been grown for thousands of years in South America. Due to its natural drought and salinity tolerance, quinoa has emerged as an agronomically important crop for production in marginal soils, in highly variable climates, and as part of diverse crop rotations. Primary areas of quinoa research have focused on improving resistance to abiotic stresses and disease, improving yields, and increasing nutrition. However, an evolving issue impacting quinoa seed end-use quality is preharvest sprouting (PHS), which is when seeds with little to no dormancy experience a rain event prior to harvest and sprout on the panicle. Far less is understood about the mechanisms that regulate quinoa seed dormancy and seed viability. This review will cover topics including seed dormancy, orthodox and unorthodox dormancy programs, desiccation sensitivity, environmental and hormonal mechanisms that regulate seed dormancy, and breeding and non-breeding strategies for enhancing resistance to PHS in quinoa.

Response of bitter and sweet Chenopodium quinoa varieties to cucumber mosaic virus: Transcriptome and small RNASeq perspective

Saponins are secondary metabolites with antiviral properties. Low saponin (sweet) varieties of quinoa (Chenopodium quinoa) have been developed because seeds high in saponins taste bitter. The aim of this study was to elucidate the role of saponin in resistance of quinoa to Cucumber mosaic virus (CMV). Differential gene expression was studied in time-series study of CMV infection. High-throughput transcriptome sequence data were obtained from 36 samples (3 varieties × +/- CMV × 1 or 4 days after inoculation × 3 replicates). Translation, lipid, nitrogen, amino acid metabolism, and mono- and sesquiterpenoid biosynthesis genes were upregulated in CMV infections. In 'Red Head' (bitter), CMV-induced systemic symptoms were concurrent with downregulation of a key saponin biosynthesis gene, TSARL1, four days after inoculation. In local lesion responses (sweet and semi-sweet), TSARL1 levels remained up-regulated. Known microRNAs (miRNA) (81) from 11 miR families and 876 predicted novel miRNAs were identified. Differentially expressed miRNA and short interfering RNA clusters (24nt) induced by CMV infection are predicted to target genomic and intergenic regions enriched in repetitive elements. This is the first report of integrated RNASeq and sRNASeq data in quinoa-virus interactions and provides comprehensive understanding of involved genes, non-coding regions, and biological pathways in virus resistance.

Analysis of Whole Transcriptome RNA-seq Data Reveals Many Alternative Splicing Events in Soybean Roots under Drought Stress Conditions

Alternative splicing (AS) is a common post-transcriptional regulatory mechanism that modulates gene expression to increase proteome diversity. Increasing evidence indicates that AS plays an important role in regulating plant stress responses. However, the mechanism by which AS coordinates with transcriptional regulation to regulate drought responses in soybean remains poorly understood. In this study, we performed a genome-wide analysis of AS events in soybean (Glycine max) roots grown under various drought conditions using the high-throughput RNA-sequencing method, identifying 385, 989, 1429, and 465 AS events that were significantly differentially spliced under very mild drought stress, mild drought stress, severe drought stress, and recovery after severe drought conditions, respectively. Among them, alternative 3' splice sites and skipped exons were the major types of AS. Overall, 2120 genes that experienced significant AS regulation were identified from these drought-treated root samples. Gene Ontology term analysis indicated that the AS regulation of binding activity has vital roles in the drought response of soybean root. Notably, the genes encoding splicing regulatory factors in the spliceosome pathway and mRNA surveillance pathway were enriched according to the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Splicing regulatory factor-related genes in soybean root also responded to drought stress and were alternatively spliced under drought conditions. Taken together, our data suggest that drought-responsive AS acts as a direct or indirect mode to regulate drought response of soybean roots. With further in-depth research of the function and mechanism of AS in the process of abiotic stress, these results will provide a new strategy for enhancing stress tolerance of plants.

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.

Coordinated mechanisms of leaves and roots in response to drought stress underlying full-length transcriptome profiling in Vicia sativa L

BACKGROUND

Common vetch (Vicia sativa L.) is an important self-pollinating annual forage legume and is of interest for drought prone regions as a protein source to feed livestock and human consumption. However, the development and production of common vetch are negatively affected by drought stress. Plants have evolved common or distinct metabolic pathways between the aboveground and underground in response to drought stress. Little is known regarding the coordinated response of aboveground and underground tissues of common vetch to drought stress.

RESULTS

Our results showed that a total of 30,427 full-length transcripts were identified in 12 samples, with an average length of 2278.89 bp. Global transcriptional profiles of the above 12 samples were then analysed via Illumina-Seq. A total of 3464 and 3062 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Gene Ontology (GO) enrichment analyses identified that the dehydrin genes and Δ¹-pyrroline-5-carboxylate synthase were induced for the biosynthesis of proline and water conservation. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis results indicated that the DEGs were significantly enriched in hormone signal transduction, starch and sucrose metabolism, and arginine and proline metabolism, and various drought response candidate genes were also identified. Abscisic acid (ABA; the AREB/ABF-SnRK2 pathway) regulates the activity of AMY3 and BAM1 to induce starch degradation in leaves and increase carbon export to roots, which may be associated with the drought stress responses in common vetch. Among the co-induced transcription factors (TFs), AREB/ABF, bHLH, MYB, WRKY, and AP2/ERF had divergent expression patterns and may be key in the crosstalk between leaves and roots during adaption to drought stress. In transgenic yeast, the overexpression of four TFs increased yeast tolerance to osmotic stresses.

CONCLUSION

The multipronged approach identified in the leaves and roots broadens our understanding of the coordinated mechanisms of drought response in common vetch, and further provides targets to improve drought resistance through genetic engineering.

The combined effect of Cr(III) and NaCl determines changes in metal uptake, nutrient content, and gene expression in quinoa (Chenopodium quinoa Willd.)

Many areas of the world are affected simultaneously by salinity and heavy metal pollution. Halophytes are considered as useful candidates in remediation of such soils due to their ability to withstand both osmotic stress and ion toxicity deriving from high salt concentrations. Quinoa (Chenopodium quinoa Willd) is a halophyte with a high resistance to abiotic stresses (drought, salinity, frost), but its capacity to cope with heavy metals has not yet been fully investigated. In this pot experiment, we investigated phytoextraction capacity, effects on nutrient levels (P and Fe), and changes in gene expression in response to application of Cr(III) in quinoa plants grown on saline or non-saline soil. Plants were exposed for three weeks to 500 mg kg^-1 soil of Cr(NO₃)₃·9H₂O either in the presence or absence of 150 mM NaCl. Results show that plants were able tolerate this soil concentration of Cr(III); the metal was mainly accumulated in roots where it reached the highest concentration (ca. 2.6 mg g^-1 DW) in the presence of NaCl. On saline soil, foliar Na concentration was significantly reduced by Cr(III). Phosphorus translocation to leaves was reduced in the presence of Cr(III), while Fe accumulation was enhanced by treatment with NaCl alone. A real-time RT-qPCR analysis was conducted on genes encoding for sulfate, iron, and phosphate transporters, a phytochelatin, a metallothionein, glutathione synthetase, a dehydrin, Hsp70, and enzymes responsible for the biosynthesis of proline (P5CS), glycine betaine (BADH), tocopherols (TAT), and phenolic compounds (PAL). Cr(III), and especially Cr(III)+NaCl, affected transcript levels of most of the investigated genes, indicating that tolerance to Cr is associated with changes in phosphorus and sulfur allocation, and activation of stress-protective molecules. Moderately saline conditions, in most cases, enhanced this response, suggesting that the halophytism of quinoa could contribute to prime the plants to respond to chromium stress.

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.

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.

Impact of heat and drought stress on peroxisome proliferation in quinoa

Although peroxisomes play a key role in plant metabolism under both normal and stressful growth conditions, the impact of drought and heat stress on the peroxisomes remains unknown. Quinoa represents an informative system for dissecting the impact of abiotic stress on peroxisome proliferation because it is adapted to marginal environments. Here we determined the correlation of peroxisome abundance with physiological responses and yield under heat, drought and heat plus drought stresses in eight genotypes of quinoa. We found that all stresses caused a reduction in stomatal conductance and yield. Furthermore, H₂ O₂ content increased under drought and heat plus drought. Principal component analysis demonstrated that peroxisome abundance correlated positively with H₂ O₂ content in leaves and correlated negatively with yield. Pearson correlation coefficient for yield and peroxisome abundance (r = -0.59) was higher than for commonly used photosynthetic efficiency (r = 0.23), but comparable to those for classical stress indicators such as soil moisture content (r = 0.51) or stomatal conductance (r = 0.62). Our work established peroxisome abundance as a cellular sensor for measuring responses to heat and drought stress in the genetically diverse populations. As heat waves threaten agricultural productivity in arid climates, our findings will facilitate identification of genetic markers for improving yield of crops under extreme weather patterns.

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.

Effect of Water Stress during Grain Filling on Yield, Quality and Physiological Traits of Illpa and Rainbow Quinoa (Chenopodium quinoa Willd.) Cultivars

The total area under quinoa (Chenopodium quinoa Willd.) cultivation and the consumption of its grain have increased in recent years because of its nutritional properties and ability to grow under adverse conditions, such as drought. Climate change scenarios predict extended periods of drought and this has emphasized the need for new crops that are tolerant to these conditions. The main goal of this work was to evaluate crop yield and quality parameters and to characterize the physiology of two varieties of quinoa grown under water deficit in greenhouse conditions. Two varieties of quinoa from the Chilean coast (Rainbow) and altiplano (Illpa) were used, grown under full irrigation or two different levels of water deficit applied during the grain filling period. There were no marked differences in yield and quality parameters between treatments, but the root biomass was higher in plants grown under severe water deficit conditions compared to control. Photosynthesis, transpiration and stomatal conductance decreased with increased water stress in both cultivars, but the coastal variety showed higher water use efficiency and less discrimination of 13C under water deficit. This response was associated with greater root development and a better stomatal opening adjustment, especially in the case of Rainbow. The capacity of Rainbow to increase its osmoregulant content (compounds such as proline, glutamine, glutamate, K and Na) could enable a potential osmotic adjustment in this variety. Moreover, the lower stomatal opening and transpiration rates were also associated with higher leaf ABA concentration values detected in Rainbow. We found negative logarithmic relationships between stomatal conductance and leaf ABA concentration in both varieties, with significant R2 values of 0.50 and 0.22 in Rainbow and Illpa, respectively. These moderate-to-medium values suggest that, in addition to ABA signaling, other causes for stomatal closure in quinoa under drought such as hydraulic regulation may play a role. In conclusion, this work showed that two quinoa cultivars use different strategies in the face of water deficit stress, and these prevent decreases in grain yield and quality under drought conditions.

Transcriptome datasets from leaves and fruits of the sweet cherry cultivars ‘Bing’, ‘Lapins’ and ‘Rainier’

Sweet cherry fruits from different cultivars have different pre- and post-harvest qualities. Here we present the transcriptome profile datasets of leaves and mature fruits of three sweet cherry cultivars (‘Bing’, ‘Lapin’ and ‘Rainier’). Using 454 GS-FLX technology (454 Life Sciences, Roche), transcriptomes of leaves and mature fruits were obtained from these cultivars. These transcriptome data sets are reported here.

Nitrogen physiology of contrasting genotypes of Chenopodium quinoa Willd. (Amaranthaceae)

Quinoa has been highlighted as a promising crop to sustain food security. The selection of physiological traits that allow identification genotypes with high Nitrogen use efficiency (NUE) is a key factor to increase Quinoa cultivation. In order to unveil the underpinning mechanisms for N-stress tolerance in Quinoa, three genotypes with similar phenology, but different NUE were developed under high (HN) or low (LN) nitrogen conditions. N metabolism processes and photosynthetic performance were studied after anthesis and in correlation with productivity to identify principal traits related to NUE. We found that protein content, net photosynthesis and leaf dry-mass were determinant attributes for yield at both HN and LN conditions. Contrastingly, the enhancement of N related metabolites (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{NH}}}_{4}^{+}$$\end{document}NH4+, proline, betacyanins) and processes related with re-assimilation of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{NH}}}_{4}^{+}$$\end{document}NH4+, including an increment of glutamine synthetase activity and up-regulation of CqAMT1,1 transporter expression in leaves, were negatively correlated with grain yield at both N conditions. Biochemical aspects of photosynthesis and root biomass were traits exclusively associated with grain yield at LN. The impact of N supply on seed quality is discussed. These results provide new insights towards the understanding the N metabolism of Quinoa.

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.

Genetic mapping of psl locus and quantitative trait loci for angular leaf spot resistance in cucumber (Cucumis sativus L.)

One of the most important cucumber diseases is bacterial angular leaf spot (ALS), whose increased occurrence in open-field production has been observed over the last years. To map ALS resistance genes, a recombinant inbred line (RIL) mapping population was developed from a narrow cross of cucumber line Gy14 carrying psl resistance gene and susceptible B10 line. Parental lines and RILs were tested under growth chamber conditions as well as in the field for angular leaf spot symptoms. Based on simple sequence repeat and DArTseq, genotyping a genetic map was constructed, which contained 717 loci in seven linkage groups, spanning 599.7 cM with 0.84 cM on average between markers. Monogenic inheritance of the lack of chlorotic halo around the lesions, which is typical for ALS resistance and related with the presence of recessive psl resistance gene, was confirmed. The psl locus was mapped on cucumber chromosome 5. Two major quantitative trait loci (QTL) psl5.1 and psl5.2 related to disease severity were found and located next to each other on chromosome 5; moreover, psl5.1 was co-located with psl locus. Identified QTL were validated in the field experiment. Constructed genetic map and markers linked to ALS resistance loci are novel resources that can contribute to cucumber breeding programs.

Nitrogen Supply Affects Photosynthesis and Photoprotective Attributes During Drought-Induced Senescence in Quinoa

Drought during senescence has become more common in Mediterranean climates in recent years. Chenopodium quinoa Willd has been identified as tolerant to poor soil conditions and drought. Previous observations have found that sufficient nitrogen (N) supply mitigates yield losses under terminal drought conditions. However, there is no understanding of the mechanisms behind this effect. We hypothesized that N up-regulates both photosynthetic and photoprotective elements during drought-induced senescence, alleviating the negative impact of drought on yield. The role of N supply and terminal drought on photoprotection was tested using three Chilean quinoa genotypes from different climatic zones: Faro, UdeC9, and BO78. Plants were grown under high nitrogen (HN) or low nitrogen (LN) conditions and subjected to terminal drought at the onset of senescence. Photosynthetic and photochemical and non-photochemical processes were evaluated at both the onset of drought and after 15 days of drought conditions. N supplementation modified most of the physiological parameters related to photochemical dissipation of energy, photosynthesis, and yield in quinoa. In contrast, water restriction did not affect photosynthesis in quinoa, and its effect on yield was dependent on the genotype. A significant interaction N × G was observed in photosynthesis, relative water content, protein content, Fv/Fm, and chlorophylls. In general, Faro was able to maintain higher levels of these attributes under LN conditions than UdeC9 and BO78. In addition, the interacting effects of N × W regulated the level of most pigments in quinoa as well as the photoprotective induction of non-photochemical quenching (NPQ) during senescence. During terminal drought at LN conditions, Faro presented a larger NPQ induction under drought conditions than UdeC9 and BO78, which was supported by a larger zeaxanthin content and de-epoxidation state of the xanthophyll pool. Interestingly, BO78 did not induce NPQ in response to drought-induced senescence but instead enhanced the content of betacyanins. This response needs to be researched in future works. Finally, we observed that LN supply reduced the correlationship between the de-epoxidation state of the xanthophyll cycle and NPQ. This could be an indication that N supply not only compromised the capacity for photosynthetic performance in quinoa plants, but also affected the plasticity of thermal dissipation, restricting further changes during drought-induced senescence.

Analysis of Stress Resistance Using Next Generation Techniques

Food security for a growing world population remains one of the most challenging tasks. Rapid climate change accelerates the loss of arable land used for crop production, while it simultaneously imposes increasing biotic and abiotic stresses on crop plants. Analysis and molecular understanding of the factors governing stress tolerance is in the focus of scientific and applied research. One plant is often mentioned in the context with stress resistance—Chenopodium quinoa. Through improved breeding strategies and the use of next generation approaches to study and understand quinoa’s salinity tolerance, an important step towards securing food supply is taken.

Genome-Wide Characterization of Heat-Shock Protein 70s from Chenopodium quinoa and Expression Analyses of Cqhsp70s in Response to Drought Stress

Heat-shock proteins (HSPs) are ubiquitous proteins with important roles in response to biotic and abiotic stress. The 70-kDa heat-shock genes (Hsp70s) encode a group of conserved chaperone proteins that play central roles in cellular networks of molecular chaperones and folding catalysts across all the studied organisms including bacteria, plants and animals. Several Hsp70s involved in drought tolerance have been well characterized in various plants, whereas no research on Chenopodium quinoa HSPs has been completed. Here, we analyzed the genome of C. quinoa and identified sixteen Hsp70 members in quinoa genome. Phylogenetic analysis revealed the independent origination of those Hsp70 members, with eight paralogous pairs comprising the Hsp70 family in quinoa. While the gene structure and motif analysis showed high conservation of those paralogous pairs, the synteny analysis of those paralogous pairs provided evidence for expansion coming from the polyploidy event. With several subcellular localization signals detected in CqHSP70 protein paralogous pairs, some of the paralogous proteins lost the localization information, indicating the diversity of both subcellular localizations and potential functionalities of those HSP70s. Further gene expression analyses revealed by quantitative polymerase chain reaction (qPCR) analysis illustrated the significant variations of Cqhsp70s in response to drought stress. In conclusion, the sixteen Cqhsp70s undergo lineage-specific expansions and might play important and varied roles in response to drought stress.