Transcriptome analysis in Concholepas concholepas (Gastropoda, Muricidae): mining and characterization of new genomic and molecular markers

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae

BACKGROUND

The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes.

RESULTS

Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs.

CONCLUSION

The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

The first use of LC-MS/MS proteomic approach in the brown mussel Perna perna after bacterial challenge: Searching for key proteins on immune response

The brown mussel Perna perna is a valuable fishing resource, primarily in tropical and subtropical coastal regions. Because of their filter-feeding habits, mussels are directly exposed to bacteria in the water column. Escherichia coli (EC) and Salmonella enterica (SE) inhabit human guts and reach the marine environment through anthropogenic sources, such as sewage. Vibrio parahaemolyticus (VP) is indigenous to coastal ecosystems but can be harmful to shellfish. In this study, we aimed to assess the protein profile of the hepatopancreas of P. perna mussel challenged by introduced - E. coli and S. enterica - and indigenous marine bacteria - V. parahaemolyticus. Bacterial-challenge groups were compared with non-injected (NC) and injected control (IC) - that consisted in mussels not challenged and mussels injected with sterile PBS-NaCl, respectively. Through LC-MS/MS proteomic analysis, 3805 proteins were found in the hepatopancreas of P. perna. From the total, 597 were significantly different among conditions. Mussels injected with VP presented 343 proteins downregulated compared with all the other conditions, suggesting that VP suppresses their immune response. Particularly, 31 altered proteins - upregulated or downregulated - for one or more challenge groups (EC, SE, and VP) compared with controls (NC and IC) are discussed in detail in the paper. For the three tested bacteria, significantly different proteins were found to perform critical roles in immune response at all levels, namely: recognition and signal transduction; transcription; RNA processing; translation and protein processing; secretion; and humoral effectors. This is the first shotgun proteomic study in P. perna mussel, therefore providing an overview of the protein profile of the mussel hepatopancreas, focused on the immune response against bacteria. Hence, it is possible to understand the immune-bacteria relationship at molecular levels better. This knowledge can support the development of strategies and tools to be applied to coastal marine resource management and contribute to the sustainability of coastal systems.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Liver transcriptome analysis of the Sparus macrocephlus in response to Vibrio parahaemolyticus infection

The black seabream (Sparus macrocephlus) is an economically pivotal aquaculture species cultured in China and Southeast Asian countries. To understand the molecular immune mechanisms underlying the response to Vibrio parahaemolyticus, a comparative gene transcription analysis were performed with utilized fresh livers of V. parahaemolyticus-immunized Sparus macrocephlus with a control group through RNA-Seq technology. A total of 256663 contigs were obtained after excluded the low-quality sequences and assembly. The average length of contigs collected from this research is 1066.93 bp. Furthermore, blast analysis indicates 30747 contigs were annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was employed to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. A total of 14470 genes were discovered in 303 KEGG pathways. RSEM and EdgeR were introduced to estimate 3841 genes significantly different expressed (False Discovery Rate<0.001) which includes 4072 up-regulated genes and 3771 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes were conducted to reveal the major immune-related pathways which refer to the toll-like receptor, complement, coagulation cascades, and chemokine signaling pathways. In addition, 92175 potential simple sequence repeats (SSRs) and 121912 candidate single nucleotide polymorphisms (SNPs) were detected and identified sequencely in the Sparus macrocephlus liver transcriptome. This research characterized a gene expression pattern for normal and the V. parahaemolyticus -immunized Sparus macrocephlus for the first time and not only sheds new light on the molecular mechanisms underlying the host-V. parahaemolyticus interaction but contribute to facilitate future studies on Sparus macrocephlus gene expression and functional genomics.

Immune-Related Functional Differential Gene Expression in Koi Carp (Cyprinus carpio) after Challenge with Aeromonas sobria

