Transcriptome sequencing of Zhikong scallop (Chlamys farreri) and comparative transcriptomic analysis with Yesso scallop (Patinopecten yessoensis)

De Novo Transcriptomic Characterization Enables Novel Microsatellite Identification and Marker Development in Betta splendens

The wild populations of the commercially valuable ornamental fish species, Betta splendens, and its germplasm resources have long been threatened by habitat degradation and contamination with artificially bred fish. Because of the lack of effective marker resources, population genetics research projects are severely hampered. To generate genetic data for developing polymorphic simple sequence repeat (SSR) markers and identifying functional genes, transcriptomic analysis was performed. Illumina paired-end sequencing yielded 105,505,486 clean reads, which were then de novo assembled into 69,836 unigenes. Of these, 35,751 were annotated in the non-redundant, EuKaryotic Orthologous Group, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes and Gene Ontology databases. A total of 12,751 SSR loci were identified from the transcripts and 7970 primer pairs were designed. One hundred primer pairs were randomly selected for PCR validation and 53 successfully generated target amplification products. Further validation demonstrated that 36% (n = 19) of the 53 amplified loci were polymorphic. These data could not only enrich the genetic information for the identification of functional genes but also effectively facilitate the development of SSR markers. Such knowledge would accelerate further studies on the genetic variation and evolution, comparative genomics, linkage mapping and molecular breeding in B. splendens.

In silico determination of nitrogen metabolism in microbes from extreme conditions using metagenomics

The acid ponds of the Danakil Depression in northern Ethiopia are polyextreme environments that exceed the normal physicochemical limits of pH, salinity, ion content, and temperature. We tested for the occurrence of DNA-based life in this environment using Metagenomic Shotgun DNA sequencing approaches. The obtained sequences were examined by the bioinformatic tools MetaSpades, DIAMOND and MEGAN 6-CE, and we were able to bin more than 90% of the metagenomics contigs of Dallol and Black Water to the Bacteria domain, and to the Proteobacteria phylum. Predictions of gene function based on SEED disclosed the presence of different nutrient cycles in the acid ponds. For this study, we focused on partial or completely sequenced genes involved in nitrogen metabolism. The KEGG nitrogen metabolism pathway mapping results for both acid ponds showed that all the predicted genes are involved directly or indirectly in the assimilation of ammonia and no dissimilation or nitrification process was identified. Furthermore, the deduced nitrogen fixation in the two acid ponds based on SEED classification indicated the presence of different sets of nitrogen fixing (nif) genes for biosynthesis and maturation of nitrogenase. Based on the in silico analysis, the predicted proteins involved in nitrogen fixation, especially the cysteine desulfurase and [4Fe-4S] ferredoxin, from both acid ponds are unique with less than 80% sequence similarity to the next closest protein sequence. Considering the extremity of the environmental conditions of the two acid ponds in the Danakil depression, this metagenomics dataset can add to the study of unique gene functions in nitrogen metabolism that enable thriving biocommunities in hypersaline and highly acidic conditions.

Hyperaccumulation of cadmium by scallop Chlamys farreri revealed by comparative transcriptome analysis

Cadmium (Cd) is a hazardous environmental contaminant, which has a serious effect on the ecosystem, food safety and human health. Scallop could accumulate high concentration of Cd from the environment and has been regarded as a Cd hyper-accumulator. In this work, we investigated the antioxidative defense, detoxification and transport of Cd in the kidneys of scallops by transcriptome analysis. A total of 598 differentially expressed genes including 387 up-regulated and 211 down-regulated ones were obtained during Cd exposure, and 46 up-regulated and 260 down-regulated ones were obtained during depuration. Cadmium exposure could cause oxidative stress in the kidneys, which was particularly shown in the pathways involved in proteasome and oxidative phosphorylation. The mRNA expression of 5 metallothionein (MT) genes were overexpressed under Cd exposure and significantly decreased during Cd depuration, which played a vital role in Cd chelation and detoxification. The expression of divalent metal transporter (DMT) genes were down-regulated insignificantly during accumulation and depuration of Cd, which suggested that the DMT played little roles in Cd transport in scallops. A positive relationship in the expression of the zinc transporter (ZIP6 and ZIP1) genes with Cd exposure and depuration was observed, which confirmed its important role for Cd uptake in the kidneys of scallops. 26S proteasome activities and MT expression were Cd-dependent. This study supplied the important reference on the hyperaccumulation of Cd by scallops and identified some effective bioindicators for the environmental risk assessment.

Histological and Expression Differences Among Different Mantle Regions of the Yesso Scallop (Patinopecten yessoensis) Provide Insights into the Molecular Mechanisms of Biomineralization and Pigmentation

