Sequencing and de novo assembly of visceral mass transcriptome of the critically endangered land snail Satsuma myomphala: Annotation and SSR discovery

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers

BACKGROUND

The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult.

OBJECTIVE

In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi.

RESULTS

After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT.

CONCLUSION

The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

A review of the endangered mollusks transcriptome under the threatened species initiative of Korea

Transcriptome studies for conservation of endangered mollusks is a proactive approach towards managing threats and uncertainties facing these species in natural environments. The population of these species is declining due to habitat destruction, illicit wildlife trade, and global climate change. These activities risk the free movement of species across the wild landscape, loss of breeding grounds, and restrictions in displaying the physiological attributes so crucial for faunal welfare. Gastropods face the most negative ecological effects and have been enlisted under Korea's protective species consortium based on their population dynamics in the last few years. Moreover, with the genetic resources restricted for such species, conservation by informed planning is not possible. This review provides insights into the activities under the threatened species initiative of Korea with special reference to the transcriptome assemblies of endangered mollusks. The gastropods such as Ellobium chinense, Aegista chejuensis, Aegista quelpartensis, Incilaria fruhstorferi, Koreanohadra kurodana, Satsuma myomphala, and Clithon retropictus have been represented. Moreover, the transcriptome summary of bivalve Cristaria plicata and Caenogastropoda Charonia lampas sauliae is also discussed. Sequencing, de novo assembly, and annotation identified transcripts or homologs for the species and, based on an understanding of the biochemical and molecular pathways, were ascribed to predictive gene function. Mining for simple sequence repeats from the transcriptome have successfully assisted genetic polymorphism studies. A comparison of the transcriptome scheme of Korean endangered mollusks with the genomic resources of other endangered mollusks have been discussed with homologies and analogies for dictating future research.

Analysis of Population Genetic Diversity and Genetic Structure of Schizothorax biddulphi Based on 20 Newly Developed SSR Markers

To protect the germplasm resources of Schizothorax biddulphi, we developed and used 20 pairs of polymorphic microsatellite primers to analyze the genetic diversity and structure of populations. A total of 126 samples were collected from the Qarqan River (CEC), Kizil River (KZL), and Aksu River (AKS) in Xinjiang, China. The results showed that 380 alleles were detected in 20 pairs of primers and the average number of alleles was 19.0. The effective allele numbers and Nei’s gene diversity ranged from 1.1499 to 1.1630 and 0.0962 to 0.1136, respectively. The Shannon index range suggested low levels of genetic diversity in all populations. The genetic distance between the CEC and AKS populations was the largest, and the genetic similarity was the smallest. There was a significant genetic differentiation between CEC and the other two populations. The UPGMA clustering tree was constructed based on population genetic distance, and the clustering tree constructed by individuals showed that the AKS population and KZL population were clustered together, and the CEC population was clustered separately. Also, the group structure analysis also got the same result. It can be seen that although the three populations of S. biddulphi do not have high genetic diversity, the differentiation between the populations was high and the gene flow was limited, especially the differentiation between the CEC population and the other two populations. This study not only provided genetic markers for the research of S. biddulphi but the results of this study also suggested the need for enhanced management of S. biddulphi populations.

Constructing a de novo transcriptome and a reference proteome for the bivalve Scrobicularia plana: Comparative analysis of different assembly strategies and proteomic analysis

Scrobicularia plana is a coastal and estuarine bivalve widely used in ecotoxicological studies. However, the underlying molecular mechanisms for S. plana pollutant responses are hardly known due to the lack of molecular databases. Thus, in this study we present a holistic approach to assess a robust reference transcriptome and proteome of this clam. A mixture of control and metal-exposed individuals was used for mRNA isolation. Four sets of high quality filtered preprocessed reads were generated (two quality scores and two sequenced lengths) and assembled with Mira, Ray and Trinity algorithms. The sixty-four generated assemblies were refined, filtered and evaluated for their proteomic quality. Eight assemblies presented top Detonate scores but one was selected due to its compactness and biological representation, which was generated: (i) from the highest quality dataset (Q20L100), (ii) using Trinity algorithm with all k-mers (AtKa), (iii) removing redundancy by CD-HIT (RR80), and (iv) filtering out poor contigs (F), that was subsequently named Q20L100AtKaRR80F. S. plana proteomic analysis revealed 10,017 peptide groups that corresponded to 2066 proteins with a wide coverage of molecular functions and biological processes, confirming the strength of the database generated.

