De novo characterisation of the greenlip abalone transcriptome (Haliotis laevigata) with a focus on the heat shock protein 70 (HSP70) family

Transcriptomic responses to thermal stress in hybrid abalone (Haliotis discus hannai ♀ × H. fulgens ♂)

China is the world’s largest abalone producing country. Currently, summer mortality caused by high temperature, is one of the biggest challenges for abalone aquaculture industry. The hybrid abalone (Haliotis discus hannai ♀ × H. fulgens ♂) was conferred on the “new variety”. It has heterosis for thermal tolerance and has been cultured at large-scale in southern China. In this study, a transcriptome analysis was performed to identify the related genes in this hybrid abalone under thermal stress and recovery stage. Compared to control group (18°C), a total of 75, 2173, 1050, 1349, 2548, 494, and 305 differentially expressed genes (DEGs) were identified at 21°C, 24°C, 27°C, 30°C, 32°C, 29°C, and 26°C, respectively. In this study, 24°C is the critical temperature at which the abalone is subjected to thermal stress. With the temperature rising, the number of stress-responsive genes increased. During the temperature recovering to the optimum, the number of stress-responsive genes decreased gradually. Thus, this hybrid abalone has a rapid response and strong adaptability to the temperature. Under the thermal stress, the abalone triggered a complicated regulatory network including degrading the misfolded proteins, activating immune systems, negative regulation of DNA replication, and activating energy production processes. The more quickly feedback regulation, more abundant energy supply and more powerful immune system might be the underlying mechanisms to fight against thermal stress in this hybrid abalone. These findings could provide clues for exploring the thermal-response mechanisms in abalone. The key genes and pathways would facilitate biomarker identification and thermal-tolerant abalone breeding studies.

Phylogenetic annotation of Caenorhabditis elegans heat shock protein 70 genes

Annotation of the 70 kDa heat shock proteins (Hsp70s) has been chaotic especially in invertebrates. In this study, we validated an emerging nomenclature of Hsp70s, which can be potentially applied to all metazoan Hsp70s, by conducting a genome-wide annotation of Caenorhabditis elegans Hsp70s. Using the phylogenetic annotation, the seven canonical C. elegans Hsp70s were successfully classified into four known lineages, cytosolic A, cytosolic B, endoplasmic reticulum, and mitochondria. Motifs specific to each lineage were all conserved in the C. elegans Hsp70s. From these results, we propose new aliases of C. elegans Hsp70s that should help future annotation of this important molecular chaperone.

Effect of low-frequency noise on the survival rate and immunity of infected Vibrio parahaemolyticus sea slug (Onchidium reevesii)

Anthropogenic noise in the marine environment has become a global environmental pollutant that affects the behavior, physiology and immunity of marine animals. However, the resistance of marine animals to pathogens while under the influence of noise is a topic that has received little attention. To assess the immune defense response of sea slugs against pathogens when exposed to low frequency noise, we performed 120 h exposure experiments on sea slugs after a Vibrio parahaemolyticus application in low frequency noise at 500 Hz and 1000 Hz. We found that after the infection with V. parahaemolyticus, the survival rate of the sea slugs decreased, the apoptosis rate and reactive oxygen species (ROS) production of hemocytes increased significantly (P < 0.05), the proliferation of hemocytes accelerated, the activities of enzymes such as superoxide dismutase (SOD), catalase (CAT), alkaline phosphatase (AKP), alanine transaminase (ALT) and lysozyme (LZM) in the hepatopancreas increased significantly, and the expression of TNF signaling pathway-related genes (TNF-α, FADD, Caspase 8, Caspase 3) and Hsp70 genes were generally upregulated. In addition, exposure of sea slug after infected with V. parahaemolyticus to low frequency noise resulted in a significant increase in both antioxidant and immune parameters, which were positively correlated with frequency. The results showed that noise frequency and exposure time had an interactive effect on the above indicators. In summary, low-frequency noise exposure increases the risk of pathogenic infections in sea slugs and exacerbates the negative effects on the antioxidant capacity and immune metabolism of the organism.