The molecular mechanisms of shell formation and pigmentation are issues of great interest in molluscan studies due to the unique physical and biological properties of shells. The Yesso scallop, Patinopecten yessoensis, is one of the most important maricultural bivalves in Asian countries, and its shell color shows polymorphism. To gain more information about the underlying mechanisms of shell formation and pigmentation, this study presents the first analyses of histological and transcriptional differences between different mantle regions of the Yesso scallop, which are thought to be responsible for the formation of different shell layers. The results showed major microstructural differences between the edge and central mantles, which were closely associated with their functions. Different biomineralization-related GO functions, which might participate in the formation of different shell layers, were significantly enriched in the different mantle regions, indicating the different molecular functions of the two mantle regions in shell formation. The melanogenesis pathway, which controls melanin biosynthesis, was the most significantly enriched pathway in the DEGs between the two mantle regions, indicating its important role in shell pigmentation. Tyr, the key and rate-limiting gene in melanogenesis, was expressed at a remarkably high level in the central mantle, while the upstream regulatory genes included in melanogenesis were mainly upregulated in the edge mantle, suggesting the different molecular functions of the two mantle regions in shell pigmentation.

Transcriptomic Analysis of Marine Gastropod Hemifusus tuba Provides Novel Insights into Conotoxin Genes

The marine gastropod Hemifusus tuba is served as a luxury food in Asian countries and used in traditional Chinese medicine to treat lumbago and deafness. The lack of genomic data on H. tuba is a barrier to aquaculture development and functional characteristics of potential bioactive molecules are poorly understood. In the present study, we used high-throughput sequencing technologies to generate the first transcriptomic database of H. tuba. A total of 41 unique conopeptides were retrieved from 44 unigenes, containing 6-cysteine frameworks belonging to four superfamilies. Duplication of mature regions and alternative splicing were also found in some of the conopeptides, and the de novo assembly identified a total of 76,306 transcripts with an average length of 824.6 nt, of which including 75,620 (99.1%) were annotated. In addition, simple sequence repeats (SSRs) detection identified 14,000 unigenes containing 20,735 SSRs, among which, 23 polymorphic SSRs were screened. Thirteen of these markers could be amplified in Hemifusus ternatanus and seven in Rapana venosa. This study provides reports of conopeptide genes in Buccinidae for the first time as well as genomic resources for further drug development, gene discovery and population resource studies of this species.

De novo assembly, characterization of tissue-specific transcriptomes and identification of immune related genes from the scallop Argopecten purpuratus

The scallop Argopecten purpuratus is one of the most economically important cultured mollusks on the coasts from Chile and Peru but its production has declined, in part, due to the emergence of mass mortality events of unknown origin. Driven by this scenario, increasing progress has been made in recent years in the comprehension of immune response mechanisms in this species. However, it is still not entirely understood how different mucosal interfaces participate and cooperate with the immune competent cells, the hemocytes, in the immune defense. Thus, in this work we aimed to characterize the transcriptome of three tissues with immune relevance from A. purpuratus by next-generation sequencing and de novo transcriptome assembly. For this, 18 cDNA libraries were constructed from digestive gland, gills and hemocytes tissues of scallops from different immune conditions and sequenced by the Illumina HiSeq4000 platform. A total of 967.964.884 raw reads were obtained and 967.432.652 clean reads were generated. The clean reads were de novo assembled into 46.601 high quality contigs and 32.299 (69.31%) contigs were subsequently annotated. In addition, three de novo specific assemblies were performed from clean reads obtained from each tissue cDNA libraries for their comparison. Gene ontology (GO) and KEGG analyses revealed that annotated sequences from digestive gland, gills and hemocytes could be classified into both general and specific subcategory terms and known biological pathways, respectively, according to the tissue nature. Finally, several immune related candidate genes were identified, and the differential expression of tissue-specific genes was established, suggesting they could display specific roles in the host defense. The data presented in this study provide the first insight into the tissue specific transcriptome profiles of A. purpuratus, which should be considered for further research on the interplay between the hemocytes and mucosal immune responses.

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.

HD-Marker: a highly multiplexed and flexible approach for targeted genotyping of more than 10,000 genes in a single-tube assay

Targeted genotyping of transcriptome-scale genetic markers is highly attractive for genetic, ecological, and evolutionary studies, but achieving this goal in a cost-effective manner remains a major challenge, especially for laboratories working on nonmodel organisms. Here, we develop a high-throughput, sequencing-based GoldenGate approach (called HD-Marker), which addresses the array-related issues of original GoldenGate methodology and allows for highly multiplexed and flexible targeted genotyping of more than 12,000 loci in a single-tube assay (in contrast to fewer than 3100 in the original GoldenGate assay). We perform extensive analyses to demonstrate the power and performance of HD-Marker on various multiplex levels (296, 795, 1293, and 12,472 genic SNPs) across two sequencing platforms in two nonmodel species (the scallops Chlamys farreri and Patinopecten yessoensis), with extremely high capture rate (98%-99%) and genotyping accuracy (97%-99%). We also demonstrate the potential of HD-Marker for high-throughput targeted genotyping of alternative marker types (e.g., microsatellites and indels). With its remarkable cost-effectiveness (as low as $0.002 per genotype) and high flexibility in choice of multiplex levels and marker types, HD-Marker provides a highly attractive tool over array-based platforms for fulfilling genome/transcriptome-wide targeted genotyping applications, especially in nonmodel organisms.

