Next Article in Journal
Pericytes in Microvessels: From “Mural” Function to Brain and Retina Regeneration
Previous Article in Journal
Sequential Transcriptome Changes in the Penumbra after Ischemic Stroke
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 24
10.3390/ijms20246350
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single-Molecule Long-Read Sequencing Reveals the Diversity of Full-Length Transcripts in Leaves of Gnetum (Gnetales)

by
Nan Deng
1,2,†,
Chen Hou
3,†,
Fengfeng Ma
1,
Caixia Liu
1,* and
Yuxin Tian
1,2,*
1
Institute of Ecology, Hunan Academy of Forestry, Changsha 410004, China
2
Hunan Cili Forest Ecosystem State Research Station, Cili 417100, China
3
School of life Sciences, Sun Yat-sen University, Guangzhou 510275, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2019, 20(24), 6350; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20246350
Submission received: 12 October 2019 / Revised: 1 December 2019 / Accepted: 11 December 2019 / Published: 17 December 2019
(This article belongs to the Section Molecular Biology)
Browse Figures
Versions Notes

Abstract

:
The limitations of RNA sequencing make it difficult to accurately predict alternative splicing (AS) and alternative polyadenylation (APA) events and long non-coding RNAs (lncRNAs), all of which reveal transcriptomic diversity and the complexity of gene regulation. Gnetum, a genus with ambiguous phylogenetic placement in seed plants, has a distinct stomatal structure and photosynthetic characteristics. In this study, a full-length transcriptome of Gnetum luofuense leaves at different developmental stages was sequenced with the latest PacBio Sequel platform. After correction by short reads generated by Illumina RNA-Seq, 80,496 full-length transcripts were obtained, of which 5269 reads were identified as isoforms of novel genes. Additionally, 1660 lncRNAs and 12,998 AS events were detected. In total, 5647 genes in the G. luofuense leaves had APA featured by at least one poly(A) site. Moreover, 67 and 30 genes from the bHLH gene family, which play an important role in stomatal development and photosynthesis, were identified from the G. luofuense genome and leaf transcripts, respectively. This leaf transcriptome supplements the reference genome of G. luofuense, and the AS events and lncRNAs detected provide valuable resources for future studies of investigating low photosynthetic capacity of Gnetum.

1. Introduction

Gnetum L., Ephedra L., and Welwitschia Hook.f. together form a monophyletic group—order Gnetales [1,2]. The phylogenetic relationship of the Gnetales with other seed plant groups has been a longstanding but controversial question of seed plant evolution [3,4,5,6,7]. On the basis of morphological and anatomical data, the Gnetales was inferred to be a sister group of angiosperms [2,8]. However, the most recent molecular phylogenies have placed the Gnetales as a sister group to Pinaceae or entire conifers [3,4,5,9]. Genus Gnetum, comprising approximately 40 species, are woody climbers, with a few exceptions being shrubs and trees [1,10,11]. Most of Gnetum species, such as the Asian lianoid Gnetum, are rich in bioactive compounds, which can be used for medicine [12,13,14].
Gnetum are characterized by their pinnate leaf veins, decussate leaves, and the presence of vessels in the stems [15]. The photosynthetic capacity of Gnetum was found to be lower than that of other seed plant groups, suggesting that this character acts as an intrinsic property and reflects a unique evolutionary history [1,16]. A very recent study uncovered the unique morphology of the Gnetum stomata [17], which to some extent explained the reasons for the low photosynthetic capacity, but underlying molecular mechanisms remain unknown. Wan et al. [3] proposed that cellulose synthase (CesA), cellulose synthase-like (Csl), and the WUSCHEL-related homeobox (WOX) family could play important roles of affecting gross morphology and leave development in Gnetum. These reported genes, however, are insufficient to explain the low photosynthetic capacity with regard to photosynthetic rate, stomatal conductance, and water transport capacity [1]. The bHLH (basic/helix-loop-helix) gene family has been reported to regulate stomatal development and photosynthesis, e.g., to initiate the asymmetrical division of protodermal cells, to terminate the meristematic stem cell identity, and to trigger flowering in response to blue light [18,19,20]. The bHLH family is characterized by bHLH (hidden Markov models, HMM accession: PF14215.6) and HLH (HMM accession: PF00010.26) domains [21]. It is the second-largest gene family in Arabidopsis (147 genes) and is also present in other plants (ranging from 1 to 553 genes). The discovery of the bHLH family promotes the understanding of molecular mechanisms that regulate low photosynthetic capacity in Gnetum.
Transcriptomes provide valuable gene resources for those gymnosperms whose reference genome is lacking. For those species that already have their whole genome sequenced, the transcriptomic data act as a good supplement to the genomic data, such as the detection of alternative splicing (AS) events and long non-coding RNAs (lncRNAs). Alternative splicing, which generates more than one transcript mediated by spliceosome (a large protein and RNA complex), is a major post-transcriptional regulatory mechanism in plants [22,23]. Alternative splicing has been proven to be involved in the regulation of growth, development, signal transduction, flowering, and responses to various environmental cues in plants [24,25,26,27]. Additionally, in recent years, lncRNAs have been reported as important regulators of gene expression, such as precursors of microRNAs (miRNAs) and miRNA target mimics [28,29].
At present, most transcriptomic data are generated using short reads sequencing, e.g., Illumina sequencing platforms. However, it is challenging to detect the presence of different isoforms, repetitive sequences, and transposable elements, because different transcripts which possess identical exons cannot be distinguished [30]. Thus, it is difficult to predict AS events and lncRNAs using short-read sequencing [31]. With the development of third-generation sequencing, the sequencing read length has grown rapidly, providing a more direct view of RNA molecules [32]. The third-generation sequencing technology overcomes the difficulties of short-read sequencing by generating the full-length sequence in a single molecule [33]. Currently, the average read length for single-molecule real-time (SMRT) sequencing is 12–15 kb, and the maximum length can reach 60 kb (PacBio Sequel platform). To date, SMRT sequencing has been applied to investigate full-length transcriptome in quite a few angiosperms species, e.g., sorghum [34], coffee [33], sugarcane [35], cotton [36], and the rubber tree [37], but barely been applied in gymnosperms.
This study performed the first SMRT transcriptome derived from Gnetum luofuense C.Y. Cheng leaves at different developmental stages using PacBio Sequel platform. To improve the sequence quality generated by PacBio Sequel platform, Illumina sequencing platforms were used to generate short reads. The annotation and structure analysis of the full-length transcripts are expected to improve the annotation of the reference genome of G. luofuense and to better understand the complexity of transcriptome in gymnosperms. Besides, the newly generated data can serve as the reference for differentially expressed analyses of G. luofuense leaves in further studies.