Analysis of Genetic Diversity of Banana Weevils (Cosmopolites sordidus) (Coleoptera: Curculionidae) Using Transcriptome-Derived Simple Sequence Repeat Markers

The banana weevil, Cosmopolites sordidus (Germar) (Coleoptera: Curculionidae) is an economically important insect pest of bananas. It causes up to 100% yield losses and substantial lifespan reduction in bananas. Advances in genomics, proteomics, and sequencing technologies have provided powerful pathways to genotyping disastrous pests such as C. sordidus. However, such technologies are often not available to the majority of rural subtropical African banana growers and pest control managers. This study was therefore motivated by the need to create cheap and easily accessible C. sordidus genotyping methods that could be deployed by banana pest control managers to the benefit of C. sordidus control programs in the tropics where such advanced technologies are not readily accessible. We used an in-house C. sordidus transcriptome from the an-ongoing study from which we mined an array of simple sequence repeat (SSR) markers. Of these, six highly polymorphic transcriptome-derived SSR markers were used to successfully genotype within and among banana weevil population genetic diversity of 12 C. sordidus populations collected from four banana-growing agro-ecological zones (AEZs) in Uganda. The developed transcriptome-derived SSR markers can be used by researchers in population genetics for characterization of the C. sordidus and identification of new genes that are linked to traits of particular interest. The significant genetic diversity revealed in C. sordidus provides pertinent information for integrated pest management strategies.

Long-read RNA sequencing of Pacific abalone Haliotis discus hannai reveals innate immune system responses to environmental stress

Haliotis discus hannai is a commercially important mollusk species, and the abalone aquaculture sector has been jeopardized by deteriorating environmental circumstances such as bacterial infection and thermal stress during the hot summers. However, due to a paucity of genetic information, such as transcriptome resources, our understanding of their stress adaptation is restricted. In this research, using single-molecule long-read (SMRT) sequencing technology, a library composed of ten tissues (i.e., haemocytes, gills, muscle, hepatopancreas, digestive tract, mantle, mucous gland, ovary, testis and head) was constructed and sequenced. In all, 41,855 high-quality unique transcripts, among which 24,778 were successfully annotated. Additionally, 13,463 SSRs, 1,169 transcription factors, and 18,124 lncRNAs were identified in H. discus hannai transcriptome. Furthermore, multiple immune-related transcripts were identified according to KEGG annotation, and a portion of these transcripts were mapped into several classical immune-related pathways, including the PI3K-AKT signaling pathway and Toll-like receptor signaling pathway. Additionally, 24 typical sequences related to the immunity pathway were detected by RT-PCR; the results showed that most of the immune-related genes showed significantly high expression at 72 h after bacterial challenges and thermal stress, especially the expression level of genes in gills was significantly higher than that in haemocytes under V. parahaemolyticus stress at 24 h. At the same time. The analysis of alternative splicing identified several innate immunity-related functions genes, including CD109 and caspase 2. These results suggest that the complex immune system, particularly the powerful innate immunity system, was crucial for H. discus hannai response to numerous environmental challenges.

Identification, classification and functional characterization of HSP70s in rainbow trout (Oncorhynchus mykiss) through multi-omics approaches

Heat shock protein 70s (HSP70s) are known to play vital biological processes in rainbow trout. However, information on the numerous roles and classification of many different HSP70s is insufficient. The purpose of this study was to investigate the characteristics of all HSP70s in rainbow trout using multi-dimensional genomic and transcriptomic analyses for inspecting HSP70 homologs, phylogenetic characteristics, DNA motifs, and transcription factor binding sites (TFBSs). Also, the transcriptomic results in conditions of acute thermal stress and Ichthyophthirius multifiliis infection were used to characterize the expression of all HSP70 homologs, and the isoforms of the most sensitive HSP70 were predicted in silico. A total of 23 HSP70s were identified, and they were divided into seven evolutionary groups (groups 1-7). Groups 1 and 2 had relatively longer phylogenetic distances compared to the other groups, which can speculate origin of groups 1 and 2 HSP70s would be different compared to others. With transcriptomic profiling, most HSPs belonging to group 3 showed highly sensitive responses to I. multifiliis infection, not thermal stress, but the group 6 HSP70s had the opposite expression tendencies. Likewise, the composition of the TFBS in each HSP70 was consistent with its group classification. Since TFBSs are widely known to influence transcriptomic expression, they could be one of the major reasons for the different patterns of expression within the HSP70 groups. Moreover, this study demonstrated several isoforms of HSP70a, by far the most sensitive HSP70s, under several stress environments such as hypoxia, thermal, and overcrowding stress. This is an important fundamental study to expand the understanding of HSP70s in rainbow trout as well as for selecting the most sensitive biomarkers for types of stress.