In order to understand the molecular basis underlying the host immune response of koi carp (Cyprinus carpio), Illumina HiSeq^TM 2000 is used to analyze the muscle and spleen transcriptome of koi carp infected with Aeromonas sobria (A. sobria). De novo assembly of paired-end reads yielded 69,480 unigenes, of which the total length, average length, N50, and GC content are 70,120,028 bp, 1037 bp, 1793 bp, and 45.77%, respectively. Annotation is performed by comparison against various databases, yielding 42,229 (non-redundant protein sequence (NR): 60.78%), 59,255 (non-redundant nucleotide (NT): 85.28%), 35,900 (Swiss-Prot: 51.67%), 11,772 (clusters of orthologous groups (COG): 16.94%), 33,057 (Kyoto Encyclopedia of Genes and Genomes (KEGG): 47.58%), 18,764 (Gene Ontology (GO): 27.01%), and 32,085 (Interpro: 46.18%) unigenes. Comparative analysis of the expression profiles between bacterial challenge fish and control fish identifies 7749 differentially expressed genes (DEGs) from the muscle and 7846 DEGs from the spleen. These DEGs are further categorized with KEGG. Enrichment analysis of the DEGs and unigenes reveals major immune-related functions, including up-regulation of genes related with Toll-like receptor signaling, complement and coagulation cascades, and antigen processing and presentation. The results from RNA-Seq data are also validated and confirmed the consistency of the expression levels of seven immune-related genes after 24 h post infection with qPCR. Microsatellites (11,534), including di-to hexa nucleotide repeat motifs, are also identified. Altogether, this work provides valuable insights into the underlying immune mechanisms elicited during bacterial infection in koi carp that may aid in the future development of disease control measures in protection against A. sobria.

Functional insights into the testis transcriptome of the edible sea urchin Loxechinus albus

The edible sea urchin Loxechinus albus (Molina, 1782) is a keystone species in the littoral benthic systems of the Pacific coast of South America. The international demand for high-quality gonads of this echinoderm has led to an extensive exploitation and decline of its natural populations. Consequently, a more thorough understanding of L. albus gonad development and gametogenesis could provide valuable resources for aquaculture applications, management, conservation and studies about the evolution of functional and structural pathways that underlie the reproductive toolkit of marine invertebrates. Using a high-throughput sequencing technology, we explored the male gonad transcriptome of this highly fecund sea urchin. Through a de novo assembly approach we obtained 42,530 transcripts of which 15,544 (36.6%) had significant alignments to known proteins in public databases. From these transcripts, approximately 73% were functionally annotated allowing the identification of several candidate genes that are likely to play a central role in developmental processes, nutrient reservoir activity, sexual reproduction, gamete generation, meiosis, sex differentiation, sperm motility, male courtship behavior and fertilization. Additionally, comparisons with the male gonad transcriptomes of other echinoderms revealed several conserved orthologous genes, suggesting that similar functional and structural pathways underlie the reproductive development in this group and other marine invertebrates.

De Novo Transcriptome Analysis of Differential Functional Gene Expression in Largemouth Bass (Micropterus salmoides) after Challenge with Nocardia seriolae

Largemouth bass (Micropterus salmoides) are common hosts of an epizootic bacterial infection by Nocardia seriolae. We conducted transcriptome profiling of M. salmoides to understand the host immune response to N. seriolae infection, using the Illumina sequencing platform. De novo assembly of paired-end reads yielded 47,881 unigenes, the total length, average length, N50, and GC content of which were 49,734,288, 1038, 1983 bp, and 45.94%, respectively. Annotation was performed by comparison against non-redundant protein sequence (NR), non-redundant nucleotide (NT), Swiss-Prot, Clusters of Orthologous Groups (COG), Kyoto Encyclopaedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Interpro databases, yielding 28,964 (NR: 60.49%), 36,686 (NT: 76.62%), 24,830 (Swissprot: 51.86%), 8913 (COG: 18.61%), 20,329 (KEGG: 42.46%), 835 (GO: 1.74%), and 22,194 (Interpro: 46.35%) unigenes. Additionally, 8913 unigenes were classified into 25 Clusters of Orthologous Groups (KOGs) categories, and 20,329 unigenes were assigned to 244 specific signalling pathways. RNA-Seq by Expectation Maximization (RSEM) and PossionDis were used to determine significantly differentially expressed genes (False Discovery Rate (FDR) < 0.05) and we found that 1384 were upregulated genes and 1542 were downregulated genes, and further confirmed their regulations using reverse transcription quantitative PCR (RT-qPCR). Altogether, these results provide information on immune mechanisms induced during bacterial infection in largemouth bass, which may facilitate the prevention of nocardiosis.

