Introduction
Soil salinization has become a major challenge for sustainable development of global agriculture. As a result, cultivation of salt-tolerant crop varieties has become a focus of plant breeding. However, development of effective breeding strategies would be significantly enhanced by improving our understanding of salt tolerance mechanisms in plants and identifying genes required for adaptation (Flowers and Colmer, ; Shabala, 2013; Ismail and Horie, ).
Na+ is not essential for plant survival; in fact, there are many halophytes whose growth is improved upon the addition of a certain degree of NaCl (Bose et al., ). As a result of this unique survival strategy, halophytes have been used to elucidate a series of mechanisms that use to cope with salt stress. In particular, succulent halophytes, a representative type of anatomical adaptation in plants by increasing cell size due to the expansion of vacuole volume under salinity, have been suggested to take advantage of tissue- and organelle-specific ion compartmentalization to maintain high K+/Na+ ratios in response to high salinity (Munns and Tester, 2008; Shabala and Mackay, 2011). Furthering our molecular understanding of the ion transport systems and compartmentalization mechanisms used by halophytes is crucial for mapping the complex salt tolerance network in plants in order to optimize breeding of salt-tolerant crops.
Halogeton glomeratus (H. glomeratus), an annual, succulent halophyte, is a typical salt-tolerant halophyte with highly succulent leaves that has been used to characterize the physiological, cellular, and molecular mechanisms involved in the response to salt stress (Wang et al., 2015a,b, 2016).
The majority of next-generation transcriptome sequencing (NGS) technologies, or RNA sequencing (RNA-Seq) are not capable of assembling full-length transcripts due to the short lengths of sequencing-reads, which are likely to produce partial reference transcriptomes (Martin and Wang, ). For the RNA-Seq of non-model organisms, one of the principal strategies is de novo assembly. In this way, we could retrieve plenty of transcripts by piecing together short, low-quality sequencing-reads in a more economical approach, and provide a new opportunity to cognize the transcriptome characteristics of non-model organisms, but the disadvantages of the de novo assembly are obvious, such as a greater sequencing depth, highly sensitive to sequencing errors, presence of chimeric molecules in the sequencing dataset, and difficult to distinguish highly similar transcripts (Martin et al., ; Martin and Wang, ). In a word, it is still challenging to construct of a comprehensive transcriptome from short sequencing-reads. However, recent development of Pacific Biosciences (PacBio) long-read technologies for transcriptome sequencing (Iso-Seq) can overcome those limitations to capture the full catalog of transcripts and their variants (Sharon et al., 2013; Minoche et al., ). Thus, for the first time, we applied the Iso-Seq protocol to analyze the H. glomeratus transcriptome. These findings serve to improve H. glomeratus transcriptome annotation and provide abundant gene resources related to improve salt tolerance in plants.
Value of the data
Halogeton glomeratus is not only highly valuable for the characterization of salt tolerance mechanisms in halophytes, but may also serve as a source of stress tolerance genes for plant breeding.
So far, there is lack of well-characterized transcriptomic profiling of H. glomeratus. Here, we provided the first comprehensive transcriptome analysis using single-molecule sequencing.
These data will help elucidate salt tolerance mechanisms of H. glomeratus and provide salt-tolerant gene resources for developing stress-tolerant plants.
Data
Transcriptome plant materials and NaCl treatments
Halogeton glomeratus seedlings were grown in plastic pots in a growth chamber and were irrigated daily with half-strength Hoagland's nutrient solution as previously described by Wang et al. (2015a). One month after germination (August 10, 2016), two independent experiments were conducted to assess different durations of salt stress. For the first experiment, seedlings were transferred to half-strength Hoagland's nutrient solution supplemented with 0, 100, 200, or 400 mM NaCl for a period of 3 days. The leaf samples were collected for each condition. For the second experiment, the seedlings were transferred to half-strength Hoagland's nutrient solution supplemented with 200 mM NaCl for 0, 2, 6, 24, or 72 h. The root samples were collected for each condition and were rinsed with distilled water. All samples were frozen in liquid N2 immediately.
Generation of a barcoding library and single-molecule sequencing
After grinding the tissue samples, total RNA was extracted from each sample using the RNeasy Plant Mini kit (Qiagen; Hilden, Germany) and quantified using an Agilent 2100 bioanalyzer. The total RNA from the 9 samples were pooled for sequencing of an H. glomeratus reference transcriptome using single-molecule long-read technology (Iso-Seq) via the PacBio RS II sequencing platform.
In brief, total RNA (10 ng) was reverse transcribed using the SMARTer PCR cDNA Synthesis Kit (Clontech; Mountain View, CA, USA). After polymerase chain reaction (PCR) optimization, large-scale PCR reactions were performed to synthesize second-strand cDNA, which was subjected to BluePippin™ size selection using the following bins: 1–2, 2–3, 3–6, and 5–10 kb. After size selection, another round of amplification was performed using 12 PCR cycles. The amplified and size-selected cDNA products were used to generate SMRTbell template libraries for sequencing based on the Iso-Seq protocol. A total of 12 SMRT cells (Table 1) were sequenced on a PacBio RS II platform (Completed on January 2, 2017). Raw sequence reads (FASTQ format) are available through the NCBI database (https://www.ncbi.nlm.nih.gov/) with the following accession information (Data was released on September 22, 2017, size is 1.6 Gb):
BioProject ID: PRJNA359784.
BioProject ID: SAMN06298282.
Table 1
| Database | Cell number | All reads | Five prime reads (%) | Three prime reads (%) | Poly-A reads (%) | Full-length non-chimeric reads (%) | Full-length non-chimeric read length (bp) |
|---|---|---|---|---|---|---|---|
| 1–2 kb | 4 | 152,748 | 68.01 | 69.89 | 69.12 | 57.96 | 1,657 |
| 2–3 kb | 4 | 176,650 | 60.10 | 62.56 | 62.02 | 50.05 | 2,848 |
| 3–6 kb | 2 | 52,148 | 72.13 | 74.09 | 73.44 | 58.28 | 3,835 |
| 5–10 kb | 2 | 51,874 | 53.04 | 61.68 | 59.97 | 35.33 | 6,221 |
Proportion of SMRT reads containing cDNA primers and poly(A) tail and their classification.
Bioinformatics analysis of isoform sequencing
Raw SMRT sequencing reads were classified as full-length (FL) non-chimeric, FL chimeric, non-FL, or short reads based on detection of 5′- and 3′-primer and polyA-tail sequences using SMRT Analysis v2.3.0 (http://smrt-analysis.readthedocs.io/en/latest/SMRT-Analysis-Software-Installation-v2.3.0/). Only reads with both 5′- and 3′-primers and polyA-tail sequences were considered to be FL. FL non-chimeric reads were clustered into consensus sequence (CS) reads using the ICE (Interative Clustering and Error Correction) algorithm and further polished using Quiver. Finally, high quality CS reads from each library were merged, and redundancy was removed to identify unique isoforms present in the H. glomeratus transcriptome using the Cd-hit program (Li and Godzik, ) for further analysis. All of the isoform sequences are publicly available on Figshare at https://figshare.com/articles/Single-molecule_long-read_transcriptome_dataset_of_halophyte_Halogeton_glomeratus/5345464.
A total of 433,420 reads of insert (ROIs) were generated. Full-length (FL) non-chimeric transcripts made up 57.96% of the 1–2 kb transcripts, 50.05% of the 2–3 kb transcripts, 58.28% of the 3–6 kb transcripts, and 35.33% of the 5–10 kb transcripts, and the libraries were characterized by average lengths of 1,657, 2,848, 3,835, and 6,221 bp, respectively (Table 1). The FL non-chimeric transcripts were processed using Quiver to yield 64,910 high quality consensus isoforms. After additional merging and discarding of redundancies, we obtained 54,835 consensus isoforms with a mean length of 2,663 nt and a N50 length of 2,937 nt (Figures 1A,B).
Figure 1
The isoform sequences were used as queries for sequence homology searches (E < 0.0001) performed with Blast, Blast2GO (Conesa et al., ), and InterProScan5 (Quevillon et al., 2005) to identify functional annotation terms from the non-redundant protein (NR), non-redundant nucleotide (NT), Gene Ontology (GO), Clusters of Orthologous Groups (COG), Kyoto Encyclopaedia of Genes and Genomes (KEGG), SwissProt, and Interpro databases. For isoform coding DNA sequence (CDS) predictions, we selected the segment of the isoform that best mapped to the functional databases as its CDS. For isoforms that were not annotated in the databases, we used ESTScan (Iseli et al., ) to predict CDSs. The predict CDSs of isoforms are publicly available on Figshare at: https://figshare.com/articles/Single-molecule_long-read_transcriptome_dataset_of_halophyte_Halogeton_glomeratus/5345563.
A total of 53,230 (97.07%) isoforms were matched to known genes, and 52,141 (95.07%) isoforms were identified to have intact coding DNA sequences (CDS) (Figures 1C–F). At the same time, our previous Illumina RNA-Seq data of H. glomeratus (Wang et al., 2015b) was re-annotated by the update versions of the database used in the present Iso-Seq transcriptome sequencing, including NR, NT, GO, COG, KEGG, SwissProt, and Interpro databases. In comparison to the results from Illumina RNA-Seq analysis of H. glomeratus, the length distribution, functional annotation, and coding sequence quantity of the Iso-Seq transcripts were significantly improved (Figure 1). In particular, with respect to the species distribution of annotation from the NR database, we found that 98.31% of the annotated isoforms showed the highest similarity to sequences from the three most prevalent species: beta vulgaris subsp. vulgaris (54.72%), spinacia oleracea (28.33%), and vitis vinifera (15.31%) (Figure 1E). Like H. glomeratus, beta vulgaris subsp. vulgaris and spinacia oleracea belongs to Chenopodiaceae. In contrast, for unigenes identified using Illumina RNA-Seq, the top species identified was vitis vinifera (16.19%), and none of the top five species belong to Chenopodiaceae (Figure 1I). Besides, to better understand the advantages of the PacBio Iso-Seq reference transcriptome, we mapped the Illumina de novo assembled transcripts (50,267) to the Iso-Seq consensus transcripts (54,835). Among Iso-Seq consensus transcripts, a total of 21,008 exhibited homology to Illumina de novo assembled 36,548 transcripts (E < 0.0001), of which 8,236 (39.20%) exhibited homology with at least two Illumina de novo assembled transcripts, including 4,668 (22.22%), 1,860 (8.85%) and 1,708 (20.74%) exhibited homology with two, three, and the more Illumina de novo assembled transcripts, respectively (Supplementary Figure 1). These results are successfully consistent with comprehensive assembly and annotation of a reference H. glomeratus transcripts using Iso-Seq technology that is substantially more contigous and contains a higher proportion of intact CDS.
Statements
Author contributions
HW designed the experiments. JW, LY, BL, YM, and XM collected the samples and performed the experiments. JW, LY, and BL analyzed RNA-Seq data. JW and LY wrote the manuscript. All authors read and approved the final manuscript.
Acknowledgments
This work was supported by National Basic Research Program of China (973 program, 2014CB160313), China Agriculture Research System (CARS-05), National Natural Science Foundation of China (No. 31460347), and Fostering Foundation for the Excellent Ph.D. Dissertation of Gansu Agricultural University (YBPY2014001).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2017.00197/full#supplementary-material
Supplementary Figure 1Statistics of the correspondence between the Illumina unigenes and PacBio reference transcripts of H. glomeratus transcriptome.
References
1
BoseJ.Rodrigo-MorenoA.ShabalaS. (2014). ROS homeostasis in halophytes in the context of salinity stress tolerance. J. Exp. Bot. 65, 2808–2837. 10.1093/jxb/ert430
2
ConesaA.GötzS.García-GómezJ. M.TerolJ.TalónM.RoblesM. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics21, 3674–3676. 10.1093/bioinformatics/bti610
3
FlowersT. J.ColmerT. D. (2008). Salinity tolerance in halophytes. New Phytol.179, 945–963. 10.1111/j.1469-8137.2008.02531.x
4
IseliC.JongeneelC. V.BucherP. (1999). ESTscan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol.99, 138–148.
5
IsmailA. M.HorieT. (2017). Genomics, physiology, and molecular breeding approaches for improving salt tolerance. Annu. Rev. Plant Biol.68, 405–434. 10.1146/annurev-arplant-042916-040936
6
LiW.GodzikA. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics22, 1658–1659. 10.1093/bioinformatics/btl158
7
MartinJ. A.WangZ. (2011). Next-generation transcriptome assembly. Nat. Rev. Genet.12, 671–682. 10.1038/nrg3068
8
MartinJ.BrunoV. M.FangZ.MengX.BlowM.ZhangT.et al. (2010). Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads. BMC Genomics11:663. 10.1186/1471-2164-11-663
9
MinocheA. E.DohmJ. C.SchneiderJ.HoltgräweD.ViehöverP.MontfortM.et al. (2015). Exploiting single-molecule transcript sequencing for eukaryotic gene prediction. Genome Biol.16:184. 10.1186/s13059-015-0729-7
10
MunnsR.TesterM. (2008). Mechanisms of salinity tolerance. Plant Biol. 59, 651–681. 10.1146/annurev.arplant.59.032607.092911
11
QuevillonE.SilventoinenV.PillaiS.HarteN.MulderN.ApweilerR.et al. (2005). InterProScan: protein domains identifier. Nucleic Acids Res.33, 116–120. 10.1093/nar/gki442
12
ShabalaS. (2013). Learning from halophytes: physiological basis and strategies to improve abiotic stress tolerance in crops. Ann. Bot.112, 65–69. 10.1093/aob/mct205
13
ShabalaS.MackayA. (2011). Ion transport in halophytes. Adv. Bot. Res.57, 151–187. 10.1016/B978-0-12-387692-8.00005-9
14
SharonD.TilgnerH.GrubertF.SnyderM. (2013). A single-molecule long-read survey of the human transcriptome. Nat. Biotechnol.31, 1009–1014. 10.1038/nbt.2705
15
WangJ.MengY.LiB.MaX.LaiY.SiE.et al. (2015a). Physiological and proteomic analyses of salt stress response in the halophyte Halogeton glomeratus. Plant Cell Environ.38:655–669. 10.1111/pce.12428
16
WangJ.LiB.MengY.MaX.LaiY.SiE.et al. (2015b). Transcriptomic profiling of the salt-stress response in the halophyte Halogeton glomeratus. BMC Genomics16:169. 10.1186/s12864-015-1373-z
17
WangJ.YaoL.LiB.MengY.MaX.LaiY.et al. (2016). Comparative proteomic analysis of cultured suspension cells of the halophyte Halogeton glomeratus by iTRAQ provides insights into response mechanisms to salt stress. Front. Plant Sci.7:110. 10.3389/fpls.2016.00110
Summary
Keywords
Halophyte, H. glomeratus, Iso-Seq, transcriptome, salt stress
Citation
Wang J, Yao L, Li B, Meng Y, Ma X and Wang H (2017) Single-Molecule Long-Read Transcriptome Dataset of Halophyte Halogeton glomeratus. Front. Genet. 8:197. doi: 10.3389/fgene.2017.00197
Received
25 August 2017
Accepted
20 November 2017
Published
01 December 2017
Volume
8 - 2017
Edited by
Graziano Pesole, Università degli studi di Bari Aldo Moro, Italy
Reviewed by
RaffaEle A. Calogero, Università degli Studi di Torino, Italy; Matteo Chiara, Università degli Studi di Milano, Italy
Updates
Copyright
© 2017 Wang, Yao, Li, Meng, Ma and Wang.
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) or licensor 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: Huajun Wang whuajun@yahoo.com
This article was submitted to RNA, a section of the journal Frontiers in Genetics
†These authors have contributed equally to this work.
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.