Investigating the biochemical effects of heat stress and sample quenching approach on the metabolic profiling of abalone (Haliotis iris)

INTRODUCTION

Ocean temperatures have been consistently increasing due to climate change, and the frequency of heatwave events on shellfish quality is a growing concern worldwide. Typically, shellfish growing areas are in remote or difficult to access locations which makes in-field sampling and sample preservation of shellfish heat stress difficult. As such, there is a need to investigate in-field sampling approaches that facilitate the study of heat stress in shellfish.

OBJECTIVES

This study aims to apply a gas chromatography-mass spectrometry (GC-MS) based metabolomics approach to examine molecular mechanisms of heat stress responses in shellfish using abalone as a model, and compare the effects of different quenching protocols on abalone metabolic profiles.

METHODS

Twenty adult Haliotis iris abalone were exposed to two temperatures (14 °C and 24 °C) for 24 h. Then, haemolymph and muscle tissues of each animal were sampled and quenched with 4 different protocols (liquid nitrogen, dry ice, cold methanol solution and normal ice) which were analyzed via GC-MS for central carbon metabolites.

RESULTS

The effects of different quenching protocols were only observed in muscle tissues in which the cold methanol solution and normal ice caused some changes in the observed metabolic profiles, compared to dry ice and liquid nitrogen. Abalone muscle tissues were less affected by thermal stress than haemolymph. There were 10 and 46 compounds significantly influenced by thermal stress in muscle and haemolymph, respectively. The changes of these metabolite signatures indicate oxidative damage, disturbance of amino acid and fatty acid metabolism, and a shift from aerobic metabolism to anaerobic pathways.

CONCLUSIONS

The study provided insights into the heat response of abalone, which could be useful for understanding the effects of marine heatwaves and summer mortality events on abalone. Dry ice appeared to be a suitable protocol, and safer in-field alternative to liquid nitrogen, for quenching of abalone tissues.

Allele-specific expression and gene regulation help explain transgressive thermal tolerance in non-native hybrids of the endangered California tiger salamander (Ambystoma californiense)

Hybridization between native and non-native species is an ongoing global conservation threat. Hybrids that exhibit traits and tolerances that surpass parental values are of particular concern, given their potential to outperform native species. Effective management of hybrid populations requires an understanding of both physiological performance and the underlying mechanisms that drive transgressive hybrid traits. Here, we explore several aspects of the hybridization between the endangered California tiger salamander (Ambystoma californiense; CTS) and the introduced barred tiger salamander (Ambystoma mavortium; BTS). We assayed critical thermal maximum (CTMax) to compare the ability of CTS, BTS and F1 hybrids to tolerate acute thermal stress, and found that hybrids exhibit a wide range of CTMax values, with 33% (4/12) able to tolerate temperatures greater than either parent. We then quantified the genomic response, measured at the RNA transcript level, of each salamander, to explore the mechanisms underlying thermal tolerance strategies. We found that CTS and BTS have strikingly different values and tissue-specific patterns of overall gene expression, with hybrids expressing intermediate values. F1 hybrids display abundant and variable degrees of allele-specific expression (ASE), likely arising from extensive compensatory evolution in gene regulatory mechanisms between CTS and BTS. We found evidence that the proportion of genes with allelic imbalance in individual hybrids correlates with their CTMax, suggesting a link between ASE and expanded thermal tolerance that may contribute to the success of hybrid salamanders in California. Future climate change may further complicate management of CTS if hybrid salamanders are better equipped to deal with rising temperatures.

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.