Scallop genome reveals molecular adaptations to semi-sessile life and neurotoxins

Bivalve molluscs are descendants of an early-Cambrian lineage superbly adapted to benthic filter feeding. Adaptations in form and behavior are well recognized, but the underlying molecular mechanisms are largely unknown. Here, we investigate the genome, various transcriptomes, and proteomes of the scallop Chlamys farreri, a semi-sessile bivalve with well-developed adductor muscle, sophisticated eyes, and remarkable neurotoxin resistance. The scallop's large striated muscle is energy-dynamic but not fully differentiated from smooth muscle. Its eyes are supported by highly diverse, intronless opsins expanded by retroposition for broadened spectral sensitivity. Rapid byssal secretion is enabled by a specialized foot and multiple proteins including expanded tyrosinases. The scallop uses hepatopancreas to accumulate neurotoxins and kidney to transform to high-toxicity forms through expanded sulfotransferases, probably as deterrence against predation, while it achieves neurotoxin resistance through point mutations in sodium channels. These findings suggest that expansion and mutation of those genes may have profound effects on scallop's phenotype and adaptation.

Piwi1 is essential for gametogenesis in mollusk Chlamys farreri

Piwi (P-element induced wimpy testis) is an important gene involved in stem cell maintenance and gametogenesis in vertebrates. However, in most invertebrates, especially mollusks, the function of Piwi during gametogenesis remains largely unclear. To further understand the function of Piwi during gametogenesis, full-length cDNA of Piwi1 from scallop Chlamys farreri (Cf-Piwi1) was characterized, which consisted of a 2,637 bp open reading frame encoding an 878-amino acid protein. Cf-Piwi1 mRNA was mainly localized in the spermatogonia, spermatocytes, oogonia, oocytes of early development and intra-gonadal somatic cells. Additionally, the knockdown of Cf-Piwi1 by injection of Cf-Piwi1-dsRNA (double-stranded RNA) into scallop adductor led to a loss of germ cells in C. farreri gonads. Apoptosis was observed mainly in spermatocytes and oocytes of early development, as well as in a small number of spermatogonia and oogonia. Our findings indicate that Cf-Piwi1 is essential for gametogenesis in the scallop C. farreri.

Construction of an Ostrea edulis database from genomic and expressed sequence tags (ESTs) obtained from Bonamia ostreae infected haemocytes: Development of an immune-enriched oligo-microarray

The flat oyster, Ostrea edulis, is one of the main farmed oysters, not only in Europe but also in the United States and Canada. Bonamiosis due to the parasite Bonamia ostreae has been associated with high mortality episodes in this species. This parasite is an intracellular protozoan that infects haemocytes, the main cells involved in oyster defence. Due to the economical and ecological importance of flat oyster, genomic data are badly needed for genetic improvement of the species, but they are still very scarce. The objective of this study is to develop a sequence database, OedulisDB, with new genomic and transcriptomic resources, providing new data and convenient tools to improve our knowledge of the oyster's immune mechanisms. Transcriptomic and genomic sequences were obtained using 454 pyrosequencing and compiled into an O. edulis database, OedulisDB, consisting of two sets of 10,318 and 7159 unique sequences that represent the oyster's genome (WG) and de novo haemocyte transcriptome (HT), respectively. The flat oyster transcriptome was obtained from two strains (naïve and tolerant) challenged with B. ostreae, and from their corresponding non-challenged controls. Approximately 78.5% of 5619 HT unique sequences were successfully annotated by Blast search using public databases. A total of 984 sequences were identified as being related to immune response and several key immune genes were identified for the first time in flat oyster. Additionally, transcriptome information was used to design and validate the first oligo-microarray in flat oyster enriched with immune sequences from haemocytes. Our transcriptomic and genomic sequencing and subsequent annotation have largely increased the scarce resources available for this economically important species and have enabled us to develop an OedulisDB database and accompanying tools for gene expression analysis. This study represents the first attempt to characterize in depth the O. edulis haemocyte transcriptome in response to B. ostreae through massively sequencing and has aided to improve our knowledge of the immune mechanisms of flat oyster. The validated oligo-microarray and the establishment of a reference transcriptome will be useful for large-scale gene expression studies in this species.

Mantle Branch-Specific RNA Sequences of Moon Scallop Amusium pleuronectes to Identify Shell Color-Associated Genes