Transcriptome of the Australian Mollusc Dicathais orbita Provides Insights into the Biosynthesis of Indoles and Choline Esters

Dicathais orbita is a mollusc of the Muricidae family and is well known for the production of the expensive dye Tyrian purple and its brominated precursors that have anticancer properties, in addition to choline esters with muscle-relaxing properties. However, the biosynthetic pathways that produce these secondary metabolites in D. orbita are not known. Illumina HiSeq 2000 transcriptome sequencing of hypobranchial glands, prostate glands, albumen glands, capsule glands, and mantle and foot tissues of D. orbita generated over 201 million high quality reads that were de novo assembled into 219,437 contigs. Annotation with reference to the Nr, Swiss-Prot and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases identified candidate-coding regions in 76,152 of these contigs, with transcripts for many enzymes in various metabolic pathways associated with secondary metabolite biosynthesis represented. This study revealed that D. orbita expresses a number of genes associated with indole, sulfur and histidine metabolism pathways that are relevant to Tyrian purple precursor biosynthesis, and many of which were not found in the fully annotated genomes of three other molluscs in the KEGG database. However, there were no matches to known bromoperoxidase enzymes within the D. orbita transcripts. These transcriptome data provide a significant molecular resource for gastropod research in general and Tyrian purple producing Muricidae in particular.

De novo Transcriptome Generation and Annotation for Two Korean Endemic Land Snails, Aegista chejuensis and Aegista quelpartensis, Using Illumina Paired-End Sequencing Technology

Aegista chejuensis and Aegista quelpartensis (Family-Bradybaenidae) are endemic to Korea, and are considered vulnerable due to declines in their population. The limited genetic resources for these species restricts the ability to prioritize conservation efforts. We sequenced the transcriptomes of these species using Illumina paired-end technology. Approximately 257 and 240 million reads were obtained and assembled into 198,531 and 230,497 unigenes for A. chejuensis and A. quelpartensis, respectively. The average and N50 unigene lengths were 735.4 and 1073 bp, respectively, for A. chejuensis, and 705.6 and 1001 bp, respectively, for A. quelpartensis. In total, 68,484 (34.5%) and 77,745 (33.73%) unigenes for A. chejuensis and A. quelpartensis, respectively, were annotated to databases. Gene Ontology terms were assigned to 23,778 (11.98%) and 26,396 (11.45) unigenes, for A. chejuensis and A. quelpartensis, respectively, while 5050 and 5838 unigenes were mapped to 117 and 124 pathways in the Kyoto Encyclopedia of Genes and Genomes database. In addition, we identified and annotated 9542 and 10,395 putative simple sequence repeats (SSRs) in unigenes from A. chejuensis and A. quelpartensis, respectively. We designed a list of PCR primers flanking the putative SSR regions. These microsatellites may be utilized for future phylogenetics and conservation initiatives.

Transcriptome analysis of the Larimichthys crocea liver in response to Cryptocaryon irritans

