Integrated Analysis of Long Non-coding RNAs (lncRNAs) and mRNAs Reveals the Regulatory Role of lncRNAs Associated With Salt Resistance in Camellia sinensis

Tea plant (Camellia sinensis), an important economic crop, is seriously affected by various abiotic stresses, including salt stress, which severely diminishes its widespread planting. However, little is known about the roles of long non-coding RNAs (lncRNAs) in transcriptional regulation under salt stress. In this study, high-throughput sequencing of tea shoots under salt-stress and control conditions was performed. Through sequencing analysis, 16,452 unique lncRNAs were identified, including 172 differentially expressed lncRNAs (DE-lncRNAs). The results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of their cis- and trans-target genes showed that these DE-lncRNAs play important roles in many pathways such as the galactinol synthase (GOLS), calcium signaling pathway, and interact with transcription factors (TFs) under salt stress. The data from the gene-specific antisense oligodeoxynucleotide-mediated reduction in the lncRNA MSTRG.139242.1 and its predicted interacting gene, TEA027212.1 (Ca2+-ATPase 13), in tea leaves revealed that MSTRG.139242.1 may function in the response of tea plants to high salinity. In addition, 12 lncRNAs were predicted to be target mimics of 17 known mature miRNAs, such as miR156, that are related to the salt-stress response in C. sinensis. Our results provide new insights into lncRNAs as ubiquitous regulators in response to salt stress in tea plants.


INTRODUCTION
Tea plant [Camellia sinensis (L.) O. Kuntze] is an important economic crop in many countries, such as China, Japan, and India. Tea made from tea plant leaves is one of the most consumed drinks worldwide. However, tea plants are seriously affected by many abiotic stresses, such as cold, salt, and drought stresses, due to their specific environmental requirements, which severely diminish the widespread planting of tea plants (Upadhyaya and Panda, 2013). Therefore, it is important for the tea industry to breed tea plant species that exhibit strong resistance and adaptability. Salt stress is a form of abiotic stress that severely affects the growth of tea plants and the quality of tea products (Upadhyaya and Panda, 2013). In our previous report, we performed an Illumina RNA sequencing (RNA-Seq) to compare the transcriptomes of tea plants treated with and without NaCl, and many differentially expressed mRNAs involved in signal transduction pathways, transcription factors (TFs), and other functional genes under salt stress were identified . Previous reports have also shown that many mRNAs and proteins in tea plants, such as CsSnRK2 , CsAQP (Yue et al., 2014), and CsVQ (Guo et al., 2018), respond to salt stress. Long non-coding RNAs (lncRNAs) are important regulatory factors that respond to various abiotic stresses in plants; however, their functions in responding to salt stress have not been studied in tea plants. Therefore, the identification of important lncRNAs that respond to salt stress in tea plants is necessary.
LncRNAs are non-coding RNAs that are longer than 200 nt and usually have low protein-coding potential (Mercer et al., 2009). According to their genomic origins, lncRNAs are broadly divided into four types: intronic lncRNAs, intergenic lncRNAs (lincRNAs), sense lncRNAs, and antisense lncRNAs (Liu et al., 2012). LncRNAs are considered "transcriptional noise" because of their poor conservation and the limited evidence for their functions (Ponjavic et al., 2007). However, increasing research has revealed that lncRNAs function in transcriptional regulation and epigenetic gene regulation, and some lncRNAs function as cis-or trans-regulators in various biological processes (Ponting et al., 2009). In mammals, many lncRNAs have been proven to be related to a wide range of diseases, especially cancers and neurodegenerative diseases (Ma et al., 2012). In recent years, a large number of lncRNAs that are important in plant growth and development (Heo and Sung, 2011;Ding et al., 2012), biotic stress responses , and abiotic stress responses (Swiezewski et al., 2009;Qin et al., 2017) have been identified. For example, the reduced expression of the rice lncRNA LDMAR can lead to male sterility under long-day conditions in rice (Ding et al., 2012). In Arabidopsis, the antisense lncRNA FLORE can regulate the circadian clock to photoperiodic flowering (Henriques et al., 2017). A nucleus-localized lncRNA, DRIR, can enhance drought-and salt-stress tolerance in Arabidopsis (Qin et al., 2017). Currently, with the development of deep-sequencing technology, genome-wide lncRNA analyses have been performed on many species, such as Arabidopsis (Liu et al., 2012), rice (Zhang Y.C. et al., 2014), Zea mays (Boerner and McGinnis, 2012), wheat (Shumayla et al., 2017), and other species (Zhang and Chen, 2013). However, in many non-model plants, studies on lncRNAs are rather limited. Recently, Varshney et al. (2019) discovered 33,400 putative lncRNAs in different tissues of tea plants, which was the first report of lncRNAs in tea plants.
In this study, to explore the early response of lncRNA under salt stress in tea plants, we reanalyzed the transcriptomes of tea plants treated with and without NaCl based on the latest tea tree genome (Wei et al., 2018). LncRNAs were identified and classified systematically. The functions of these lncRNAs were predicted, and differentially expressed lncRNAs (DE-lncRNAs) that responded to salt stress were identified. In addition, the cis-targeting and trans-targeting relationship between lncRNAs and mRNAs and the possible interactions between lncRNAs and abiotic stress-related miRNAs were found to be involved in salt stress in tea plants. Notably, in our previous report, the Ca 2+ -transporting ATPase was shown to be involved in the salt-stress response in tea plants , and in this study, the differentially expressed Ca 2+ -transporting ATPase 13 (TEA027212.1) was proven to be co-expressed with the lncRNA MSTRG.139242.1, which indicates that this lncRNA may participate in Ca 2+ signal transduction in response to salt stress. In addition, six DE-lncRNAs were randomly selected for quantitative real-time PCR (qRT-PCR) validation to confirm the reliability of the expression levels obtained from the RNA-Seq transcriptome. Overall, the results obtained in our study provide a valuable resource for studying lncRNAs involved in salt stress in tea plants and will enhance our understanding of the putative regulatory functions of lncRNAs in plants.