Amusium pleuronectes (Linnaeus) that secretes red- and white-colored valves in two branches of mantle tissues is an excellent model for shell color research. High-throughput transcriptome sequencing and profiling were applied in this project to reveal the detailed molecular mechanism of this phenotype differentiation. In this study, 50,796,780 and 54,361,178 clean reads were generated from the left branch (secreting red valve, RS) and right branch (secreting white valve, WS) using the Illumina Hiseq 2000 platform. De novo assembly generated 149,375 and 176,652 unigenes with an average length of 764 bp and 698 bp in RS and WS, respectively. Kyoto encyclopedia of genes and genomes (KEGG) metabolic pathway analysis indicated that the differentially expressed genes were involved in 228 signaling pathways, and 43 genes were significantly enriched (P<0.01). Nineteen of 20 differentially expressed vitellogenin genes showed significantly high expression in RS, which suggested that they probably played a crucial role in organic pigment assembly and transportation of the shell. Moreover, 687 crystal formation-related (or biomineralization-related) genes were detected in A. pleuronectes, among which 144 genes exhibited significant difference between the two branches. Those genes could be classified into shell matrix framework participants, crystal nucleation and growth-related elements, upstream regulation factors, Ca level regulators, and other classifications. We also identified putative SNP and SSR markers from these samples which provided the markers for genetic diversity analysis, genetic linkage, QTL analysis. These results provide insight into the complexity of shell color differentiation in A. pleuronectes so as valuable resources for further research.

Transcriptomics Analysis of Crassostrea hongkongensis for the Discovery of Reproduction-Related Genes

BACKGROUND

The reproductive mechanisms of mollusk species have been interesting targets in biological research because of the diverse reproductive strategies observed in this phylum. These species have also been studied for the development of fishery technologies in molluscan aquaculture. Although the molecular mechanisms underlying the reproductive process have been well studied in animal models, the relevant information from mollusks remains limited, particularly in species of great commercial interest. Crassostrea hongkongensis is the dominant oyster species that is distributed along the coast of the South China Sea and little genomic information on this species is available. Currently, high-throughput sequencing techniques have been widely used for investigating the basis of physiological processes and facilitating the establishment of adequate genetic selection programs.

RESULTS

The C.hongkongensis transcriptome included a total of 1,595,855 reads, which were generated by 454 sequencing and were assembled into 41,472 contigs using de novo methods. Contigs were clustered into 33,920 isotigs and further grouped into 22,829 isogroups. Approximately 77.6% of the isogroups were successfully annotated by the Nr database. More than 1,910 genes were identified as being related to reproduction. Some key genes involved in germline development, sex determination and differentiation were identified for the first time in C.hongkongensis (nanos, piwi, ATRX, FoxL2, β-catenin, etc.). Gene expression analysis indicated that vasa, nanos, piwi, ATRX, FoxL2, β-catenin and SRD5A1 were highly or specifically expressed in C.hongkongensis gonads. Additionally, 94,056 single nucleotide polymorphisms (SNPs) and 1,699 simple sequence repeats (SSRs) were compiled.

CONCLUSIONS

Our study significantly increased C.hongkongensis genomic information based on transcriptomics analysis. The group of reproduction-related genes identified in the present study constitutes a new tool for research on bivalve reproduction processes. The large group of molecular markers discovered in this study will be useful for population screening and marker assisted selection programs in C.hongkongensis aquaculture.

The genome-wide identification of mitogen-activated protein kinase kinase (MKK) genes in Yesso scallop Patinopecten yessoensis and their expression responses to bacteria challenges

Mitogen-activated protein kinase kinases (MKK) are the essential components of the evolutionarily conserved MAPK signaling cascade, which regulates a variety of cellular activities and innate immune responses. Although MKK genes have been extensively studied in various vertebrate and invertebrate species, they have not been systematically characterized in bivalves. In this study, we identified and characterized five MKK genes (PyMKK1/2, PyMKK4, PyMKK5, PyMKK3/6 and PyMKK7) in the Yesso scallop (Patinopecten yessoensis). Phylogenetic and protein structural analyses were conducted to determine their identities and evolutionary relationships. To gain insights into the possible roles of MKK genes during scallop innate immune responses, quantitative real-time PCR (qRT-PCR) was used to investigate their expression profiles during different developmental stages in samples taken from healthy adult tissues and hemocytes after Micrococcus luteus and Vibrio anguillarum bacterial infections. The Yesso scallop MKKs (PyMKKs) were found to have highly conserved structural features compared to the MKK genes from other invertebrate species. Using qRT-PCR analysis, three distinct expression patterns were detected among the PyMKKs over the course of ten different developmental stages. In adult scallops, the majority of the PyMKKs were highly expressed in mantle, gill, muscle and hemocytes. The differential expression patterns of the five PyMKKs after M. luteus (Gram-positive) and V. anguillarum (Gram-negative) bacterial infections suggested their possible involvement in the innate immune response and provide the foundation and resource for the further study on innate immune response of MAPK signal pathway in mollusk.

Differential gene expression analysis of benzo(a)pyrene toxicity in the clam, Ruditapes philippinarum

