Abstract
Tea (Camellia sinensis) is one of the significant cash crops in China. As a leaf crop, nitrogen supply can not only increase the number of new shoots and leaves but also improve the tenderness of the former. However, a conundrum remains in science, which is the molecular mechanism of nitrogen use efficiency, especially long non-coding RNA (lncRNA). In this study, a total of 16,452 lncRNAs were identified through high-throughput sequencing analysis of lateral roots under nitrogen stress and control conditions, of which 9,451 were differentially expressed lncRNAs (DE-lncRNAs). To figure out the potential function of nitrogen-responsive lncRNAs, co-expression clustering was employed between lncRNAs and coding genes. KEGG enrichment analysis revealed nitrogen-responsive lncRNAs may involve in many biological processes such as plant hormone signal transduction, nitrogen metabolism and protein processing in endoplasmic reticulum. The expression abundance of 12 DE-lncRNAs were further verified by RT-PCR, and their expression trends were consistent with the results of RNA-seq. This study expands the research on lncRNAs in tea plants, provides a novel perspective for the potential regulation of lncRNAs on nitrogen stress, and valuable resources for further improving the nitrogen use efficiency of tea plants.
Introduction
As a crucial element in plant growth, nitrogen is an essential component of biological compounds such as nucleic acid, protein, chlorophyll, and plant hormones. Nitrogen also regulates plants’ absorption and utilization of other elements, such as phosphorus and potassium (Novoa and Loomis, 1981; Xu et al., 2012). In agricultural production, a large amount of synthetic nitrogen fertilizer has been used in crop fields to promote the yield of crops (; Zhang et al., 2012). Nitrogen fertilizer abuse increases agricultural costs and becomes one of the chief culprits of environmental pollution. (). Therefore, reducing the application of nitrogen fertilizer and improving the nitrogen use efficiency of crops are momentous to agriculture development. To this end, more and more researchers have carried out studies on crops with low nitrogen tolerance, providing new solutions for the underlying molecular mechanisms of nitrogen regulation and the improvement of nitrogen use efficiency (Kant et al., 2011; Krapp, 2015; ). Plant lateral roots (LRs) can absorb nitrogen in the form of organic (amino acids and peptides) and inorganic (nitrate and ammonium nitrogen), and different nitrogen sources can regulate the growth of LRs (). The number and location of lateral roots determine the spatial configuration of roots, which means that their formation and development will directly affect plant growth (). Nitrate use was reported to be positively associated with the development of LRs (Lynch, 2013). In other words, reducing the application of nitrate fertilizer can promote the growth of LRs, thereby improving the plant’s nutrient uptake capacity and nitrogen use efficiency ().
Long non-coding RNA (lncRNA) is an RNA transcript longer than 200 nucleotides with no or limited protein-coding capacity (Kung et al., 2013; Kopp and Mendell, 2018). According to their genomic location, they are usually divided into i (intron lncRNAs), o (overlapping lncRNAs), u (intergenic lncRNAs), and x (antisense lncRNA) types (Roberts et al., 2011). Much evidence, including epigenetics, transcriptional and post-transcriptional regulation in the form of RNA, can prove that lncRNAs have regulatory functions in gene expression (Sunkar et al., 2007; ). Through next-generation sequencing technology and bioinformatics methods, many lncRNAs have been found in Arabidopsis, wheat, maize, rice, and other model plants. It is significant in flowering regulation, photomorphogenesis, stress response, and other growth pathways. (Wang et al., 2014; ; Wang and Chekanova, 2017; Zhao et al., 2018; Liang et al., 2022). For example, through genome-wide analysis of Arabidopsis thaliana full-length cDNA database, 76 ncRNAs were identified, including 5 small interfering RNAs (siRNA) precursors and 14 natural antisense transcripts of protein-coding genes. A set of 127 RNA sequencing samples including total RNAseq datasets and PacBio fl-cDNA datasets in maize was used for identifying 1,077 differentially multiple abiotic stress-responsive TE-lncRNAs, and 39 are hubs in co-expression networks, including a small number that are evolutionary conserved (Lv et al., 2019). In the study of low nitrogen stress in barley, 498 lncRNAs were identified, of which 487 were newly discovered, and 56 lncRNAs responsive to low nitrogen stress were identified (). 637 lncRNAs responsive to nitrogen and 664 lncRNAs responsive to drought were identified in maize seedlings (Zhang et al., 2014; Lv et al., 2016). However, studies of lncRNAs in many non-model plants are relatively limited, and correspondingly, studies on the genome-wide identification and analysis of tea plants lncRNAs are relatively few. Besides, there have been no published studies on the molecular mechanism of lncRNA response to nitrogen in tea plants.
In China, tea (Camellia sinensis) is one of the vital cash crops and stands out among traditional industries. As a leaf crop, tea plants consume a lot of nitrogen nutrition and have high requirements for soil nitrogen (Venkatesan et al., 2004; Zhu et al., 2014). The results showed that nitrogen nutrition could promote the germination and elongation of tea shoots, increase the number, weight and area of new shoots and leaves, and advance the tenderness of new shoots (Ruan et al., 2019). Nitrogen can benefit tea plants’ vegetative growth and thwart their reproductive growth, thus increasing yield (Ruan et al., 2019). Within a reasonable dosage range, applying nitrogen in tea gardens can enrich the types of tea aroma substances, improve the freshness of tea leaves, and increase the content of amino acids, tea polyphenols, catechins and chlorophyll, especially free amino acids. (Ruan et al., 2010; Lin et al., 2021).Therefore, improving nitrogen absorption and utilization efficiency has been a focal point of tea nutrition research in recent years. Screening and breeding tea varieties with high nitrogen utilization is paramount to the tea industry (). Most previous studies focused on revealing the coding genes regulated by nitrogen (Xia et al., 2017; Wei et al., 2018). Thanks to the release of the tea genome, we used High-Throughput Sequencing technology to analyze the lncRNAs expression profile of tea roots under high and low nitrogen conditions, analyzed these lncRNAs regulatory coding genes, and identified the function of lncRNAs involved in nitrogen metabolism. In conclusion, our research results provide many valuable references for clarifying the response mechanism of lncRNAs to nitrogen in LRs of tea plants and expand a new path for improving the nitrogen use efficiency of tea plants.
Materials and methods
Plant materials and nitrogen treatments
Purebred tea seeds (Camellia sinensis cv. Fuding dabai) were germinated in perlite and then cultured in nutrient solution (0.75mM (NH4)2SO4, 0.25 mM Ca(NO3)2•4(H2O)3, 0.05 mM KH2PO4, 0.35 mM K2SO4, 0.395 mM CaCl2, 0.21 mM MgSO4, 35.0 μM NaFeEDTA, 46.1 μM H3BO3, 2.0 μM MnSO4, 0.3 μM CuSO4, 2.0 μM ZnSO4 and 0.5 μM Na2MoSO4). Nitrogen was divided into three concentrations: 0.25 mM (labeled as low nitrogen, LN), 1 mM (labelled as control, CK), and 2.5 mM (labelled as high nitrogen, HN). The control experiment (1 mM) was supplemented with 0.75 mM ammonium and 0.25 mM nitrate using (NH4)2SO4 and Ca (NO3)2•4(H2O), respectively, which was the best combination of N concentration for seedlings growth. The growth conditions of seedlings in the light culture box were as follows: 28/25°C (day/night), 75% relative humidity, 16/8 h (light/darkness) photoperiod, and 300 μmol−2 s−1 light intensity. The liquid culture medium was changed every 5 days. The lateral roots of the seedlings were sampled after 10 weeks of seedling growth. The control and nitrogen treatments were repeated three times (CK-1, CK-2, CK-3; LN-1, LN-2, LN-3, HN-1, HN-2, HN-3). The samples were immediately soaked in liquid nitrogen and stored in a refrigerator of -80°C for RNA-seq analysis and qRT-PCR verification.
RNA isolation, library construction and RNA sequencing
The total RNA of LH, CK, HN was isolated with the plant RNA extraction kit with DNase (TIANDZ, Inc., Beijing, China) according to the manufacturer’s protocol. RNA was purified and concentrated using NanoDrop2000 Spectrophotometer (Thermo Fisher Scientific, USA), 1.2% agarose gel electrophoresis, and Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA). High-quality RNA samples are used for library construction, and then the total RNA-seq library is constructed and sequenced using the IlluminaHiSeq platform. All sequencing data were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (accession number PRJNA595712).
Transcriptome assembly
The SolexaQA++ v3.1 program was applied to execute quality trimming using the Q30 value (). After removing rRNA, low-quality reads, aptamer sequences, and contaminating reads, the remaining clean reads were aligned with the reference genome of C. sinensis var. sinensis by Hierarchical Indexing of Spliced Transcript Alignment (HISat) software (Pertea et al., 2016; Xia et al., 2017). Use the gffCompare program to annotate the assembled transcript and the unknown transcript and then filter out possible lncRNAs by pfam databases (Kong et al., 2007; ; ).
Identification of LncRNAs
In this study, a strict calculation method was used to determine the lncRNA of tea plants (; Xiao et al., 2015). The primary screening of transcripts should meet the following conditions: with a class code of “i”, “x”, “u”, “o”, and “e”, a length ≥ 200 bp, and fragments per kilobase of transcript per million mapped reads (FPKM) value ≥ 0.1(Kelley and Rinn, 2012; Luo et al., 2022). Subsequently, transcripts were compared with uniref90 and Pfam protein databases using the CPC2 program to evaluate their protein-coding potential (Kang et al., 2017). Transcripts that meet these conditions are eventually considered as candidates for lncRNAs for further analysis: Non-coding transcripts larger than 200 bp; FPKM > 1; a CPC score < -1.
Differential expression analysis
The expression level of lncRNAs was quantitatively detected by FPKM value by StringTie software. DESeq2 software package (1.10.1) was used to analyze the differential expression between nitrogen concentration and control treatments. The P value of the result is adjusted by the methods of Benjamini and Hochberg to control the false discovery rate. The lncRNA with the adjustment value P < 0.01 and the absolute value |log2(FPKM) ratio| ≥ 1 by DESeq is considered to be differentially expressed (Wan et al., 2020).
KEGG pathway enrichment analysis
To elucidate the potential functions of differentially expressed lncRNAs (DE-lncRNAs), co-expression analysis between Nitrogen-responsive lncRNAs and coding genes was employed for generating clusters with different expression patterns. Coding genes of clusters were used for the KEGG pathway enrichment by clusterProfiler (Yu et al., 2012). In the KEGG enrichment analysis, a false discovery rate ≤ 0.05 was used as a criterion to identify significantly enriched pathways ().
Validation and quantification of LncRNAs
Total RNA was extracted from the nine samples (CK-1, CK-2, CK-3; LN-1, LN-2, LN-3, HN-1, HN-2, HN-3) using Trizol reagent (Invitrogen). To validate the lncRNAs, Real-time PCR was performed Bio-Rad Real-time thermal cycler CFX96 with SYBR Premix ExTaq™ Kit (Takara Co. Ltd., Japan). The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and β-actin were used as controls. The 2-ΔΔCt method was used to determine the relative expression levels (Livak and Schmittgen, 2001). Three biological replicates per sample. All qRT-PCR primers were designed with NCBI primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/). Differences between groups were analyzed by one-way analysis of variance and Duncan’s test, and P < 0.05 was considered to be significantly different between treatment groups.
Result
Phenotypic response of C.sinensis lateral roots to nitrogen treatments
In different nitrogen treatments, the growth of the underground part of the seedlings changed significantly in the tenth week (Figure 1A). The main results were as follows: with the increase in nitrogen concentration, the number and the length of lateral roots decreased significantly (Figure 1B, C). This showed that low nitrogen conditions could promote the growth of lateral roots of seedlings. On the contrary, it was inhibited under high nitrogen conditions.
Figure 1
Genome−wide prediction of LncRNAs candidates in C. sinensis
To investigate the lncRNAs in tea plants under nitrogen treatments, transcriptome sequencing analysis of three biological repeats of HN, CK and HN was carried out using the Illumina HiSeq2500 platform. Each sample had an average base of 9GB, and the lowest value of Q30 was 92.83% (Table S1). All clean reads were compared with the tea plants reference genome. A total of 86814 transcripts were identified, and the results showed that 38% (33515) of transcripts matched exactly and entirely with the intron chain and 30% (26279) were annotated as multi-exon with at least one junction match, while 7% (6049) transcripts were annotated as containment of reference (reverse containment) (Figure 2A). Total RNA-seq technology could recover transcripts located in gene intergenic or intronic regions. After filtering the data through CPC2 and Pfam, a total of 16,452 candidate transcripts were obtained.
Figure 2
DE-lncRNAs under nitrogen treatments
To identify nitrogen-responsive lncRNAs in LRs of tea plants, the expression of lncRNAs was compared between the LN, CK and HN treatments. Based on the standardized FPKM value, 9,451 DE-lncRNAs were revealed in LN, CK, and HN treatments (Figure 2B). Moreover, 357 DE-lncRNAs overlapped in those treatments (Figure 2C). Furthermore, in the comparison of LN and CK treatments, there were 6,455 DE-lncRNAs, of which 3,805 were up-regulated and 2,650 were down-regulated (Figure 3A); in the comparison of HN and CK treatments, there were 4,844 DE-lncRNAs, of which 2,879 were up-regulated and 1,965 were down-regulated (Figure 3B); and 3,431 DE-lncRNAs were identified in the comparison of HN and LN, of which 1,780 were up-regulated and 1,651 were down-regulated (Figure 3C).
Figure 3
Enrichment analysis of nitrogen-responsive coding genes
The nitrogen-responsive coding genes were submitted into KEGG pathway enrichment, and the top 20 enriched pathways were shown in (Figure 3D-F). In the LN vs. CK comparisons, 2,161 coding genes were annotated in 128 KEGG pathways. Among them, the most enriched pathway was metabolic pathways (447), followed by biosynthesis of secondary metabolites (248), photosynthesis-antenna proteins (61) and starch and sucrose metabolism (35) (Figure 3D, Table S2). In addition, 1,514 coding genes were enriched in 116 KEGG pathways between HN vs. CK comparisons. Coding genes identified in metabolic pathways (327) were the most abundant, followed by biosynthesis of secondary metabolites (197), protein processing in endoplasmic reticulum (42) and amino sugar and nucleotide sugar metabolism (28) (Figure 3E, Table S3). Among the LN vs. HN comparisons, 1,140 coding genes were verified in 113 KEGG pathways. Moreover, the most significant pathway was metabolic pathways (242), followed by Biosynthesis of secondary metabolites (151), plant-pathogen interaction (36) and glutathione metabolism (20) (Figure 3F, Table S4).
Co-expression clustering analysis of nitrogen-responsive LncRNAs and coding genes in lateral roots development
In different nitrogen treatments, the lncRNAs with the same trend of coding genes expression were selected and annotated by KEGG (Figure 4). In LN treatment, compared with CK and HN treatments, the up-regulated and down-regulated DE-lncRNAs were clustered into cluster1 and cluster2, respectively, and then these up-regulated and down-regulated DE-lncRNAs were annotated to cluster 1 and cluster 2 KEGG pathways, severally. The down-regulated DE-lncRNAs compared with CK and HN was shown in cluster1. The cluster 1 KEGG pathway identified that 720 DE-lncRNAs were annotated in 104 KEGG pathways. The pathway with the most aggregation was ribosome (65), followed by protein processing in endoplasmic reticulum (26) and RNA transport (20) (Table S5); the up-regulated DE-lncRNAs were shown in cluster 2, and KEGG pathway (cluster 2) identified that 1311 DE-lncRNAs were annotated in 122 KEGG pathways, such as metabolic pathways (275), biosynthesis of secondary metabolites (145) and plant-pathogen interaction (32) (Table S6). In CK treatment, compared with LN and HN treatments, the up-regulated and down-regulated DE-lncRNAs were clustered into cluster 3 and 4, respectively, and then these up and down-regulated DE-lncRNAs were annotated to cluster 3 and cluster 4 KEGG pathways, severally. The down-regulated DE-lncRNAs were shown in cluster 3. The KEGG pathway (cluster 3) identified that 971 DE-lncRNAs were annotated in 121 KEGG pathways, such as metabolic pathways (202), biosynthesis of secondary metabolites (114) and carbon metabolism (38) (Table S7); the up-regulated DE-lncRNAs were shown in cluster 4, and KEGG pathway (cluster 4) identified that 815 DE-lncRNAs were annotated in 116 KEGG pathways, such as biosynthesis of secondary metabolites (99), biosynthesis of amino acids (25), plant hormone signal transduction (24) (Table S8). In HN treatment, compared with LN and CK treatments, the up-regulated and down-regulated DE-lncRNAs were clustered into cluster 5 and cluster 6, respectively, and then these up and down-regulated DE-lncRNAs were annotated to cluster 5 and cluster 6 KEGG pathways, severally. The up-regulated DE-lncRNAs were shown in cluster 5, and KEGG pathway (cluster 5) identified that 993 DE-lncRNAs were annotated in 114 KEGG pathways, such as protein processing in endoplasmic reticulum (45), ribosome (30) and RNA transport (25) (Table S9); the down-regulated DE-lncRNAs were shown in cluster 6, and KEGG pathway (cluster 6) identified that 556 DE-lncRNAs were annotated in 102 KEGG pathways, such as metabolic pathways (114), biosynthesis of secondary metabolites (74), and plant hormone signal transduction (25) (Table S10). The clustering results showed that nitrogen treatment could affect the expression of lncRNAs, associated with the nitrate transporter and controlling the growth of tea plants LRs. In addition, nitrogen treatment regulated the hormone signals of tea plants LRs, associated with auxin synthesis and auxin signals genes.
Figure 4
Validation of lncRNAs expression using qRT-PCR
To determine the reliability of RNA-seq transcriptome results, 12 DE-lncRNAs were selected for qRT-PCR verification. The result showed that the expression level of these DE-lncRNAs was in line with that of transcripts estimated from sequence data, suggesting the repeatability and accuracy of RNA-seq data (Figure 5). The identified DE-lncRNAs sequences were recorded in Table S11.
Figure 5
Discussion
Improving nitrogen use efficiency is of great significance in agricultural production and environmental protection.To promote nitrogen use efficiency and adapt to barren land, improving the resistance mechanism to low nitrogen or nitrogen deficiency is advisable. (). Previous studies have shown that the early nodulin gene was found in rice, which was identified by transcriptional analysis under low nitrogen stress and was related to increased nitrogen use efficiency (). At the same time, barley is becoming a model crop for similar research because of its adaptability to barren land (). Due to the particularity system of tea plants, nitrogen use efficiency is critical in recent research (Zhang et al., 2020). It includes nitrogen absorption efficiency, utilization efficiency and transport efficiency.
During past years, research on improving the nitrogen use efficiency of tea plants has been the top priority. In addition, nitrogen can also regulate the synthesis and transport of auxin, gibberellin, and other plant hormones through nitrogen signals, thus affecting the growth and development of plant LRs (Sun et al., 2017). In previous studies, nitrogen deficiency induced the expression of high-affinity nitrate transporters NRT1, NRT2.4 and NRT2.5 in Arabidopsis thaliana roots (Lezhneva et al., 2014). Moreover, NRT1 and NRT2 can transport not only but also auxin. So, under a low nitrogen environment, NRT1 and NRT2 are highly expressed and regulate the formation of LRs by inducing auxin accumulation and transport (Krouk et al., 2010). In this study, low nitrogen could also promote the growth of tea plants’ LRs. However, accumulating evidence shows that lncRNAs are momentous in various plants’ biological processes, but their characteristics, expression patterns and potential functions in response to nitrogen in tea plants have not been reported. In this study, we identified the lncRNAs of tea plants according to the latest genome of tea plants, and confirmed the lncRNAs related to nitrogen response and LRs development. This study has gained discoveries of the regulatory function of lncRNAs on nitrogen in tea plants and provides a resource for further study of the molecular mechanism of nitrogen response.
In previous studies, 7,245 lncRNA were identified in RNA-seq analysis under nitrogen-deficiency treatment in maize, of which 673 were corresponding to low nitrogen deficiency (Lv et al., 2016). 918 differentially expressed lncRNAs was found in the study of nitrogen deficiency in rice (Shin et al., 2018). In the research on nitrogen deficiency in Populus, 126 differentially expressed lncRNAs were identified in N-starved (). In this study, 16452 candidates lncRNAs revealed in LN and CK and HN treatments were selected from 9 LRs samples of tea plants. The results showed that 9451 lncRNAs with different functions were differentially expressed under nitrogen treatments. Therefore, these lncRNAs are involved in the response of tea plants to nitrogen stress, and the results can provide candidate genes for the study of nitrogen utilization in tea plants.
In this study, the results of KEGG pathway analysis showed that the DE-lncRNAs of the three treatments existed in multiple KEGG pathways, and the top three were “metabolic pathways”, “biosynthesis of secondary metabolites”, and “glutathione metabolism” (Figure 4). The results not only revealed that these pathways were related to nitrogen stress response, but also indicated that some components were involved in nitrogen stress response processes simultaneously. In addition, each group has specific DEGS and enrichment pathways, which means that there are different response mechanisms under different N stress. In addition, the up-regulated and down-regulated DE-lncRNAs of each group compared with the other groups were analyzed by KEGG pathway. The results demonstrated response mechanism varied with different N stresses.
In LN treatment, compared to the other two groups (CK, HN), up-regulated DE-lncRNAs were mainly annotated in “metabolic pathways” and “biosynthesis of secondary metabolites”. This indicated that nitrogen deficiency could promote the expression of coding genes in these two pathways. Considering the involvement of stress response: It has been reported that for many plant species, an increase in plant secondary metabolism caused by biotic or abiotic stress is part of the stress response (Sathiyabama et al., 2016). In line with this, previous studies manifested that the production of terpenoids and cannabinoids in plants is stimulated under the stress of nitrogen and phosphorus nutrient deficiency. The results showed a negative correlation between secondary metabolism and plant stress. That is, the concentration of secondary metabolites decreased with the increase of nitrogen supply, and the optimal nitrogen condition and high nitrogen condition in plants have nothing to do with secondary metabolism (Shiponi and Bernstein, 2021a; Shiponi and Bernstein, 2021b). The result is supported by many other plant studies. Furthermore, the carbon nutrient balance hypothesis points out that the production of nitrogen-rich primary metabolites is limited under low nitrogen conditions. Plant metabolism and energy consumption changed from the production of nitrogen-containing metabolites to the production of nitrogen-free metabolites, promoting the accumulation of secondary metabolites (Massad et al., 2012; ). This result is also consistent with the results of DE-lncRNAs analysis in our study.
Low-affinity transport system (LATS) and high-affinity transport system (HATS) are two NO3 uptake systems in plants They are mediated by nitrate transporter 1 (NRT1) and nitrate transporter 2 (NRT2), respectively (Tong et al., 2005; Wang et al., 2012). The results showed that with the decrease in nitrogen use efficiency, high-affinity nitrate transporters (NRT2) activity replaced low-affinity nitrate transporters (NRT1) activity, so NRT2 may be the key target to improve nitrogen use efficiency, especially under the condition of low nitrogen condition (). Moreover, studies have shown that under low nitrogen conditions, NRT2 not only transports NO3, but also participates in the transport of auxin, thus promoting plant root growth (Krouk et al., 2010). In this study, the up-regulated DE-lncRNA in low nitrogen treatment were clustered in the “nitrogen metabolism” pathway, including LJ.29731 and LJ.29732 corresponding to coding genes NRT2.3 and NRT2.4. This indicates that lncRNAs cis-regulate the expression of corresponding NRT2 coding genes under low nitrogen treatment and participate in the growth of lateral roots of tea plants (Figure 1, Figure 4, cluster2).
Plant hormones can regulate plant growth through different signal pathways and contribute to cell elongation, seed germination, flower formation, leaf morphogenesis, fruit development and so on. The same is true for resistance to biotic and abiotic stress (Qaadt, 2020). In this study, the normal nitrogen (CK), compared with the other two groups (LN and HN), up-regulated DE-lncRNAs were clustered in the “plant hormone signal transduction pathway”. Auxin is a crucial signal of plant growth and the higher the concentration, the more inhibited it is. At the same time, the precise distribution of auxin in plants determines the top-down organ morphology of plants (). AUX/IAA family transcription factors and auxin response factor (ARF) are two main transcription factors in the signal transduction pathway. In the auxin signal transduction pathway, ARF binds to the auxin response element in the auxin response gene promoter to activate the expression of the auxin response genes (). At low auxin concentration, AUX/IAA binds to ARF to form a heterodimer, which inhibits the transcription of ARFs, and then the expression of ARF-regulated auxin-responsive genes was inhibited. When auxin rises to a certain concentration, auxin binds to its receptor SCF-TIRl (). By interacting with AUX/IAA, AUX/IAA proteins were ubiquitinated and ARFs were released. Thus, auxin response genes are expressed and plants undergo a series of growth responses (). In this study, compared with CK and HN treatments, 14 DE-lncRNAs coding to ARF genes were up-regulated in LN treatment, which indicated that lncRNAs regulated the expression of ARF genes in low nitrogen treatment, thus increasing the transcriptional activity of ARFs (Table S12). Jasmonic acid (JA) and their derivatives methyl jasmonate (MeJA) play an essential role in plant morphological, physiological and biochemical processes in response to drought, cold and salt stress (; Zhang et al., 2021). Different nitrogen concentrations are also involved in the metabolism and synthesis of jasmonic acid. In previous studies, JAR1 (jasmonate resistant 1)-activated JA and ethylene signaling pathways could inhibit roots growth in Boron deficiency, and low nitrogen could inhibit the content of jasmonic acid (Li et al., 2020; ). In this study, three DE-lncRNAs (LJ.13289, LJ.13289, LJ.33982) corresponding to jasmonate resistant 4 were up-regulated in LN treatment. SAURs were discovered and named because of their response to the rapid induction of auxin. The expression of SAURs was also regulated by many internal and external factors (Wiesel et al., 2015). In addition to being induced by auxin, a small part of the expression of SAURs showed auxin inhibition or non-response (Ren and Gray, 2015). Many studies have shown that SAURs can also be induced by other plant hormones, such as brassinolide, gibberellin and cytokinin, while abscisic acid and jasmonic acid inhibit it (). In addition, light, temperature, moisture and other environmental factors also affect the expression of SAURs (Ouyang et al., 2015). In this study, 7 DE-lncRNAs coding to SAUR genes were down-regulated in HN treatment compared with LN and HN treatments (LJ.1615, LJ.12746, LJ.11860, LJ.5020, LJ.12505, LJ.23544, LJ.36187). To sum up, DE-lncRNAs are involved in the synthesis and metabolism of hormones in tea plants under different nitrogen conditions.
Glutathione is a reductant that regulates signal molecules and scavenges free radicals and ROS in the redox system of plant cells. Glutathione participates in the metabolic process of all eukaryotes to improve plants’ tolerance to abiotic stresses such as salt, temperature, drought, heavy metals and autotoxicity (; Koprivova et al., 2010; Niu et al., 2017). Glutathione adjusts continuously in plants, giving play to different metabolic processes and extremely significant antioxidant advantages, which can effectively scavenge free radicals in organisms, which cannot be reflected by glutamine, lycopene and other oxidants (Nakamura et al., 2019). Glutathione can oxidize GSH by reactive oxygen species and regulate the reduction of GSSG by glutathione reductase, thus regulating the balance between GSH (reduced) and GSSG (oxidized) (). Exogenous GSSG could not induce root growth under normal conditions but could promote root development under auxin treatment (Tyburski and Tretyn, 2010). Therefore, the interaction between auxin and GSSG can regulate plant root development. It is reported that plants can inhibit the growth of lateral roots by regulating the reduction ratio of GSH/GSSG in the presence of auxin (Wei et al., 2013). Glutathione S-transferase (GSTs) converts GSH to GSSG, while glutathione reductase (GR) induces the reduction of GSSG to GSH, so glutathione S-transferase plays an essential role in plant growth and development (Jiang et al., 2010). In this study, cluster analysis of different treatments showed that a total of 82 De-lncRNAs coding glutathione S-transferase genes were distributed in six clusters, indicating that in nitrogen treatment, lncRNAs could regulate the coding glutathione S-transferase genes and participated in the growth and development of tea plants (Figure 4, Table S13).
Conclusion
Through RNA-seq and bioinformatics, this study identified 16,452 LncRNAs in the whole genome of tea plants, among which 9,451 DE-lncRNAs were differentially expressed under nitrogen stress. The KEGG pathway further analyzed the biological and regulation function of these lncRNAs. To sum up, this study expands the cognition of lncRNAs in tea plants and explores a path for further study of tea plants response to nitrogen and improvement of nitrogen use efficiency of woody plants.
Statements
Data availability statement
The data presented in the study are deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive repository, accession number PRJNA595712.
Author contributions
SH: Project administration. Methodology, Formal analysis, Investigation, Data curation, Writing-original draft. XL: Methodology, Formal analysis, Investigation, Data curation, Writing-original draft. XC: Formal analysis, Investigation. WX: Formal analysis, Data curation, Writing-original draft. YH: Investigation, Visualization. JL: Investigation. HM: Conceptualization. ZZ: Formal analysis. YZ: Investigation. AJ: Visualization. RH: Formal analysis. All authors contributed to the article and approved the submitted version.
Acknowledgments
This study was financially supported by the Science and technology projects of Nanjing (202210013), Expert Workstation of Yunnan Province (202105AF150045), National Natural Science Foundation of China (32172628), Top-notch Academic Programs Project of Jiangsu Higher Education Institutions (2022) and China Agriculture Research System of MOF and MARA (CARS-19). We thank Dr. Yuehua Ma (Central laboratory of College of Horticulture, Nanjing Agricultural University) for using Bio-Rad Real-time thermal cycler CFX96.
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.2023.1080427/full#supplementary-material
References
1
AlbornozF. (2016). Crop responses to nitrogen overfertilization: a review. Sci. Hortic.205, 79–83. doi: 10.1016/j.scienta.2016.04.026
2
AliM.BaekK.-H. (2020). Jasmonic acid signaling pathway in response to abiotic stresses in plants. Int. J. Mol. Sci.21, 621. doi: 10.3390/ijms21020621
3
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215, 403–410. doi: 10.1016/S0022-2836(05)80360-2
4
BaiM. Y.ShangJ. X.OhE.FanM.BaiY.ZentellaR.et al. (2012). Brassinosteroid, gibberellin and phytochrome impinge on a common transcription module in arabidopsis. Nat. Cell Biol.14 (8), 810–817. doi: 10.1038/ncb2546
5
BashandyT.GuilleminotJ.VernouxT.Caparros-RuizD.LjungK.MeyerY.et al. (2010). Interplay between the NADP-linked thioredoxin and glutathione systems in arabidopsis auxin signaling. Plant Cell.22 (2), 376–391. doi: 10.1105/tpc.109.071225
6
BiY. M.KantS.ClarkJ.GiddaS.MingF.XuJ.et al. (2009). Increased nitrogen-use efficiency in transgenic rice plants over-expressing a nitrogen-responsive early nodulin gene identified from rice expression profiling. Plant Cell Environ.32, 1749–1760. doi: 10.1111/j.1365-3040.2009.02032.x
7
CaleyD.PinkR.TrujillanoD.CarterD. (2010). Long noncoding RNAs, chromatin, and development. ScientificWorldJournal10, 90–102. doi: 10.1100/tsw.2010.7
8
ChekanovaJ. A. (2015). Long non-coding RNAs and their functions in plants. Curr. Opin. Plant Biol.27, 207–216. doi: 10.1016/j.pbi.2015.08.003
9
ChenZ.JiangQ.JiangP.ZhangW.HuangJ.LiuC.et al. (2020). Novel low-nitrogen stress-responsive long non-coding RNAs (lncRNA) in barley landrace B968 (Liuzhutouzidamai) at seedling stage. BMC Plant Biol.20 (1), 142. doi: 10.1186/s12870-020-02350-2
10
ChenZ.LiY.LiuC.WangY.HeT.GuoG.et al. (2018). Reference gene selection for quantitative RT-PCR normalisation in barley under low-nitrogen stress, based on RNAseq data. J. Cereal Sci.82, 213–215. doi: 10.1016/j.jcs.2018.06.014
11
ChenM.WangC.BaoH.ChenH.WangY. (2016). Genome-wide identification and characterization of novel lncRNAs in populus under nitrogen deficiency. Mol. Genet. Genom.291, 1663–1680. doi: 10.1007/s00438-016-1210-3
12
CoxM. P.PetersonD. A.BiggsP. J. (2010). SolexaQA: At-a-glance quality assessment of illumina second-generation sequencing data. BMC Bioinf.11 (1), 485. doi: 10.1186/1471-2105-11-485
13
De SmetI.LauS.MayerU.JürgensG. (2010). Embryogenesis-the humble beginnings of plant life. Plant J.61 (6), 959–970. doi: 10.1111/j.1365-313X.2010.04143.x
14
DuanL. N.DietrichD.NgC. H.ChanP. M. Y.BhaleraoR.BennettM. J.et al. (2013). Endodermal ABA signaling promotes lateral root quiescence during salt stress in arabidopsis seedlings. Plant Cell.25 (1), 324–341. doi: 10.1105/tpc.112.107227
15
El-GebaliS.MistryJ.BatemanA.EddyS. R.LucianiA.PotterS. C.et al. (2018). The pfam protein families database in 2019. Nucleic Acids Res.47, D427–D432. doi: 10.1093/nar/gky995
16
FigueiredoM. R. A.StraderL. C. (2022). Intrinsic and extrinsic regulators of Aux/IAA protein degradation dynamics. Trends Biochem. Sci.47 (10), 865–874. doi: 10.1016/j.tibs.2022.06.004
17
GoodA. G.BeattyP. H. (2011). Fertilizing nature: a tragedy of excess in the commons. PloS Bio l.9, e1001124. doi: 10.1371/journal.pbio.1001124
18
GruffmanL.JamtgardS.NasholmT. (2014). Plant nitrogen status and cooccurrence of organic and inorganic nitrogen sources influence root uptake by scots pine seedlings. Tree Physiol.34 (2), 205–213. doi: 10.1093/treephys/tpt121
19
HanM.OkamotoM.BeattyP. H.RothsteinS. J.GoodA. G. (2015). The genetics of nitrogen use efficiency in crop plants. Annu. Rev. Genet.49, 1–9. doi: 10.1146/annurev-genet-112414-055037
20
HanH.RakusováH.VerstraetenI.ZhangY.FrimlJ. (2020). SCFTIR1/AFB auxin signaling for bending termination during shoot gravitropism. Plant Physiol.183 (1), 37–40. doi: 10.1104/pp.20.00212
21
HanZ. X.RanaM. M.LiuG. F.GaoM. J.LiD. X.WuF. G.et al. (2016a). Green tea flavour determinants and their changes over manufacturing processes. Food Chem.212, 739–748. doi: 10.1016/j.foodchem.2016.06.049
22
HodgesD. M.AndrewsC. J.JohnsonD. A.HamiltonR. I. (1996). Antioxidant compound responses to chilling stress in differentially sensitive inbred maize lines. Physiol. Plant98, 685–692. doi: 10.1111/j.1399-3054.1996.tb06672.x
23
HuangN. C.LiuK. H.LoH. J.TsayY. F. (1999). Cloning and functional characterization of an arabidopsis nitrate transporter gene that encodes a constitutive component of low-affinity uptake. Plant Cell.11 (8), 1381–1392. doi: 10.1105/tpc.11.8.1381
24
HuangY.WangS.ShiL.XuF. (2021). JASMONATE RESISTANT 1 negatively regulates root growth under boron deficiency in arabidopsis. J. Exp. Bot.72 (8), 3108–3121. doi: 10.1093/jxb/erab041
25
HuangN. C.LiuK. H.LoH. J.TsayY. F. (1999). Cloning and functional characterization of an Arabidopsis nitrate transporter gene that encodes a constitutive component of low-affinity uptake. Plant Cell11 (8), 1381–1392.
26
HuS.ZhangM.YangY.XuanW.ZouZ.ArkorfulE.et al. (2020). A novel insight into nitrogen and auxin signaling in lateral root formation in tea plant [Camellia sinensis (L.) o. kuntze]. BMC Plant Biol.20 (1), 232. doi: 10.1186/s12870-020-02448-7
27
IyerM. K.NiknafsY. S.MalikR.SinghalU.SahuA.HosonoY.et al. (2015). The landscape of long noncoding RNAs in the human transcriptome. Nat. Genet.47 (3), 199–208. doi: 10.1038/ng.3192
28
JiangH. W.LiuM. J.ChenI. C.HuangC. H.ChaoL. Y.HsiehH. L. (2010). A glutathione s-transferase regulated by light and hormones participates in the modulation of arabidopsis seedling development. Plant Physiol.154 (4), 1646–1658. doi: 10.1104/pp.110.159152
29
KangY.YangD.KongL.HouM.MengY.WeiL.et al. (2017). CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features. Nucleic Acids Res.45, W12–W16. doi: 10.1093/nar/gkx428
30
KantS.BiY.RothsteinS. J. (2011). Understanding plant response to nitrogen limitation for the improvement of crop nitrogen use efficiency. J. Exp. Bot.62 (4), 1499–1509. doi: 10.1093/jxb/erq297
31
KelleyD.RinnJ. (2012). Transposable elements reveal a stem cell-specific class of long noncoding RNAs. Genom. Biol.13, R107. doi: 10.1186/gb-2012-13-11-r107
32
KongL.ZhangY.YeZ.LiuX.ZhaoS.WeiL.et al. (2007). CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res.35 (suppl2), W345–9.
33
KoppF.MendellJ. T. (2018). Functional classification and experimental dissection of long noncoding RNAs. Cell172, 393–407. doi: 10.1016/j.cell.2018.01.011
34
KoprivovaA.MugfordS. T.KoprivaS. (2010). Arabidopsis root growth dependence on glutathione is linked to auxin transport. Plant Cell29, 1157–1167. doi: 10.1007/s00299-010-0902-0
35
KrappA. (2015). Plant nitrogen assimilation and its regulation: a complex puzzle with missing pieces. Curr. Opin. Plant Biol.25, 115–122. doi: 10.1016/j.pbi.2015.05.010
36
KroukG.LacombeB.BielachA.Perrine-WalkerF.MalinskaK.MounierE.et al. (2010). Nitrate-regulated auxin transport by NRT1.1 defines a mechanism for nutrient sensing in plants. Dev. Cell18 (6), 927–937. doi: 10.1016/j.devcel
37
KungJ. T. Y.ColognoriD.LeeJ. T. (2013). Long noncoding RNAs: past, present, and future. Genetics193 (3), 651–669. doi: 10.1534/genetics.112.146704
38
LezhnevaL.KibaT.Feria-BourrellierA. B.LafougeF.Boutet-MerceyS.ZoufanP.et al. (2014). The arabidopsis nitrate transporter NRT2.5 plays a role in nitrate acquisition and remobilization in nitrogen-starved plants. Plant J.80 (2), 230–241. doi: 10.1111/tpj.12626
39
LiangW.HuF.QiW.ZhaoC.ChenT.WangC.et al. (2022). Comprehensive transcriptome analysis of GS3 near-isogenic lines during panicle development in rice (Oryza sativa l.). Front. Genet.13. doi: 10.3389/fgene.2022.857143
40
LiS.HamaniA. K. M.SiZ.LiangY.GaoY.DuanA. (2020). Leaf gas exchange of tomato depends on abscisic acid and jasmonic acid in response to neighboring plants under different soil nitrogen regimes. Plants (Basel).9 (12), 1674. doi: 10.3390/plants9121674
41
LinZ. H.ChenC. S.ZhongQ. S.RuanQ. C.ChenZ. H.YouX. M.et al. (2021). The GC-TOF/MS-based metabolomic analysis reveals altered metabolic profiles in nitrogen-deficient leaves and roots of tea plants (Camellia sinensis). BMC Plant Biol.21 (1), 506. doi: 10.1186/s12870-021-03285-y
42
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(–delta delta c) method. Methods25 (4), 402–408. doi: 10.1006/meth.2001.1262
43
LuoC.HeB.ShiP.XiJ.GuiH.PangB.et al. (2022). Transcriptome dynamics uncovers long non-coding RNAs response to salinity stress in chenopodium quinoa. Front. Plant Sci.13. doi: 10.3389/fpls.2022.988845
44
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(1), 864. doi: 10.1186/s12864-019-6245-5
45
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. doi: 10.1186/s12864-016-2650-1
46
LynchJ. P. (2013). Steep, cheap and deep: an ideotype to optimize water and n acquisition by maize root systems. Ann. Bot-London.112 (2), 347–357. doi: 10.1093/aob/mcs293
47
MassadT. J.DyerL. A.VegaC. G. (2012). Costs of defense and a test of the carbon-nutrient balance and growth-differentiation balance hypotheses for two co-occurring classes of plant defense. PloS One7, e47554. doi: 10.1371/journal.pone.0047554
48
NakamuraS. I.WongkaewA.NakaiY.RaiH.Ohkama-OhtsuN. (2019). Foliar-applied glutathione activates zinc transport from roots to shoots in oilseed rape. Plant Sci., 424–434. doi: 10.1016/j.plantsci.2018.10.018
49
NiuX.SongL.XiaoY.GeW. (2017). Drought-tolerant plant growth-promoting rhizobacteria associated with foxtail millet in a semi-arid agroecosytem and their potential in alleviating drought stress. Front. Microbiol.8 (1), 2580. doi: 10.3389/fmicb.2017.02580
50
NovoaR.LoomisR. S. (1981). Nitrogen and plant production. Plant Soil.58 (1-3), 177–204. doi: 10.1007/BF02180053
51
OuyangF.MaoJ. F.WangJ.ZhangS.LiY. (2015). Transcriptome analysis reveals that red and blue light regulate growth and phytohormone metabolism in Norway spruce [Picea abies (L.) karst]. PloS One10 (8), e0127896. doi: 10.1371/journal.pone.0127896
52
PerteaM.KimD.PerteaG. M.LeekJ. T.SalzbergS. L. (2016). Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and ballgown. Nat. Protoc.11, 1650–1667. doi: 10.1038/nprot.2016.095
53
QaadtR. (2020). Phytohormone signaling mechanisms and genetic methods for their modulation and detection. Curr. Opin. Plant Biol.57, 31–40. doi: 10.1016/j.pbi.2020.05.011
54
RenH.GrayW. M. (2015). SAUR proteins as effectors of hormonal and environmental signals in plant growth. Mol. Plant8 (8), 1153–1164. doi: 10.1016/j.molp.2015.05.003
55
RobertsA.PimentelH.TrapnellC.PachterL. (2011). Identification of novel transcripts in annotated genomes using RNA-seq. Bioinformatics27, 2325–2329. doi: 10.1093/bioinformatics/btr355
56
RuanJ.HaerdterR.GerendásJ. (2010). Impact of nitrogen supply on carbon/nitrogen allocation: a case study on amino acids and catechins in green tea [Camellia sinensis (L.) o. kuntze] plants. Plant Biol.12, 724–734. doi: 10.1111/j.1438-8677.2009.00288.x
57
RuanL.WeiK.WangL. Y.ChengH.WuL.LiH. (2019). Characteristics of free amino acids (the quality chemical components of tea) under spatial heterogeneity of different nitrogen forms in tea (Camellia sinensis) plants. Molecules24, 415. doi: 10.3390/molecules24030415
58
SathiyabamaM.BernsteinN.AnusuyaS. (2016). Chitosan elicitation for increased curcumin production and stimulation of defence response in turmeric (Curcuma longa l.). Ind. Crops Prod.89, 87–94. doi: 10.1016/j.indcrop.2016.05.007
59
ShinS. Y.JeongJ. S.LimJ. Y.KimT.ParkJ. H.KimJ. K.et al. (2018). Transcriptomic analyses of rice (Oryza sativa) genes and non-coding RNAs under nitrogen starvation using multiple omics technologies. BMC Genom.19, 532. doi: 10.1186/s12864-018-4897-1
60
ShiponiS.BernsteinN. (2021a). Response of medical cannabis (Cannabis sativa l.) genotypes to p supply under long photoperiod: functional phenotyping and the ionome. Ind. Crops Prod.161, 113154. doi: 10.1016/j.indcrop.2020.113154
61
ShiponiS.BernsteinN. (2021b). The highs and lows of p supply in medical cannabis: effects on cannabinoids, the ionome, and morpho-physiology. Front. Plant Sci.12, 657323. doi: 10.3389/fpls.2021.657323
62
SunkarR.ChinnusamyV.ZhuJ.ZhuJ. K. (2007). Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci.12 (7), 301–309. doi: 10.1016/j.tplants.2007.05.001
63
SunC. H.YuJ. Q.HuD. G. (2017). Nitrate: a crucial signal during lateral roots development. Front. Plant Sci.8. doi: 10.3389/fpls.2017.00485
64
TongY.ZhouJ. J.LiZ.MillerA. J. (2005). A two-component high-affinity nitrate uptake system in barley. Plant J.41 (3), 442–450. doi: 10.1111/j.1365-313X.2004.02310.x
65
TyburskiJ.TretynA. (2010). Glutathione and glutathione disulfide affect adventitious root formation and growth in tomato seedling cuttings. Acta Physiol. Plant32 (2), 411–417. doi: 10.1007/s11738-009-0418-9
66
VenkatesanS.MurugesanS.GanapathyM. N.VermaD. P. (2004). Long−term impact of nitrogen and potassium fertilizers on yield, soil nutrients and biochemical parameters of tea. J. Sci. Food Agr.84, 1939–1944. doi: 10.1002/jsfa.1897
67
WangH. V.ChekanovaJ. A. (2017). Long noncoding RNAs in plants. Adv. Exp. Med. Biol.1008, 133–154. doi: 10.1007/978-981-10-5203-3_5
68
WangY.FanX.LinF.HeG.TerzaghiW.ZhuD.et al. (2014). Arabidopsis noncoding RNA mediates control of photomorphogenesis by red light. Proc. Natl. Acad. Sci. U S A.111 (28), 10359–10364. doi: 10.1073/pnas.1409457111
69
WangY. Y.HsuP. K.TsayY. F. (2012). Uptake, allocation and signaling of nitrate. Trends Plant Sci.17 (8), 458–467. doi: 10.1016/j.tplants.2012.04.006
70
WanS.ZhangY.DuanM.HuangL.WangW.XuQ.et al. (2020). Integrated analysis of long non-coding RNAs (lncRNAs) and mRNAs reveals the regulatory role of lncRNAs associated with salt resistance in camellia sinensis. Front. Plant Sci.11, 218. doi: 10.3389/fpls.2020.00218
71
WeiK.WangL. Y.ChengH.ZhangC. C.MaC. L.ZhangL. Q.et al. (2013). Identification of genes involved in indole-3-butyric acid-induced adventitious root formation in nodal cuttings of camellia sinensis (L.) by suppression subtractive hybridization. Gene514 (2), 91–98. doi: 10.1016/j.gene.2012.11.008
72
WeiC. L.YangH.WangS. B.ZhaoJ.LiuC.GaoL. P.et al. (2018). Draft genome sequence of camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. Proc. Natl. Acad. Sci. U. S. A.115, E4151–E4158. doi: 10.1073/pnas.1719622115
73
WieselL.DavisJ. L.MilneL.et al. (2015). A transcriptional reference map of defence hormone responses in potato[J]. Sci. Rep.5, 15229. doi: 10.1038/srep15229
74
XiaoH. M.YuanZ. T.GuoD. H.HouB. F.YinC. L.ZhangW. Q.et al. (2015). Genome-wide identification of long noncoding RNA genes and their potential association with fecundity and virulence in rice brown planthopper, nilaparvata lugens. BMC Genomics16 (1), 749. doi: 10.1186/s12864-015-1953-y
75
XiaE. H.ZhangH. B.ShengJ.LiK.ZhangQ. J.KimC.et al. (2017). The tea tree genome provides insights into tea flavor and independent evolution of caffeine biosynthesis. Mol. Plant10, 866–877. doi: 10.1016/j.molp.2017.04.002
76
XuG.FanX.MillerA. J. (2012). Plant nitrogen assimilation and use efficiency. Annu. Rev. Plant Biol.63, 153–182. doi: 10.1146/annurev-arplant-042811-105532
77
YuG.WangL. G.HanY.HeQ.-Y. (2012). clusterProfiler: an r package for comparing biological themes among gene clusters. OMICS: A J. Integr. Biol.16, 284–287. doi: 10.1089/omi.2011.0118
78
ZhangW.HanZ.GuoQ.LiuY.ZhengY.WuF.et al. (2014). Identification of maize long non-coding RNAs responsive to drought stress. PloS One9, e98958. doi: 10.1371/journal.pone.0098958
79
ZhangQ.HeK.HuoH. (2012). Cleaning china’s air. Nature484, 161–162. doi: 10.1038/484161a
80
ZhangX.LiuH.Pilon-SmitsE.HuangW.WangP.WangM.et al. (2020). Transcriptome-wide analysis of nitrogen-regulated genes in tea plant (Camellia sinensis l. o. kuntze) and characterization of amino acid transporter CsCAT9.1. Plants (Basel).9 (9), 1218. doi: 10.3390/plants9091218
81
ZhangS.ZhuL.ShenC.JiZ.ZhangH.ZhangT.et al. (2021). Natural allelic variation in a modulator of auxin homeostasis improves grain yield and nitrogen use efficiency in rice. Plant Cell.33 (3), 566–580. doi: 10.1093/plcell/koaa037
82
ZhaoX.LiJ.LianB.GuH.LiY.QiY. (2018). Global identification of arabidopsis lncRNAs reveals the regulation of MAF4 by a natural antisense RNA. Nat. Commun.9 (1), 5056. doi: 10.1038/s41467-018-07500-7
83
ZhuT.ZhangJ.MengT.ZhangY.YangJ.MüllerC.et al. (2014). Tea plantation destroys soil retention of NO3– and increases N2O emissions in subtropical China. Soil Biol. Biochem.73, 106–114. doi: 10.1016/j.soilbio.2014.02.016
Summary
Keywords
long non-coding RNAs, nitrogen stress, RNA-seq, Camellia sinensis, lateral roots
Citation
Hu S, Hu Y, Mei H, Li J, Xuan W, Jeyaraj A, Zhao Z, Zhao Y, Han R, Chen X and Li X (2023) Genome-wide analysis of long non-coding RNAs (lncRNAs) in tea plants (Camellia sinensis) lateral roots in response to nitrogen application. Front. Plant Sci. 14:1080427. doi: 10.3389/fpls.2023.1080427
Received
26 October 2022
Accepted
23 January 2023
Published
24 February 2023
Volume
14 - 2023
Edited by
Weicong Qi, Jiangsu Academy of Agricultural Sciences (JAAS), China
Reviewed by
Kang Wei, Tea Research Institute, Chinese Academy of Agricultural Sciences, China; Lin Meng, Tobacco Research Institute, Chinese Academy of Agricultural Sciences, China
Updates
Copyright
© 2023 Hu, Hu, Mei, Li, Xuan, Jeyaraj, Zhao, Zhao, Han, Chen and Li.
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: Xinghui Li, lxh@njau.edu.cn
†These authors share first authorship
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.