Plant Materials and NaCl Treatments
One-year-old "Pingyangtezao" tea plant cutting seedlings with consistent growth were pre-incubated in nutrient solution under standard growth conditions (temperature: 25 ± 3 • C, air relative humidity: 60-70%, photoperiod: 12 h light/12 h dark) for 2 months. These tea plants were used for stress assays under the same growth conditions. In order to analyze the early response of tea plant to high salinity and evaluate the early salt responsive genes in the tea plant leaves, salt treatment was conducted by soaking the roots in nutrient solution that contained 250 mM NaCl (salt treatment) or standard nutrient solution (control), and fresh leaf samples (first and second leaves) were randomly collected at 4 h when tea plants began to show the symptoms of salt stress and wilt slightly after treatment. Control and NaCl treatment were both repeated three times (Control1,Control2,and Control3 and Salt1,Salt2,and Salt3).The samples were immediately immersed in liquid nitrogen and stored at −80 • C until RNA-Seq, qRT-PCR validation, and further analysis.

Analysis of Transcriptomic Data Based on the Tea Tree Reference Genome
In our previous report, we sequenced six tea plant samples using RNA-Seq technology. A total of 79.40 Gb of clean data was obtained, and the Q30 base percentage of each sample was greater than 92.93% . The identification of lncRNAs was executed according to the pipeline shown in Figure 1A. All RNA-Seq datasets were aligned to the reference genome of C. sinensis var. sinensis (Wei et al., 2018) using the HISAT2 system 1 (Kim et al., 2015) to reconstruct the transcriptome. After the alignment, the StringTie software 2 (Pertea et al., 2016) was used to assemble reads into transcripts and for quantification. The assembled transcripts were annotated using the gffcompare program, 3 and the unknown transcripts were used to screen for putative lncRNAs.

