PacBio Single-Molecule Long-Read Sequencing Reveals Genes Tolerating Manganese Stress in Schima superba Saplings

Abstract

Schima superba (Theaceae) is a subtropical evergreen tree and is used widely for forest firebreaks and gardening. It is a plant that tolerates salt and typically accumulates elevated amounts of manganese in the leaves. With large ecological amplitude, this tree species grows quickly. Due to its substantial biomass, it has a great potential for soil remediation. To evaluate the thorough framework of the mRNA, we employed PacBio sequencing technology for the first time to generate S. Superba transcriptome. In this analysis, overall, 511,759 full length non-chimeric reads were acquired, and 163,834 high-quality full-length reads were obtained. Overall, 93,362 open reading frames were obtained, of which 78,255 were complete. In gene annotation analyses, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Genes (COG), Gene Ontology (GO), and Non-Redundant (Nr) databases were allocated 91,082, 71,839, 38,914, and 38,376 transcripts, respectively. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT), and Coding Non-Coding Index (CNCI) databases and observed 8,551, 9,174, 20,720, and 18,669 lncRNAs, respectively. Moreover, nine genes were randomly selected for the expression analysis, which showed the highest expression of Gene 6 (Na_Ca_ex gene), and CAX (CAX-interacting protein 4) was higher in manganese (Mn)-treated group. This work provided significant number of full-length transcripts and refined the annotation of the reference genome, which will ease advanced genetic analyses of S. superba.

Introduction

Schima superba is naturally grown in mining areas and abundantly absorbs different metals, particularly manganese (Mn). Due to this absorbance ability, this tree is being used as phytoremediation agent in metal toxic soils (Yang et al., 2008). This tree belongs to the family Theaceae, and it is commonly used for gardening and forest firebreak. Characteristically, the leaves of this tree contain an unusual high concentration of Mn (Yang et al., 2021).

Hyper-accumulation of harmful metals in different plants helps in cleaning or decontaminating soils. Contamination of metals in the soil can be reduced through plants’ absorbance and transportation, which is a kind of green environmental remediation technology (Padmavathiamma and Li, 2007). Woody plants have high biomass and long-life cycle, and those are unique advantages in the phytoremediation of heavy metal. But the competency of woody trees to accumulate heavy metals is usually very low, and most of heavy metals are accumulated in the roots of woody trees (Capuana, 2020). As reported by Baker and Brooks (1989), more than 10,000 mg/kg concentrations of Mn in above ground tissues of plants are considered as the hyper-accumulation of Mn. Along with having a potential to accumulate Mn, S. superba has vast ecological amplitude, dominant canopy species (Tang et al., 2020) and has considerable biomass, which suggests its potential use in the soil remediation (Luo et al., 2006).

Although, researchers have reported the hyper-accumulation of Mn in leaves, roots, and shoots, no study has performed the transcriptomic analysis of S. superba till date. Since the establishment of large-scale sequencing technology, transcript sequences have become a significant source of studying gene expression and regulation (Levy and Myers, 2016). The restriction of sequencing read length in the second-generation, the full-length record acquired by splicing is not finished, and the third-generation sequencing technology signified by Pacific Biosciences (PacBio) has successfully overcome this issue (Jiao and Schneeberger, 2017; Slatko et al., 2018). Based on its long-read length, high-quality full-length record can be legitimately acquired to distinguish simple sequence repeat (SSR), fulfill functional annotation of transcripts, and lncRNA (Wang et al., 2015; Xie et al., 2020). Transcriptomics provides information of the complete set of RNA transcripts in a specific tissue or cell under different developmental or physiological stages, involving ribosomal RNA, transfer RNA, messenger RNA, and non-coding RNAs. Moreover, it concentrates on the gene expression, gene structure, and function at RNA level that help to reveal the molecular processes involved in various biological mechanisms (Dong and Chen, 2013).