Polycyclic aromatic hydrocarbons (PAHs) are known for their carcinogenic, teratogenic and mutagenic properties. Benzo(a)pyrene (BaP) possesses the greatest carcinogenic potential among the various PAHs. In this study, digital gene expression (DGE) was performed to investigate the digestive gland transcriptome profile of the clam Ruditapes philippinarum exposed to BaP. A total of 10,508,312 and 11,414,297 clean reads were generated respectively, from control and BaP exposure DGE libraries. One hundred and forty-five differentially expressed genes were detected after comparing two libraries with 58 up-regulated and 87 down-regulated genes. GO annotation and KEGG pathway analyses were performed on all genes to understand their biological functions and processes. The results showed that numerous enriched differentially expressed genes are related to growth and development, antioxidant metabolism, apoptosis and detoxification metabolism. Quantitative real-time PCR was performed to verify the expressed genes of DGE. Our results provide evidences that RNA-seq is a powerful tool for toxicology and capable of generating novel and valuable information at the transcriptome level for characterizing deleterious effects caused by BaP.

Transcriptomic response to estrogen exposure in the male Zhikong scallop, Chlamys farreri

The transcriptomes of Zhikong scallop exposed to 17β-estradiol were determined using the Roche/454. A total of 51,997 unigenes, representing 45,030 contigs and 6967 singlets were obtained. And 14,028, 19,798 and 14,981 of these unigenes were annotated from the non-redundant nucleic acid database, non-redundant protein database and Swiss protein database, respectively. A total of 10,699 unigenes were further annotated to biological processes (9080), molecular functions (8692) and cellular components (7829) using the GO, and 8945 unigenes were mapped to biological pathways including the metabolism (2862) and genetic information processing (2263). Most importantly, 16,692 unigenes and 18,686 unigenes in testis, and 10,492 unigenes and 13,186 unigenes in digestive gland were up-regulated significantly after exposure to 50 and 500 ng E2/L; while 10,212 unigenes and 9409 unigenes in testis and 10,629 unigenes and 9463 unigenes in digestive gland were down-regulated. These valuable information provides insights into the mechanisms in invertebrate exposure to EDCs.

Comparative transcriptome analysis of two oysters, Crassostrea gigas and Crassostrea hongkongensis provides insights into adaptation to hypo-osmotic conditions

Environmental salinity creates a key barrier to limit the distribution of most aquatic organisms. Adaptation to osmotic fluctuation is believed to be a factor facilitating species diversification. Adaptive evolution often involves beneficial mutations at more than one locus. Bivalves hold great interest, with numerous species living in waters, as osmoconformers, who maintain the osmotic pressure balance mostly by free amino acids. In this study, 107,076,589 reads from two groups of Crassostrea hongkongensis were produced and the assembled into 130,629 contigs. Transcripts putatively involved in stress-response, innate immunity and cell processes were identified according to Gene ontology and KEGG pathway analyses. Comparing with the transcriptome of C. gigas to characterize the diversity of transcripts between species with osmotic divergence, we identified 182,806 high-quality single nucleotide polymorphisms (SNPs) for C. hongkongensis, and 196,779 SNPs for C. gigas. Comparison of 11,602 pairs of putative orthologs allowed for identification of 14 protein-coding genes that experienced strong positive selection (Ka/Ks>1). In addition, 45 genes that may show signs of moderate positive selection (1 ≥ Ka/Ks>0.5) were also identified. Based on Ks ratios and divergence time between the two species published previously, we estimated a neutral transcriptome-wide substitution mutation rate of 1.39 × 10(-9) per site per year. Several genes were differentially expressed across the control and treated groups of each species. This is the first time to sequence the transcriptome of C. hongkongensis and provide the most comprehensive transcriptomic resource available for it. The increasing amount of transcriptome data on Crassostrea provides an excellent resource for phylogenetic analysis. A large number of SNPs identified in this work are expected to provide valuable resources for future marker and genotyping assay development. The analysis of natural selection provides an innovative view on the adaptation within species and sets the basis for future genetic and evolutionary studies.

Gene discovery through transcriptome sequencing for the invasive mussel Limnoperna fortunei

The success of the Asian bivalve Limnoperna fortunei as an invader in South America is related to its high acclimation capability. It can inhabit waters with a wide range of temperatures and salinity and handle long-term periods of air exposure. We describe the transcriptome of L. fortunei aiming to give a first insight into the phenotypic plasticity that allows non-native taxa to become established and widespread. We sequenced 95,219 reads from five main tissues of the mussel L. fortunei using Roche’s 454 and assembled them to form a set of 84,063 unigenes (contigs and singletons) representing partial or complete gene sequences. We annotated 24,816 unigenes using a BLAST sequence similarity search against a NCBI nr database. Unigenes were divided into 20 eggNOG functional categories and 292 KEGG metabolic pathways. From the total unigenes, 1,351 represented putative full-length genes of which 73.2% were functionally annotated. We described the first partial and complete gene sequences in order to start understanding bivalve invasiveness. An expansion of the hsp70 gene family, seen also in other bivalves, is present in L. fortunei and could be involved in its adaptation to extreme environments, e.g. during intertidal periods. The presence of toll-like receptors gives a first insight into an immune system that could be more complex than previously assumed and may be involved in the prevention of disease and extinction when population densities are high. Finally, the apparent lack of special adaptations to extremely low O₂ levels is a target worth pursuing for the development of a molecular control approach.