2. Results

2.1. Transcriptome from PacBio Sequel Sequencing

To reveal the complexity of the transcriptome in G. luofuense leaves, a pooled sample at different developmental stages was sequenced with the PacBio Sequel platform. In total, 9.98 GB of raw data was generated, made up of 3,689,825 subreads with an average length of 2705 bp and N50 of 3168 bp. A total of 185,089 circular consensus sequences (CCSs) were obtained, of which the full-length reads were 143,578 (Table 1). The full-length non-chimeric (flnc) reads were characterized by possessing complete 3′/5′ terminal primers and poly(A) tails. A total of 139,488 flnc reads were generated from G. luofuense leaf transcriptome, accounting for 75% of the CCSs (Table 1). After error correction, 80,496 polished consensus reads were obtained from Quiver, ranging from 167 bp to 14,735 bp (N50 of 3614 bp). The length distribution of the subreads, CCSs, and flnc reads was shown in Figure S1. To decrease PacBio Sequel sequencing errors, all PacBio sequenced transcripts were improved by the comparisons to the Illumina sequenced reads, whereupon the total numbers of nucleotides, mean lengths, and N50 were increased (Table 2). The lengths of de novo assembled unigenes were much shorter with the mean length of 972 bp than those generated from PacBio Sequel sequencing.

2.2. Transcript Structure Analysis

After error correction with the Illumina-sequenced data, 77,380 (96.31%) reads were mapped to the reference genome of G. luofuense using GMAP, leaving 3116 reads unmapped. Among the mapped reads, 6678 (8.30%) were mapped to multiple locations on the reference genome, 70,702 (87.83%) were uniquely mapped, and 43,299 (53.79%) and 27,403 (34.04%) were mapped to the positive strands and negative strands, respectively (Table 3, Figure 1A). Of the mapped reads, high-quality reads (i.e., with coverage and identity values over 98%) accounted for over 80% (Figure 1B), and the curve of the corrected isoform numbers reached a saturation level (Figure 1C). Overall, 93.01% of the de novo reads were successfully mapped to the reference genome (Table 3), where 5665 polished consensus reads (16.48%) were classified as isoforms from known genes, 23,443 (68.19%) were classified as novel isoforms from known genes, and 5269 (15.33%) were novel isoforms from novel genes (Figure 1D, File S1). The mean length and N50 of all isoforms were 3024 bp and 3394 bp, respectively. Figure S2A shows the density and number of the mapped reads on the reference genome of G. montanum (= G. luofuense). For 96.31% of all generated reads, scaffold498063 (the longest scaffold) had the most abundant mapped reads, whereas scaffold809851 and scaffold761035 had the lowest (Figure S2B).

2.3. Transcript Annotation and Classification

All the consensus isoforms were annotated by querying seven databases, i.e., gene ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), Eukaryotic Orthologous Groups/Clusters of Orthologous Groups (KOG/COG), Protein Family (Pfam), non-redundant protein sequences (NR), nucleotide sequence (NT), and Swiss-Prot. A total of 8980 isoforms had hits on all seven databases, and 34,667 had hits on at least one database. The NR database had the highest number of isoform annotations (34,170), followed by KEGG (33,813), whereas the NR database (15,467) had the lowest (Table 4). Furthermore, 5269 isoforms of novel genes were annotated against at least one database, and 116 isoforms had at least one hit in all databases. The NR database had the highest number of hits (3782). Additionally, 23,443 novel isoforms of known genes were searched against at least one database, and 6548 annotations of known genes had at least one significant hit on all databases, with the most being in NR database (22,167). All NR annotations were distributed among 256 species, where most of the consensus isoforms were homologous to those from Picea sitchensis (Bong.) Carrière (8319), Amborella trichopoda Baill. (4153), and Nelumbo nucifera Gaertn. (2606) (Figure 2).
All the GO annotations were assigned to 51 GO categories (GO level 2, Figure 3). Among them, ‘binding’ (GO:0005488) represented the largest group (4056, 18.2%), followed by ‘metabolic process’ (GO:0008152) (3248, 14.6%) and ‘catalytic activity’ (3240, 14.5%). The 4006 novel genes were assigned to 43 GO categories, with the top three categories, i.e., ‘binding’, ‘catalytic activity’, and ‘metabolic process’. Besides, a total of 288 KEGG pathways were identified, of which novel genes were involved in 225 KEGG pathways. Among these KEGG pathways, the pathway ‘metabolism’ had the most abundant annotations of all genes and novel genes, followed by the orthologous term ‘signal transduction’ (908 for all genes and 146 for novel genes) (Figure 4).

2.4. Identification of LncRNAs and Fusion Genes

In total, 8890, 3955, 8871, and 5540 lncRNAs were identified by the Coding-Non-Coding-Index (CNCI), Coding Potential Calculator (CPC), Pfam, and PLEK databases, respectively (Figure 5A). Additionally, 1660 lncRNAs were identified by the four methods, with full lengths mainly ranging from 1000 bp to 5000 bp and lengths on average ranging from 203 bp to 8106 bp (File S2, Figure 5B). Among them, only 10 lncRNAs were identified as known transcripts whereas the others were identified as novel transcripts (Figure S3). All detected lncRNAs were subdivided into the following four types: lincRNAs (long intergenic non-coding RNAs, 556, 33.49%), sense intronic lncRNAs (525, 31.63%), antisense lncRNAs (487, 29.34%), and sense overlapping lncRNAs (92, 5.54%) (Figure 5C). Most of the identified lncRNAs had one or two exon(s), which was different from the exon distribution of regular mRNAs (Figure 5D). Additionally, 1174 fusion gene events were identified, of which 533 could be identified among the known genes, whereas the remaining fusion genes could not be linked to any known genes, probably there are still many undetected genes in the reference genome of G. luofuense (Figure S3).

2.5. Transcription Factor, Alternative Polyadenylation, and Alternative Splicing Analyses