In this research, a single molecule real-time sequencing (SMRT) has been completed to get the full-length transcriptome of S. superba. Utilizing the transcriptome data functional annotations of transcripts, transcriptional factors (TFs) analysis, prediction of lnRNA, and basic sequence repeat (SSR) analysis were performed. In addition, the expression pattern of various genes provides basic understanding of transcriptional regulation of S. superba. This study will not only assist future functional genomics research but also provide an avenue for further selective genetic engineering and breeding projects of S. superba.

Materials and Methods

Pot Experiment and Sample Collection

A pot test was carried out in the glasshouse of Shanghai Jiao Tong University, China, by using 1-year-old S. superba saplings. For this purpose, 18 saplings, with same growth potential (about 35 cm high), were taken and randomly allocated into two groups: control (CK) and Mn concentration 100 mM (WT), in three replicates. Plant tissues were collected after 1, 5, and 10 days of transplantation in quarter-strength of Hoagland solution. The tissues were collected in liquid nitrogen and processed through the isolation of RNA at −80°C.

RNA Extraction

RNA extraction from the collected plant samples (0.2 g leaves) was performed with the RNeasy plus Mini Kit (Qiagen, Valencia, United States). Using the Nano drop ND-1000 (Nano drop Technologies, Rockland, DE, United States) spectrophotometer, its quality and quantity were calculated after observing RNA on agarose gel. By using Qubit^® RNA Assay Kit and RNA Nano 6000 Assay Kit, RNA quality has been identified clearly.

Library Construction and SMRT Sequencing

For sequencing, C2 sequencing reagents were used in Pacific Biosciences (PacBio) real-time sequencer. Purified RNA was utilized for synthesizing cDNA by SMRT PCR cDNA Synthesis Kit (Clontech, CA, United States). Full-length cDNA of different sizes was selected and cDNA libraries were constructed using the BluePippin^® (SageScience, Beverly, MA, United States). After BluePippin screening, the fragments were pacified to large-scale PCR to get sufficient total cDNA and evaluate using Qubit fluorometer (Life Technologies, Carlsbad, CA, United States). The libraries’ uniqueness was maintained by using Agilent Bio analyzer 2100 system and SMRT sequencing was achieved.

Error Correction and Quality Filtering

Sequence statistics were obtained by using the SMRT link 5.1 software. Relevant parameters such as no-polish TRUE, min-length 50, max-drop-fraction 0.8, min-predicted-accuracy 0.8, min-passes 2, min-z score -9999.0, and max-length 15000 were followed to produce CCS by subreading BAM files. The output was CCS.BAM files, which, using pbclassify.py, were categorized into full length and non-full length reads. Full and non-full length FASTA files were fed into the cluster stage, and additional nucleotide errors were corrected by LoRDEC software in the consensus readings. Unnecessary data in the corrected consensus reads were separated using CD-HIT (-c 0.95 -T 6 -G 0 -aL 0.00 -aS 0.99) to acquire final transcripts for the subsequent investigation.

Functional Annotation of Transcripts

We recognized functional annotations matching by searching to the non-redundant nucleotide database (Nr), Swiss-Prot, Pfam, non-supervised orthologous groups egg (NOG), clusters of orthologous groups (COG), eukaryotic orthologous groups (KOG), Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and BLAST alignment databases, and all transcript sequences were examined for homology (Bairoch and Apweiler, 2000; Li et al., 2002; Tatusov et al., 2003; Kanehisa et al., 2004).

Identification of TFs, lncRNAs, and SSR

Plant transcription factors were predicted using iTAK v1.7a (Zheng et al., 2016). Four tools, CNCIv2 (Sun et al., 2013), CPCvcpc-0.9-r2 with e-value “1e-10” (Kong et al., 2007), Pfam-scan (E 0.001 -domE 0.001) (Finn et al., 2016) and PLEKv1.2 with min length 200 (Li et al., 2002) were chosen to predict candidate lncRNAs. Transcripts predicted with coding capacity by all of the mentioned above tools were separated out, and those without coding prospective were considered as candidate set of lncRNAs.

Development of SSR Markers