Gonad transcriptome analysis of pearl oyster Pinctada margaritifera: identification of potential sex differentiation and sex determining genes

BACKGROUND

Black pearl farming is based on culture of the blacklip pearl oyster Pinctada margaritifera (Mollusca, lophotrochozoa), a protandrous hermaphrodite species. At first maturation, all individuals are males. The female sex appears progressively from two years old, which represents a limitation for broodstock conditioning for aquaculture production. In marine mollusks displaying hermaphroditic features, data on sexual determinism and differentiation, including the molecular sex determining cascade, are scarce. To increase genomic resources and identify the molecular mechanisms whereby gene expression may act in the sexual dimorphism of P. margaritifera, we performed gonad transcriptome analysis.

RESULTS

The gonad transcriptome of P. margaritifera was sequenced from several gonadic samples of males and females at different development stages, using a Next-Generation-Sequencing method and RNAseq technology. After Illumina sequencing, assembly and annotation, we obtained 70,147 contigs of which 62.2% shared homologies with existing protein sequences, and 9% showed functional annotation with Gene Ontology terms. Differential expression analysis identified 1,993 differentially expressed contigs between the different categories of gonads. Clustering methods of samples revealed that the sex explained most of the variation in gonad gene expression. K-means clustering of differentially expressed contigs showed 815 and 574 contigs were more expressed in male and female gonads, respectively. The analysis of these contigs revealed the presence of known specific genes coding for proteins involved in sex determinism and/or differentiation, such as dmrt and fem-1 like for males, or foxl2 and vitellogenin for females. The specific gene expression profiles of pmarg-fem1-like, pmarg-dmrt and pmarg-foxl2 in different reproductive stages (undetermined, sexual inversion and regression) suggest that these three genes are potentially involved in the sperm-oocyte switch in P. margaritifera.

CONCLUSIONS

The study provides a new transcriptomic tool to study reproduction in hermaphroditic marine mollusks. It identifies sex differentiation and potential sex determining genes in P. margaritifera, a protandrous hermaphrodite species.

Conservation and divergence in the frog immunome: pyrosequencing and de novo assembly of immune tissue transcriptomes

BACKGROUND

Frogs are a diverse group of vertebrates for which limited genomic resources are available. Natural frog populations face a multitude of threats, including habitat degradation, infectious disease, and environmental change. Characterizing the functional genomics of anuran tissues in general - and the immune system in particular - will enhance our knowledge of genetic and epigenetic responses to environmental threats and inform conservation and recovery efforts.

RESULTS

To increase the number of species with genomic datasets and characterize gene expression in immune-related tissues, we sequenced the transcriptomes of three tissues from two frogs (Espadarana prosoblepon and Lithobates yavapaiensis) on the Roche 454 GS FLX platform. Our sequencing produced 8881 E. prosoblepon and 5428 L. yavapaiensis annotated gene products after de novo assembly and Gene Ontology classification. Transcripts of the innate and acquired immune system were expressed in all three tissues. Inflammatory response and acquired immunity transcripts were significantly more diverged between E. prosoblepon and L. yavapaiensis compared to innate immunity and immune system development transcripts. Immune-related transcripts did not show an overall elevated rate of functional evolution, with the exception of glycosyl proteases, which include lysozymes, central bacterial and fungal-killing enzymes of the innate immune system.

CONCLUSIONS

The three frog transcriptomes provide more than 600 Mbp of new genomic data, and will serve as a valuable framework for future comparative studies of non-model anurans. Additionally, we show that immune gene divergence varies by functional group and that transcriptome studies can be useful in comparing rates of evolutionary change across gene families.

Deep transcriptome sequencing of Pecten maximus hemocytes: a genomic resource for bivalve immunology