The large yellow croaker (Larimichthys crocea) is an economically important marine fish cultured in China and East Asian countries and is facing a serious threat from Cryptocaryon irritans, which is a protozoan ectoparasite that infects most reared marine fish species. To understand the molecular immune mechanisms underlying the response to C. irritans, we first performed a comparative gene transcription analysis using livers from C. irritans-immunized L. croceas and from a control group through RNA-Seq technology. After the removal of low-quality sequences and assembly, 51360 contigs were obtained, with an average length of 1066.93 bp. Further, a blast analysis indicates that 30747 contigs can be annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was used to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. Moreover, 14470 genes were found in 303 KEGG pathways. We used RSEM and EdgeR to determine that 3841 genes were significantly differentially expressed (FDR < 0.001), including 2129 up-regulated genes and 1712 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes revealed major immune-related pathways, including the toll-like receptor, complement and coagulation cascades, and chemokine signaling pathways. In addition, 28748 potential simple sequence repeats (SSRs) were detected from 12776 transcripts, and 62992 candidate single nucleotide polymorphisms (SNPs) were identified in the L. croceas liver transcriptome. This study characterized a gene expression pattern for normal and C. irritans-immunized L. croceas for the first time and not only sheds new light on the molecular mechanisms underlying the host-C. irritans interaction but also facilitates future studies on L. croceas gene expression and functional genomics.

Transcriptome profiles of wild and cultured South African abalone, Haliotis midae

This report describes the use of pyrosequencing technologies to generate the first comparative analysis of de novo assembled transcriptome data from cultured and wild specimens of the South African abalone. The transcriptome data and database described here provide a significant genomic resource for abalone research. The data set annotated 11,240 genes, which matched genes with known functions in other species. A large number of transmembrane protein domains (4087) that may indicate a high portion of undiscovered gene receptors were identified. Further, we detected an interesting set of transcription factors (516) that are valuable candidates for participating in regulatory events in developmental (such as cell proliferation and differentiation) and reproductive processes.

Concholepas concholepas Ferritin H-like subunit (CcFer): Molecular characterization and single nucleotide polymorphism associated to innate immune response

Ferritin has been identified as the principal protein of iron storage and iron detoxification, playing a pivotal role for the cellular homeostasis in living organisms. However, recent studies in marine invertebrates have suggested its association with innate immune system. In the present study, one Ferritin subunit was identified from the gastropod Concholepas concholepas (CcFer), which was fully characterized by Rapid Amplification of cDNA Ends technique. Simultaneously, a challenge test was performed to evaluate the immune response against Vibrio anguillarum. The full length of cDNA Ccfer was 1030 bp, containing 513 bp of open reading frame that encodes to 170 amino acid peptide, which was similar to the Ferritin H subunit described in vertebrates. Untranslated Regions (UTRs) were identified with a 5'UTR of 244 bp that contains iron responsive element (IRE), and a 3'UTR of 273 bp. The predicted molecular mass of deduced amino acid of CcFer was 19.66 kDa and isoelectric point of 4.92. Gene transcription analysis revealed that CcFer increases against infections with V. anguillarum, showing a peak expression at 6 h post-infection. Moreover, a single nucleotide polymorphism was detected at -64 downstream 5'UTR sequence (SNP-64). Quantitative real time analysis showed that homozygous mutant allele (TT) was significantly associated with higher expression levels of the challenged group compared to wild (CC) and heterozygous (CT) variants. Our findings suggest that CcFer is associated to innate immune response in C. concholepas and that the presence of SNPs may involve differential transcriptional expression of CcFer.

Natural product research in the Australian marine invertebrate Dicathais orbita

The predatory marine gastropod Dicathais orbita has been the subject of a significant amount of biological and chemical research over the past five decades. Natural products research on D. orbita includes the isolation and identification of brominated indoles and choline esters as precursors of Tyrian purple, as well as the synthesis of structural analogues, bioactivity testing, biodistributional and biosynthetic studies. Here I also report on how well these compounds conform to Lipinski’s rule of five for druglikeness and their predicted receptor binding and enzyme inhibitor activity. The composition of mycosporine-like amino acids, fatty acids and sterols has also been described in the egg masses of D. orbita. The combination of bioactive compounds produced by D. orbita is of interest for further studies in chemical ecology, as well as for future nutraceutical development. Biological insights into the life history of this species, as well as ongoing research on the gene expression, microbial symbionts and biosynthetic capabilities, should facilitate sustainable production of the bioactive compounds. Knowledge of the phylogeny of D. orbita provides an excellent platform for novel research into the evolution of brominated secondary metabolites in marine molluscs. The range of polarities in the brominated indoles produced by D. orbita has also provided an effective model system used to develop a new method for biodistributional studies. The well characterized suite of chemical reactions that generate Tyrian purple, coupled with an in depth knowledge of the ecology, anatomy and genetics of D. orbita provide a good foundation for ongoing natural products research.