A total of 1974 transcription factors (TFs) were identified by iTAK (Ithaca, NY, USA), of which 1824 TFs were annotated to 82 families (Figure 6A). Sucrose nonfermenting 2 (SNF2) accounted for the greatest proportion in the known TF families, followed by coumarate-3-hydroxylase (C3H) (Figure 6A). In addition, 5647 genes in the Gnetum leaves were found to have at least one supported poly(A) site (Figure 6B). The transcripts with one poly(A) site were in the highest proportion (50.40%), followed by genes with two poly(A) sites (24.14%). The results of APA analysis were presented in File S3.
A total of 12,998 AS events (assigned to 4459 genes) were detected, comprising 5,400 (41.54%) retained introns, 2,876 (22.13%) alternative 3′ splice sites, 2043 (15.27%) alternative 5′ splice sites, 1290 (9.92%) skipped exons, 1077 (8.29%) alternative first exons, 276 (2.12%) alternative last exons, and 36 (0.28%) mutually exclusive exons (Figure 6C). Additionally, an alternative 3′ splice site was detected in all the transcripts from novel genes.

2.6. Phylogenetic Analysis of bHLH Genes in G. luofuense

A total of 67 bHLH genes were detected by searching against the reference genome of G. luofuense with regard to the bHLH or HLH domain. The information of the bHLH gene family members in Gnetum was presented in File S4. Besides the bHLH and HLH domains, the KIX_2 (HMM accession: PF16987.5), AAA_33 (HMM accession: PF13671.6), Macro (HMM accession: PF01661.21), DcpS_C (HMM accession: PF11969.8), and zf-C2HE (HMM accession: PF16278.5) domains were also identified in the Gnetum bHLH gene family. A rooted neighbor-joining tree of the 67 bHLH proteins from G. luofuense was shown in Figure S4. A further concentrated phylogeny of 30 bHLH genes (including 15 novel genes) detected in the full-length transcriptome was reconstructed (Figure 7), four subfamilies were separated from one another, however some deep nodes had low statistical support. In addition, almost all the detected bHLH proteins contained conserved motifs 1 and 2 and the two motifs were close to each other, suggesting the presence of the two conserved domains. In contrast, motif 5 had the second widest distribution, whereas motif 3 and 4 were only found in one gene (TnS0004498063t28).

3. Discussion

In the present study, the third-generation sequencing technology was applied to generate full-length transcriptome of G. luofuense leaves. 80,496 polished consensus reads were obtained with an average length of 3223 bp. The PacBio Sequel platforms gained considerably longer transcripts than those generated in Illumina platforms, providing a superior overview of the Gnetum transcriptome. After mapping the reads of PacBio and Illumina sequencing against the reference genome, respectively, we were able to visualize the coverages and AS events (Figure 8). Fifteen isoforms with a high coverage were detected (with lengths of over 5000 bp) from the bHLH gene family, providing a huge advantage over the short reads which needed assembly. PacBio Sequel sequencing is subject to higher rates of errors [38], and a previous study reported error rates could reach to 11–14% [39]. Therefore, in this study, we applied Illumina-sequenced data to correct errors generated from Pacbio Sequel sequencing and made the complete assessment of our assembled full-length transcriptome.
LncRNAs, which are operationally defined as RNA genes with a length of over 200 bp, have elusive functions as they are not responsible for protein coding [40]. Although lncRNAs play an important role in gene regulation in plants [28], their numbers, characteristics, and genetic patterns remain unclear [40]. So far, lncRNAs have been identified in many angiosperms, such as A. thaliana [41], Oryza sativa [42], and Morus notabilis [43], very few lncRNAs have been identified in gymnosperms. The lncRNAs in plants have been known to participate in root, stem, and leaf development [44]. For example, lncRNA-HID1 acted as a factor of promoting photomorphogenesis in light [45]. Another case shows that the overexpression of lncRNA-npc48 resulted in A. thaliana leaf serration and delaying flowering time [46]. Moreover, lncRNA (TWISTED LEAF) played an important cis-regulatory role to regulate the expression of gene OsMYB60 during leaf development [47]. Function prediction of lncRNAs is challenging owing to the lack of the homology between closely related species [48]. In this study, we used four methods to identify lncRNAs in Gnetum leaves, whereupon 1660 transcript sequences were identified as putative lncRNAs and most of them contained one exon or two exons. These results provide a base to better understand the varied roles of lncRNAs in gymnosperms.
A large number (5269, 15.33% of total transcripts) of isoforms from novel genes were detected but with a lack of annotation. The newly detected genes are deemed to enrich the knowledge of the reference genome of G. luofuense, although our sampling for long-read transcriptome sequencing is restricted to leaves. The predicted novel genes and novel isoforms may provide valuable opportunities to figure out the gene functions involved in leaf development of G. luofuense. In addition, 23,443 novel isoforms of known genes were detected, the result suggests that AS events frequently occur in Gnetum. It is more confident to predict AS events based on full-length transcriptome than those assembled by short reads because of the complexity of the AS mechanism in eukaryotic cells [35]. Our results indicate that the retention of introns was the main AS type found in the Gnetum leaves, congruent with sugarcane [35] and strawberry [49] as reported in previous studies.
Members of bHLH gene family diverge largely between gymnosperms and angiosperms, the later possesses an even larger range of gene numbers from 150 to 553. To date, Picea abies is known to have the most abundant bHLH genes (107) among the gymnosperms, followed by G. montanum (64). In this study, 30 bHLH genes including 15 novel genes were identified from the leaf transcriptome data, but the deep divergence of the phylogenetic tree received poor support. This result is also found in other organisms, it is probably because the conserved domains of bHLH genes are short and the remaining regions are highly divergent [50,51]. According to the neighbor-joining tree, the bHLH genes identified in G. luofuense leaves were resolved to four clades, but phylogenetic relationships between the different clades were vague because of the poor support. The bHLH genes may have recent common origins derived from several genomic duplication events, and the sequence divergence outside the domains may have experienced extensive shuffling events afterward [52].

4. Materials and Methods

4.1. Samples Selection and RNA Extraction