Immunity-related genes and signaling pathways under hypoxic stresses in Haliotis diversicolor: a transcriptome analysis

Due to increased temperatures and aquaculture density, thermal and hypoxia stresses have become serious problems for the aquaculture of abalone Haliotis diversicolor. Stresses lead to immunosuppression, which can cause severe negative impacts on aquaculture farms. To study the mechanism of immunosuppression after hypoxia stress and bacterial challenge, transcriptomes of H. diversicolor hemocytes involved in immunity were profiled. A total of 307,395,572 clean reads were generated and assembled into 99,774 unigenes. KEGG analysis indicated that 225 unigenes with immunologic function were mapped into immune-related pathways. Expression of 41 unigenes measured by quantitative real-time PCR (qRT-PCR) showed consistent results with that of transcriptome analysis. When exposure challenge of Vibrio parahaemolyticus, it is indicated that the PI3K-AKT, MAPK, NF-κB and P53 signal pathways were involved in the hypoxia-induced immunosuppression of H. diversicolor. Furthermore, when the AKT gene (HdAKT) was inhibited by double-stranded RNA (dsRNA), expression levels of HdAKT was lower than the blank and control group in hemocytes at 4 h, 12 h and 24 h (p < 0.05).

Expression of Heat Shock Proteins in Thermally Challenged Pacific Abalone Haliotis discus hannai

Summer mortality, caused by thermal conditions, is the biggest threat to abalone aquaculture production industries. Various measures have been taken to mitigate this issue by adjusting the environment; however, the cellular processes of Pacific abalone (Haliotis discus hannai) have been overlooked due to the paucity of genetic information. The draft genome of H. discus hannai has recently been reported, prompting exploration of the genes responsible for thermal regulation in Pacific abalone. In this study, 413 proteins were systematically annotated as members of the heat shock protein (HSP) super families, and among them 26 HSP genes from four Pacific abalone tissues (hemocytes, gill, mantle, and muscle) were differentially expressed under cold and heat stress conditions. The co-expression network revealed that HSP expression patterns were tissue-specific and similar to those of other shellfish inhabiting intertidal zones. Finally, representative HSPs were selected at random and their expression patterns were identified by RNA sequencing and validated by qRT-PCR to assess expression significance. The HSPs expressed in hemocytes were highly similar in both analyses, suggesting that hemocytes could be more reliable samples for validating thermal condition markers compared to other tissues.

RNA Sequencing Analysis to Capture the Transcriptome Landscape during Tenderization in Sea Cucumber Apostichopus japonicus

Sea cucumber (Apostichopus japonicus) is an economically significant species in China having great commercial value. It is challenging to maintain the textural properties during thermal processing due to the distinctive physiochemical structure of the A. japonicus body wall (AJBW). In this study, the gene expression profiles associated with tenderization in AJBW were determined at 0 h (CON), 1 h (T_1h), and 3 h (T_3h) after treatment at 37 °C using Illumina HiSeq™ 4000 platform. Seven-hundred-and-twenty-one and 806 differentially expressed genes (DEGs) were identified in comparisons of T_1h vs. CON and T_3h vs. CON, respectively. Among these DEGs, we found that two endogenous proteases-72 kDa type IV collagenase and matrix metalloproteinase 16 precursor-were significantly upregulated that could directly affect the tenderness of AJBW. In addition, 92 genes controlled four types of physiological and biochemical processes such as oxidative stress response (3), immune system process (55), apoptosis (4), and reorganization of the cytoskeleton and extracellular matrix (30). Further, the RT-qPCR results confirmed the accuracy of RNA-sequencing analysis. Our results showed the dynamic changes in global gene expression during tenderization and provided a series of candidate genes that contributed to tenderization in AJBW. This can help further studies on the genetics/molecular mechanisms associated with tenderization.

Transcriptomic Analysis Reveals Insights on Male Infertility in Octopus maya Under Chronic Thermal Stress