Identification of LncRNAs in Tea Plants
The prediction of lncRNAs consisted of two parts: basic filtering and potential coding capability screening. For basic filtering, transcripts with a class code of "i, " "x, " "u, " "o, " and "e" (Lv et al., 2013), a length ≥ 200 bp, an exon number ≥ 2, and a fragments per kilobase of transcript per million mapped reads (FPKM) value ≥ 0.1 (Kelley and Rinn, 2012) were selected. Subsequently, four methods, including Coding Potential Calculator (CPC) 4 (Kong et al., 2007), Coding-Non-Coding Index (CNCI) 5 , Coding Potential Assessment Tool (CPAT) 6 , and Protein family (Pfam) 7 (Finn et al., 2014) protein domain analysis were used to evaluate the coding capabilities of the transcripts, and potential lncRNAs were filtered out by combining the four results. Next, the identified lncRNAs were classified into four categories: lincRNA, antisense-lncRNA, intronic-lncRNA, and sense-lncRNA using cuffcompare. 8 In addition, lncRNAs and mRNAs were comparatively analyzed according to their transcript length and exon number.

Prediction of Potential Target Genes of LncRNAs
The prediction of potential target genes of lncRNAs was based on the two interaction modes between lncRNAs and mRNAs. First, because lncRNA can regulate the expression of nearby genes, we searched coding genes 100 kb upstream and downstream of lncRNAs for cis-target genes. The other interaction mode of lncRNA and mRNA is due to the complementary pairing of bases. The LncTar  tool was used for target gene prediction of lncRNAs. The free energy and standard free energy of paired sites were calculated, and the target genes with standard free energy threshold < −0.1 were considered transtarget genes of lncRNAs.

Identification and Functional Analysis of DE-LncRNAs
The expression levels of lncRNAs were quantified using FPKM values with the StringTie software. Differential expression analysis of the NaCl and Control groups was performed using the DESeq R package (1.10.1). 9 The resulting P-values were adjusted using Benjamini and Hochberg's approach for controlling the false discovery rate. LncRNAs with an adjusted value of P < 0.01 and an absolute value of log2 (fold change) > 1 found by DESeq were considered as differentially expressed.
Functional analyses of DE-lncRNAs were performed by annotation and classification of their cis-and trans-target genes. BLASTX alignment (value of E < 10 −5 ) (Altschul et al., 1997) between target genes and public databases was performed, which included the Clusters of Orthologous Groups (COG) of proteins (Tatusov et al., 2000) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases (Kanehisa et al., 2004). Gene Ontology (GO) classifications were conducted using Blast2GO and WEGO software (Ashburner et al., 2000).

Prediction of LncRNAs as Endogenous Target Mimics for miRNAs Under Salt Stress
The 172 DE-lncRNAs and 539 known conserved mature miRNAs of C. sinensis (Zhu and Luo, 2013;Jeyaraj et al., 2017a,b;Sun et al., 2018) were submitted to psRNATarget 10 using a maximum expectation of 2.0 Dai et al., 2018) for miRNA-lncRNA interaction prediction. Less than three mismatches and G/U pairs were allowed within the lncRNA and miRNA pairing regions. The coexpression network was established using Cytoscape, 11 based on the DE-lncRNAs, target mRNAs of DE-lncRNAs, and miRNAs interacting with DE-lncRNAs.

LncRNA and Target mRNA Suppression in Tea Plants Using Antisense Oligonucleotides (AsODNs)
A gene-suppression assay in tea plants using antisense oligonucleotides (AsODNs) was conducted according to Liu et al. (2018). ODN sequences were selected using the Soligo software. 12 The sequences of input and selected ODNs are shown in Supplementary Data Sheet S1: Table S1. To increase the efficiency of gene suppression, two or three independent ODNs were synthesized of each gene and mixed in equal proportions during treatment. The ODNs were purified by HPLC, and the concentration was diluted to 10 µM. Fresh tea shoots with a bud and two leaves were excised and inserted into centrifuge tubes containing 1 ml of mixed ODNs for treatment. All leaves from the shoots were harvested after 24 h of treatment, and each treatment had three biological repeats. At the same time, a random non-sense control was used (Supplementary Data Sheet S1: Table S1). A nucleotide BLAST search of the tea tree genome CDS and lncRNA sequences with this sequence was used to verify that the sequences did not overlap with other sequences.