Pecten maximus, the king scallop, is a bivalve species with important commercial value for both fisheries and aquaculture, traditionally consumed in several European countries. Major problems in larval rearing, however, still limit hatchery-based seed production. High mortalities during early larval stages, likely related to bacterial pathogens, represent the most relevant bottleneck. To address this issue, understanding host defense mechanisms against microbes is extremely important. In this study next-generation RNA-sequencing was carried on scallop hemocytes. To enrich for immune-related transcripts, cDNA libraries from hemocytes challenged in vivo with inactivated-Vibrio anguillarum and in vitro with pathogen-associated molecular patterns, as well as unchallenged controls, were sequenced yielding 216,444,674 sequence reads. De novo assembly of the scallop hemocyte transcriptome consisted of 73,732 contigs (31% annotated). A total of 934 contigs encoded proteins with a known immune function, grouped into several functional categories. Particular attention was reserved to Toll-like receptors (TLRs), a family of pattern recognition receptors (PRRs) involved in non-self recognition. Through mining the scallop hemocyte transcriptome, at least four TLRs could be identified. The organization of canonical TLR domains demonstrated that single cysteine cluster and multiple cysteine cluster TLRs co-exist in this species. In addition, preliminary data concerning their mRNA level following bacterial challenge suggested that different members of this family could exhibit opposite responses to pathogenic stimuli. Finally, a global analysis of differential expression comparing gene-expression levels in in vitro and in vivo stimulated hemocytes against controls provided evidence on a large set of transcripts involved in the great scallop immune response.

A scallop IGF binding protein gene: molecular characterization and association of variants with growth traits

Background

Scallops represent economically important aquaculture shellfish. The identification of genes and genetic variants related to scallop growth could benefit high-yielding scallop breeding. The insulin-like growth factor (IGF) system is essential for growth and development, with IGF binding proteins (IGFBPs) serving as the major regulators of IGF actions. Although an effect of IGF on growth was detected in bivalve, IGFBP has not been reported, and members of the IGF system have not been characterized in scallop.

Results

We cloned and characterized an IGFBP (PyIGFBP) gene from the aquaculture bivalve species, Yesso scallop (Patinopecten yessoensis, Jay, 1857). Its full-length cDNA sequence was 1,445 bp, with an open reading frame of 378 bp, encoding 125 amino acids, and its genomic sequence was 10,193 bp, consisting of three exons and two introns. The amino acid sequence exhibited the characteristics of IGFBPs, including multiple cysteine residues and relatively conserved motifs in the N-terminal and C-terminal domains. Expression analysis indicated that PyIGFBP was expressed in all the tissues and developmental stages examined, with a significantly higher level in the mantle than in other tissues and a significantly higher level in gastrulae and trochophore larvae than in other stages. Furthermore, three single nucleotide polymorphisms (SNPs) were identified in this gene. SNP c.1054A>G was significantly associated with both shell and soft body traits in two populations, with the highest trait values in GG type scallops and lowest in AG type ones.

Conclusion

We cloned and characterized an IGFBP gene in a bivalve, and this report also represents the first characterizing an IGF system gene in scallops. A SNP associated with scallop growth for both the shell and soft body was identified in this gene. In addition to providing a candidate marker for scallop breeding, our results also suggest the role of PyIGFBP in scallop growth.

Large-scale development of gene-associated single-nucleotide polymorphism markers for molluscan population genomic, comparative genomic, and genome-wide association studies

Mollusca is the second most diverse group of animals in the world. Despite their perceived importance, omics-level studies have seldom been applied to this group of animals largely due to a paucity of genomic resources. Here, we report the first large-scale gene-associated marker development and evaluation for a bivalve mollusc, Chlamys farreri. More than 21,000 putative single-nucleotide polymorphisms (SNPs) were identified from the C. farreri transcriptome. Primers and probes were designed and synthesized for 4500 SNPs, and 1492 polymorphic markers were successfully developed using a high-resolution melting genotyping platform. These markers are particularly suitable for population genomic analysis due to high polymorphism within and across populations, a low frequency of null alleles, and conformation to neutral expectations. Unexpectedly, high cross-species transferability was observed, suggesting that the transferable SNPs may largely represent ancestral genetic variations that have been preserved differentially among subfamilies of Pectinidae. Gene annotations were available for 73% of the markers, and 65% could be anchored to the recently released Pacific oyster genome. Large-scale association analysis revealed key candidate genes responsible for scallop growth regulation, and provided markers for further genetic improvement of C. farreri in breeding programmes.

Sequencing-based gene network analysis provides a core set of gene resource for understanding thermal adaptation in Zhikong scallop Chlamys farreri

Marine organisms are commonly exposed to variable environmental conditions, and many of them are under threat from increased sea temperatures caused by global climate change. Generating transcriptomic resources under different stress conditions are crucial for understanding molecular mechanisms underlying thermal adaptation. In this study, we conducted transcriptome-wide gene expression profiling of the scallop Chlamys farreri challenged by acute and chronic heat stress. Of the 13 953 unique tags, more than 850 were significantly differentially expressed at each time point after acute heat stress, which was more than the number of tags differentially expressed (320-350) under chronic heat stress. To obtain a systemic view of gene expression alterations during thermal stress, a weighted gene coexpression network was constructed. Six modules were identified as acute heat stress-responsive modules. Among them, four modules involved in apoptosis regulation, mRNA binding, mitochondrial envelope formation and oxidation reduction were downregulated. The remaining two modules were upregulated. One was enriched with chaperone and the other with microsatellite sequences, whose coexpression may originate from a transcription factor binding site. These results indicated that C. farreri triggered several cellular processes to acclimate to elevated temperature. No modules responded to chronic heat stress, suggesting that the scallops might have acclimated to elevated temperature within 3 days. This study represents the first sequencing-based gene network analysis in a nonmodel aquatic species and provides valuable gene resources for the study of thermal adaptation, which should assist in the development of heat-tolerant scallop lines for aquaculture.