Simple sequence repeats were identified by MISA v1.0¹ (Beier et al., 2017), with default parameters. MISA can recognize seven types of SSR, namely mono, di, tri, tetra, penta, and hexa nucleotide by analyzing transcript sequences. By using tool Batch primer 3, SSR primers were designed. DNA was extracted from S. superba for PCR amplification, and the product was extracted in 8% polyacrylamide gel.

qRT-PCR Analysis

All 36 saplings of control group (CK) and Mn-treated group (WT) were utilized by RNA extraction via RNeasy plus Mini Kit (Qiagen, Valencia, United States). The cDNA was synthesized using SMRT PCR cDNA synthesis kit (Clontech, United States). The qRT PCR was performed to see the expression of nine different transporter genes (Supplementary Table S5). The reaction mixture consisted of 7.8 μl of cDNA and 10 μl of 2× qPCR. The process involve initial denaturation at 95°C for 10 min, which is followed by 45 cycles of denaturation at 95°C for 15 s, while annealing and extension at 60°C for 60 1 min.

Results

SMRT Sequencing Data of S. superba

Plant tissues were used for RNA extraction and library construction of cDNA. After removing adaptor sequences, a total of 1,311,552,27 low-quality sequences were obtained. It comprised of 1–6 kb data. After self-correction of subread sequences (with minPasses = 39), overall 588,479 circular consensus sequence (CCS) were obtained. Moreover, 511,759 full-length non-chimeric (FLNC) picks were also identified. Iterative sequence clustering revealed 167,782 consensus isoforms, with a mean read length of 2187 bp. Together with non-FL sequence, the Quiver program corrected the consistent sequence in each cluster, yielding 163834 high-quality transcripts (high- quality isoforms) with more than 97.65% accuracy. In addition, the low-quality record was corrected with the corresponding sequencing data to enhance the accuracy of the isoforms. Finally, using CD HIT (3.3) tools, 97,520 bar non-redundant transcript sequences were obtained (Tables 1, 2 and Figure 1).

FIGURE 1

Open Reading Frame and AS Event Prediction

About 93,362 ORFs were predicted using Trans Decoder software. Overall, 78,255 complete ORFs were recognized, and the complete ORFs were successfully analyzed for length distribution (Figure 2); 7,489 AS events were revealed out of all transcripts obtained by SMRT sequencing (Supplementary Table S1). Owing to the lack of an unavailable S. Superba reference genome, the classification of the forms of AS events in future research would be justified (Figure 2).

FIGURE 2

Long Non-coding RNA Identification

These are a group of non-coding poly-A RNAs that are involved in growth and stress response of the plant. In this work, we utilize four computational methods to recognize lncRNAs, associated with Pfam, CPC, CPAT, and CNCI databases. Overall, 8,551, 9,174, 20,720, and 18,669 lncRNAs were recognized in the Pfam, CPC, CPAT, and CNCI, respectively (Supplementary Table S2). By filtering the record of less than 300 bp, 2,011 transcripts were examined as lncRNAs by all four mechanisms (Figure 3).

FIGURE 3

Transcription Factor Prediction

These are the key regulators of gene expression and play a significant role in plant growth and development. In this study, TFs are divided into 64 families, and 9423 putative TFs were identified (Supplementary Table S3). TFs in S. superba transcriptome were mostly related to the RLK-Pelle-DLS (473,5.02%), bHLH (278,2.95%), MYB-related (230,2.44%), CAMK-CDPK (218,2.31%), C2H2 (218,2.31%), C3H (202,2.14%), RWP-RK (189,2.01%), bZIP (177,1.88%), TKL-PI-4 (159,1.69%), and NAC (152,1.61%) families (Figure 4).

FIGURE 4

Functional Annotation of Transcripts

All 16,3834 unique SMRT transcripts were practically annotated by seven data storage, such as GO, KOG, Pfam, Swiss-Prot, COG, Nr, and KEGG using BLAST (7) software [version 2.2.26] (Table 3). By comparing the transcript sequence to NR with homologous species, 68,630 genes were annotated (Figure 5).