qRT-PCR
qRT-PCR was performed to validate the expression patterns of salt-responsive lncRNAs. Total RNA was isolated from the six tea samples (Salt1, Salt2, Salt3, Control1, Control2, and Control3) using TRNzol reagent (TIANGEN, Beijing, China). First-strand cDNA was synthesized using the PrimeScript RT reagent kit with gDNA Eraser (TaKaRa, Dalian, China), after which qRT-PCR detection was completed using the EvaGreen qPCR MasterMix-No Dye kit (ABM, Richmond, BC, Canada) on a StepOne Plus PCR instrument (ABI, United States) following the manufacturers' protocols. The tea plant Csβ-actin gene was chosen as the reference gene in accordance with previous methods , and the primers used in this assay are shown in Supplementary Data sheet S1: Table S2. The relative expression levels of the genes were calculated using the 2 − CT method (Livak and Schmittgen, 2001), and the data are presented as the mean ± SD from three independent biological replicates. Group differences were tested using one-way ANOVA and Duncan's test, and significant differences among various treatment groups are represented by different letters (P < 0.05). The SPSS 20.0 software was used to determine significant differences between the Control and 250 mM NaCl treatment data.

Identification of LncRNAs in C. sinensis
To identify salt-responsive RNAs in C. sinensis, six cDNA libraries from three controls (Control1, Control2, Control3) and three salt-treated (Salt1, Salt2, Salt3) leaf samples that were treated with 250 mM NaCl for 4 h were constructed and sequenced using an Illumina HiSeq2500 platform in our previous report . In this report, the RNA-Seq results were aligned to the latest tea tree genome, and the alignment results are shown in Table 1. The mapped ratio of each sample ranged from 88.99 to 91.10%. After the data were filtered with CPC, CNCI, Pfam, and CPAT, 16,452 candidate lncRNAs were obtained from the six samples ( Figure 1B). Then, the identified lncRNAs were classified into four categories: 11,010 lincRNAs (66.9%), 2,375 antisense-lncRNAs (14.4%), 1,863 intronic-lncRNAs (11.3%), and 1,204 sense-lncRNAs (7.3%) ( Figure 1C).
To understand the differences between lncRNAs and mRNAs in sequence and structure, lncRNAs were analyzed and compared with mRNAs, according to the transcript length and number of exons. The average length of lncRNAs was 822 nucleotides, and the average length of mRNAs was 1,348 nucleotides. As shown in Figure 2A, 75.19% of lncRNAs were < 1,000 nucleotides, while only 54.28% of mRNAs were < 1,000 nucleotides. In addition, 94.44% of lncRNAs had two (83.15%) or three (11.29%) exons, while 10.05% of mRNAs had two exons and 46.99% of mRNAs had only one exon ( Figure 2B).

DE-lncRNAs Involved in NaCl Stress
In this study, 172 DE-lncRNAs in tea plants were identified by DESeq using the FPKM value, including 101 upregulated and 71 downregulated lncRNAs (Figure 3 and Supplementary Table S1), suggesting that these lncRNAs may function in response to salt stress.
To analyze the potential functions of these lncRNAs, the predicted cis-target genes encoding genes 100 kb upstream or downstream of these DE-lncRNAs and trans-target genes predicted by complementary pairing of bases were searched and annotated. In total, 250 cis-and 421 trans-target genes of these DE-lncRNAs were predicted (Supplementary Tables S2, S3). Furthermore, GO and KEGG enrichment analyses were used to investigate the potential functions of the cis-and transtarget genes, respectively. As shown in Figures 4A,B, 116 cistarget and 181 trans-target genes were annotated in the GO   database. For both cis-and trans-target genes, "Cell part" (GO: 0005623) was the largest subcategory in the cellular component category; "Catalytic activity" (GO: 0003824) in the molecular function category, and "Metabolic process" (GO: 0008152) in the biological process category were the most abundant terms. In addition, 50 cis-and 93 trans-target genes were assigned to 50 and 63 KEGG pathways, respectively (Figures 5A,B), which mainly included "Galactose metabolism" (ko00052), "Biosynthesis of amino acids" (ko01230), and "Ribosome" (ko03010). The results showed that these pathways may play important roles in the salt-stress response of tea plants.