Suppressive subtractive hybridisation transcriptomics provides a novel insight into the functional role of the hypobranchial gland in a marine mollusc

The hypobranchial gland present in gastropods is an organ whose function is not clearly understood. Involved in mucus production, within members of the family Muricidae it is also the source of the ancient dye Tyrian purple and its bioactive precursors. To gain further insights into hypobranchial gland biology, suppressive subtractive hybridisation was performed on hypobranchial gland and mantle tissue from the marine snail Dicathais orbita creating a differentially expressed cDNA library. 437 clones were randomly sequenced, analysed and annotated and 110 sequences had their functions putatively identified. Importantly this approach identified a putative gene involved in Tyrian purple biosynthesis, an arylsulphatase gene. Confirmation of the upregulation of arylsulphatase in the hypobranchial gland compared to the mantle was demonstrated using quantitative real-time PCR. Other genes identified as playing an important role in the hypobranchial gland were those involved in mucus protein synthesis, choline ester regulation, protein and energy production. This study confirms that the hypobranchial gland is involved in the production of mucus secretion and also identifies it as a site of chemical interaction and biosynthesis. This study lays the foundation for a better understanding of the enzymatic production of Tyrian purple precursors within the gland.

Complete mitochondrial genome of Concholepas concholepas inferred by 454 pyrosequencing and mtDNA expression in two mollusc populations

Despite the great relevance of mitochondrial genome analysis in evolutionary studies, there is scarce information on how the transcripts associated with the mitogenome are expressed and their role in the genetic structuring of populations. This work reports the complete mitochondrial genome of the marine gastropod Concholepas concholepas, obtained by 454 pryosequencing, and an analysis of mitochondrial transcripts of two populations 1000 km apart along the Chilean coast. The mitochondrion of C. concholepas is 15,495 base pairs (bp) in size and contains the 37 subunits characteristic of metazoans, as well as a non-coding region of 330 bp. In silico analysis of mitochondrial gene variability showed significant differences among populations. In terms of levels of relative abundance of transcripts associated with mitochondrion in the two populations (assessed by qPCR), the genes associated with complexes III and IV of the mitochondrial genome had the highest levels of expression in the northern population while transcripts associated with the ATP synthase complex had the highest levels of expression in the southern population. Moreover, fifteen polymorphic SNPs were identified in silico between the mitogenomes of the two populations. Four of these markers implied different amino acid substitutions (non-synonymous SNPs). This work contributes novel information regarding the mitochondrial genome structure and mRNA expression levels of C. concholepas.

Ubiquitin-conjugating enzyme E2-like gene associated to pathogen response in Concholepas concholepas: SNP identification and transcription expression

Ubiquitin-conjugated E2 enzyme (UBE2) is one of the main components of the proteasome degradation cascade. Previous studies have shown an increase of expression levels in individuals challenged to some pathogen organism such as virus and bacteria. The study was to characterize the immune response of UBE2 gene in the gastropod Concholepas concholepas through expression analysis and single nucleotide polymorphisms (SNP) discovery. Hence, UBE2 was identified from a cDNA library by 454 pyrosequencing, while SNP identification and validation were performed using De novo assembly and high resolution melting analysis. Challenge trials with Vibrio anguillarum was carried out to evaluate the relative transcript abundance of UBE2 gene from two to thirty-three hours post-treatment. The results showed a partial UBE2 sequence of 889 base pair (bp) with a partial coding region of 291 bp. SNP variation (A/C) was observed at the 546th position. Individuals challenged by V. anguillarum showed an overexpression of the UBE2 gene, the expression being significantly higher in homozygous individuals (AA) than (CC) or heterozygous individuals (A/C). This study contributes useful information relating to the UBE2 gene and its association with innate immune response in marine invertebrates.