FIGURE 5

GO Classified Transcripts

Transcripts GO classification statistics demonstrated 71,839 unique genes, which were enriched in major categories of molecular function, cellular component, biological component, and catalytic activity (Figures 3, 5). This analysis also helped us to get transcripts COG classification statistics (Figure 6).

FIGURE 6

COG Classification

To examine the functional annotation and classification of S. Superba for further study, all the transcripts were checked against the COG database clusters². This investigation indicated that 38,914 transcripts were allocated into 24 categories (Figure 7). The highest group was general function prediction only (4,690, 11.1%), followed by signal transduction (4,291, 10.15%) whereas carbohydrate transport and metabolism was 3,981, 9.42%. Our results revealed six groups percentage were less than 1.00%, i.e., nuclear structure and modification, chromatin structure and dynamics including extracellular structure (Figure 7).

FIGURE 7

KEGG Annotated Transcripts

KEGG data storage interpreted a total of 67,426 sequences and plotted 367 operative categories in S. superba. Among them, metabolism was the largest category. The functional annotations of all 78,559 unique transcripts were detected in this work (Supplementary Table S4). A significant number of genes, especially interrelated in salt-tolerance and fatty acid component of S. superba were annotated, such as oxidative phosphorylation (1073), plant hormone signal transduction (506), fatty acid biosynthesis (246), biosynthesis of unsaturated fatty acids (94), and α-linolenic acid metabolism (199). We also recognized matches to our unique transcripts in clusters of orthologous class of proteins (COG) (44,376, 56.49%), Pfam database (41,535, 52.87%), and Swiss-Prot (58,535, 74.51%) (Table 4).

Analysis of Gene Expression in qRT PCR

For the expression analysis of potential high expression genes, nine transporter genes were selected randomly for qRT PCR. Results of qRT PCR analysis showed the highest expression by Gene 6 (Na_Ca_ex), in Mn-treated group (WT) on day 1, 5, and 10, while the expression of CAX (CAX-interacting protein 4) was higher in Mn-treated group (WT) at day 5 (Figure 8).

FIGURE 8

SSR Identification

After the screening of 95,258 obtained transcripts, about 58,396 possible SSRs were recognized from 29,075 transcripts. Among them 26,312 composed one SSR, and 19,007 contained two loci or more. In addition, 29,075 and 18,324 transcripts contained one SSR and at least two SSRs, respectively. Whereas observed compound formations were 18,324 SSRs. As shown in Table 5, the number of mono-, di-, tri-, tetra-, penta-, and hexa-nucleotide repeats were 44,045, 51,716, 14,521, 1,675, 753, and 988, accordingly. SSRs with 10 repeat units (16949, 14.9%) were the most abundant, followed by those with 6 (12563, 11%), 11 (10564, 9.2%), and 7 (7995, 87.03%) (Table 6).

Development of SSR Markers

By using Primer 3.0 software, 95,374 pair of primers were designed, and among them, 200 were randomly selected for PCR (Supplementary Table S6). The product of 180 PCR pairs was examined successfully, with an efficiency rate of 90%. However, the remaining 20 pairs of primer failed to amplify at various annealing temperatures.

Discussion

Since the development of high-throughput sequencing technology, transcriptomic analysis has become a valuable technology to study gene expression and regulation (Wang et al., 2015; Dharshini et al., 2016). However, due to read length limitation of the second-generation sequencing in different organisms, the full-length transcript obtained by splicing is not complete (Rhoads and Au, 2015). SMRT sequencing technology effectively solves this problem. This study is the first transcriptome analysis of S. superba using SMRT sequencing technology. The sequencing peaks were gained from leaf samples using the Pacific Biosciences Iso-Seq platform. After sequencing, 511,759 full-length non-chimerics (FLNC) were identified. Previous studies also reported similar results that newly discovered transcripts by SMRT in red clover and wheat were 45 bp longer than the known transcripts (Dong et al., 2015; Chao et al., 2018). Subsequently, our results also showed that the SMRT technology is very efficient for full-length cDNA sequencing and provided a rich resource for further functional genomics analysis in S. superba (Workman et al., 2018; Zuo et al., 2018).