High-resolution linkage and quantitative trait locus mapping aided by genome survey sequencing: building up an integrative genomic framework for a bivalve mollusc

Genetic linkage maps are indispensable tools in genetic and genomic studies. Recent development of genotyping-by-sequencing (GBS) methods holds great promise for constructing high-resolution linkage maps in organisms lacking extensive genomic resources. In the present study, linkage mapping was conducted for a bivalve mollusc (Chlamys farreri) using a newly developed GBS method-2b-restriction site-associated DNA (2b-RAD). Genome survey sequencing was performed to generate a preliminary reference genome that was utilized to facilitate linkage and quantitative trait locus (QTL) mapping in C. farreri. A high-resolution linkage map was constructed with a marker density (3806) that has, to our knowledge, never been achieved in any other molluscs. The linkage map covered nearly the whole genome (99.5%) with a resolution of 0.41 cM. QTL mapping and association analysis congruously revealed two growth-related QTLs and one potential sex-determination region. An important candidate QTL gene named PROP1, which functions in the regulation of growth hormone production in vertebrates, was identified from the growth-related QTL region detected on the linkage group LG3. We demonstrate that this linkage map can serve as an important platform for improving genome assembly and unifying multiple genomic resources. Our study, therefore, exemplifies how to build up an integrative genomic framework in a non-model organism.

Identification of reference genes for qRT-PCR analysis in Yesso scallop Patinopecten yessoensis

BACKGROUND

Bivalves comprise around 30,000 extant species and have received much attention for their importance in ecosystems, aquaculture and evolutionary studies. Despite the increasing application of real-time quantitative reverse transcription PCR (qRT-PCR) in gene expression studies on bivalve species, little research has been conducted on reference gene selection which is critical for reliable and accurate qRT-PCR analysis. For scallops, systematic evaluation of reference genes that can be used among tissues or embryo/larva stages is lacking, and β-actin (ACT) is most frequently used as qRT-PCR reference gene without validation.

RESULTS

In this study, 12 commonly used candidate reference genes were selected from the transcriptome data of Yesso scallop (Patinopectenyessoensis) for suitable qRT-PCR reference genes identification. The expression of these genes in 36 tissue samples and 15 embryo/larva samples under normal physiological conditions was examined by qRT-PCR, and their expression stabilities were evaluated using three statistic algorithms, geNorm, NormFinder, and comparative ∆Ct method. Similar results were obtained by the three approaches for the most and the least stably expressed genes. Final comprehensive ranking for the 12 genes combing the results from the three programs showed that, for different tissues, DEAD-box RNA helicase (HELI), ubiquitin (UBQ), and 60S ribosomal protein L16 (RPL16) were the optimal reference genes combination, while for different embryo/larva stages, gene set containing Cytochrome B (CB), Cytochrome C (CC), Histone H3.3 (His3.3), and Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were recommended for qRT-PCR normalization. ACT was among the least stable genes for both adult tissues and embryos/larvae.

CONCLUSIONS

This work constitutes the first systematic analysis on reference genes selection for qRT-PCR normalization in scallop under normal conditions. The suitable reference genes we recommended will be useful for the identification of genes related to biological processes in Yesso scallop, and also in the reference gene selection for other scallop or bivalve species.

Development of a rapid and efficient method for non-lethal DNA sampling and genotyping in scallops

Non-lethal DNA sampling has long appealed to researchers studying population and conservation genetics, as it does not necessitate removing individuals permanently from their natural environment or destroying valuable samples. However, such an approach has not yet been well established in bivalves. In this study, we demonstrate that the gill represents a good source of tissue for non-lethal sampling in scallops. Removal of a few gill filaments caused no noticeable behavioral abnormalities or increased mortality rates in Zhikong scallop (Chlamys farreri) during a three-month period of observation. To facilitate rapid gill-based DNA extraction, six methods (MA-MF) were designed and evaluated, each requiring less than one hour of processing time. The optimal method was identified as MF, in terms of maintaining DNA integrity and genotyping accuracy. Further optimization of MF method by orthogonal experimental design suggested that the utilization of gills could be limited to 2 mg of sample, which is sufficient for performing up to 20,000 PCR reactions. We also demonstrate the excellent cross-species utility of MF in two additional scallop species, Yesso scallop (Patinopecten yessoensis) and bay scallop (Argopecten irradians). Taken together, our study provides a rapid and efficient approach for applying non-lethal DNA sampling in bivalve species, which would serve as a valuable tool for maintaining bivalve populations and conservation genetics, as well as in breeding studies.