Octopus maya endemic to the Yucatan Peninsula, Mexico, is an ectotherm organism particularly temperature-sensitive. Studies in O. maya females show that temperatures above 27°C reduce the number of eggs per spawn, fertilization rate and the viability of embryos. High temperatures also reduce the male reproductive performance and success. However, the molecular mechanisms are still unknown. The transcriptomic profiles of testes from thermally stressed (30°C) and not stressed (24°C) adult male octopuses were compared, before and after mating to understand the molecular bases involved in the low reproductive performance at high temperature. The testis paired-end cDNA libraries were sequenced using the Illumina MiSeq platform. Then, the transcriptome was assembled de novo using Trinity software. A total of 53,214,611 high-quality paired reads were used to reconstruct 85,249 transcripts and 77,661 unigenes with an N50 of 889 bp length. Later, 13,154 transcripts were annotated implementing Blastx searches in the UniProt database. Differential expression analysis revealed 1,881 transcripts with significant difference among treatments. Functional annotation and pathway mapping of differential expressed transcripts revealed significant enrichment for biological processes involved in spermatogenesis, gamete generation, germ cell development, spermatid development and differentiation, response to stress, inflammatory response and apoptosis. Remarkably, the transcripts encoding genes such as ZMYND15, KLHL10, TDRD1, TSSK2 and DNAJB13, which are linked to male infertility in other species, were differentially expressed among the treatments. The expression levels of these key genes, involved in sperm motility and spermatogenesis were validated by quantitative real-time PCR. The results suggest that the reduction in male fertility at high temperature can be related to alterations in spermatozoa development and motility.

Different Transcriptomic Responses to Thermal Stress in Heat-Tolerant and Heat-Sensitive Pacific Abalones Indicated by Cardiac Performance

The Pacific abalone Haliotis discus hannai is one of the most economically important mollusks in China. Even though it has been farmed in southern China for almost 20 years, summer mortality remains the most challengeable problem for Pacific abalone aquaculture recently. Here, we determined the different heat tolerance ability for five selective lines of H. discus hannai by measuring the cardiac performance and Arrhenius breakpoint temperature (ABT). The Red line (RL) and Yangxia line (YL) were determined as the most heat-sensitive and most heat-tolerant line, respectively. Heart rates for RL were significantly lower than those of the YL at the same temperature (p < 0.05). The differentially expressed genes (DEGs), which were enriched in several pathways including cardiac muscle contraction, glutathione metabolism and oxidative phosphorylation, were identified between RL and YL at control temperature (20°C) and heat stress temperature (28.5°C, the ABT of the RL) by RNA-seq method. In the RL, 3370 DEGs were identified between the control and the heat-stress temperature, while only 1351 DEGs were identified in YL between these two temperature tests. Most of these DEGs were enriched in the pathways such as protein processing in endoplasmic reticulum, nucleotide binding and oligomerization domain (NOD) like receptor signaling, and ubiquitin mediated proteolysis. Notably, the most heat-tolerant line YL used an effective heat-protection strategy based on moderate transcriptional changes and regulation on the expression of key genes.

Integrative transcriptome analysis and discovery of genes involving in immune response of hypoxia/thermal challenges in the small abalone Haliotis diversicolor

In recent years, the abalone aquaculture industry has been threatened by the deteriorating environmental conditions, such as hypoxia and thermal stress in the hot summers. It is necessary to investigate the molecular mechanism in response to these environmental challenges, and subsequently understand the immune defense system. In this study, the transcriptome profiles by RNA-seq of hemocytes from the small abalone Haliotis diversicolor after exposure to hypoxia, thermal stress, and hypoxia plus thermal stress were established. A total of 103,703,074 clean reads were obtained and 99,774 unigenes were assembled. Of the 99,774 unigenes, 47,154 and 20,455 had homologous sequences in the Nr and Swiss-Prot protein databases, while 16,944 and 10,840 unigenes could be classified by COG or KEGG databases, respectively. RNAseq analysis revealed that the differentially expressed genes (DEGs) after challenges of hypoxia, thermal stress, or hypoxia plus thermal stress were 24,189, 29,165 and 23,665, among which more than 3000 genes involved in at least 230 pathways, including several classical immune-related pathways. The genes and pathways that were involved in immune response to hypoxia/thermal challenges were identified by transcriptome analysis and further validated by quantitative real-time PCR and RNAi technology. The findings in this study can provide information on H. diversicolor innate immunity to improve the abalone aquaculture industry, and the analysis of the potential immune-related genes in innate immunity signaling pathways and the obtained transcriptome data can provide an invaluable genetic resource for the study of the genome and functional genes.