Sequencing of Transcriptome and Small RNA Revealed the Xenobiotic Metabolism-Related Genes and Potential Regulatory miRNA in Asian Tramp Snail

The Asian tramp snail, Bradybaena similaris (Ferusssac), is an invasive land snail species and has been a rising agricultural pest in south of China. As a pest, it also plays a role in transmission of Angiostrongylus cantonensis. However, present studies on this species are rare and the molecular information is limited. For this purpose, we sequenced the transcriptome and small RNA of B. similaris collected from citrus orchards. In total, 89,747 unigenes with an N50 size of 1287 bp and an average length of 817 bp were generated from ∼8.9 Gb transcriptome and 31 Mb clean reads were generated from ∼36 Mb small RNA library. To demonstrate the usefulness of these two datasets, we analyzed a series of genes associated with xenobiotic metabolism and core RNAi machinery. Analysis of the transcripts resulted in annotation of 126 putative genes encoding cytochrome P450 monooxygenases (CYP, 45), carboxyl/cholinesterases (CCE, 13), glutathione-S-transferases (GST, 24), and ATP-binding cassette transporters (ABC, 44). Analysis of the small RNA detected 42 miRNAs. In addition, four genes involved in small RNA pathways (miRNA, piRNA, and siRNA) were identified, and a total of 430 genes that can be targeted by miRNAs were predicted. Moreover, we found that a few miRNAs could target certain genes involved in xenobiotic metabolism. Therefore, we believe that these two datasets and the characterization of the identified/predicted genes will facilitate the molecular study of this species as well as other land snails with agricultural importance.

De novo transcriptome analysis of Lantana camara L. revealed candidate genes involved in phenylpropanoid biosynthesis pathway

Lantana camara L. is an economically important essential oil producing plant belonging to family Verbenaceae. It is used in medication for treating various diseases like cancer, ulcers, tumor, asthma and fever. The plant is a useful source of essential bioactive compounds such as steroids, flavonoids and phenylpropanoid glycosides etc. Nonetheless, very little is known about the genomic or transcriptomic resources of L. camara, and this might be the reason of hindering molecular studies leading to identification of improved lines. Here we used Illumina sequencing platform and performed the L. camara leaf (LCL) and root (LCR) de novo transcriptome analyses. A total of 70,155,594 and 84,263,224 clean reads were obtained and de novo assembly generated 72,877 and 513,985 unigenes from leaf (LCL) and root (LCR) respectively. Furthermore, the pathway analysis revealed the presence of 229 and 943 genes involved in the phenylpropanoid biosynthesis in leaf and root tissues respectively. Similarity search was performed against publically available genome databases and best matches were found with Sesamum indicum (67.5%) that were much higher than that of Arabidopsis thaliana (3.9%). To the best of our knowledge, this is the first comprehensive transcriptomic analysis of leaf and root tissues of this non-model plant from family Verbenaceae and may serve as a baseline for further molecular studies.

The transcriptome analysis of the whole-body of the gastropod mollusk Limax flavus and screening of putative antimicrobial peptide and protein genes

The gastropod mollusk Limax flavus, one of the most widespread pests in China, is used to treat infectious diseases in traditional Chinese medicine. However, little genomic information is available for this non-model species. In this study, the whole-body transcriptome of L. flavus was sequenced using next generation sequencing technology. A total of 6.81 Gb clean reads were obtained, which were assembled into 150,766 transcripts with 132,206 annotated unigenes. Functionally classification assigned 30,542 unigenes to 56 Gene Ontology terms, 16,745 unigenes were divided into 26 euKaryotic Ortholog Groups of proteins categories, and 13,854 unigenes were assigned to 230 Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, we identified 17,251 simple sequence repeats and several kinds of antimicrobial peptide and protein (AMPs) genes. The transcriptome data of L. flavus will provide a valuable genomic resource for further studies on this species, and the AMPs identified in L. flavus will support its medical potential.

De novo transcriptome of the whole-body of the gastropod mollusk Philomycus bilineatus, a pest with medical potential in China

Philomycus bilineatus is a highly common gastropod mollusk pest in China and is also utilized to treat infectious diseases. However, no genomic resources are available for this non-model species. In the present study, the transcriptomic analysis of P. bilineatus was completed. After sequencing using the next generation sequencing technology, 9.11 Gb of clean reads were obtained, which led to the assembly and annotation of 145,523 transcripts and 125,690 unigenes. Unigenes were functionally classified using Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG). A total of 27,554 unigenes were assigned into 55 GO terms, 13,989 unigenes were differentiated into 26 KOG categories, and 16,368 unigenes were assigned to 229 KEGG pathways. Furthermore, 16,614 simple sequence repeats (SSRs), 38 olfactory genes, and 40 antimicrobial peptide/protein genes were identified. The transcriptome profile of P. bilineatus will provide a valuable genomic resource for further study, will promote the development of new pest management strategies through interference of chemosensory communication, and will support potential medicinal uses of this species.

Transcriptome analysis of air-breathing land slug, Incilaria fruhstorferi reveals functional insights into growth, immunity, and reproduction

Background

Incilaria (= Meghimatium) fruhstorferi is an air-breathing land slug found in restricted habitats of Japan, Taiwan and selected provinces of South Korea (Jeju, Chuncheon, Busan, and Deokjeokdo). The species is on a decline due to depletion of forest cover, predation by natural enemies, and collection. To facilitate the conservation of the species, it is important to decide on a number of traits related to growth, immunity and reproduction addressing fitness advantage of the species.

Results

The visceral mass transcriptome of I. fruhstorferi was enabled using the Illumina HiSeq 4000 sequencing platform. According to BUSCO (Benchmarking Universal Single-Copy Orthologs) method, the transcriptome was considered complete with 91.8% of ortholog genes present (Single: 70.7%; Duplicated: 21.1%). A total of 96.79% of the raw read sequences were processed as clean reads. TransDecoder identified 197,271 contigs that contained candidate-coding regions. Of a total of 50,230 unigenes, 34,470 (68.62% of the total unigenes) annotated to homologous proteins in the Protostome database (PANM-DB). The GO term and KEGG pathway analysis indicated genes involved in metabolism, phosphatidylinositol signalling system, aminobenzoate degradation, and T-cell receptor signalling pathway. Many genes associated with molluscan innate immunity were categorized under pathogen recognition receptor, TLR signalling pathway, MyD88 dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related. The reproduction-associated unigenes showed homology to protein fem-1, spermatogenesis-associated protein, sperm associated antigen, and testis expressed sequences, among others. In addition, we identified key growth-related genes categorized under somatotrophic axis, muscle growth, chitinases and collagens. A total of 4822 Simple Sequence Repeats (SSRs) were also identified from the unigene sequences of I. fruhstorferi.

Conclusions

This is the first available genomic information for non-model land slug, I. fruhstorferi focusing on genes related to growth, immunity, and reproduction, with additional focus on microsatellites and repeating elements. The transcriptome provides access to greater number of traits of unknown relevance in the species that could be exploited for in-depth analyses of evolutionary plasticity and making informed choices during conservation planning. This would be appropriate for understanding the dynamics of the species on a priority basis considering the ecological, health, and social benefits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5526-3) contains supplementary material, which is available to authorized users.