Leaves of G. luofuense were collected from a bamboo garden at Sun Yat-sen University, Guangzhou, China on 8–9 May 2018. To obtain a good representation of the Gnetum transcriptome, leaves at different developmental stages (from young to old) were collected to cover the entire transcriptome. All the leaf samples for RNA isolation were stored in RNA protection reagent (Qiagen, Hilden, Germany) at −20 ℃.
The total RNA of the sample was extracted using a TRIzol kit (Invitrogen, Carlsbad, CA, USA) according to a previously described two-step protocol [53]. The RNA quality, integrity, and quantity were determined using the following four methods: 1) the RNA samples were examined by agarose gel electrophoreses to verify whether the RNA had degraded or not; 2) the RNA purity was tested using a NanoDrop spectrophotometer (ThermoFisher Scientific, Wilmington, DE, USA); 3) the RNA was quantified using Qubit; and 4) the RNA integrity was examined using the Agilent Bioanalyzer 2100 system (Beijing, China). The concentration of total RNA extracted was over 300 μg/μL, with a mass yield of at least 5 μg. The OD260/280 value for each sample was 2.0–2.2, and the OD260/230 was 1.8–2.1. For PacBio Sequel sequencing, two rounds of sample pooling were carried out. Firstly, we combined the leaf samples of G. luofuense at three developmental stages for cDNA library construction. Secondly, leave samples were pooled for the Illumina RNA-Seq, and the methods of RNA extraction, sequencing, and assembly referred to our previous study [13].

4.2. Sample Preparation and PacBio Sequel Sequencing

PacBio Sequel sequencing was performed with the Clontech SMARTer PCR cDNA Synthesis Kit and BluePippin Size-Selection System (Sage Science, Beverly, MA, USA). One pooled sample was subjected to purification and size selection according to the PacBio Sequel protocol: (1) The cDNA was synthesized using a Clontech SMARTer PCR cDNA Synthesis kit (Clontech, Takara Bio Inc., Shiga, Japan) and amplified using a KAPA HIFI PCR kit (Kapa Biosystems, Boston, MA, USA); (2) the cDNA was purified using a QIAquick PCR Purification kit (Qiagen, Hilden, Germany) and then precipitated and normalized with a Trimmer-2 cDNA Normalization kit (Evrogen, Moscow, Russia). Full-length cDNA damage/terminal repair and SMRTbell template preparation were then carried out. The size-selection protocol was applied because the smaller cDNAs were abundant and should be, therefore, preferentially sequenced. The RNA-seq datasets generated are available from the NCBI Sequence Read Archive database (SRA) under BioProject number accession PRJNA572572.

4.3. PacBio Sequencing Data Processing

The sequencing data were processed using PacBio SMRTlink (5.1) software. Firstly, circular consensus sequences (CCSs) were generated from the subreads.bam file adapters (effective insert of single molecules), with the following parameter settings: min_length, 200; max_drop_fraction, 0.8; no_polish, true; min_zscore, −9999; min_passes, 1; min_predicted_accuracy, 0.8; and max_length, 18,000. All the CCSs were classified as full-length reads and non-full-length reads with the following two parameters: ignoring poly(A), false and minSeqLength, 200. The full-length reads were identified considering the presence of 5′ adapter sequences, 3′ adapter sequences, and poly(A) tails. The non-full-length and full-length isoforms were clustered to generate the consensus using an isoform-level clustering algorithm. To obtain full-length polished consensus sequences, Quiver software was used to correct PacBio sequencing errors with the following parameters: hq_quiver_min_accuracy, 0.99; bin_by_primer, false; bin_size_kb, 1; qv_trim_5p, 100; and qv_trim_3p, 30. To reduce the error rates in transcript prediction and ensure transcriptome completeness, a homology search against the Pfam database was conducted. Additional nucleotide errors in consensus reads were corrected using the Illumina RNA-Seq data with the software LoRDEC (Helsinki, Finland) [54]. Additionally, consensus sequences were mapped to the reference genome of G. montanum (= G. luofuense) [3] using GMAP [55] with the following parameter settings: -f samse -n 0; min-trimmed-coverage, 0.85, and min-identity, 0.9.

4.4. Structure Analysis and Annotation

Gene structure analysis was performed using the TAPIS pipeline, the GMAP output files in bam and gff/gtf formats were used for gene and transcript determination, read clusters that overlapped non-annotated genes were classified as novel genes, and AS and alternative polyadenylation (APA) events were analyzed. Additionally, fusion transcripts were determined as the transcripts mapping to two or more long-distance-range genes and were validated by at least two Illumina reads. The ANGEL pipeline, a long-read implementation of ANGLE [56], was used to determine the protein-coding sequences from the cDNAs. The transcription factors (TFs) were predicted using software iTAK [24].
LncRNAs, which do not encode proteins, are a type of RNA with a length above 200 nt. In this study, lncRNAs were detected using the following databases: Coding-Non-Coding-Index (CNCI) [57], Coding Potential Calculator (CPC) [58], PLEK [59], and Pfam. CNCI, which was set to default parameters in this study, profiles adjoining nucleotide triplets to distinguish protein-coding from non-coding sequences. CPC was used to assess the integrity of the open reading frame in a transcript, whereas Pfam Scan was used to identify known protein family domains with the parameters (-E 0.001 and -domE 0.001). Any transcript with a Pfam hit was excluded. Default parameters were used for the Pfam searches. The NCBI’s eukaryotic protein database was used to clarify the coding and non-coding transcripts with e-value < 10−10.
Gene Ontology (GO) enrichment analysis of the full-length transcriptome was conducted using the GOseq package [60] implemented in R (R core team 2018). All full-length transcripts were submitted to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/) to obtain KEGG Orthology assignments. In addition, assembled transcripts were searched against the NCBI’s non-redundant protein sequences (NR), nucleotide sequence (NT) databases, the Swiss-Prot (a manually annotated and reviewed protein sequence database) database with e-value < 10−5, and the Pfam (database of protein families and domains) database with e-value < 10−10. KOBAS software [61] was used to test the statistical enrichment of the transcriptome in the KEGG pathways. Picard tools and SAM tools [62] were, respectively, used to mark and sort duplicated reads and to re-order the alignment results for each sample and to detect single nucleotide polymorphism.

4.5. Phylogenetic Analysis