Gonadal Transcriptome Analysis of Pacific Abalone Haliotis discus discus: Identification of Genes Involved in Germ Cell Development

Little is known about the molecular mechanisms governing gonadal developmental processes in abalones. Here, we conducted transcriptome analysis of Pacific abalone Haliotis discus discus for gene discovery in the brain, ovary, testis, and unfertilized eggs. Among the annotated unigenes, 48.6% of unigenes were identified by Venn diagram analysis as having universal or tissue-specific expression. Twenty-three genes with gonad-biased gene ontology (GO) terms were first obtained. Secondly, 36 genes were found by screening known gene names related to germ cell development. Finally, 17 genes were obtained by querying the annotated unigene database for zygotically expressed gonadal genes (ovary and testis) and maternally expressed gonadal genes (ovary, testis, and unfertilized eggs) using keywords related to reproduction. To further verify tissue distribution pattern and subcellular localization of these genes, RT-PCR and in situ hybridization were performed using a unigene encoding a germ cell marker, vasa, as control. The results showed that vasa was expressed mainly in the early developmental stages of germ cells in both sexes. One of the candidate genes, vitelline envelope zona pellucida domain protein 12 (ZP12), was expressed in the primordial germ cells of immature gonad and early developmental stages of germ cells of the adult female. The results obtained from the present study suggest that vasa and ZP12 are involved in germ cell development of Pacific abalone and that ZP12 is an especially useful germ cell-specific marker in immature adults. The current gonadal transcriptome profile is an extensive resource for future reproductive molecular biology studies of this species.

In-depth proteomic analyses of Haliotis laevigata (greenlip abalone) nacre and prismatic organic shell matrix

Background

The shells of various Haliotis species have served as models of invertebrate biomineralization and physical shell properties for more than 20 years. A focus of this research has been the nacreous inner layer of the shell with its conspicuous arrangement of aragonite platelets, resembling in cross-section a brick-and-mortar wall. In comparison, the outer, less stable, calcitic prismatic layer has received much less attention. One of the first molluscan shell proteins to be characterized at the molecular level was Lustrin A, a component of the nacreous organic matrix of Haliotis rufescens. This was soon followed by the C-type lectin perlucin and the growth factor-binding perlustrin, both isolated from H. laevigata nacre, and the crystal growth-modulating AP7 and AP24, isolated from H. rufescens nacre. Mass spectrometry-based proteomics was subsequently applied to to Haliotis biomineralization research with the analysis of the H. asinina shell matrix and yielded 14 different shell-associated proteins. That study was the most comprehensive for a Haliotis species to date.

Methods

The shell proteomes of nacre and prismatic layer of the marine gastropod Haliotis laevigata were analyzed combining mass spectrometry-based proteomics and next generation sequencing.

Results

We identified 297 proteins from the nacreous shell layer and 350 proteins from the prismatic shell layer from the green lip abalone H. laevigata. Considering the overlap between the two sets we identified a total of 448 proteins. Fifty-one nacre proteins and 43 prismatic layer proteins were defined as major proteins based on their abundance at more than 0.2% of the total. The remaining proteins occurred at low abundance and may not play any significant role in shell fabrication. The overlap of major proteins between the two shell layers was 17, amounting to a total of 77 major proteins.

Conclusions