In past years, PacBio platform was used to sequence the full-length transcriptional group of tea (Camellia sinensis) (Qiao et al., 2019). Finally, 213,389 polished consensus sequences were obtained, 223,120 CDS sequences were predicted, 195,062 SSR loci were detected, and 5,785 transcription factors belonging to 60 transcription factor families were predicted. PacBio RS II platform was used to analyze the full-length transcription group of sunflowers, 10.43 Gb clean data were obtained, and 38,302 de-redundant sequences were obtained (Chen et al., 2018). Previous research on S. superba was mainly focused on genetic mapping, gene expression, and proteomic analysis (Chen et al., 2014; Zhang et al., 2019, Yang et al., 2021). However, these techniques could not provide and assembled full-length transcripts. On the other hand, PacBio SMRT single-molecule long-read sequencing revealed a better capacity in capturing transcript sequences, mainly long transcript sequences. The transcriptome sequencing of 18 samples was completed, and a total of 124.45 G of clean data were obtained, with Q30 reaching 85%. These data were corrected by Illumina clean reads, after the integration and quality control of data from the two platforms. After sequencing, 588,479 circular consensus (CCS) reads were obtained, including 511,759 full-length reads non-chimeric (FLNC) sequences. Previous studies showed that the length of the FLNC sequence emulates the length of the cDNA sequence in whole sequence library, whereas the library can be examined by the length of FLNC sequence (Chao et al., 2018). The full-length non-chimeric sequences were clustered to obtain 167,782 consensus sequences, and 163,834 high-quality consensus sequences were obtained by polishing the consensus sequences. The low-quality consensus sequences were corrected with the second-generation transcriptome data, and they were consistent with high-quality. The sequences were merged and subjected to de-redundancy analysis to obtain 41,112 transcript sequences. Each sample obtained 7,489 alternative splicing, and a total of 113,698 SSR and 78,255 complete CDS regions were obtained. Previous studies showed variations in the percentage of tetra, penta, and hexa nucleotide repeats observed in sorghum (5.4%), rice (2.54%), Populus (1.66%), and Medicago (0.94%) (Ref 4). Our findings of sequence library provide basic information for improving the S. superba draft genome annotation and complete characterization of this species transcriptome.

Results of qRT PCR analysis indicated the highest expression of Gene 6 (CAX-interacting protein 4), in Mn-treated group (WT) on day 1, 5, and 10. Our results are similar with the study of Pittman and Hirschi (2016). Previous studies indicated that CAX shows high expression under salt stress. Moreover, the higher expression of CAX (Na_Ca_ex) in Mn-treated group (WT) was also observed. The CAX proteins are one of the five transporter families that constitute Ca²⁺/cation antiporters (CaCA) superfamily (Shigaki et al., 2006; Emery et al., 2012). Although the CAX family members were at first identified as Ca²⁺ transporters, related work revealed their ability to transport a wide array of ions (Pittman and Hirschi, 2016). In plants, the CAXs mediate efflux of ions present into the vacuole. This study has provided enriched information about expression and regulation of genes in stress conditions of S. superba. In addition, this study provided for the first time a full-length transcriptome of the S. superba using the SMRT sequencing method. The transcriptome design in this research will assist future research on functional genomics and facilitate support for advance genetic engineering of S. superba.

References

• 1

  BairochA.ApweilerR. (2000). The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000.Nucleic Acids Res.2845–48. 10.1093/nar/28.1.45

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  BakerA. J.BrooksR. (1989). Terrestrial higher plants which hyperaccumulate metallic elements. a review of their distribution, ecology and phytochemistry.Biorecovery1281–126.

  • Google Scholar

• 3

  BeierS.ThielT.MünchT.ScholzU.MascherM. (2017). MISA-web: a web server for microsatellite prediction.Bioinformatics332583–2585. 10.1093/bioinformatics/btx198

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  CapuanaM. (2020). A review of the performance of woody and herbaceous ornamental plants for phytoremediation in urban areas.IForest13:139. 10.3832/ifor3242-013

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  ChaoY.YuanJ.LiS.JiaS.HanL.XuL. (2018). Analysis of transcripts and splice isoforms in red clover (Trifolium pratense L.) by single-molecule long-read sequencing.BMC Plant Biol.18:300. 10.1186/s12870-018-1534-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  ChenJ.HuW. J.WangC.LiuT. W.ChalifourA.ShenZ. J.et al (2014). Proteomic analysis reveals differences in tolerance to acid rain in two broad-leaf tree species, Liquidambar formosana and Schima superba.PLoS One9:e102532. 10.1371/journal.pone.0102532

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  ChenJ.TangX.RenC.WeiB.WuY.WuQ.et al (2018). Full-length transcriptome sequences and the identification of putative genes for flavonoid biosynthesis in safflower.BMC Genomics19:548. 10.1186/s12864-018-4946-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  DharshiniS.ChakravarthiM.ManojV. M.NaveenaraniM.KumarR.MeenaM.et al (2016). De novo sequencing and transcriptome analysis of a low temperature tolerant Saccharum spontaneum clone IND 00-1037.J. Biotechnol.231280–294. 10.1016/j.jbiotec.2016.05.036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  DongL.LiuH.ZhangJ.YangS.KongG.ChuJ. S.et al (2015). Single-molecule real-time transcript sequencing facilitates common wheat genome annotation and grain transcriptome research.BMC Genomics16:1039. 10.1186/s12864-015-2257-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  DongZ.ChenY. (2013). Transcriptomics: advances and approaches.Sci. China Life Sci.56960–967. 10.1007/s11427-013-4557-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  EmeryL.WhelanS.HirschiK. D.PittmanJ. K. (2012). Protein phylogenetic analysis of Ca2+/cation antiporters and insights into their evolution in plants.Front. Plant Sci.3:1. 10.3389/fpls.2012.00001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  FinnR. D.CoggillP.EberhardtR. Y.EddyS. R.MistryJ.MitchellA. L.et al (2016). The Pfam protein families database: towards a more sustainable future.Nucleic Acids Res.44D279–D285.

  • Google Scholar

• 13

  JiaoW. B.SchneebergerK. (2017). The impact of third generation genomic technologies on plant genome assembly.Curr. Opin. Plant Biol.3664–70. 10.1016/j.pbi.2017.02.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  KanehisaM.GotoS.KawashimaS.OkunoY.HattoriM. (2004). The KEGG resource for deciphering the genome.Nucleic Acids Res.32D277–D280.

  • Google Scholar

• 15

  KongL.ZhangY.YeZ. Q.LiuX. Q.ZhaoS. Q.WeiL.et al (2007). CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine.Nucleic Acids Res.35W345–W349.

  • Google Scholar

• 16

  LevyS. E.MyersR. M. (2016). Advancements in next-generation sequencing.Annu. Rev. Genomics Hum. Genet.1795–115.

  • Google Scholar

• 17

  LiW.JaroszewskiL.GodzikA. (2002). Tolerating some redundancy significantly speeds up clustering of large protein databases.Bioinformatics1877–82. 10.1093/bioinformatics/18.1.77

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  LuoC.ShenZ.LiX.BakerA. J. (2006). Enhanced phytoextraction of Pb and other metals from artificially contaminated soils through the combined application of EDTA and EDDS.Chemosphere631773–1784. 10.1016/j.chemosphere.2005.09.050

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  PadmavathiammaP. K.LiL. Y. (2007). Phytoremediation technology: hyper-accumulation metals in plants.Water Air Soil Pollut.184105–126. 10.1007/s11270-007-9401-5

  • CrossRef
  • Google Scholar

