Kangqi Zhou†
Xianhui Pan
Yong Lin- Guangxi Key Laboratory for Aquatic Genetic Breeding and Healthy Aquaculture, Guangxi Academy of Fishery Sciences, Nanning, China
Cipangopaludina chinensis is an economically important aquatic snail with high medicinal value. However, molecular biology research on C. chinensis is limited by the lack of a reference genome, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. chinensis. Herein, we conducted the first full-length transcriptome analysis of C. chinensis using PacBio single-molecule real-time (SMRT) sequencing technology. We identified a total of 26,312 unigenes with an average length of 2,572 bp, of which the largest number of zf-c2h2 transcription factor families (120,18.24%) were found, and also observed that the majority of the 8,058 SSRs contained 4-7 repeat units, which provided data for subsequent work on snail genetics Subsequently, 91.86% (24,169) of the genes were successfully annotated to the four major databases, while the highest homology was observed with Pomacea canaliculata. Functional annotation revealed that the majority of transcripts were enriched in metabolism, signal transduction and Immune-related pathways, and several candidate genes involved in drug metabolism and immune response were identified (e.g., CYP1A1, CYP2J, CYP2U1, GST, ,PIK3, PDE3A, PRKAG). This study lays a foundation for future molecular biology research and provides a reference for studying genes associated with the medicinal value of C. chinensis.
1 Introduction
Cipangopaludina chinensis is a common, large freshwater snail, that is, widely distributed in lakes, rice fields, ponds and other waters around the world. As a primary consumer, it plays an important role in energy transfer and water purification in the ecosystem. The meat and shell of the snail have edible and medicinal value, can be used as feed for poultry and aquatic animals, and are favored by consumers and farmers (Chen L. et al., 2019; Luo et al., 2021; Zhou et al., 2021). In recent years, with the rise of the Guangxi snail rice noodle industry in China, the demand for snails has increased, which has promoted the development of the snail industry from relying on natural fishing to artificial breeding (Qin., 2020). As a result, scholars are paying increasing attention to this snail farming model (He et al., 2018; Lei., 2020; Qiu., 2020; Jie., 2021; Wu et al., 2021). Some researchers have explored the potential of C. chinensis in medicine. In an alcohol-induced liver injury model, the significant antioxidant and hepatoprotective activity of C. chinensis polysaccharide (CCPS) was demonstrated (Jiang et al., 2013; Fan et al., 2014). And also found that CCPS can inhibit neovascularization of endothelial cells by mediating their “proliferation, migration and lumen formation” (Xiong et al., 2017; Xiong et al., 2020). Meanwhile, an in vitro model using a 2.2.15 cell line of human hepatocellular carcinoma cells (HepG2) revealed the significant anti-hepatitis B virus effect of CCPS (Liu et al., 2013). However, fewer studies have focused on the transcriptome level from the perspective of C. chinensis molecular biology and breeding so far.
Currently, single-molecule real-time (SMRT) long-read sequencing technology based on PacBio platform is the most reliable full-length cDNA molecular sequencing method, which can truly reflect the transcriptome information of the sequenced species and is conducive to the study of variable splicing, gene fusion, alleles and other mRNA structures, and is now widely used in the whole transcriptome analysis of plant and animal species (Korlach., 2013; Roberts et al., 2013; Lin et al., 2017; Ardui et al., 2018; Li et al., 2018; Chen S. et al., 2019; Yang et al., 2020). Here, based on the Sequel platform, we used PacBio’s SMRT to sequence the full-length transcriptome of C. chinensis. Then, we conducted a systematic analysis of all sequences, including functional annotation and identification of simple sequence repeats (SSRs), coding DNA sequences (CDSs), non-coding RNAs (ncRNAs), transcription factors (TFs), and alternative splicing (AS) events. We aimed to obtain a large amount of useful sequence information. This study obtained a full-length transcriptome of this snail for the first time, which will provide valuable data for subsequent gene function research, molecularly assisted breeding, and biomedical functions.
2 Materials and Methods
2.1 Sample Collection and RNA Extraction
A healthy individual of C. chinensis (weight 20.30 g, shell height 46.05 mm, shell width 33.27 mm) was collected from the rice snail breeding demonstration base in Ligao Village (23.37°N, 111.29°E), Liuzhou, Guangxi, China. Eight tissues were collected, namely the foot, intestine, hepatopancreas, gill, gonad, esophagus, blood, and lip. The samples were frozen with liquid nitrogen and stored at −80°C for later use. The trial was performed in accordance with the Institutional Animal Care and Use Committees of Guangxi Academy (CGA-00927).
The total RNA from these tissues was extracted separately, and equal amounts of RNA were mixed to construct a library after quality testing. Briefly, total RNA was extracted by grinding each tissue in TRIzol reagent (Life technologies, United States) on dry ice and processed following the manufacturer’s protocol. The integrity of the RNA was determined with an Agilent 2,100 Bioanalyzer and agarose gel electrophoresis. The purity and concentration of the RNA were determined with a Nanodrop micro-spectrophotometer (Thermo Fisher).
2.2 Library Construction and Single-Molecule Real-Time Sequencing
Referring to the methods of Huang et al., (2021). Briefly, mRNA was enriched using Oligo (dT) magnetic beads (PacBio biosciences, United States). The enriched mRNA (2 μg) was reverse transcribed into cDNA using a Clontech SMARTer PCR cDNA Synthesis Kit (Clontech, United States). PCR cycle optimization was used to determine the optimal amplification cycle number for the downstream large-scale PCR reactions. The optimized cycle number was then used to generate double-stranded cDNA. In addition, >5-kb size selection was performed using the BluePippin™ Size-Selection System and the resulting cDNA was mixed equally with non-size-selected cDNA. Large-scale PCR was performed for subsequent SMRTbell library construction. The cDNA samples were barcoded, then pooled and constructed as a “single” sample in the SMRTbell library, with different barcodes distinguishing the data from the different samples. The cDNAs were DNA damage repaired, end repaired, and ligated to sequencing adapters. The SMRTbell template was annealed to the sequencing primer and bound to polymerase, and sequenced on the PacBio Sequel II platform by Gene Denovo Biotechnology Co. (Guangzhou, China). A total of 57 Gb reads were generated from the SMRT library.
2.3 Data Processing
The raw sequencing reads of the cDNA libraries were analyzed using the unigene sequencing (Iso-Seq) pipeline supported by Pacific Biosciences (Gordon et al., 2015). First, high-quality circular consensus sequences (CCSs) were extracted from the subreads BAM file. The integrity of the transcripts was judged on the basis of whether the CCS reads contained 5′ primer, 3′ primer, and poly A structures. The sequences containing all three structures were called full-length sequences (FL reads). Subsequently, primers, barcodes, poly A tails, and concatemers of full passes were removed to obtain full-length non-chimeric (FLNC) reads. The FLNC reads were clustered to generate complete unigenes. Similar FLNC reads were clustered hierarchically using minimap2 to obtain a consistent sequence (unpolished consensus unigenes) (Li., 2018). The quiver algorithm was then used to further correct the consistent sequences. Two strategies were used to improve the accuracy of PacBio reads. First, we used non-full-length reads to polish the obtained cluster consensus Unigenes with the Quiver software to produce FL, polished, high-quality consensus sequences (accuracy of ≥99%). Second, Illumina short reads were obtained from the same samples with the LoRDEC to further correct low-quality isoforms (Salmela & Rival, 2014). On the basis of the results, both high-quality unigenes (prediction accuracy was ≥ 0.99) were combined for subsequent analysis.
2.4 Basic Annotation of Unigenes
To annotate them, the unigenes were BLAST analyzed against the NCBI non-redundant protein (Nr) database (https://www.ncbi.nlm.nih.gov), the Swiss-Prot protein database (https://www.expasy.ch/sprot), the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Kanehisa & Goto, 2000), and the COG/KOG database (https://www.ncbi.nlm.nih.gov/COG) with the BLASTx program (https://www.ncbi.nlm.nih.gov/BLAST/) at an E-value threshold of 1e-5 to evaluate sequence similarity to genes of other species. Gene Ontology (GO) annotation was performed using the Blast2GO software with the Nr annotation results for the unigenes (Conesa et al., 2005). Unigenes with the top 20 highest scores and no shorter than 33 HSP (High-scoring Segment Pair) hits were selected for Blast2GO analysis (Table 1). Then, functional classification of the unigenes was performed using the WEGO software (Ye et al., 2006).
TABLE 1. Description of full-length sequencing in Cipangopaludina chinensis.
2.5 Structure Analysis
Open reading frames (ORFs) were detected using the ANGEL software with the unigene sequences to obtain CDSs, protein sequences, and UTR sequences (Shimizu et al., 2006). Protein-coding sequences of the unigenes were aligned with hmmscan to Animal TFdb (https://bioinfo.life.hust.edu.cn/AnimalTFDB/#!/) to predict TF families. CNCI (version 2) was used to assess the protein-coding potential of transcripts without annotations with default parameters for potential ncRNAs (Sun et al., 2013). ncRNAs were classified according to their secondary structures and sequence conservation.
2.6 Simple Sequence Repeats Prediction and Primer Design
The MIcroSAtellite (MISA, https://pgrc.ipk-gatersleben.de/misa/) tool was employed for microsatellite mining across the whole transcriptome. Initially SSRs of 1–6 nucleotide units were identified with the minimum repeat unit number defined as 15 for mono-nucleotides, 6 for di-nucleotides, 5 for tri-nucleotides, and 4 each for tetra-, penta-, and hexa-nucleotides. If the distance between two SSRs was shorter than 100 bp, they were considered one SSR.
2.7 Alternative Splicing Detection
To analyze AS events in the transcript unigenes, COding GENome reconstruction Tool (Cogent) was first used to partition transcripts into gene families based on k-mer similarity and reconstruct each family into a coding reference genome based on De Bruijn graph methods (Li et al., 2017). Then, the SUPPA tool was used to analyze AS events in the transcript unigenes (Alamancos et al., 2015).
3 Results
3.1 Single-Molecule Real-Time Sequencing Data Analysis
Using the PacBio Sequel II sequencing platform, we obtained 57,428,331,313 bp of original transcript data (about 57 Gb) and a total of 26,658,437 subreads with an average subread length of 2,154 bp and an N50 length of 2,466 bp. After CCS analysis, 818,895 CCS reads were obtained with a mean CCS read length and pass number of 2,428 bp and 29, respectively. By clustering and correction, 697,861 FLNC reads were identified with an average length of 2,315 bp and an N50 length of 2,530 bp. Finally, the redundant reads were removed by CD-HIT and 26,312 unigenes with a mean length of 2,572 bp were obtained (Table 1). The length distribution of unigenes is shown in Figure 1. The majority of unigenes were between 1,000 and 5,000 bp, and the longest unigene was about 14,000 bp.
FIGURE 1. Length distribution of unigenes obtained from the Cipangopaludina chinensis library.
3.2 Functional Annotation of C. chinensis Sequences
The 26,312 unigenes were annotated using four databases, with 24,151 (91.79%), 22,756 (86.49%), 13,934 (52.96%), and 16,423 (62.42%) unigenes matched to the Nr, KEGG, KOG, and Swiss-Prot databases, respectively (Figure 2A). A total of 24,169 (91,86%) unigenes were annotated in all the databases (Figure 2B). After aligning the unigene sequences with the Nr database using BlastX, the sequence with the best comparison result for each unigene was taken as the corresponding homologous sequence, the species to which the homologous sequence belonged was determined, and the homologous sequences of each species were counted. The results showed that the species with the highest homology was Pomacea canaliculata (17,089 unigenes), followed by Mizuhopecten yessoensis (1,199) and Aplysia californica (969) (Figure 2C).
FIGURE 2. Functional annotation of Cipangopaludina chinensis. (A) Statistics of the transcripts annotated in different databases. (B) Venn diagram of annotations in the NR, GO, KEGG, KOG, and Swiss-Prot databases. (C) Distribution of the top 10 species with matched transcripts in the NR database. (D) Distribution of GO terms for all annotated transcripts in the biological process, cellular component, and molecular function ontologies. (E) KEGG pathways enriched in the transcripts. (F) COG categories of the transcripts.
A total of 13,214 sequences were annotated using the GO function classification library under the three broad categories: molecular function (17,422 unigenes), cellular component (48,494), and biological process (80,752). “Binding” (8,587, 49.29%), “cell” (8,587, 17.71%), and “cellular process” (10,215, 12.65%) were the most common terms among the annotated unigenes in the three categories mentioned above (Figure 2D).
In KEGG pathway analysis, the unigenes were assigned to five main categories: cellular processes (4,595 unigenes), environmental information processing (4,728), genetic information processing (1,797), metabolism (8,156), and organismal systems (10,363). The signal transduction (4,297, 16.33%) category had the largest number of unigenes. Next was infectious diseases (4,230, 16.08%) and cancers (3,068, 11.66%) (Figure 2E). Based on the KEGG annotation information we were able to further obtain Pathway annotations for unigene and focus on some of the more applicable pathways including PI3K-Akt signaling pathway (436, 5.34%), cGMP-PKG signaling pathway (389, 4.77%), Drug metabolism-cytochrome P450 (91, 1.12%), Platinum drug resistance (127, 1.56%), AMPK signaling pathway (267,3.27%), NF-kappa B signaling pathway (119, 1.46%), T cell receptor signaling pathway (71, 0.87%), Toll and Imd signaling pathway (80, 0.98%) (Table 2). From these, we also identified some candidate genes involved in drug metabolism and immune related, such as cytochrome P450 family (CYP1A1, CYP2J and CYP2U1), phosphoinositide-3-kinase regulatory subunit (PIK3), cGMP-inhibited 3′,5′- cyclic phosphodiesterase A (PDE3A), glutathione S-transferase (GST), 5′-AMP-activated protein kinase (PRKAG), NF-kappa-B inhibitor alpha (NFKBIA), cell division control protein 42 (CDC42), baculoviral IAP repeat-containing protein 2/3 (BIRC2_3), and others (Table 2). These provide data for the subsequent development of biologic drugs in snail.
TABLE 2. Regional distribution of some SSRs in the full-length transcripts of Cipangopaludina chinensis.
All unigenes were classified into 25 categories in the KOG functional classification. (Figure 2F). The highest number of unigenes was annotated to the general function prediction only category (2,807, 10.67%), followed by the signal transduction mechanisms (2,563, 9.74%) and posttranslational modification, protein turnover, and chaperones (2,196, 8.35%) categories.
3.3 Structure Analysis of C. chinensis
3.3.1 Coding DNA Sequences Prediction
The numbers and lengths of 5′ UTRs, 3′ UTRs, and CDSs were identified using the ANGEL software. 5′ UTRs were observed in 23,559 unigenes and the lengths were mostly less than 1,000 nt, while 3′ UTRs were observed in 23,772 unigenes and the lengths were mainly less than 3,000 nt. We predicted that there were 24,991 CDSs in the library. Most of the CDSs (20,491 unigenes, 81.99%) were less than 2,000 nt in length, while 4,396 unigenes (17.59%) were in the range of 2,000–5,000 nt and only 104 unigenes (representing 0.42%) exceeded 5,000 nt (Figure 3A).
FIGURE 3. Sequence structure analysis of Cipangopaludina chinensis. (A) Length distribution of CDSs. (B) Venn diagram of lncRNAs identified by the CNCI and CPC methods. (C) Length density distribution of the lncRNA and mRNA transcripts. (D) Families and numbers of the top 10 TFs produced by SMRT. (E) Summary of SSR types in the full-length transcripts of C. chinensis. (F) Alternative splicing events predicted by SUPPA.
3.3.2 Identification of ncRNAs and Transcription Factors
ncRNA analysis was performed on 2,143 full-length transcribed sequences that were not annotated to the four major databases. We identified 1,287 and 1,952 ncRNAs using CPC and CNCI, respectively (Figure 3B). Integrating the results of these two methods, a total of 1,277 ncRNA transcripts were predicted, accounting for only 4.85% of the total data. Additionally, comparing the length distribution density of mRNAs with the predicted ncRNAs showed that the proportion of ncRNAs in the 64–2,000 nt range was greater than that of the mRNAs (Figure 3C).
TFs are protein molecules with a specific structure and function in regulating gene expression in aquatic animals. We compared the predicted protein sequences with the Animal TFdb database using the hmmscan software and identified 658 TFs from 48 TF families. In Figure 3D, the top 10 TF families are listed, among which zf-c2h2 family (120, 18.24%) has the largest number, followed by the IF-bZIP family (94, 14.29%) (Figure 3D).
3.3.3 Simple Sequence Repeats Analysis
The MISA software was used to search all unigenes of the transcriptome and a total of 8,058 SSRs were found. The SSR statistics showed that the most abundant type was di-nucleotide repeat sequences (4,647), accounting for 57.67%, followed by tri-nucleotide (1,274, 15.81%), tetra-nucleotide (1,239, 15.37%), and mono-nucleotide repeats (816, 10.13%). However, the proportions of penta-nucleotide and hexa-nucleotide repetitive sequences were very small, only 0.48 and 0.53%, respectively. In addition, most SSRs had 4–7 repeat units (Figure 3E). Because CDSs were not predicted for some transcripts, the total count of SSRs in different regions was 6,037. Among these regions, 3′ UTRs had the most (4,667), followed by 5′ UTRs and CDSs with the same number, 685 (Table 2).
3.3.4 Alternative Splicing Prediction
Based on the Cogent software assembly sequence as a reference, the SUPPA software was used for AS event prediction. In our study, a total of 1,085 unigenes were assigned to six AS events (alternative 3′ splice sites, alternative 5′ splice sites, alternative first exons, mutually exclusive exons, retained introns, and exon skipping). Retained introns (RIs) accounted for the highest proportion of all events, at 66.45% (721). Next were alternative 3′ splice sites (A3) and alternative 5′ splice sites (A5), with 180 (16.59%) and 162 (14.93%) events, respectively. The other three events, alternative first exons (15, 1.39%), mutually exclusive exons (6, 0.55%), and exon skipping (1, 0.09%), accounted for only 2.03% (Figure 3F).
Discussion
SMRT sequencing has the characteristics of high accuracy, reliability, long reading, and sequencing speed, and it enables animal genomic research without a reference genome (Hestand & Ameur, 2019; Jayakumar & Sakakibara, 2019; Zhang et al., 2019; Xiao & Zhou, 2020; Deng et al., 2021). Since the C. chinensis is an economically important snail with high medicinal value, and it has vital scientific research potential (Zhou et al., 2022). The current study was the first to obtain full-length transcriptome data of the snail using third-generation sequencing technology. By data analysis, 26,312 unigenes with a mean length of 2,572 bp were obtained in the study. Deng et al. used the same technology to sequence the no-reference-genome species Coelomactra antiquata and found a total of 39,209 unigenes with an average length of 2,732 bp (Deng et al., 2021). Similarly, Yang et al. found 63,801 full-length transcripts in the three developmental stages of Rhynchophorus ferrugineus, with an average length of 2,964 bp (Yang et al., 2020). Compared with the above organisms, the transcript data volume and sequence length of C. chinensis were relatively small.
We performed functional annotations on all the obtained transcripts in the four major databases, of which 91.86% of unigenes were successfully annotated and only 8.14% were unannotated. The percentage of annotated unigenes is very high, possibly because there have been many molecular biology studies on gastropod mollusks (e.g., P. canaliculata, Cipangopaludina cathayensis, and Tegillarca granosa) in the past, and the data collected in these databases are relatively complete; thus, the genomic information of C. chinensis can be compared with existing data (Chen et al., 2017). According to the species distribution statistics of the NR database, P. canaliculata has the most homologous sequences, revealing that P. canaliculata is closely related to C. chinensis. This coincides with the findings of Mu et al., who performed an evolutionary relationship analysis based on the mitochondrial genome, which suggested that P. canaliculata has a close evolutionary relationship with C. cathayensis (homologous to C. chinensis) (Mu et al., 2015). In terms of the annotation results from the GO, KEGG, and KOG databases, the number of transcripts identified in KEGG was highest, among which the numbers of transcripts involved in intracellular signal transduction, infectious diseases, and cancer were the largest. In-depth analysis also observed a number of pathways and candidate genes involved in drug metabolism and immune response (Table 2). Previously, some scholars found that C. chinensis flesh is rich in biopolysaccharides and confirmed that they play an active role in anti-tumor activity, lowering blood sugar, and anti-hepatitis B virus activity (Jiang et al., 2013; Shi et al., 2016; Xiong et al., 2020). This implies that C. chinensis genes may have many biological activities, participate in material metabolism and immune processes in animals, and have important biomedical value.
The structure of the non-redundant unigenes was also analyzed in our study. A total of 1,277 ncRNAs and 8,058 SSRs were identified. ncRNAs have complex and precise regulatory function in cell differentiation, ontogeny, and gene expression (Wu et al., 2014; Quan et al., 2018). The transcripts produced by 4–9% of mammalian genome sequences are ncRNAs (Lu and Huang, 2017; Han et al., 2018). In recent years, research on ncRNAs has progressed rapidly, but the functions of most ncRNAs are still unclear, especially in lower animals (Niu et al., 2019). In contrast, research on SSRs has been relatively comprehensive and SSRs have been widely used in practical applications, such as genetic diversity analysis, genetic map construction, and molecular marker development (Ali et al., 2019; Sun et al., 2019; Cregan et al., 2020) It was reported that 16,717 and 20,106 SSRs were found in P. canaliculata and C. antiquata, respectively, based on transcriptome sequencing (Mu et al., 2015). This shows that the number of SSRs found in C. chinensis in this study was relatively low. Furthermore, we also predicted the TFs of C. chinensis and obtained detailed TF numbers and family classifications, among which the zf-c2h2 family was the most abundant. In plants, zf-c2h2 TFs participate in all aspects of growth and development (especially root and floral development), as well as biotic and abiotic stress responses (Agarwal et al., 2007; Jiang & Pan, 2012; Liu et al., 2015). In animals, zf-c2h2 TFs are more involved in tumors, cancers, and related gene regulation (Hu et al., 2019; Hu et al., 2021). However, the role of zf-c2h2 TFs in the lower animal C. chinensis needs further research.
Conclusion
C. chinensis is an aquatic animal product with high economic value and has potential for medicinal development. More importantly, it plays an important role as a carbon sink in the ecosystem. In this study, PacBio SMRT sequencing technology was used to successfully construct a high-quality transcript library of C. chinensis (26,312 unigenes with an average length of 2,572 bp), and the sequences were functionally annotated (including 24,151 Nr, 22,756 KEGG, 13,934 KOG, and 16,423 Swiss-Prot sequences) and structurally analyzed (including 24,991 CDSs, 1,277 ncRNAs, 6 AS events, 658 TFs, and 8,058 SSRs) using bioinformatics software. These results provide comprehensive scientific data for future research on functional gene mining, molecular breeding, and biomedicine research of C. chinensis.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA817831/.
Ethics Statement
The animal study was reviewed and approved by the Institutional Animal Care and Use Committees of Guangxi Academy (CGA-00927).
Author Contributions
YL, KZ, and XP together designed this experiment, analyzed data, and wrote the paper. LW and XD taken part in doing this experiment and proofing the manuscript. YH, ZC, and JQ performed data collection and part experiment.
Funding
This study was supported by the Natural Science Foundation of Guangxi Province (2020GXNSFBA297067), Specific Research Project of Guangxi for Bases and Talents (AD21220010), Guangxi Innovation-Driven Development Project (AA20302019-5) and The National Modern Agricultural Technology System (nycytxgxcxtd-2021-08).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Agarwal, P., Arora, R., Ray, S., Singh, A. K., Singh, V. P., Takatsuji, H., et al. (2007). Genome-wide Identification of C2H2 Zinc-Finger Gene Family in Rice and Their Phylogeny and Expression Analysis. Plant Mol. Biol. 65 (4), 467–485. doi:10.1007/s11103-007-9199-y
Alamancos, G. P., Pagès, A., Trincado, J. L., Bellora, N., and Eyras, E. (2015). Leveraging Transcript Quantification for Fast Computation of Alternative Splicing Profiles. Rna 21 (9), 1521–1531. doi:10.1261/rna.051557.115
Ali, A., Pan, Y. B., Wang, Q. N., Wang, J. D., Chen, J. L., and Gao, S. J. (2019). Genetic Diversity and Population Structure Analysis of Saccharum and Erianthus Genera Using Microsatellite (SSR) Markers. Sci. Rep. 9 (1), 395. doi:10.1038/s41598-018-36630-7
Ardui, S., Ameur, A., Vermeesch, J. R., and Hestand, M. S. (2018). Single Molecule Real-Time (SMRT) Sequencing Comes of Age: Applications and Utilities for Medical Diagnostics. Nucleic acids Res. 46 (5), 2159–2168. doi:10.1093/nar/gky066
Chen, H., Xiao, G., Chai, X., Lin, X., Fang, J., and Teng, S. (2017). Transcriptome Analysis of Sex-Related Genes in the Blood Clam Tegillarca Granosa. PLoS One 12 (9), e0184584. doi:10.1371/journal.pone.0184584
Chen, L., Du, X., Wen, Y., Huang, J., Luo, F., Fan, H., et al. (2019a). Analysis of Meat Content and Nutrient Composition of Three Kinds of Snails in Guangxi Area. Meat Ind. 10, 20–22.
Chen, S., Qiu, G., and Yang, M. (2019b). SMRT Sequencing of Full-Length Transcriptome of Seagrasses Zostera Japonica. Sci. Rep. 9 (1), 14537. doi:10.1038/s41598-019-51176-y
Conesa, A., Gotz, S., Garcia-Gomez, J. M., Terol, J., Talon, M., and Robles, M. (2005). Blast2GO: a Universal Tool for Annotation, Visualization and Analysis in Functional Genomics Research. Bioinformatics 21 (18), 3674–3676. doi:10.1093/bioinformatics/bti610
Cregan, P. B., Akkaya, M. S., Bhagwat, A. A., Lavi, U., and Rongwen, J. (2020). “Length Polymorphisms of Simple Sequence Repeat (SSR) DNA as Molecular Markers in Plants,” in Plant Genome Analysis (Boca Raton: CRC Press), 47–56. doi:10.1201/9781003068907-5
Deng, A., Li, J., Yao, Z., Afriyie, G., Chen, Z., Guo, Y., et al. (2021). SMRT Sequencing of the Full-Length Transcriptome of the Coelomactra Antiquata. Front. Genet. 12, 741243. doi:10.3389/fgene.2021.741243
Fan, Q., Li, C., and Wang, K. (2014). Protective Effect of Cipangopaludina Chinensis Polysaccharide on Alcoholic Hepatic Injury in Mice. J. Pathogen Biol. 9 (2), 105–108.
Gordon, S. P., Tseng, E., Salamov, A., Zhang, J., Meng, X., Zhao, Z., et al. (2015). Widespread Polycistronic Transcripts in Fungi Revealed by Single-Molecule Mrna Sequencing. Plos One 10 (7), e0132628. doi:10.1371/journal.pone.0132628
Han, C., Hu, J., Hu, S., Zhu, R., Bai, Y., and Zhu, H. (2018). Research Progress of Long Non-coding RNA. Lett. Biotechnol. 29 (1), 123–130.
He, J., Wu, H., Feng, P., Chen, Z., Gan, X., Li, H., et al. (2018). Comprehensive Planting and Breeding Technology of Rice Snails. Henan Nongye 23, 50–51.
Hestand, M. S., and Ameur, A. (2019). The Versatility of SMRT Sequencing. Genes 10 (1), 24. doi:10.3390/genes10010024
Hu, H., Zhang, Q., Hu, F. F., Liu, C. J., and Guo, A. Y. (2021). A Comprehensive Survey for Human Transcription Factors on Expression, Regulation, Interaction, Phenotype and Cancer Survival. Brief. Bioinform 22 (5), bbab002. doi:10.1093/bib/bbab002
Hu, H., Miao, Y.-R., Jia, L.-H., Yu, Q.-Y., Zhang, Q., and Guo, A.-Y. (2019). AnimalTFDB 3.0: a Comprehensive Resource for Annotation and Prediction of Animal Transcription Factors. Nucleic acids Res. 47 (D1), D33–D38. doi:10.1093/nar/gky822
Huang, Y., Zhang, L., Huang, S., and Wang, G. (2021). Full-length Transcriptome Sequencing of Heliocidaris Crassispina Using PacBio Single-Molecule Real-Time. Fish Shellfish Immunol. 120, 507–514. doi:10.1016/j.fsi.2021.12.014
Jayakumar, V., and Sakakibara, Y. (2019). Comprehensive Evaluation of Non-hybrid Genome Assembly Tools for Third-Generation PacBio Long-Read Sequence Data. Briefings Bioinforma. 20 (3), 866–876. doi:10.1093/bib/bbx147
Jiang, C., Jiao, Y., Chen, X., Li, X., Yan, W., Yu, B., et al. (2013). Preliminary Characterization and Potential Hepatoprotective Effect of Polysaccharides from Cipangopaludina Chinensis. Food Chem. Toxicol. 59, 18–25. doi:10.1016/j.fct.2013.05.036
Jiang, L., and Pan, L.-j. (2012). Identification and Expression of C2H2 Transcription Factor Genes in Carica Papaya under Abiotic and Biotic Stresses. Mol. Biol. Rep. 39, 7105–7115. doi:10.1007/s11033-012-1542-y
Jie, Q. (2021). Research Progress on Rice Field Culture of Cipangopaludina Cahayensis. Rural Econ. Science-Technology 21, 56–57.
Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucl. Acids Res. 28, 27–30. doi:10.1093/nar/28.1.27
Lei, D. (2020). Analysis of Challenges and Countermeasures Facing Liuzhou Snail Powder Industry under the Background of "Internet +. Enterp. Sci. Technol. Dev. 7, 23–24.
Li, H. (2018). Minimap2: Pairwise Alignment for Nucleotide Sequences. Bioinformatics 34 (18), 3094–3100. doi:10.1093/bioinformatics/bty191
Li, J., Harata-Lee, Y., Denton, M. D., Feng, Q., Rathjen, J. R., Qu, Z., et al. (2017). Long Read Reference Genome-free Reconstruction of a Full-Length Transcriptome from Astragalus Membranaceus Reveals Transcript Variants Involved in Bioactive Compound Biosynthesis. Cell. Discov. 3, 17031. doi:10.1038/celldisc.2017.31
Li, Y., Fang, C., Fu, Y., Hu, A., Li, C., Zou, C., et al. (2018). A Survey of Transcriptome Complexity in Sus scrofa Using Single-Molecule Long-Read Sequencing. DNA Res. 25 (4), 421–437. doi:10.1093/dnares/dsy014
Lin, H., Lin, X., Zhu, J., Yu, X. Q., Xia, X., Yao, F., et al. (2017). Characterization and Expression Profiling of Serine Protease Inhibitors in the Diamondback Moth, Plutella Xylostella (Lepidoptera: Plutellidae). BMC genomics 18 (1), 162. doi:10.1186/s12864-017-3583-z
Liu, Q., Wang, Z., Xu, X., Zhang, H., Li, C., and Prasad, M. (2015). Genome-Wide Analysis of C2H2 Zinc-Finger Family Transcription Factors and Their Responses to Abiotic Stresses in Poplar (Populus trichocarpa). Plos One 10 (8), e0134753. doi:10.1371/journal.pone.0134753
Liu, X. Y., Li, C. P., and Wang, K. X. (2013). Experimental Study on Polysaccharide of Cipangopaludina Chinensis against HBV In Vitro. Zhongguo Zhong Yao Za Zhi 38 (6), 879–883.
Lu, C., and Huang, Y. H. (2017). Progress in Long Non-coding RNAs in Animals. Yi Chuan 39 (11), 1054–1065. doi:10.16288/j.yczz.17-120
Luo, H., Chen, L., Jing, T., Sun, W., Li, Z., Zhou, M., et al. (2021). Muscle Main Nutrients of Four Snails in the Viviparidae. J. Fish. China 46 (10), 1–8.
Mu, X., Hou, G., Song, H., Xu, P., Luo, D., Gu, D., et al. (2015). Transcriptome Analysis between Invasive Pomacea Canaliculata and Indigenous Cipangopaludina Cahayensis Reveals Genomic Divergence and Diagnostic Microsatellite/SSR Markers. BMC Genet. 16 (1), 12–13. doi:10.1186/s12863-015-0175-2
Niu, C., Xue, L., Ji, H., and Li, S. (2019). Progress in Research on Biological Function of lncRNA. Chin. J. Biol. 2, 228–232.
Quan, J., Pan, X., Zhao, L., Li, Z., Dai, K., Yan, F., et al. (2018). LncRNA as a Diagnostic and Prognostic Biomarker in Bladder Cancer: a Systematic Review and Meta-Analysis. Ott 11, 6415–6424. doi:10.2147/ott.s167853
Roberts, R. J., Carneiro, M. O., and Schatz, M. C. (2013). The Advantages of SMRT Sequencing. Genome Biol. 14 (6), 405. doi:10.1186/gb-2013-14-6-405
Salmela, L., and Rivals, E. (2014). LoRDEC: Accurate and Efficient Long Read Error Correction. Bioinformatics 30 (24), 3506–3514. doi:10.1093/bioinformatics/btu538
Shi, Y., Xiong, Q., Wang, X., Li, X., Yu, C., Wu, J., et al. (2016). Characterization of a Novel Purified Polysaccharide from the Flesh of Cipangopaludina Chinensis. Carbohydr. Polym. 136, 875–883. doi:10.1016/j.carbpol.2015.09.062
Shimizu, K., Adachi, J., and Muraoka, Y. (2006). ANGLE: a Sequencing Errors Resistant Program for Predicting Protein Coding Regions in Unfinished cDNA. J. Bioinform. Comput. Biol. 04 (03), 649–664. doi:10.1142/s0219720006002260
Sun, L., Luo, H., Bu, D., Zhao, G., Yu, K., Zhang, C., et al. (2013). Utilizing Sequence Intrinsic Composition to Classify Protein-Coding and Long Non-coding Transcripts. Nucleic Acids Res. 41, e166. doi:10.1093/nar/gkt646
Sun, L., Zhou, Q., Huang, M., Zhou, C., Xu, S., Yang, S., et al. (2019). Using the SSR Fluorescent Labling to Establish SSR Fingerprint for 131 Cultivars in Bougainvillea (Nyctaginaceae). J. Nanjing For. Univ. Nat. Sci. Ed. 43 (6), 60–66. doi:10.3969/j.issn.1000-2006.201812045
Wu, W., Ding, X., Li, W., Liu, C., Xie, S., and Wu, Z. (2021). Cultivation and Breeding Techniques of Alpine Thin Shell Snails and Rice Fields. Special Econ. Animal Plant 24 (10), 4.
Wu, Z., Liu, X., Liu, L., Deng, H., Zhang, J., Xu, Q., et al. (2014). Regulation of lncRNA Expression. Cell. Mol. Biol. Lett. 19 (4), 561–575. doi:10.2478/s11658-014-0212-6
Xiao, T., and Zhou, W. (2020). The Third Generation Sequencing: the Advanced Approach to Genetic Diseases. Transl. Pediatr. 9 (2), 163–173. doi:10.21037/tp.2020.03.06
Xiong, Q., Hao, H., He, L., Jing, Y., Xu, T., Chen, J., et al. (2017). Anti-inflammatory and Anti-angiogenic Activities of a Purified Polysaccharide from Flesh of Cipangopaludina Chinensis. Carbohydr. Polym. 176 (1), 152–159. doi:10.1016/j.carbpol.2017.08.073
Xiong, Q., Li, H., Zhou, L., Liang, J., Zhang, Z., Han, Y., et al. (2020). A Sulfated Polysaccharide from the Edible Flesh of Cipangopaludina Chinensis Inhibits Angiogenesis to Enhance Atherosclerotic Plaque Stability. J. Funct. Foods 66, 103800. doi:10.1016/j.jff.2020.103800
Yang, H., Xu, D., Zhuo, Z., Hu, J., and Lu, B. (2020). SMRT Sequencing of the Full-Length Transcriptome of the Rhynchophorus Ferrugineus (Coleoptera: Curculionidae). PeerJ 8, e9133. doi:10.7717/peerj.9133
Ye, J., Fang, L., Zheng, H., Zhang, Y., Chen, J., Zhang, Z., et al. (2006). WEGO: a Web Tool for Plotting GO Annotations. Nucleic Acids Res. 34 (2), W293–W297. doi:10.1093/nar/gkl031
Zhang, W., Jia, B., and Wei, C. (2019). PaSS: a Sequencing Simulator for PacBio Sequencing. BMC Bioinforma. 20 (1), 352–357. doi:10.1186/s12859-019-2901-7
Zhou, K., Lin, Y., Pang, H., Wei, Z., Qin, J., Huang, Y., et al. (2021). Comparison Analysis of Muscles Nutrient Composition and Digestive Enzyme in Cipangopaludina Chinensis. Freshw. Fish. 6 (12), 13–18.
Keywords: Cipangopaludina chinensis, full-length transcriptome, SMRT sequencing, functional annotation, structure prediction
Citation: Zhou K, Chen Z, Du X, Huang Y, Qin J, Wen L, Pan X and Lin Y (2022) SMRT Sequencing Reveals Candidate Genes and Pathways With Medicinal Value in Cipangopaludina chinensis. Front. Genet. 13:881952. doi: 10.3389/fgene.2022.881952
Received: 23 February 2022; Accepted: 26 May 2022;
Published: 16 June 2022.
Edited by:
Huabing Li, Shanghai Jiao Tong University, ChinaReviewed by:
Cheng Qin, Zunyi Vocational and Technical College, ChinaTsukasa Fukunaga, Waseda University, Japan
Copyright © 2022 Zhou, Chen, Du, Huang, Qin, Wen, Pan and Lin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xianhui Pan, MTUwMDIzODEyNjFAMTYzLmNvbQ==; Yong Lin, bGlubm4yMDA1QDEyNi5jb20=
†These authors have contributed equally to this work and share first authorship