LncRNAs Participate in Galactinol
In addition, many TF families play vital roles in regulating plant resistance mechanisms under abiotic stress. In the present study, 23 DE-lncRNA target genes were identified, including HSF, GATA, MADS-box, bHLH, ERF, WRKY, ICE, LHY, and MYC TFs. Except all of the bHLH, LHY and MYC TFs and some of the TFs in GATA and MADS families were downregulated; the rest of the TFs were upregulated ( Table 3).
LncRNA MSTRG.139242.1 May Interact With Its Nearby mRNATEA027212.1 Among all of the target genes of the DE-lncRNAs, we found an upregulated coding gene, TEA027212.1 (Scaffold765 1270655-1272793), which was annotated as Ca 2+ -transporting ATPase 13, that was located only 5,911 bp adjacent to the upregulated lncRNA MSTRG.139242.1 (Scaffold765 1278704-1412939) according to the tea tree genome. To verify their expression relationship, a gene suppression assay in tea  Figure 7A). However, gene suppression of TEA027212.1 from 0 to 6 h significantly reduced MSTRG.139242.1 transcript levels from 6 to 12 h ( Figure 7B). These results indicate that the lncRNA MSTRG.139242.1 may be regulated by its nearby coding gene, TEA027212.1, and may be involved in Ca 2+ transport in response to salt stress in tea plants.

Data Deposition
The sequencing data were deposited in the NCBI Sequence Read Archive (SRA) database 13 under accession number 13 http://www.ncbi.nlm.nih.gov/Traces/sra