A total of 225 bHLH protein sequences in Arabidopsis thaliana were downloaded from PlantTFDB (http://planttfdb.cbi.pku.edu.cn/) and were searched for homology against the reference genome of G. montanum (= G. luofuense) [3] and the full-length transcriptome of G. luofuense generated in this study. The domain architectures of putative bHLH transcription factors were annotated using the Pfam HMMs: bHLH (PF14215.6) and HLH (PF00010.26). All bHLH proteins were aligned using ClustalW, and neighbor-joining trees were constructed with 1000 bootstrap iteration using software MEGA (Hachioji, Japan) [63]. The motif was predicted using MEME version 5.0.5 (Reno, NV, USA), where the number of motifs was set as 6 and the motif site distribution was set as 0 or 1 occurrence per sequence.

5. Conclusions

The present study generated 80,496 full-length transcripts with a long N50 length to reveal the complexity of full-length transcriptome during leaf development of G. luofuense. To overcome the defects of PacBio Sequel sequencing, Illumina sequenced data of the same sample were applied to improve the quality of the consensus reads. It is noteworthy that a large number of novel genes and novel isoforms were detected in the present study, improving our understanding of the reference genome of G. luofuense. The full-length transcriptome could serve as a reference to further detect differentially expressed genes and isoforms and to seek for genes of interest in order to uncover gene functions during leaf development of G. luofuense. Moreover, the AS events and lncRNAs detected in the leaf transcriptome also provide additional resources for the study of the low photosynthetic characteristic of this genus.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/24/6350/s1.

Author Contributions

Investigation, N.D. and F.M.; methodology, N.D. and C.H.; writing—original draft, N.D. and C.H.; writing—review and editing, N.D. and C.H.; resources, C.H.; software, F.M.; funding acquisition, Y.T.; supervision, Y.T.; validation, C.L.; project administration, C.L.

Funding

This research was funded by the Special Fund for Hunan Forestry Science and Technology (XLK201812).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Deng, N.; Hou, C.; Liu, C.; Li, M.; Bartish, I.; Tian, Y.; Chen, W.; Du, C.; Jiang, Z.; Shi, S. Significance of photosynthetic characters in the evolution of Asian Gnetum (Gnetales). Front. Plant Sci. 2019, 10, 30. [Google Scholar] [CrossRef]
  2. Doyle, J.A. Molecular and fossil evidence on the origin of angiosperms. Rev. Earth Planet. Sci. 2012, 40, 301–326. [Google Scholar] [CrossRef]
  3. Wan, T.; Liu, Z.M.; Li, L.F.; Leitch, A.R.; Leitch, I.J.; Lohaus, R.; Liu, Z.J.; Xin, H.P.; Gong, Y.B.; Liu, Y. A genome for gnetophytes and early evolution of seed plants. Nat. Plants 2018, 4, 82. [Google Scholar] [CrossRef] [PubMed]
  4. Ickert-Bond, S.M.; Renner, S.S. The Gnetales: Recent insights on their morphology, reproductive biology, chromosome numbers, biogeography, and divergence times. J. Syst. Evol. 2016, 54, 1–16. [Google Scholar] [CrossRef] [Green Version]
  5. Gong, Y.B.; Yang, M.; Vamosi, J.C.; Yang, H.M.; Mu, W.X.; Li, J.K.; Wan, T. Wind or insect pollination? Ambophily in a subtropical gymnosperm Gnetum parvifolium (Gnetales). Plant Species Biol. 2016, 31, 272–279. [Google Scholar] [CrossRef]
  6. Sarah, M. Phylogenetic relationships among seed plants: Persistent questions and the limits of molecular data. Am. J. Bot. 2009, 96, 228–236. [Google Scholar]
  7. Price, R.A. Systematics of the Gnetales: A review of morphological and molecular evidence. Int. J. Plant Sci. 1996, 157, S40–S49. [Google Scholar] [CrossRef]
  8. Crane, P.R. Phylogenetic analysis of seed plants and the origin of angiosperms. Ann. Mo. Bot. Gard. 1985, 72, 716–793. [Google Scholar] [CrossRef]
  9. Bojian, Z.; Takahiro, Y.; Yang, Z.; Masami, H. The position of gnetales among seed plants: Overcoming pitfalls of chloroplast phylogenomics. Mol. Biol. Evol. 2010, 27, 2855–2863. [Google Scholar]
  10. Chen, H.; Wikström, N.; Strijk, J.S.; Rydin, C. Resolving phylogenetic relationships and species delimitations in closely related gymnosperms using high-throughput NGS, Sanger sequencing and morphology. Plant Syst. Evol. 2016, 302, 1345–1365. [Google Scholar]
  11. Hou, C.; Humphreys, A.M.; Thureborn, O.; Rydin, C. New insights into the evolutionary history of Gnetum (Gnetales). Taxon 2015, 64, 239–253. [Google Scholar] [CrossRef]
  12. Deng, N.; Liu, C.; Chang, E.; Ji, J.; Yao, X.; Yue, J.; Bartish, I.; Chen, L.; Jiang, Z.; Shi, S. High temperature and UV-C treatments affect stilbenoid accumulation and related gene expression levels in Gnetum parvifolium. Electron. J. Biotechnol. 2016, 25, 43–49. [Google Scholar] [CrossRef]
  13. Deng, N.; Chang, E.; Li, M.; Ji, J.; Yao, X.; Banish, I.V.; Liu, J.; Ma, J.; Chen, L.; Jiang, Z.; et al. Transcriptome characterization of Gnetum parvifolium reveals candidate genes involved in important secondary metabolic pathways of flavonoids and stilbenoids. Front. Plant Sci. 2016, 7, 174. [Google Scholar] [CrossRef] [PubMed]
  14. Tanaka, T.; Iliya, I.; Ito, T.; Furusawa, M.; Nakaya, K.I.; Iinuma, M.; Shirataki, Y.; Matsuura, N.; Ubukata, M.; Murata, J. Stilbenoids in lianas of Gnetum parvifolium. Chem. Pharm. Bull. 2001, 39, 858–862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Markgraf, F. Monographie der Gattung Gnetum. Bull. Jar. Bot. Buitenz. 1930, 10, 407–511. [Google Scholar]
  16. Feild, T.S.; Lawong, B. Xylem hydraulic and photosynthetic function of Gnetum (Gnetales) species from Papua New Guinea. New Phytol. 2008, 177, 665–675. [Google Scholar] [CrossRef]
  17. Rudall, P.J.; Rice, C.L. Epidermal patterning and stomatal development in Gnetales. Ann. Bot. 2019, 124, 149–164. [Google Scholar] [CrossRef]
  18. Yawen, L.; Xu, L.; Kunwu, L.; Hongtao, L.; Chentao, L. Multiple bHLH proteins form heterodimers to mediate CRY2-dependent regulation of flowering-time in Arabidopsis. PLoS Genet. 2013, 9, e1003861. [Google Scholar]
  19. Lynn Jo, P.; Sloan, D.B.; Bogenschutz, N.L.; Torii, K.U. Termination of asymmetric cell division and differentiation of stomata. Nature 2007, 445, 501–505. [Google Scholar]
  20. Kyoko, O.I.; Bergmann, D.C. Arabidopsis FAMA controls the final proliferation/differentiation switch during stomatal development. Plant Cell 2006, 18, 2493–2505. [Google Scholar]
  21. Gabriela, T.O.; Enamul, H.; Quail, P.H. The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell 2003, 15, 1749–1770. [Google Scholar]
  22. Thatcher, S.R.; Zhou, W.; Leonard, A.; Wang, B.B.; Beatty, M.; Zastrow-Hayes, G.; Zhao, X.; Baumgarten, A.; Li, B. Genome-wide analysis of alternative splicing in Zea mays: Landscape and genetic regulation. Plant Cell 2014, 26, 3472. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wahl, M.C.; Will, C.L.; Reinhard, L. The spliceosome: Design principles of a dynamic RNP machine. Cell 2009, 136, 701–718. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhang, Q.; Zhang, X.; Wang, S.; Tan, C.; Zhou, G.; Li, C. Involvement of alternative splicing in barley seed germination. PLoS ONE 2016, 11, e0152824. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Yang, S.; Tang, F.; Zhu, H. Alternative splicing in plant immunity. Int. J. Mol. Sci. 2014, 15, 10424–10445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Estelle, R.; Cabrito, T.R.; Batista, R.A.; Hussein, M.A.M.; Teixeira, M.C.; Alekos, A.; Isabel, S.C.; Paula, D. Intron retention in the 5′UTR of the novel ZIF2 transporter enhances translation to promote zinc tolerance in Arabidopsis. PLoS Genet. 2014, 10, e1004375. [Google Scholar]
  27. Blencowe, B.J. Alternative splicing: New insights from global analyses. Cell 2006, 126, 37–47. [Google Scholar] [CrossRef] [Green Version]
  28. Jun, L.; Huan, W.; Nam-Hai, C. Long noncoding RNA transcriptome of plants. Plant. Biotechnol. J. 2015, 13, 319–328. [Google Scholar]
  29. Ariel, F.; Romero-Barrios, N.; Jégu, T.; Benhamed, M.; Crespi, M. Battles and hijacks: Noncoding transcription in plants. Trends Plant Sci. 2015, 20, 362–371. [Google Scholar] [CrossRef]
  30. Bernard, E.; Jacob, L.; Mairal, J.; Vert, J.P. Efficient RNA isoform identification and quantification from RNA-Seq data with network flows. Bioinformatics 2013, 30, 2447–2455. [Google Scholar] [CrossRef] [Green Version]
  31. Tilgner, H.; Raha, D.; Habegger, L.; Mohiuddin, M.; Gerstein, M.; Snyder, M. Accurate identification and analysis of human mRNA isoforms using deep long read sequencing. G3-Genes Genom. Genet. 2013, 3, 387–397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Donald, S.; Hagen, T.; Fabian, G.; Michael, S. A single-molecule long-read survey of the human transcriptome. Nat. Biotechnol. 2013, 31, 1009–1014. [Google Scholar]
  33. Cheng, B.; Furtado, A.; Henry, R.J. Long-read sequencing of the coffee bean transcriptome reveals the diversity of full-length transcripts. Gigascience 2017, 6, 1–13. [Google Scholar] [CrossRef] [PubMed]
  34. Abdelghany, S.E.; Hamilton, M.; Jacobi, J.L.; Ngam, P.; Devitt, N.; Schilkey, F.; Benhur, A.; Reddy, A.S.N. A survey of the sorghum transcriptome using single-molecule long reads. Nat. Biotechnol. 2016, 7, 11706. [Google Scholar]
  35. Hoang, N.V.; Furtado, A.; Mason, P.J.; Marquardt, A.; Kasirajan, L.; Thirugnanasambandam, P.P.; Botha, F.C.; Henry, R.J. A survey of the complex transcriptome from the highly polyploid sugarcane genome using full-length isoform sequencing and de novo assembly from short read sequencing. BMC Genom. 2017, 18, 395. [Google Scholar] [CrossRef]
  36. Wang, M.; Wang, P.; Liang, F.; Ye, Z.; Li, J.; Shen, C.; Pei, L.; Wang, F.; Hu, J.; Tu, L. A global survey of alternative splicing in allopolyploid cotton: Landscape, complexity and regulation. New Phytol. 2017, 217, 163. [Google Scholar] [CrossRef] [Green Version]
  37. Makita, Y.; Kawashima, M.; Lau, N.S.; Othman, A.S.; Matsui, M. Construction of Pará rubber tree genome and multi-transcriptome database accelerates rubber researches. BMC Genom. 2018, 19, 922. [Google Scholar] [CrossRef] [Green Version]
  38. Kin Fai, A.; Underwood, J.G.; Lawrence, L.; Wing Hung, W. Improving PacBio long read accuracy by short read alignment. PLoS ONE 2012, 7, e46679. [Google Scholar]
  39. Roberts, R.J.; Carneiro, M.O.; Schatz, M.C. The advantages of SMRT sequencing. Genome Biol. 2013, 14, 405. [Google Scholar] [CrossRef]
  40. Rinn, J.L.; Chang, H.Y. Genome regulation by long noncoding RNAs. Annu. Rev. Biochem. 2012, 81, 145–166. [Google Scholar] [CrossRef] [Green Version]
  41. Jun, L.; Choonkyun, J.; Jun, X.; Huan, W.; Shulin, D.; Lucia, B.; Catalina, A.H.; Nam-Hai, C. Genome-wide analysis uncovers regulation of long intergenic noncoding RNAs in Arabidopsis. Plant Cell 2012, 24, 4333–4345. [Google Scholar]
  42. He, F.; Liu, Q.; Zheng, L.; Cui, Y.; Shen, Z.; Zheng, L. RNA-Seq analysis of rice roots rveals the involvement of post-transcriptional regulation in response to cadmium stress. Front. Plant Sci. 2015, 6, 1136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Song, X.; Sun, L.; Luo, H.; Ma, Q.; Zhao, Y.; Pei, D. Genome-wide identification and characterization of long non-coding RNAs from mulberry (Morus notabilis) RNA-seq Data. Genes 2016, 7, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kwenda, S.; Birch, P.R.J.; Moleleki, L.N. Genome-wide identification of potato long intergenic noncoding RNAs responsive to Pectobacterium carotovorum subspecies brasiliense infection. BMC Genom. 2016, 17, 614. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Yuqiu, W.; Xiuduo, F.; Fang, L.; Guangming, H.; William, T.; Danmeng, Z.; Xing Wang, D. Arabidopsis noncoding RNA mediates control of photomorphogenesis by red light. Proc. Natl. Acad. Sci. USA 2014, 111, 10359. [Google Scholar]
  46. Amor, B.B.; Wirth, S.; Merchan, F.; Laporte, P.; D’Aubentoncarafa, Y.; Hirsch, J.; Maizel, A.; Mallory, A.; Lucas, A.; Deragon, J.M. Novel long non-protein coding RNAs involved in Arabidopsis differentiation and stress responses. Genome Res. 2009, 19, 57. [Google Scholar] [CrossRef] [Green Version]
  47. Liu, X.; Li, D.; Zhang, D.; Yin, D.; Zhao, Y.; Ji, C.; Zhao, X.; Li, X.; He, Q.; Chen, R. A novel antisense long noncoding RNA, twisted leaf, maintains leaf blade flattening by regulating its associated sense R2R3-MYB gene in rice. New Phytol. 2018, 218, 774–788. [Google Scholar] [CrossRef] [Green Version]
  48. Yandell, M.; Ence, D. A beginner’s guide to eukaryotic genome annotation. Nat. Rev. Genet. 2012, 13, 329–342. [Google Scholar] [CrossRef]
  49. Li, Y.; Dai, C.; Hu, C.; Liu, Z.; Kang, C. Global identification of alternative splicing via comparative analysis of SMRT- and Illumina-based RNA-seq in strawberry. Plant J. 2017, 90, 164. [Google Scholar] [CrossRef] [Green Version]
  50. Grandori, C.; Cowley, S.M.; James, L.P.; Eisenman, R.N. The Myc/Max/Mad network and the transcriptional control of cell behavior. Annu. Rev. Cell Dev. Biol. 2000, 16, 653. [Google Scholar] [CrossRef]
  51. Atchley, W.R.; Fitch, W.M. A natural classification of the basic helix-loop-helix class of transcription factors. Proc. Natl. Acad. Sci. USA 1997, 94, 5172–5176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Morgenstern, B.; Atchley, W.R. Evolution of bHLH transcription factors: Modular evolution by domain shuffling? Mol. Biol. Evol. 1999, 16, 1654–1663. [Google Scholar] [CrossRef] [PubMed]
  53. Furtado, A. RNA extraction from developing or mature wheat seeds. Methods Mol. Biol. 2014, 1099, 23. [Google Scholar] [PubMed]
  54. Leena, S.; Eric, R. LoRDEC: Accurate and efficient long read error correction. Bioinformatics 2014, 30, 3506–3514. [Google Scholar]
  55. Wu, T.; Watanabe, C. GMAP: A genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 2005, 21, 1859. [Google Scholar] [CrossRef] [Green Version]
  56. Shimizu, K.; Adachi, J.; Muraoka, Y. ANGLE: A sequencing errors resistant program for predicting protein coding regions in unfinished cDNA. J. Bioinform. Comput. Biol. 2011, 4, 649–664. [Google Scholar] [CrossRef]
  57. Liang, S.; Haitao, L.; Dechao, B.; Guoguang, Z.; Kuntao, Y.; Changhai, Z.; Yuanning, L.; Runsheng, C.; Yi, Z. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013, 41, e166. [Google Scholar]
  58. Lei, K.; Yong, Z.; Zhi-Qiang, Y.; Xiao-Qiao, L.; Shu-Qi, Z.; Liping, W.; Ge, G. CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007, 35, W345–W349. [Google Scholar]
  59. Li, A.; Zhang, J.; Zhou, Z. PLEK: A tool for predicting long non-coding RNAs and messenger RNAs based on an improved k-mer scheme. BMC Bioinform. 2014, 15, 311. [Google Scholar] [CrossRef] [Green Version]
  60. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11. [Google Scholar] [CrossRef] [Green Version]
  61. Chen, X.; Xizeng, M.; Jiaju, H.; Yang, D.; Jianmin, W.; Shan, D.; Lei, K.; Ge, G.; Chuan-Yun, L.; Liping, W. KOBAS 2.0: A web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011, 39, 316–322. [Google Scholar]
  62. Heng, L. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 2011, 27, 2987–2993. [Google Scholar]
  63. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
Ijms 20 06350 g001 550
Figure 1. Results of mapping to reference genome. (A) Pie chart of mapped reads; (B) Range of mapping coverage and identity, x-axis represents the scale ranges, y-axis represents the percentages; (C) Saturation curve of consensus reads, x-axis represents numbers of full-length non-chimera (flnc) reads, y-axis represents numbers of genes; (D) Classification of flnc reads after mapping to the reference genome of G. montanum (= G. luofuense).
Figure 1. Results of mapping to reference genome. (A) Pie chart of mapped reads; (B) Range of mapping coverage and identity, x-axis represents the scale ranges, y-axis represents the percentages; (C) Saturation curve of consensus reads, x-axis represents numbers of full-length non-chimera (flnc) reads, y-axis represents numbers of genes; (D) Classification of flnc reads after mapping to the reference genome of G. montanum (= G. luofuense).
Ijms 20 06350 g001
Ijms 20 06350 g002 550
Figure 2. Nucleotide sequence (NR) (NCBI nonredundant protein sequences) homologous species distribution diagram of all genes and novel genes. Bars represent numbers of annotated consensus isoforms homologous to varied land plant species.
Figure 2. Nucleotide sequence (NR) (NCBI nonredundant protein sequences) homologous species distribution diagram of all genes and novel genes. Bars represent numbers of annotated consensus isoforms homologous to varied land plant species.
Ijms 20 06350 g002
Ijms 20 06350 g003 550
Figure 3. Functional annotation of full-length transcripts based on gene ontology (GO) categories. Bars represent the numbers of assignments proteins with BLASTx matches to each GO term. Upper bars represent GO annotation of all genes, and lower bars represent GO annotation of novel genes.
Figure 3. Functional annotation of full-length transcripts based on gene ontology (GO) categories. Bars represent the numbers of assignments proteins with BLASTx matches to each GO term. Upper bars represent GO annotation of all genes, and lower bars represent GO annotation of novel genes.
Ijms 20 06350 g003
Ijms 20 06350 g004 550
Figure 4. Functional annotation of full-length transcripts based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) categories. Bars represent the numbers of assignments proteins with BLASTx matches to each KEGG term.
Figure 4. Functional annotation of full-length transcripts based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) categories. Bars represent the numbers of assignments proteins with BLASTx matches to each KEGG term.
Ijms 20 06350 g004
Ijms 20 06350 g005 550
Figure 5. Analysis of long non-coding RNAs (lncRNAs) in G. luofuense leaves. (A) A Venn diagram of lncRNAs predicted by 4 methods; (B) The length density distribution of predicted lncRNAs and mRNAs in G. luofuense leaves, x-axis represents the length, y-axis represents the density; (C) Classification of predicted lncRNAs; (D) The distribution of exons numbers of mRNAs and predicted lncRNAs, x-axis represent numbers of exons, y-axis represent the percentages.
Figure 5. Analysis of long non-coding RNAs (lncRNAs) in G. luofuense leaves. (A) A Venn diagram of lncRNAs predicted by 4 methods; (B) The length density distribution of predicted lncRNAs and mRNAs in G. luofuense leaves, x-axis represents the length, y-axis represents the density; (C) Classification of predicted lncRNAs; (D) The distribution of exons numbers of mRNAs and predicted lncRNAs, x-axis represent numbers of exons, y-axis represent the percentages.
Ijms 20 06350 g005
Ijms 20 06350 g006 550
Figure 6. (A) Number of transcript factors identified in the present study; (B) Varied types of alternative polyadenylation and corresponding gene numbers; (C) A schematic graph that illustrates different types of alternative splicing.
Figure 6. (A) Number of transcript factors identified in the present study; (B) Varied types of alternative polyadenylation and corresponding gene numbers; (C) A schematic graph that illustrates different types of alternative splicing.
Ijms 20 06350 g006
Ijms 20 06350 g007 550
Figure 7. Gene structure, motif, and phylogenetic analysis of 30 bHLH genes. (A) A phylogenetic tree of bHLH transcription factors using a neighbor-joining method; (B) Conserved motifs in bHLH genes were marked in different colors; (C) Gene structures of bHLH genes in Gnetum.
Figure 7. Gene structure, motif, and phylogenetic analysis of 30 bHLH genes. (A) A phylogenetic tree of bHLH transcription factors using a neighbor-joining method; (B) Conserved motifs in bHLH genes were marked in different colors; (C) Gene structures of bHLH genes in Gnetum.
Ijms 20 06350 g007
Ijms 20 06350 g008 550
Figure 8. Different transcript isoforms of one bHLH gene detected in the leaf transcriptome mapped to the reference genome of G. luofuense.
Figure 8. Different transcript isoforms of one bHLH gene detected in the leaf transcriptome mapped to the reference genome of G. luofuense.
Ijms 20 06350 g008
Table
Table 1. Summary of PacBio sequencing data in Gnetum. luofuense.
Table 1. Summary of PacBio sequencing data in Gnetum. luofuense.
TermsAmount
Subreads bases9.98 G
Number of subreads3,689,825
Average length of subreads2750 bp
N50 of subreads3168 bp
Number of CCSs185,089
Number of sequences with 5′ terminal primers167,590
Number of sequences with 3′ terminal primers166,862
Number of sequences with poly(A) tails156,554
Number of full-length sequences143,578
Number of full-length non-chimeric (flnc) reads139,488
Average length of flnc reads3065 bp
Percentage of flnc reads75%
Number of polished consensus reads80,496
Minimum length of consensus reads167 bp
Maximum length of consensus reads14,735 bp
Average length of consensus reads3223 bp
N50 of consensus reads3614 bp
Table
Table 2. Summary of data correction of PacBio sequencing data in comparison with Illumina sequenced short reads.
Table 2. Summary of data correction of PacBio sequencing data in comparison with Illumina sequenced short reads.
TypeBefore Correction
(Pacbio Sequel)		After Correction
(Pacbio Sequel)		Short Reads
(Illumina)
Total nucleotide259,381,401260,514,86744,269,498
Total sequence80,49680,49645,566
Mean length3223 bp3237 bp972 bp
Minimum length167 bp167 bp201 bp
Maximum length14,735 bp14,734 bp12,325 bp
N503614 bp3629 bp2030 bp
N902102 bp2102 bp323 bp
Table
Table 3. Comparisons of Pacbio and Illumina sequenced data with regard to read mapping.
Table 3. Comparisons of Pacbio and Illumina sequenced data with regard to read mapping.
TermsPacBio Sequenced DataIllumina Sequenced Data
Number of ReadsPercentageNumber of ReadsPercentage
Total mapped77,38096.13%54,289,03893.01%
Unmapped31163.87%4,079,9956.99%
Multiple mapped66788.30%1,691,8372.90%
Uniquely mapped70,70287.83%52,597,20190.11%
Uniquely mappedto positive strands43,29953.79%26,287,94945.04%
Uniquely mappedto negative strands27,40334.04%26,309,25245.07%
Table
Table 4. Numbers of annotated consensus isoforms by the search against different databases.
Table 4. Numbers of annotated consensus isoforms by the search against different databases.
DatabaseTotal Number of Annotated IsoformsNumber of Novel GenesNumber of Novel Isoforms of Known GenesNumber of Isoforms of Known Genes
NR34,170378222,1678221
SwissProt28,954274119,0237190
KEGG33,813365922,0058149
KOG21,723192514,6585140
GO22,297104715,1806070
NT15,46793610,6803851
Pfam22,297104715,1806070
In all databases898011665488281
At least in one database34,667526923,4432316

Share and Cite

MDPI and ACS Style

Deng, N.; Hou, C.; Ma, F.; Liu, C.; Tian, Y. Single-Molecule Long-Read Sequencing Reveals the Diversity of Full-Length Transcripts in Leaves of Gnetum (Gnetales). Int. J. Mol. Sci. 2019, 20, 6350. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20246350

AMA Style

Deng N, Hou C, Ma F, Liu C, Tian Y. Single-Molecule Long-Read Sequencing Reveals the Diversity of Full-Length Transcripts in Leaves of Gnetum (Gnetales). International Journal of Molecular Sciences. 2019; 20(24):6350. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20246350

Chicago/Turabian Style

Deng, Nan, Chen Hou, Fengfeng Ma, Caixia Liu, and Yuxin Tian. 2019. "Single-Molecule Long-Read Sequencing Reveals the Diversity of Full-Length Transcripts in Leaves of Gnetum (Gnetales)" International Journal of Molecular Sciences 20, no. 24: 6350. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20246350

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop