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

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.

Characterization of TRAF genes and their responses to Vibrio anguillarum challenge in Argopecten scallops

The tumor necrosis factor receptor-related factor (TRAF) family has been reported to be involved in many immune pathways, such as TNFR, TLR, NLR, and RLR in animals. However, little is known about the roles of TRAF genes in the innate immune of Argopecten scallops. In this study, we first identified five TRAF genes, including TRAF₂, TRAF₃, TRAF₄, TRAF₆ and TRAF₇, but not TRAF₁ and TRAF₅, from both the bay scallop A. irradians (Air) and the Peruvian scallop A. purpuratus (Apu). The phylogenetic analysis showed that the TRAF genes in Argopecten scallops (AiTRAF) belong to the branch of molluscan TRAF family, which lacks TRAF₁ and TRAF₅. Since TRAF₆ is a key bridge factor in the tumor necrosis factor superfamily and plays an important role in innate and adaptive immunity, we cloned the ORFs of the TRAF₆ gene in both A. irradians and A. purpuratus, as well as in two reciprocal hybrids (Aip for the hybrid Air × Apu and Api for the hybrid Apu × Air). Differences in conformational and post-translational modification resulted from the variation in amino acid sequences may cause differences in activity among them. Analysis of conserved motifs and protein structural domains revealed that AiTRAF contains typical structural domains similar to those of other mollusks and has the same conserved motifs. Tissue expression of TRAF in Argopecten scallops challenged by Vibrio anguillarum was examined by qRT-PCR. The results showed that AiTRAF were higher in gill and hepatopancreas. When challenged by Vibrio anguillarum, the expression of AiTRAF was significantly increased compared with the control group, indicating that AiTRAF may play an important role in the immunity of scallops. In addition, the expression of TRAF was higher in Api and Aip than in Air when challenged by Vibrio anguillarum, suggesting that TRAF may have contributed to the high resistance of Api and Aip to Vibrio anguillarum. The results of this study may provide new insights into the evolution and function of TRAF genes in bivalves and ultimately benefit scallop breeding.

Transcriptome analysis of the endangered dung beetle Copris tripartitus (Coleoptera: Scarabaeidae) and characterization of genes associated to immunity, growth, and reproduction

BACKGROUND

Dung beetles recycle organic matter through the decomposition of feces and support ecological balance. However, these insects are threatened by the indiscriminate use of agrochemicals and habitat destruction. Copris tripartitus Waterhouse (Coleoptera: Scarabaeidae), a dung beetle, is listed as a class-II Korean endangered species. Although the genetic diversity of C. tripartitus populations has been investigated through analysis of mitochondrial genes, genomic resources for this species remain limited. In this study, we analyzed the transcriptome of C. tripartitus to elucidate functions related to growth, immunity and reproduction for the purpose of informed conservation planning.

RESULTS

The transcriptome of C. tripartitus was generated using next-generation Illumina sequencing and assembled de novo using a Trinity-based platform. In total, 98.59% of the raw sequence reads were processed as clean reads. These reads were assembled into 151,177 contigs, 101,352 transcripts, and 25,106 unigenes. A total of 23,450 unigenes (93.40%) were annotated to at least one database. The largest proportion of unigenes (92.76%) were annotated to the locally curated PANM-DB. A maximum of 5,512 unigenes had homologous sequences in Tribolium castaneum. Gene Ontology (GO) analysis revealed a maximum of 5,174 unigenes in the Molecular function category. Further, in Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, a total of 462 enzymes were associated with established biological pathways. Based on sequence homology to known proteins in PANM-DB, representative immunity, growth, and reproduction-related genes were screened. Potential immunity-related genes were categorized into pattern recognition receptors (PRRs), the Toll-like receptor signaling pathway, the MyD88- dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related transcripts. Among PRRs, we conducted detailed in silico characterization of TLR-2, CTL, and PGRP_SC2-like. Repetitive elements such as long terminal repeats, short interspersed nuclear elements, long interspersed nuclear elements and DNA elements were enriched in the unigene sequences. A total of 1,493 SSRs were identified among all unigenes of C. tripartitus.

CONCLUSIONS

This study provides a comprehensive resource for analysis of the genomic topography of the beetle C. tripartitus. The data presented here clarify the fitness phenotypes of this species in the wild and provide insight to support informed conservation planning.

Toll like receptors and their evolution in the lymnaeid freshwater snail species Radix auricularia and Lymnaea stagnalis, key intermediate hosts for zoonotic trematodes

One of the major evolutionarily conserved pathways in innate immunity of invertebrates is the toll-like receptor (TLR) pathway. However, little is known of the TLR protein family in gastropod molluscs despite their role in the transmission of human diseases, especially the common lymnaeid freshwater snail species Radix auricularia and Lymnaea stagnalis, key intermediate hosts of zoonotic trematodes. Using comparative genomics and gene prediction approaches utilising the freshwater snail Biomphalaria glabrata genome as a reference ten putative TLR proteins were identified in both R. auricularia and L. stagnalis. Phylogenetic analyses revealed that unlike other molluscs the lymnaeid species also possessed class 1 TLRs, previously thought to be unique to B. glabrata. Gene duplication events were also seen across the TLR classes in the lymnaeids with several of the genes appearing to exist as potential tandem elements in R. auricularia. Each predicted TLR was shown to possess the typical the leucine-rich repeat extracellular and TIR intracellular domains and both single cysteine clusters and multiple cysteine clusters TLRs were identified in both lymnaeid species. Principle component analyses of 3D models of the predicted TLRs showed that class 1 and 5 proteins did not cluster based on similarity of structure, suggested to be potential adaptation to a range of pathogens. This study provides the first detailed account of TLRs in lymnaeids and affords a platform for further research into the role of these proteins into susceptibility and compatibility of these snails with trematodes and their role in transmission.

Identification of peptidoglycan recognition proteins in hemocytes and kidney of common periwinkle Littorinalittorea

Peptidoglycan Recognition Proteins (PGRPs) are a diverse group of proteins involved in innate immunity. In particular, PGRPs have been shown to participate in immune pattern recognition in various mollusks. However, they have not been described in Caenogastropoda, a large molluscan group comprising sea, freshwater and land snails. In this study, four short PGRPs with molecular weights ranging from 21 to 34 kDa and their isoforms were identified and structurally characterized in the kidney and hemocytic transcriptomes of a caenogastropod mollusk Littorina littorea. All of them (LlPGRP1-4) are secretory, possess a signal peptide and a characteristic N-terminal N-acetylmuramoyl-l-alanine amidase (Ami) domain with conserved Zn²⁺ binding- and amidase catalytic sites. The shortest proteins, LlPGRP1 and LlPGRP2, have no additional conserved motifs on the N-terminus. In longer and most abundantly expressed LlPGRP3 and LlPGRP4 the Ami-domain is combined with an N-terminal SH3-domain and a cysteine-rich motif, respectively. Expression analysis showed that LlPGRPs of the common periwinkle were uninvolved in the immune response to infection with trematode Himasthla elongata though they might act in antibacterial defense.

Comparative analysis of the miRNA-mRNA regulation networks in turbot (Scophthalmus maximus L.) following Vibrio anguillarum infection

MicroRNAs could not only regulate posttranscriptional silencing of target genes in eukaryotic organisms, but also have positive effect on their target genes as well. These microRNAs have been reported to be involved in mucosal immune responses to pathogen infection in teleost. Therefore, we constructed the immune-related miRNA-mRNA networks in turbot intestine following Vibrio anguillarum infection. In our results, 1550 differentially expressed (DE) genes and 167 DE miRNAs were identified. 113 DE miRNAs targeting 89 DE mRNAs related to immune response were used to construct miRNA-mRNA interaction networks. Functional analysis showed that target genes were associated with synthesis and degradation of ketone bodies, mucin type O-Glycan biosynthesis, homologous recombination, biotin metabolism, and intestinal immune network for IgA production that were equivalent to the function of IgT and IgM in fish intestine. Finally, 10 DE miRNAs and 7 DE mRNAs were selected for validating the accuracy of high-throughput sequencing results by qRT-PCR. The results of this study will provide valuable information for the elucidation of the regulation mechanisms of miRNA-mRNA interactions involved in disease resistance in teleost mucosal immune system.

Dynamics of MiRNA Transcriptome in Turbot (Scophthalmus maximus L.) Intestine Following Vibrio anguillarum Infection

MicroRNAs (miRNAs) are a group of small non-coding RNAs, which could bind to the 3'-untranslated regions of their target mRNAs to regulate gene expression in various biological processes, including immune-regulated signaling pathways. Turbot (Scophthalmus maximus L.), an important commercial fish species in China, has been suffering with Vibrio anguillarum infection resulted in dramatic economic loss. Therefore, we investigated the expression profiles of miRNAs, as well as the immune-related miRNA-mRNA pairs in turbot intestine at 1 h, 4 h, and 12 h following V. anguillarum infection. As a result, 266 predicted novel miRNAs and 283 conserved miRNAs belonging to 92 miRNA families were detected. A total of 44 miRNAs were differentially expressed in the intestine following V. anguillarum infection. Following prediction, the potential target genes of differentially expressed miRNAs were grouped into a wide range of functional categories, including immune defense/evasion, inflammatory responses, RIG-I signaling pathway, and Toll-like receptor signaling pathway. Moreover, we selected 15 differentially expressed immune genes and their related differentially expressed miRNAs to construct an interaction network for V. anguillarum infection in turbot. These results suggested that in teleost, as in higher vertebrates, miRNAs prominently contribute to immune responses, protecting the host against infection. In addition, this is the first report of comprehensive identification of turbot miRNAs being differentially regulated in the intestine related to V. anguillarum infection. Our results provided an opportunity for further understanding of the molecular mechanisms of miRNA regulation in turbot host-pathogen interactions.