• 20

  PittmanJ. K.HirschiK. D. (2016). CAX-ing a wide net: cation/H+ transporters in metal remediation and abiotic stress signalling.Plant Biol.18741–749. 10.1111/plb.12460

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  QiaoD.YangC.ChenJ.GuoY.LiY.NiuS.et al (2019). Comprehensive identification of the full-length transcripts and alternative splicing related to the secondary metabolism pathways in the tea plant (Camellia sinensis).Sci. Rep.9:2709.

  • Google Scholar

• 22

  RhoadsA.AuK. F. (2015). PacBio sequencing and its applications.Genomics Proteomics Bioinformatics13278–289. 10.1016/j.gpb.2015.08.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  ShigakiT.ReesI.NakhlehL.HirschiK. D. (2006). Identification of three distinct phylogenetic groups of CAX cation/proton antiporters.J. Mol. Evol.63815–825. 10.1007/s00239-006-0048-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  SlatkoB. E.GardnerA. F.AusubelF. M. (2018). Overview of next-generation sequencing technologies.Curr. Protoc. Mol. Biol.122:e59.

  • Google Scholar

• 25

  SunL.LuoH.BuD.ZhaoG.YuK.ZhangC.et al (2013). Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts.Nucleic Acids Res.41:e166. 10.1093/nar/gkt646

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  TangC. Q.HanP. B.LiS.ShenL. Q.HuangD. S.LiY. F.et al (2020). Species richness, forest types and regeneration of Schima in the subtropical forest ecosystem of Yunnan, southwestern China.For. Ecosyst.71–19.

  • Google Scholar

• 27

  TatusovR. L.FedorovaN. D.JacksonJ. D.JacobsA. R.KiryutinB.KooninE. V.et al (2003). The COG database: an updated version includes eukaryotes.BMC Bioinformatics4:41. 10.1186/1471-2105-4-41

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  WangY.WangJ.ChengJ.XuD.JiangL. (2015). Genetic interactions between the Golgi Ca^{2 +} /H⁺ exchanger Gdt1 and the plasma membrane calcium channel Cch1/Mid1 in the regulation of calcium homeostasis, stress response and virulence in Candida albicans.FEMS Yeast Res.15:fov069. 10.1093/femsyr/fov069

  • CrossRef
  • Google Scholar

• 29

  WorkmanR. E.MyrkaA. M.WongG. W.TsengE.WelchK. C.Jr.TimpW. (2018). Single-molecule, full-length transcript sequencing provides insight into the extreme metabolism of the ruby-throated hummingbird Archilochus colubris.GigaScience7giy009.

  • Google Scholar

• 30

  XieJ.ZuoJ.HuangY.LiC.ChenY. (2020). The origin and germplasm collection for cultivated Dendrobium officinale K. Kimura & Migo individuals revealed by EST-SSR markers.Genet. Resour. Crop Evol.671209–1219. 10.1007/s10722-020-00906-7

  • CrossRef
  • Google Scholar

• 31

  YangS. X.DengH.LiM. S. (2008). Manganese uptake and accumulation in a woody hyperaccumulator Schima superba.Plant Soil Environ.54441–446. 10.17221/401-pse

  • CrossRef
  • Google Scholar

• 32

  YangZ.ZhangR.ZhouZ. (2021). Selection and Validation of Appropriate Reference Genes for Gene Expression Studies in Schima Superba.

  • Google Scholar

• 33

  ZhangR.YangH.ZhouZ.ShenB.XiaoJ.WangB. (2019). A high-density genetic map of Schima superba based on its chromosomal characteristics.BMC Plant Biol.19:41. 10.1186/s12870-019-1655-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  ZhengY.JiaoC.SunH.RosliH. G.PomboM. A.ZhangP.et al (2016). iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases.Mol. Plant91667–1670.

  • Google Scholar

• 35

  ZuoC.BlowM.SreedasyamA.KuoR. C.RamamoorthyG. K.Torres-JerezI.et al (2018). Revealing the transcriptomic complexity of switchgrass by PacBio long-read sequencing.Biotechnol. Biofuels11:170.

  • Google Scholar