Abstract
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.
Introduction
Quinoa (Chenopodium quinoa Willd.), an ancient annual dicotyledonous crop in the Chenopodiaceae family with approximately 7,000 years of cultivation, originated from the Andes mountains of South America (, ; ). Quinoa, like cereal crops such as rice, maize, and wheat, is most commonly consumed as seeds (; ). Quinoa was a major food crop for the Indians, for instance, the Aztec and Inca civilizations. However, the cultivation process was interrupted after the arrival of Spanish immigrants because the conquerors prohibited local quinoa cultivation. Quinoa was rediscovered hundreds of years later, particularly as the twentieth century began, by developed countries due to its comprehensive nutrition. This orphan crop that belonged to poor Andes local farmers suddenly gained attention and became popular worldwide. Compared with cereal crops such as rice, quinoa offers an excellent balance between protein, oil and carbohydrate. The gluten-free starch of quinoa is suitable for celeriac patients (; ; ; Zou et al., 2017). The protein content of quinoa is as high as 15%, with an excellent balance in amino acids, and comparable to cheese and even better than beef (). The lipid composition of quinoa is also superior due to its high concentration of polyunsaturated fatty acids like omega-3 fatty acid and docosahexaenoic acid (DHA), both of which are essential for human physiological demands and fit to a healthy diet. Finally, although the starch content in quinoa seed is relatively low compared to rice and wheat, it can simultaneously fulfil people’s daily energy consumption with a lower glycemic index. Considering its excellent qualities, quinoa grain is a unique plant product covering all the nutrients for human survival and is recommended by the United Nations (FAO) and World Health Organization (WHO) for potential worldwide consumption. Furthermore, the year 2013 was announced as “The International Year of Quinoa” by the United Nations, which implied that the potential of this emerging crop was becoming more widely recognized (; ).
Quinoa is not only nutritious but also has many excellent agronomic characteristics. Quinoa is highly resistant to various abiotic stresses, including drought, salinity, cold, and soil nutrient defection, which makes it a robust candidate for agricultural development in marginal lands with poor soil conditions (; ; ). Tidal flats and marshlands are potential resources of new farmland worldwide. In such areas, salinity is the primary limiting factor in local agriculture. High salinity can damage cells by causing ionic, osmotic, nutrient, and oxidative stresses, resulting in plants’ growth inhibition and even death (Zhu, 2002; ; Zhao et al., 2021). As a salt-tolerant crop with high economic and nutrient values, quinoa is undoubtedly an attractive choice. Therefore, salt tolerance is a remarkable agronomy trait of quinoa, and many studies have been conducted, yielding promising results. For example, it has been found that quinoa has a unique tissue (epidermal bladder cell) for salt storage to prevent the somatic cell from salt stress (Zou et al., 2017). completed the first chromosome-level quinoa genome with a size as large as 1.4 Gb, which provides an excellent resource for quinoa molecular and genetic study.
Long non-coding RNAs (LncRNAs) are often recognized as transcripts larger than 200 nt in length but have no apparent protein-coding potential (; ). LncRNAs have been discovered to be involved in a variety of biological regulatory processes, and their expression patterns are more tissue-specific than mRNA (; ). According to their relative positions to nearby protein-coding genes in the genome, LncRNAs are generally classified into five types: sense, antisense, bidirectional, intergenic, and intronic (; ; ). Previous studies have indicated that plant LncRNAs play functional roles in signal pathway transmission and molecular regulation under abiotic stresses such as salt stress. LncRNA973, for example, modulated the expression of a number of salt stress-related genes to positively regulate the response to salt stress in cotton (Zhang et al., 2019). In upland cotton, a competing endogenous RNA of miR160b regulated ARF genes in response to salt via a Long non-coding RNA-lncRNA354 (Zhang et al., 2021). A nucleus-localized drought-induced LncRNA (DRIR), which functioned in water transport and ABA signaling, could enhance the tolerance to drought and salt stress in Arabidopsis (). Currently, stress-responsive LncRNAs have been identified in many species, such as maize (, ; ), rice (Zhang et al., 2014), Arabidopsis (), pistachio (), chickpea (), and rapeseed (). However, no systematic study on salt-responsive LncRNAs in quinoa has been reported. This study focused on quinoa LncRNAs and their dynamic responses to salt stress. The findings will provide a massive amount of salt-responsive LncRNAs in quinoa and enlighten the potential patterns of LncRNAs incorporation with the coding genes.
Results
Transcriptome assembly and long non-coding RNAs identification
To profile the LncRNA transcripts in response to salt stress, we performed the time-series dynamic analysis of RNA-sequencing on quinoa roots exposed to high salinity conditions (Supplementary Table 1). In total, 464 million raw reads were generated, and 430 million clean reads were obtained after cleaning. Approximately 85.30% (367 million) of the clean reads were mapped to the quinoa reference genome and assembled into transcripts. Cleaned data (∼30 Gb) was first mapped to the reference genome, and then alignments were assembled into transcripts in each sample (Supplementary Table 2). Transcripts from all samples were merged to form a unified set of transcripts, which consisted of 188,663 genes and 234,387 transcripts. There were 146,440 transcripts (62.5% of the predicted transcripts) with only one exon. There were 1.24 transcripts per gene and 3.15 exons in a transcript.
Coding potential ability was analyzed by CPC2, from which there were 153,751 non-coding LncRNAs identified (Supplementary Table 3). A genome location study showed these LncRNAs were intensely distributed in 18 chromosomes of C. quinoa, and the distribution intensity in the two ends of each chromosome was higher than in other parts (Figure 1A). According to the relationship between a transcript and the closest reference transcript, LncRNAs were grouped into different classes represented by characters and symbols by GffCompare software (Figure 1B). “u” was the most abundant type (72.8%), showing that most LncRNAs came from the intergenic region. The second most abundant class was the “x” type (12.1%), followed by “i” (6.7%) and “p”(3.3%). These results indicated that LncRNAs seldom overlapped with reference transcripts. As shown in Figures 1C,D, the length of LncRNAs was generally much shorter than the coding RNAs, and LncRNA had less exon contained than the coding ones. The difference in expression pattern between coding RNAs and LncRNAs was measured statistically using a two-tailed Mann-Whitney U-Test (Figure 1E). The expression pattern of LncRNAs was significantly different from that of the coding genes. The overall expression level of LncRNA was significantly lower than that of coding RNAs at 0, 0.5, and 2 h of -salinity treatment (p < 2.2e-16). However, at 24 h of treatment, the overall expression level of LncRNAs increased sharply and was more significant than that of coding RNAs. The expression level of coding RNAs was relatively stable within 24 h of treatment. For further validating the reliability of LncRNAs, 50 randomly selected LncRNAs were aligned against Pacbio full-length cDNA datasets from . A 47 of 50 LncRNAs were successfully aligned into full-length transcripts. Among them, 46 full-length cDNAs corresponding to LncRNAs were also no coding potential by CPC2 program prediction. These cross-validation results revealed that a set of high-confidence LncRNAs was obtained in the study (Supplementary Table 8).
FIGURE 1
Identification of differentially expressed long non-coding RNAs in response to salt stress
Differentially expressed transcripts, including coding RNAs and LncRNAs, were then identified by DESeq2 (Supplementary Table 4). We identified 4,460 DE-LncRNAs, of which 214, 1,731, and 3,102 were identified at 0.5, 2, and 24 h of treatment (Figure 2A). Only 54 LncRNAs were differentially expressed at all the three-time points, occupying 1.2% of total DE-LncRNAs (4,460) (Figure 2B). Totally 6,791 DE-coding RNAs were also identified, 104 (1.5%) differentially expressed at three-time points (Figure 2B). At 0.5 h of treatment 75% (161) DE-LncRNAs and 77% (237) DE-coding RNAs were upregulated, while at 2 h of treatment 48% (831) DE-LncRNAs and 56% (2,197) DE-coding RNAs were upregulated. At 24 h of treatment, as high as 81% (2,513) DE-LncRNAs were upregulated, whereas only 39% (1,462) DE-coding RNAs were upregulated (Figure 2C and Supplementary Table 2).
FIGURE 2
Previous studies have suggested lncRNAs may play cis regulation role against neighboring genes. To investigate this possibility, we further measured the expression correlation between salt-responsive LncRNAs and their closest neighboring gene in either the 5′ or 3′ direction, yielding a dataset of 1,740 LncRNA-Coding Genes pairs (differentially expressed LncRNAs and their neighboring genes). The correlation analysis showed lncRNAs were strongly and highly significantly correlated with the expression of their closest neighboring gene (r = 0.346, p-value < 2.2e-16) (Figure 2D). Suggesting that lncRNAs may either be involved in cis-acting regulation or are subject to some of the same cis-acting regulatory features as their neighboring genes.
Construction of gene co-expression network and analysis of salinity responsive modules
To infer the potential biological functions of the LncRNAs, a weighted gene co-expression network consisting of both LncRNAs and coding RNAs based on expression profiles was constructed by the WGCNA program (). The soft-thresholding power was predicted and defined as 7, as it was the lowest power for which the scale-free topology fit index reached 0.90 (Figures 3A,B). There were finally 36 modules (Figure 3C) generated, and they were named from M1 to M36. The relationship between modules and salinity treatment was calculated, of which seven modules were highly relevant to salinity treatment (Figure 3D). M30 (r = 0.93, p = 1.57e-05) and M12 (r = 0.99, p = 9.12e-10) were upregulated significantly at 0.5 and 2 h, respectively; M17 (r = 0.90, p = 5.48e-05) and M32 (r = −0.99, p = 4.71e-09) modules upregulated and downregulated at 24 h, respectively; M13 (r = 0.90, p = 5.53e-05) upregulated at 0.5 and 2 h both; M11 (r = 0.90, p = 6.35e-05) upregulated at both 2 and 24 h, while M33 (r = 0.90, p = 6.35e-05) downregulated at the same time. No module was significantly regulated at all the three-time points. The percentage of LncRNAs in salinity-responsive modules ranged from 20 to 40%. The genes with the highest connectivity in each module were selected as hub genes. Totally, 210 hub genes were identified from these seven salinity-responsive modules listed above, which constituted a subnetwork (Figure 3E and Supplementary Table 6). The hub genes of M30 included both TF genes and LncRNAs, and so did M13, M12, and M11 modules. This implied that LncRNAs within these modules might interact with transcript factors and their role in salinity response. However, TF genes were not included in the hub genes of M17 and M32 modules. Instead, more than half of the hub genes were LncRNAs in them. In M17 module as high as 23 LncRNAs were at the hub position (Supplementary Table 6).
FIGURE 3
Gene ontology enrichment analysis of salinity-responsive modules
Gene ontology analysis of coding RNAs was performed to predict the function of the LncRNAs within the same module (Supplementary Table 5). The most representative GO terms of high-salinity responsive modules are shown in Figure 4A. A large part of the GO terms was enriched in all the seven high-salinity responsive modules. They were related to various critical biological processes, including biological regulation (GO:0065007), regulation of gene expression (GO:0010468), response to stress (GO:0006950), regulation of metabolic process (GO:0019222), transcription (GO:0006350) and developmental process (GO:0032502). Transport (GO:0006810) and metabolic process (GO:0008152) were enriched in six salinity-responsive modules. A small part of the GO terms was enriched in only two or three modules such as cell cycle (GO:0007049). Within these modules, there were several genes directly responding to salt stress. For instance, in M11, a transcript encoding a salt tolerance zinc finger (C2H2 type) which is highly homologous to the Arabidopsis gene (AT1G27730.1); in M12, there is a transcript encoding calcineurin B-like protein 1, which is highly homologous to Arabidopsis gene (AT4G17615.1), it might function as a positive regulator of salt and drought responses and as a negative regulator of cold response, and mediates the activation of AKT1 by CIPK proteins (CIPK6, CIPK16, and CIPK23) in response to low potassium conditions; In M13, one transcript homologous to Arabidopsis gene AT5G12010 was included, which might respond to salt stress and function in ABA-activated signaling pathway.
FIGURE 4
Validation of salt-responsive long non-coding RNAs by quantitative real-time PCR analysis
To validate the reliability of salt-responsive LncRNAs, 15 salt-responsive LncRNAs were randomly selected and then subjected the 0 and 2 h salt-treated samples to quantitative real-time PCR (qRT-PCR) to compare expression changes between replicated control and salt-treated. As a result, 11 of 15 salt-responsive LncRNAs were successfully detected and showed a high degree of consistency (r = 0.877, p = 0.000392) between RNA-seq and qRT-PCR (Figure 4B and Supplementary Table 7).
Discussion
Because of its remarkable tolerance to soil salt, drought and infertility, high-quality and well-balanced nutrient, quinoa has become popular and well-studied for its unique nature (; ; Zou et al., 2017). Quinoa is not only a species with high tolerance to salt stress, but also a plant preferring sodium (). Hence, it is acknowledged as halophytic plant species by some studies. Up to date, there have been dozens of studies carried out to dissect the characteristic. It has been demonstrated that the quinoa plant has epidermal bladder cells on its leaves and that these could pump extra sodium chloride into it (). Sodium cation could be sequestrated in leaf-cell vacuole (). The cellular potassium retention ability has been enhanced; on the other hand, cellular sodium exclusion and xylem loading in quinoa is also superior (). However, the salt tolerance trait of quinoa has a complicated underlying mechanism.
Plants have developed a series of physiological functions to alleviate stress or adapt to different environmental changes. The signal pathways and regulation systems lying behind are being dissected more and more. Non-coding RNAs have been illustrated to play important roles in the biological processes. These non-coding RNAs mainly include small non-coding RNA (sncRNA) and long non-coding (LncRNA). miRNA is one of the major parts of sncRNAs, as well as the siRNA; and they mediate gene silencing and hence regulate the expression levels of the genes (). Mul-miR3954, for example, was discovered in Arabidopsis to improve the salt-tolerance level of transgenic plants (); miR398b and miR298 could regulate Cu/Zn-SOD expression in response to ROS levels induced by salt stress (). TCONS_00009717 is a miRNA found in soybean that may be induced by salt stress, its potential target is cytochrome P450 (). Summarily, LncRNAs were also discovered and found to be involved in response to several abiotic stresses in plant species such as Arabidopsis thaliana, Zea may, and Nicotiana tabacum (; , ; Yu et al., 2019; Zhang et al., 2019, 2021; ). However, as a remarkable salt tolerant crop, the LncRNAs of quinoa responding to salt stress have not been reported yet.
In the study, we performed a time-course dynamic transcriptome analysis on quinoa under salinity treatment. Under salinity treatments, we extracted 234,387 transcripts from quinoa seedling roots. 153,751 LncRNAs and 51,667 genes were discovered among them. LncRNAs were generally shorter in length than coding RNA, which was likely due to the fact that nearly 90% of LncRNAs had only one exon (Figures 1C,D). In Figure 2B, time-course dynamic transcriptome analysis revealed a total of 4,460 identified DE-LncRNAs and 6,791 DE-coding RNAs and the DE-LncRNAs and DE-coding RNAs showed a time-dependent pattern, respectively. Meanwhile, as shown in Figure 2A, most of DE-LncRNAs were upregulated, which was probably why the overall expression level of LncRNAs increased enormously at 24 h of treatment (Supplementary Table 2).
Previous studies have suggested lncRNAs may play cis regulation role against neighboring genes. A dataset of 1,740 LncRNA-Coding Genes pairs (salt-responsive LncRNAs and their closest neighboring gene in either the 5′ or 3′ direction) was generated. Pearson correlation analysis showed highly significantly correlation level (r = 0.346, p-value < 2.2e-16) (Figure 2D), which illustrated that LncRNAs may play a cis-acting regulation role on their neighboring genes under salt stress.
The DE-LncRNAs and DE-coding RNAs were further examined, yielding 36 weighted gene co-expression network modules, seven of which showed responses to salinity stress. Highly connected hub genes in a module most likely played important roles in the same biological processes. The hub genes from each of the seven modules were shown as a subnetwork in Figure 3E. Some DE-coding RNAs are also classified as transcription factors (TF). For example, in module 11, one transcript encodes a homologous of Arabidopsis gene (ATG27730.1), which encodes a zinc finger TF involved in salt tolerance. It is also a hub gene in that module with high connectivity with others. In Arabidopsis, the homologous gene could respond to various kinds of stresses such as salt, cold, drought and oxidative. It is closely associated with the signal pathway of stress tolerance in plants. In modules 17 and 32, there were no TFs involved, which indicated that some LncRNAs might also work independently of TFs in response to salinity stress in quinoa. Another 48 TFs belonging to coding RNAs and 62 LncRNAs in this module implied that LncRNAs within these modules might interact with transcript factors and play a hub role in salinity response.
Conclusion
As a well-known salt-tolerant crop, few reports have identified long non-coding RNAs under salt stress in quinoa. Hence, we provided a bulk of LncRNAs in quinoa roots on a large scale and identified those induced or suppressed by salt treatment. Furthermore, we also predicted the potential gene-expression modules in which the LncRNAs might be included and function together with genes. Hopefully, these findings will serve as a dataset resource for further research on quinoa salt tolerance and provide a reference for quinoa breeding work.
Materials and methods
Plant material and salt treatment
Chenopodium quinoa cultivar QQ056 were acquired from the USDA-ARS National Plant Germplasm System (NPGS) with permission for scientific research. Detailed information on the variety could be found at npgsweb.ars-grin.gov (accession: PI 584524). The experiments were performed in a phytotron at 28°C under 16 h/8 h (day/night) photoperiod. Seeds of quinoa were surface sterilized, germinated for 7 days, and then moved into a half-strength Hoagland solution. Seedlings of 28-day-old were transferred to half-strength MS with 300 mM NaCl for salt treatment. Roots of quinoa were gathered at 0, 0.5, 2, and 24 h of treatment and stored at –80°C for further investigation. Three biological repeats were used at each time point.
Library preparation and RNA-seq sequencing
Total RNA was isolated, purified and concentrated with an RNAprep Pure Plant Kit (Tiangen, China). A Thermo 2000 Bioanalyzer evaluated the concentration and quality of RNA with an RNA NanoDrop (Thermo Scientific, United States). cDNA libraries were conducted using TruSeq mRNA Sample Prep Kit (Illumina, United States). RNA sequencing (RNA-seq) was performed on the Illumina NovaSeq 6000 platform (Illumina, United States) at Annoroad Gene Technology (Beijing) Co., Ltd., China.
RNA-seq data analysis
Raw data were trimmed with low-quality bases and short reads (< 50 bp) using Fastp v0.23.2 (). Cleaned reads were then mapped into the quinoa reference genome () by splice-aware alignment method using STAR v2.7.10a with two-pass mode (). The mapped reads of each sample were separately assembled by the reference annotation-based transcript (RABT) assembly algorithm and then combined with known transcript annotation () into an updated GTF file using StringTie v2.2.1 (). Finally, the abundance of transcripts was quantified and normalized with HTseq () and DESeq2 (). The transcripts with fragments per kilobase of transcript per million fragments mapped (FPKM) < 1 in more than three samples were excepted for downstream analysis.
Coding potential prediction of long non-coding RNAs
Computational prediction of quinoa LncRNAs was followed as described by with some custom modification. Sequences of transcripts were firstly retrieved using Gffread v0.12.2 (). Then transcripts were evaluated in their coding potential using CPC v2.0 (). Default parameters were used. Finally, non-coding transcripts longer than 200 bp were considered as LncRNAs. The updated GTF file was also compared with the reference GTF file using gffCompare () to generate class codes representing the position information between the updated transcript and the closest reference transcript.
Differential expression analysis
For differential expression analysis, we compared every time point (0.5, 2, and 24 h of treatment) with CK (0 h) using DESeq2 () based on the negative binomial distribution. Differentially expressed should be fulfilled the following criteria: (I) fold change > = 1; (II) adjusted p value < 0.05.
Gene co-expression network analysis
A weighted co-expression network was further constructed for linking coding RNAs and LncRNAs using the WGCNA program (Zhang and Horvath, 2005; ). The correlation between module eigengenes and salt treatment was then calculated. Modules that showed a significant correlation (| r| > 0.9, p < 0.001) with a specific time of treatment were recognized as salt-responsive modules (Xue et al., 2013). The hub genes of each module were worked out based on the Topological overlap matrix (TOM).
Gene ontology enrichment analysis
Gene Ontology (GO) enrichment analysis was performed with agriGO v2.0 toolkit (). Fisher’s exact test was applied for the enrichment analysis, and the false discovery rate (FDR) was assessed using the Yekutieli method. GO terms with an FDR less than 0.05 were considered significant.
Quantitative real-time PCR validation
To validate the reliability and accuracy of the LncRNA analysis, ten salt-responsive LncRNAs were chosen and measured by quantitative real-time PCR (qRT-PCR). The first-strand cDNA of 0 and 2 h samples were synthesized using the PrimeScript TM RT-PCR Kit (Takara, Japan). We conducted the qRT-PCR on an ABI StepOnePlus Real-Time PCR System (Applied Biosystems, CA, United States) with FastStart Universal SYBR Green Master (Roche, Germany) according to the manufacturer’s protocol. The primers were designed using Primer Premier v.5.0 software (Premier Biosoft International, CA, United States). CqEF1a (AUR62020767) was used as an internal standard to normalize the relative expression level and determine expression values based on the 2–ΔΔCt method. The primers for qRT-PCR were presented in Supplementary Table 7.
Statements
Data availability statement
The datasets presented in this study were deposited in National Center for Biotechnology Information (NCBI) BioProject database under accession number: PRJNA636120. All data can be found in the article/Supplementary material.
Author contributions
YL and XC conceived and designed the experiments. CL, BH, PS, JX, HG, BP, JC, and FH participated in experiments and data analyses. CL, BH, and PS wrote the manuscript with inputs and guidance from YL. All authors have read and approved the final manuscript.
Funding
This work was supported by the Jiangsu Agricultural Science and Technology Independent Innovation Fund (JASTIF) [CX(21)3127] and National Natural Science Foundation of China (31601818).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.988845/full#supplementary-material
References
1
BöhmJ.MessererM.MüllerH. M.Scholz-StarkeJ.GradognaA.ScherzerS.et al (2018). Understanding the molecular basis of salt sequestration in epidermal bladder cells of Chenopodium quinoa.Curr. Biol.283075.e–3085.e. 10.1016/j.cub.2018.08.004
2
ChenR.LiM.ZhangH.DuanL.SunX.JiangQ.et al (2019). Continuous salt stress-induced long non-coding RNAs and DNA methylation patterns in soybean roots.BMC Genomics20:730. 10.1186/s12864-019-6101-7
3
ChenS.ZhouY.ChenY.GuJ. (2018). fastp: An ultra-fast all-in-one FASTQ preprocessor.Bioinformatics34i884–i890. 10.1093/bioinformatics/bty560
4
ChenX.MengL.HeB.QiW.JiaL.XuN.et al (2022). Comprehensive transcriptome analysis uncovers hub long non-coding RNAs regulating potassium use efficiency in Nicotiana tabacum.Front. Plant Sci.13:777308. 10.3389/fpls.2022.777308
5
DengF.ZhangX.WangW.YuanR.ShenF. (2018). Identification of Gossypium hirsutum long non-coding RNAs (lncRNAs) under salt stress.BMC Plant Biol.18:23. 10.1186/s12870-018-1238-0
6
DobinA.DavisC. A.SchlesingerF.DrenkowJ.ZaleskiC.JhaS.et al (2013). STAR: Ultrafast universal RNA-seq aligner.Bioinformatics2915–21. 10.1093/bioinformatics/bts635
7
FengX.-M.QiaoY.MaoK.HaoY.-J. (2010). Ectopic overexpression of AtmiR398b gene in tobacco influences seed germination and seedling growth.Plant Cell Tissue Organ. Cult.10253–59. 10.1007/s11240-010-9704-x
8
GaiY.-P.YuanS.-S.ZhaoY.-N.ZhaoH.-N.ZhangH.-L.JiX.-L. (2018). A novel LncRNA, MuLnc1, associated with environmental stress in mulberry (Morus multicaulis).Front. Plant Sci9:669. 10.3389/fpls.2018.00669
9
JacobsenS. E.MujicaA.JensenC. R. (2003). The resistance of quinoa (Chenopodium quinoa Willd.) to adverse abiotic factors.Food Rev. Int.1999–109. 10.1081/FRI-120018872
10
JacobsenS.-E.MonterosC.ChristiansenJ. L.BravoL. A.CorcueraL. J.MujicaA. (2005). Plant responses of quinoa (Chenopodium quinoa Willd.) to frost at various phenological stages.Eur. J. Agron.22131–139. 10.1016/j.eja.2004.01.003
11
JannesarM.SeyediS. M.Moazzam JaziM.NiknamV.EbrahimzadehH.BotangaC. (2020). A genome-wide identification, characterization and functional analysis of salt-related long non-coding RNAs in non-model plant Pistacia vera L. using transcriptome high throughput sequencing.Sci. Rep.10:5585. 10.1038/s41598-020-62108-6
12
JarvisD. E.HoY. S.LightfootD. J.SchmöckelS. M.LiB.BormT. J. A.et al (2017). The genome of Chenopodium quinoa.Nature542307–312. 10.1038/nature21370
13
KangY.-J.YangD.-C.KongL.HouM.MengY.-Q.WeiL.et al (2017). CPC2: A fast and accurate coding potential calculator based on sequence intrinsic features.Nucleic Acids Res.45W12–W16. 10.1093/nar/gkx428
14
KoziołM. J. (1992). Chemical composition and nutritional evaluation of quinoa (Chenopodium quinoa Willd.).J. Food Compos. Anal.535–68. 10.1016/0889-1575(92)90006-6
15
KumarN.BharadwajC.SahuS.ShivA.ShrivastavaA. K.ReddyS. P. P.et al (2021). Genome-wide identification and functional prediction of salt- stress related long non-coding RNAs (lncRNAs) in chickpea (Cicer arietinum L.).Physiol. Mol. Biol. Plants272605–2619. 10.1007/s12298-021-01093-0
16
KungJ. T. Y.ColognoriD.LeeJ. T. (2013). Long noncoding RNAs: Past, present, and future.Genetics193651–669. 10.1534/genetics.112.146704
17
LangfelderP.HorvathS. (2008). WGCNA: An R package for weighted correlation network analysis.BMC Bioinform.9:559.
18
LiuJ.JungC.XuJ.WangH.DengS.BernadL.et al (2012). Genome-wide analysis uncovers regulation of long intergenic noncoding RNAs in Arabidopsis.Plant Cell244333–4345. 10.1105/tpc.112.102855
19
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.Genome Biol.15:550. 10.1186/s13059-014-0550-8
20
LvY.HuF.ZhouY.WuF.GautB. S. (2019). Maize transposable elements contribute to long non-coding RNAs that are regulatory hubs for abiotic stress response.BMC Genomics20:864. 10.1186/s12864-019-6245-5
21
LvY.LiangZ.GeM.QiW.ZhangT.LinF.et al (2016). Genome-wide identification and functional prediction of nitrogen-responsive intergenic and intronic long non-coding RNAs in maize (Zea mays L.).BMC Genomics17:350. 10.1186/s12864-016-2650-1
22
LvY.ZhaoL.XuX.WangL.WangC.ZhangT.et al (2013). Characterization of expressed sequence tags from developing fibers of Gossypium barbadense and evaluation of insertion-deletion variation in tetraploid cultivated cotton species.BMC Genomics14:170. 10.1186/1471-2164-14-170
23
MaQ.SuC.DongC.-H. (2021). Genome-wide transcriptomic and proteomic exploration of molecular regulations in quinoa responses to ethylene and salt stress.Plants10:2281. 10.3390/plants10112281
24
MaughanP. J.TurnerT. B.ColemanC. E.ElzingaD. B.JellenE. N.MoralesJ. A.et al (2009). Characterization of salt overly sensitive 1 (SOS1) gene homoeologs in quinoa (Chenopodium quinoa Willd.).Genome52647–657. 10.1139/G09-041
25
MoralesA.Zurita-SilvaA.MaldonadoJ.SilvaH. (2017). Transcriptional responses of Chilean quinoa (Chenopodium quinoa Willd.) under water deficit conditions uncovers ABA-independent expression patterns.Front. Plant Sci8:216. 10.3389/fpls.2017.00216
26
NegiP.MishraS.GanapathiT. R.SrivastavaA. K. (2021). Regulatory short RNAs?: A decade’s tale for manipulating salt tolerance in plants.Physiol. Plant1731535–1555. 10.1111/ppl.13492
27
PerteaG.PerteaM. (2020). GFF utilities: GffRead and GffCompare.F1000Res.9:304. 10.12688/f1000research.23297.1
28
PerteaM.KimD.PerteaG. M.LeekJ. T.SalzbergS. L. (2016). Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown.Nat. Protoc.111650–1667. 10.1038/nprot.2016.095
29
PontingC. P.OliverP. L.ReikW. (2009). Evolution and functions of long noncoding RNAs.Cell136629–641. 10.1016/j.cell.2009.02.006
30
PutriG. H.AndersS.PylP. T.PimandaJ. E.ZaniniF. (2022). Analysing high-throughput sequencing data in Python with HTSeq 2.0.Bioinformatics382943–2945. 10.1093/bioinformatics/btac166
31
QinT.ZhaoH.CuiP.AlbesherN.XiongL. (2017). A nucleus-localized long non-coding RNA enhances drought and salt stress tolerance.Plant Physiol.1751321–1336. 10.1104/pp.17.00574
32
RazzaghiF.AhmadiS. H.AdolfV. I.JensenC. R.JacobsenS.-E.AndersenM. N. (2011). Water relations and transpiration of quinoa (Chenopodium quinoa Willd.) Under Salinity and Soil Drying.J. Agron. Crop Sci.197348–360. 10.1111/j.1439-037X.2011.00473.x
33
RinnJ. L.ChangH. Y. (2012). Genome regulation by long noncoding RNAs.Annu. Rev. Biochem.81145–166. 10.1146/annurev-biochem-051410-092902
34
SchmöckelS. M.LightfootD. J.RazaliR.TesterM.JarvisD. E. (2017). Identification of putative transmembrane proteins involved in salinity tolerance in Chenopodium quinoa by integrating physiological data, RNAseq, and SNP analyses.Front. Plant Sci.8:1023. 10.3389/fpls.2017.01023
35
ShabalaS.MackayA. (2011). Ion transport in halophytes.Adv. Bot. Res.57151–199. 10.1016/B978-0-12-387692-8.00005-9
36
ShiP.GuM. (2020). Transcriptome analysis and differential gene expression profiling of two contrasting quinoa genotypes in response to salt stress.BMC Plant Biol.20:568. 10.1186/s12870-020-02753-1
37
TanX.LiS.HuL.ZhangC. (2020). Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting rapeseed (Brassica napus L.) genotypes subjected to drought stress and re-watering.BMC Plant Biol.20:81. 10.1186/s12870-020-2286-9
38
TianT.LiuY.YanH.YouQ.YiX.DuZ.et al (2017). agriGO v2.0: A GO analysis toolkit for the agricultural community, 2017 update.Nucleic Acids Res.45W122–W129. 10.1093/nar/gkx382
39
TrapnellC.WilliamsB. A.PerteaG.MortazaviA.KwanG.van BarenM. J.et al (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation.Nat. Biotechnol.28511–515. 10.1038/nbt.1621
40
van ZelmE.ZhangY.TesterinkC. (2020). Salt tolerance mechanisms of plants.Annu. Rev. Plant Biol.71403–433. 10.1146/annurev-arplant-050718-100005
41
Vega-GálvezA.MirandaM.VergaraJ.UribeE.PuenteL.MartínezE. A. (2010). Nutrition facts and functional potential of quinoa (Chenopodium quinoa willd.), an ancient Andean grain: A review.J. Sci. Food Agric.902541–2547. 10.1002/jsfa.4158
42
WuG.PetersonA. J.MorrisC. F.MurphyK. M. (2016). quinoa seed quality response to sodium chloride and sodium sulfate salinity.Front. Plant Sci.7:790. 10.3389/fpls.2016.00790
43
XueZ.HuangK.CaiC.CaiL.JiangC. Y.FengY.et al (2013). Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing.Nature500593–597. 10.1038/nature12364
44
YuY.ZhangY.ChenX.ChenY. (2019). Plant noncoding RNAs: Hidden players in development and stress responses.Annu. Rev. Cell Dev. Biol.35407–431. 10.1146/annurev-cellbio-100818-125218
45
ZhangB.HorvathS. (2005). A general framework for weighted gene co-expression network analysis.Stat. Appl. Genet. Mol. Biol4:17. 10.2202/1544-6115.1128
46
ZhangX.DongJ.DengF.WangW.ChengY.SongL.et al (2019). The long non-coding RNA lncRNA973 is involved in cotton response to salt stress.BMC Plant Biol.19:459. 10.1186/s12870-019-2088-0
47
ZhangX.ShenJ.XuQ.DongJ.SongL.WangW.et al (2021). Long noncoding RNA lncRNA354 functions as a competing endogenous RNA of miR160b to regulate ARF genes in response to salt stress in upland cotton.Plant Cell Environ.443302–3321. 10.1111/pce.14133
48
ZhangY.-C.LiaoJ.-Y.LiZ.-Y.YuY.ZhangJ.-P.LiQ.-F.et al (2014). Genome-wide screening and functional analysis identify a large number of long noncoding RNAs involved in the sexual reproduction of rice.Genome Biol.15:512. 10.1186/s13059-014-0512-1
49
ZhaoS.ZhangQ.LiuM.ZhouH.MaC.WangP. (2021). Regulation of plant responses to salt stress.Int. J. Mol. Sci.22:4609. 10.3390/ijms22094609
50
ZhuJ.-K. (2002). Salt and drought stress signal transduction in plants.Annu. Rev. Plant Biol.53247–273. 10.1146/annurev.arplant.53.091401.143329
51
ZouC.ChenA.XiaoL.MullerH. M.AcheP.HabererG.et al (2017). A high-quality genome assembly of quinoa provides insights into the molecular basis of salt bladder-based salinity tolerance and the exceptional nutritional value.Cell Res.271327–1340. 10.1038/cr.2017.124
Summary
Keywords
quinoa, LncRNA, salt stress, RNA-seq, gene co-expression
Citation
Luo C, He B, Shi P, Xi J, Gui H, Pang B, Cheng J, Hu F, Chen X and Lv Y (2022) Transcriptome dynamics uncovers long non-coding RNAs response to salinity stress in Chenopodium quinoa. Front. Plant Sci. 13:988845. doi: 10.3389/fpls.2022.988845
Received
07 July 2022
Accepted
24 August 2022
Published
20 September 2022
Volume
13 - 2022
Edited by
Jian Chen, Jiangsu University, China
Reviewed by
Sidra Tul Muntha, Kohsar University Murree, Pakistan; Bo Wei, Ghent University, Belgium
Updates
Copyright
© 2022 Luo, He, Shi, Xi, Gui, Pang, Cheng, Hu, Chen and Lv.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xi Chen, chenxi_201402@163.comYuanda Lv, Lyd0527@126.com
†These authors have contributed equally to this work
This article was submitted to Plant Bioinformatics, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.