DISCUSSION
Salt stress is an important factor that affects the growth of plants, including tea plants, under natural conditions, and the mechanism of the response to salt stress and salt stress-resistance breeding have been studied for a long period. Over the past decade, a number of genes involved in salt stress have been identified and verified, in addition to their regulatory pathways that mediate the transduction of stress signals and the process of the salt-stress response, such as the Ca 2+ signaling pathway, ABA pathway, and MAPK cascade (Boudsocq and Sheen, 2009;Cheong and Kim, 2010;Yoshida et al., 2014;Abbasi et al., 2016;Zhu, 2016). In tea plants, many coding genes, such as AQP, VQ, and SnRK2, participate in the salt-stress response (Yue et al., 2014;Guo et al., 2018;Zhang et al., 2018). However, a few studies have been performed on the identification and mechanism of salt stress-related non-coding RNAs, especially lncRNAs, although lncRNAs have been proven to be involved in many biological processes in plants (Ponjavic et al., 2007;Ponting et al., 2009). In this study, we systematically identified tea plant lncRNAs based on the latest tea tree genome to find lncRNAs associated with salt stress and preliminarily analyzed their possible interactions with mRNAs or miRNAs. A number of lncRNAs were identified for the first time to be specifically expressed under high concentrations of NaCl and involved in the stress response. A total of 16,452 candidate lncRNAs from six tea leaf samples were identified in this study. Similar to the lncRNAs identified in other species (Liu et al., 2012;Shuai et al., 2014;Zhu et al., 2015;Tian et al., 2016;Ma et al., 2018), the average length of lncRNAs in tea plants was much shorter than that of mRNAs, and most lncRNAs (66.9%) were lincRNAs. In some plants, functional lncRNAs that were differentially expressed under specific conditions were identified using transcriptome sequencing, and the functions were validated. For example, Zhu et al. (2015) identified 490 lncRNAs that were significantly upregulated in tomato ripening mutants using RNA-Seq, and two novel lncRNAs, lncRNA1459 and lncRNA1840, were proven to function in the delay of fruit ripening by gene silencing and mutagenesis (Li et al., 2018). Ding et al. (2014) conducted RNA-seq of salt-treated Arabidopsis seedlings, and Qin et al. (2017) obtained a putative lncRNA named DRIR that was later proven to increase tolerance to drought and salt stress. Similarly, 172 lncRNAs with various functions were found to be differentially expressed under NaCl stress in our results. Thus, these lncRNAs are involved in the salt-stress response of tea plants, and the results can provide candidate genes for salt tolerance studies in tea plants.
Previous studies have reported that lncRNAs play important roles in a variety of biological processes in response to stress by acting directly or indirectly on functional genes in plants Fatica and Bozzoni, 2014;Wang et al., 2017). Currently, the interactions between lncRNAs, mRNAs, and miRNAs are the main research hotspots. Thus, in this study, the interactions between lncRNAs and mRNAs, and between lncRNAs and miRNAs were predicted.
To the best of our knowledge, lncRNAs can negatively or positively regulate the expression of protein-coding genes by acting in cis or in trans; lncRNAs work in cis when their target genes are on the same chromosome and within a close distance, and lncRNAs work in trans when they affect genes on other chromosomes (Kornienko et al., 2013). In addition, recent studies found that some cis-acting lncRNAs also have the ability to act in trans (Martianov et al., 2007;Schmitz et al., 2010). In this study, possible cis-and trans-target genes were predicted and annotated to explore the possible functions of the DE-lncRNAs. In total, 250 cis-and 421 trans-target genes of these DE-lncRNAs were predicted. GO and KEGG analyses indicated that these lncRNAs participated in various pathways, such as catalytic activity, galactose metabolism, and biosynthesis of amino acids in response to salt stress. Notably, 42 cis-and 67 trans-target genes were differentially expressed under salt stress. Among these target genes, Ca 2+ -transporting ATPase (TEA027212.1) was reported to be an important gene in Ca 2+ signaling under abiotic stress (Wilkins et al., 2016), which suggests that its corresponding lncRNA, MSTRG.139242.1, may interact with Ca 2+ -transporting ATPase and may be involved in this process. Thus, a gene suppression assay of MSTRG.139242.1 and TEA027212.1 (Ca 2+ -transporting ATPase 13) using AsODNs was conducted and suggested that lncRNA MSTRG.139242.1 may be regulated by its nearby coding gene, TEA027212.1. In a previous study, the tomato lncRNA lncRNA33732 activated by WRKY1 induces RBOH expression and conferred resistance to Phytophthora infestans infection (Cui et al., 2018). In this study, we can speculate that lncRNA MSTRG.139242.1 may interact with its nearby mRNA TEA027212.1 in response to salt stress. In addition, an auxin-responsive protein (TEA005327.1) and ethylene-responsive TF ABR1 (TEA015017.1), the target genes of the lncRNAs MSTRG.49718.1 and MSTRG.54711.1, respectively, have been reported to be involved in the auxin, ethylene, and ABA signal transduction pathways in response to abiotic stresses in other plants (Pandey et al., 2005;Jain and Khurana, 2009;Feng et al., 2015), indicating that these two lncRNAs may participate in hormone signal transduction in response to salt stress in tea plants. We speculate that the expression of these target mRNAs may be regulated by their upstream lncRNAs in response to salt stress.
In addition, it is worth noting that many lncRNA target genes were classified as GOLS, calcium signaling, and some crucial TFs. GOLS was reported to be a key enzyme in the synthesis of raffinose family oligosaccharides (RFOs) and catalyzes the condensation of UDP-galactose with myo-inositol to produce galactinol as the sole donor for the synthesis of RFOs. In addition, RFOs are used for the transport and storage of carbohydrates and as compatible solutes for protection against abiotic and biotic stresses (Peters and Keller, 2009;Philippe et al., 2010). Sun Z.B. et al. (2013) have reported that TsGOLS2 enhances tolerance to high salinity and osmotic stresses in Arabidopsis thaliana. Thus, in the present study, lncRNAs targeting the eight upregulated GOLS genes may participate in RFO synthesis in response to salt stress. The Ca 2+ signaling pathway was also reported to mediate plant response to salt stress, which starts with Ca 2+ transporters such as Ca 2+ -ATPases, Ca 2+ sensors, and relay proteins such as CaM, CMLs, CDPKs, CBLs, and CIPKs (Wilkins et al., 2016). LncRNAs targeting calcium-transporting ATPase, calmodulin-interacting protein, CDPKs and CIPKs in this study were predicted to be involved in salt-stress response through the Ca 2+ signaling pathway. TFs have long been an important topic in plant resistance research, and some TF families such as WRKY, bHLH, ERF, GATA, and MYC play significant roles in salt stress adaptation (Lindemose et al., 2013). Here, the results provide many candidate lncRNAs, which may interact with vital TFs and response to salt stress.
On the other hand, lncRNAs can also interact with miRNAs. Some miRNAs have been found to respond to different environmental stresses (Zhang, 2015). Some lncRNAs can bind to miRNAs, thus inhibiting the effect of miRNAs on downstream genes. This type of lncRNA is termed an endogenous target mimic (eTM) . Jiang et al. (2019) reported that the tomato lncRNA23468 functions as an eTM to modulate NBS-LRR genes by mimicking miR482b in the tomato-P. infestans interaction. In this study, 12 lncRNAs were predicted to be target mimics of 17 known mature miRNAs. Interestingly, five lncRNAs, MSTRG.56302.5, MSTRG.116911.1, MSTRG.92784.12, MSTRG.36615.10, and MSTRG.55773.2, were predicted to be eTMs of miR156. According to previous reports, miR156 is induced by drought, and salinity stresses in several different plant species specifically regulates downstream SPL TFs and affects plant morphology, phase change, and seed germination (Stief et al., 2014). Thus, the five lncRNAs may influence the effect of miRNA156 on SPL in response to salt stress. These results provide candidate lncRNAs for exploring the interactions between lncRNAs and miRNAs in response to salt stress.