The H. laevigata shell proteome shares moderate sequence similarity at the protein level with other gastropod, bivalve and more distantly related invertebrate biomineralising proteomes. Features conserved in H. laevigata and other molluscan shell proteomes include short repetitive sequences of low complexity predicted to lack intrinsic three-dimensional structure, and domains such as tyrosinase, chitin-binding, and carbonic anhydrase. This catalogue of H. laevigata shell proteins represents the most comprehensive for a haliotid and should support future efforts to elucidate the molecular mechanisms of shell assembly.

Electronic supplementary material

The online version of this article (10.1186/s12953-018-0139-3) contains supplementary material, which is available to authorized users.

De novo assembly, gene annotation, and marker development using Illumina paired-end transcriptome sequencing in the Crassadoma gigantea

Crassadoma gigantea is an important commercial marine bivalve species in Baja California and Mexico. In this study, we have applied RNA-Seq technology to profile the transcriptome of the C. gigantea for the first time. A total of 80,832,518 raw reads were produced from a Illumina HiSeq4000 platform, and 77,306,198 (95.64%) clean reads were generated after trimming the adaptor sequences. The transcriptome assembled into 158,855 transcripts with an N50 size of 1995 bp and an average size of 1008 bp. A number of DNA repair related genes, such as MSH3, EGF, TGF, IGF, FGF, encoding different groups of growth factors were found in the transcriptome data of C. gigantean. In addition, immune related genes Toll-like receptor (TLR) including TLR1, TLR2, TLR3, TLR4, TLR5, TLR6, TLR7, TLR8, and TLR9 was also observed in C. gigantean. A set of 12 polymorphic microsatellite loci were firstly developed and characterized in C. gigantea. The results show that the number of alleles and expected heterozygosity ranged from 3 to 9 and from 0.254 to 0.820, respectively. The average polymorphic information content was 0.790. These microsatellite loci will facilitate future studies of population structure and conservation genetics in this species.

Epipodial Tentacle Gene Expression and Predetermined Resilience to Summer Mortality in the Commercially Important Greenlip Abalone, Haliotis laevigata

"Summer mortality" is a phenomenon that occurs during warm water temperature spikes that results in the mass mortality of many ecologically and economically important mollusks such as abalone. This study aimed to determine whether the baseline gene expression of abalone before a laboratory-induced summer mortality event was associated with resilience to summer mortality. Tentacle transcriptomes of 35 greenlip abalone (Haliotis laevigata) were sequenced prior to the animals being exposed to an increase in water temperature-simulating conditions which have previously resulted in summer mortality. Abalone derived from three source locations with different environmental conditions were categorized as susceptible or resistant to summer mortality depending on whether they died or survived after the water temperature was increased. We detected two genes showing significantly higher expression in resilient abalone relative to susceptible abalone prior to the laboratory-induced summer mortality event. One of these genes was annotated through the NCBI non-redundant protein database using BLASTX to an anemone (Exaiptasia pallida) Transposon Ty3-G Gag Pol polyprotein. Distinct gene expression signatures were also found between resilient and susceptible abalone depending on the population origin, which may suggest divergence in local adaptation mechanisms for resilience. Many of these genes have been suggested to be involved in antioxidant and immune-related functions. The identification of these genes and their functional roles have enhanced our understanding of processes that may contribute to summer mortality in abalone. Our study supports the hypothesis that prestress gene expression signatures are indicative of the likelihood of summer mortality.

Alternative Splicing Profile and Sex-Preferential Gene Expression in the Female and Male Pacific Abalone Haliotis discus hannai

In order to characterize the female or male transcriptome of the Pacific abalone and further increase genomic resources, we sequenced the mRNA of full-length complementary DNA (cDNA) libraries derived from pooled tissues of female and male Haliotis discus hannai by employing the Iso-Seq protocol of the PacBio RSII platform. We successfully assembled whole full-length cDNA sequences and constructed a transcriptome database that included isoform information. After clustering, a total of 15,110 and 12,145 genes that coded for proteins were identified in female and male abalones, respectively. A total of 13,057 putative orthologs were retained from each transcriptome in abalones. Overall Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyzed in each database showed a similar composition between sexes. In addition, a total of 519 and 391 isoforms were genome-widely identified with at least two isoforms from female and male transcriptome databases. We found that the number of isoforms and their alternatively spliced patterns are variable and sex-dependent. This information represents the first significant contribution to sex-preferential genomic resources of the Pacific abalone. The availability of whole female and male transcriptome database and their isoform information will be useful to improve our understanding of molecular responses and also for the analysis of population dynamics in the Pacific abalone.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