CONCLUSION
In the present study, genome-wide identification of lncRNAs was performed using high-throughput RNA-Seq and bioinformatics analysis. lncRNAs (16,452) were identified in tea plants, among which, 172 lncRNAs were differentially expressed and responsive to salt stress by cis-or trans-interaction with important coding genes. These lncRNAs may participate in GOLS, calcium signaling pathway, and interact with TFs in response to salt stress. Notably, 35 DE-lncRNAs were predicted to interact with 42 differentially expressed coding genes, which may participate in pathways such as the auxin response, ABA, and Ca 2+ signal transduction pathways under salt stress. AsODN suppression of the lncRNA MSTRG.139242.1 and its predicted interacting gene, TEA027212.1 (Ca 2+ -ATPase 13), indicated that MSTRG.139242.1 interacts with Ca 2+ -ATPase 13 in the Ca 2+ -transport pathway in response to high salinity in tea plants. In addition, 12 lncRNAs were predicted to be target mimics of 17 known mature miRNAs in C. sinensis, thereby affecting the expression of downstream functional genes. This study can provide a source of lncRNAs and benefit an in-depth understanding of the function and regulatory mechanisms in tea plant response to salt stress.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
YYu and YYa designed the experiments. SW and YZ performed the experiments and data analysis. SW and YZ wrote the manuscript. WW and QX provided valuable advice on the manuscript. LH and MD revised the manuscript. All authors discussed the results and contributed to the manuscript.