De novo assembly and annotation of the European abalone Haliotis tuberculata transcriptome

The European abalone Haliotis tuberculata is a delicacy and consequently a commercially valuable gastropod species. Aquaculture production and wild populations are subjected to multiple climate-associated stressors and anthropogenic pressures, including rising sea-surface temperatures, ocean acidification and an emerging pathogenic Vibrio infection. Transcript expression data provides a valuable resource for understanding abalone responses to variation in the biotic and abiotic environment. To generate an extensive transcriptome, we performed next-generation sequencing of RNA on larvae exposed to temperature and pH variation and on haemolymph of adults from two wild populations after experimental infection with Vibrio harveyi. We obtained more than 1.5 billion raw paired-end reads, which were assembled into 328,519 contigs. Filtration and clustering produced a transcriptome of 41,099 transcripts, of which 10,626 (25.85%) were annotated with Blast hits, and 7380 of these were annotated with Gene Ontology (GO) terms in Blast2Go. A differential expression analysis comparing all samples from the two life stages identified 5690 and 10,759 transcripts with significantly higher expression in larvae and adult haemolymph respectively. This is the greatest sequencing effort yet in the Haliotis genus, and provides the first high-throughput transcriptomic resource for H. tuberculata.

Transcriptome Analysis Revealed Changes of Multiple Genes Involved in Haliotis discus hannai Innate Immunity during Vibrio parahemolyticus Infection

Abalone (Haliotis discus hannai) is one of the most valuable marine aquatic species in Korea, Japan and China. Tremendous exposure to bacterial infection is common in aquaculture environment, especially by Vibrio sp. infections. It’s therefore necessary and urgent to understand the mechanism of H. discus hannai host defense against Vibrio parahemolyticus infection. However studies on its immune system are hindered by the lack of genomic resources. In the present study, we sequenced the transcriptome of control and bacterial challenged H. discus hannai tissues. Totally, 138 MB of reference transcriptome were obtained from de novo assembly of 34 GB clean bases from ten different libraries and annotated with the biological terms (GO and KEGG). A total of 10,575 transcripts exhibiting the differentially expression at least one pair of comparison and the functional annotations highlight genes related to immune response, cell adhesion, immune regulators, redox molecules and mitochondrial coding genes. Mostly, these groups of genes were dominated in hemocytes compared to other tissues. This work is a prerequisite for the identification of those physiological traits controlling H. discus hannai ability to survive against Vibrio infection.

De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile

Paphia textile is an important, aquaculture bivalve clam species distributed mainly in China, Philippines, and Malaysia. Recent studies of P. textile have focused mainly on artificial breeding and nutrition analysis, and the transcriptome and genome of P. textile have rarely been reported. In this work, the transcriptome of P. textile foot tissue was sequenced on an Illumina HiSeq™ 2000 platform. A total of 20,219,795 reads were generated, resulting in 4.08 Gb of raw data. The raw reads were cleaned and assembled into 54,852 unigenes with an N50 length of 829 bp. Of these unigenes, 38.92% were successfully annotated based on their matches to sequences in seven public databases. Among the annotated unigenes, 14,571 were assigned Gene Ontology terms, 5448 were classified to Clusters of Orthologous Groups categories, and 6738 were mapped to 228 pathways in the Kyoto Encyclopedia of Genes and Genomes database. For functional marker development, 5605 candidate simple sequence repeats were identified in the transcriptome and 80 primer pairs were selected randomly and amplified in a wild population of P. textile. A total of 36 loci that exhibited obvious repeat length polymorphisms were detected. The transcriptomic data and microsatellite markers will provide valuable resources for future functional gene analyses, genetic map construction, and quantitative trait loci mapping in P. textile.