Disk abalone (Haliotis discus discus) expresses a novel antistasin-like serine protease inhibitor: Molecular cloning and immune response against bacterial infection

Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica)

Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.

Make it double: identification and characterization of a Tandem-Hirudin from the Asian medicinal leech Hirudinaria manillensis

Haematophagous leeches express a broad variety of secretory proteins in their salivary glands, among them are hirudins and hirudin-like factors. Here, we describe the identification, molecular and initial functional characterization of Tandem-Hirudin (TH), a novel salivary gland derived factor identified in the Asian medicinal leech, Hirudinaria manillensis. In contrast to the typical structure of hirudins, TH comprises two globular domains arranged in a tandem-like orientation and lacks the elongated C-terminal tail. Similar structures of thrombin inhibitors have so far been identified only in kissing bugs and ticks. Expression of TH was performed in both cell-based and cell-free bacterial systems. A subsequent functional characterization revealed no evidence for a thrombin-inhibitory potency of TH.

In the footsteps of sea stars: deciphering the catalogue of proteins involved in underwater temporary adhesion

Sea stars adhere strongly but temporarily to underwater substrata via the secretion of a blend of proteins, forming an adhesive footprint that they leave on the surface after detachment. Their tube feet enclose a duo-gland adhesive system comprising two types of adhesive cells, contributing different layers of the footprint and de-adhesive cells. In this study, we characterized the catalogue of sea star footprint proteins (Sfps) in the species Asterias rubens to gain insights in their potential function. We identified 16 Sfps and mapped their expression to type 1 and/or type 2 adhesive cells or to de-adhesive cells by double fluorescent in situ hybridization. Based on their cellular expression pattern and their conserved functional domains, we propose that the identified Sfps serve different functions during attachment, with two Sfps coupling to the surface, six providing cohesive strength and the rest forming a binding matrix. Immunolabelling of footprints with antibodies directed against one protein of each category confirmed these roles. A de-adhesive gland cell-specific astacin-like proteinase presumably weakens the bond between the adhesive material and the tube foot surface during detachment. Overall, we provide a model for temporary adhesion in sea stars, including a comprehensive list of the proteins involved.

The Origin and Evolution of Antistasin-like Proteins in Leeches (Hirudinida, Clitellata)

Bloodfeeding is employed by many parasitic animals and requires specific innovations for efficient feeding. Some of these innovations are molecular features that are related to the inhibition of hemostasis. For example, bloodfeeding insects, bats, and leeches release proteins with anticoagulatory activity through their salivary secretions. The antistasin-like protein family, composed of serine protease inhibitors with one or more antistasin-like domains, is tightly linked to inhibition of hemostasis in leeches. However, this protein family has been recorded also in non-bloodfeeding invertebrates, such as cnidarians, mollusks, polychaetes, and oligochaetes. The present study aims to 1) root the antistasin-like gene tree and delimit the major orthologous groups, 2) identify potential independent origins of salivary proteins secreted by leeches, and 3) identify major changes in domain and/or motif structure within each orthologous group. Five clades containing leech antistasin-like proteins are distinguishable through rigorous phylogenetic analyses based on nine new transcriptomes and a diverse set of comparative data: the trypsin + leukocyte elastase inhibitors clade, the antistasin clade, the therostasin clade, and two additional, unnamed clades. The antistasin-like gene tree supports multiple origins of leech antistasin-like proteins due to the presence of both leech and non-leech sequences in one of the unnamed clades, but a single origin of factor Xa and trypsin + leukocyte elastase inhibitors. This is further supported by three sequence motifs that are exclusive to antistasins, the trypsin + leukocyte elastase inhibitor clade, and the therostasin clade, respectively. We discuss the implications of our findings for the evolution of this diverse family of leech anticoagulants.

Spatiotemporal Expression of Anticoagulation Factor Antistasin in Freshwater Leeches

Antistasin, which was originally discovered in the salivary glands of the Mexican leech Haementeria officinalis, was newly isolated from Helobdella austinensis. To confirm the temporal expression of antistasin during embryogenesis, we carried out semi-quantitative RT-PCR. Hau-antistasin1 was uniquely expressed at stage 4 of the cleavage and was strongly expressed in the late stages of organogenesis, as were other antistasin members. In order to confirm the spatial expression of antistasin, we performed fluorescence in situ hybridization in the late stages of organogenesis. The expression of each antistasin in the proboscis showed a similar pattern and varied in expression in the body. In addition, the spatial expression of antistasin orthologs in different leeches showed the possibility of different function across leech species. Hau-antistasin1 was expressed in the same region as hedgehog, which is a known mediator of signal transduction pathway. Hau-antistasin1 is probably a downstream target of Hedgehog signaling, involved in segment polarity signal pathway.

Isolation and Characterization of Poecistasin, an Anti-Thrombotic Antistasin-Type Serine Protease Inhibitor from Leech Poecilobdella manillensis

Antistasin, first identified as a potent inhibitor of the blood coagulation factor Xa, is a novel family of serine protease inhibitors. In this study, we purified a novel antistasin-type inhibitor from leech Poecilobdella manillensis called poecistasin. Amino acid sequencing of this 48-amino-acid protein revealed that poecistasin was an antistasin-type inhibitor known to consist of only one domain. Poecistasin inhibited factor XIIa, kallikrein, trypsin, and elastase, but had no inhibitory effect on factor Xa and thrombin. Poecistasin showed anticoagulant activities. It prolonged the activated partial thromboplastin time and inhibited FeCl₃-induced carotid artery thrombus formation, implying its potent function in helping Poecilobdella manillensis to take a blood meal from the host by inhibiting coagulation. Poecistasin also suppressed ischemic stroke symptoms in transient middle cerebral artery occlusion mice model. Our results suggest that poecistasin from the leech Poecilobdella manillensis plays a crucial role in blood-sucking and may be an excellent candidate for the development of clinical anti-thrombosis and anti-ischemic stroke medicines.

Genetic and molecular basis of the immune system in the brachiopod Lingula anatina

The extension of comparative immunology to non-model systems, such as mollusks and annelids, has revealed an unexpected diversity in the complement of immune receptors and effectors among evolutionary lineages. However, several lophotrochozoan phyla remain unexplored mainly due to the lack of genomic resources. The increasing accessibility of high-throughput sequencing technologies offers unique opportunities for extending genome-wide studies to non-model systems. As a result, the genome-based study of the immune system in brachiopods allows a better understanding of the alternative survival strategies developed by these immunologically neglected phyla. Here we present a detailed overview of the molecular components of the immune system identified in the genome of the brachiopod Lingula anatina. Our findings reveal conserved intracellular signaling pathways as well as unique strategies for pathogen detection and killing in brachiopods.

The purplish bifurcate mussel Mytilisepta virgata gene expression atlas reveals a remarkable tissue functional specialization

BACKGROUND

Mytilisepta virgata is a marine mussel commonly found along the coasts of Japan. Although this species has been the subject of occasional studies concerning its ecological role, growth and reproduction, it has been so far almost completely neglected from a genetic and molecular point of view. In the present study we present a high quality de novo assembled transcriptome of the Japanese purplish mussel, which represents the first publicly available collection of expressed sequences for this species.

RESULTS

The assembled transcriptome comprises almost 50,000 contigs, with a N50 statistics of ~1 kilobase and a high estimated completeness based on the rate of BUSCOs identified, standing as one of the most exhaustive sequence resources available for mytiloid bivalves to date. Overall this data, accompanied by gene expression profiles from gills, digestive gland, mantle rim, foot and posterior adductor muscle, presents an accurate snapshot of the great functional specialization of these five tissues in adult mussels.

CONCLUSIONS

We highlight that one of the most striking features of the M. virgata transcriptome is the high abundance and diversification of lectin-like transcripts, which pertain to different gene families and appear to be expressed in particular in the digestive gland and in the gills. Therefore, these two tissues might be selected as preferential targets for the isolation of molecules with interesting carbohydrate-binding properties. In addition, by molecular phylogenomics, we provide solid evidence in support of the classification of M. virgata within the Brachidontinae subfamily. This result is in agreement with the previously proposed hypothesis that the morphological features traditionally used to group Mytilisepta spp. and Septifer spp. within the same clade are inappropriate due to homoplasy.

Identification and functional characterization of a novel antistasin/WAP-like serine protease inhibitor from the tropical sea cucumber, Stichopus monotuberculatus

A novel antistasin/WAP-like serine protease inhibitor, named as StmAW-SPI, was identified from sea cucumber (Stichopus monotuberculatus) and functionally characterized in this study. The full-length cDNA of StmAW-SPI is 1917 bp in length with a 72 bp 5'-untranslated region (UTR), a 294 bp 3'-UTR and a 1551 bp open reading frame (ORF) encoding a protein of 516 amino acids with a deduced molecular weight of 54.56 kDa. The StmAW-SPI protein has 5-fold internal repeats (IRs) of antistasin domain and 6-fold IRs of WAP domain. For the gene structure, StmAW-SPI contains 10 exons separated by 9 introns. The StmAW-SPI mRNA expression pattern was determined using quantitative real-time PCR. The highest level of StmAW-SPI was found in the intestine, followed by coelomocytes, gonad, body wall and respiratory tree. The StmAW-SPI expressions were significantly up-regulated after polyriboinosinic polyribocytidylic acid [Poly (I:C)] or lipopolysaccharides (LPS) challenge in in vitro experiments performed in primary coelomocytes. In addition, the serine protease inhibitory activity and bacterial protease inhibitory activity of StmAW-SPI were examined, and the antibacterial activity was also demonstrated in this study. Our study, as a whole, suggested that StmAW-SPI might play a critical role in the innate immune defense of sea cucumber against microbial infections, by not only inactivating the serine protease but also inhibiting the growth of pathogens.

Identification of genes associated with shell color in the black-lipped pearl oyster, Pinctada margaritifera

Background

Color polymorphism in the nacre of pteriomorphian bivalves is of great interest for the pearl culture industry. The nacreous layer of the Polynesian black-lipped pearl oyster Pinctada margaritifera exhibits a large array of color variation among individuals including reflections of blue, green, yellow and pink in all possible gradients. Although the heritability of nacre color variation patterns has been demonstrated by experimental crossing, little is known about the genes involved in these patterns. In this study, we identify a set of genes differentially expressed among extreme color phenotypes of P. margaritifera using a suppressive and subtractive hybridization (SSH) method comparing black phenotypes with full and half albino individuals.

Results

Out of the 358 and 346 expressed sequence tags (ESTs) obtained by conducting two SSH libraries respectively, the expression patterns of 37 genes were tested with a real-time quantitative PCR (RT-qPCR) approach by pooling five individuals of each phenotype. The expression of 11 genes was subsequently estimated for each individual in order to detect inter-individual variation. Our results suggest that the color of the nacre is partially under the influence of genes involved in the biomineralization of the calcitic layer. A few genes involved in the formation of the aragonite tablets of the nacre layer and in the biosynthesis chain of melanin also showed differential expression patterns. Finally, high variability in gene expression levels were observed within the black phenotypes.

Conclusions

Our results revealed that three main genetic processes were involved in color polymorphisms: the biomineralization of the nacreous and calcitic layers and the synthesis of pigments such as melanin, suggesting that color polymorphism takes place at different levels in the shell structure. The high variability of gene expression found within black phenotypes suggests that the present work should serve as a basis for future studies exploring more thoroughly the expression patterns of candidate genes within black phenotypes with different dominant iridescent colors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1776-x) contains supplementary material, which is available to authorized users.

Three novel clade B serine protease inhibitors from disk abalone, Haliotis discus discus: Molecular perspectives and responses to immune challenges and tissue injury

Serine protease inhibitors (SERPINs) control cellular protease activity in order to maintain cellular homeostasis. The immune and inflammatory responses of invertebrate clade B SERPINs have not been widely reported. In the present study, three proteins with high similarity to clade B SERPINs, referred to as AbSERPIN-1, AbSERPIN-2 and AbSERPIN-3, were identified from disk abalone (Haliotis discus discus). While AbSERPIN-1 (399 aa) was of a typical size for this protein class, AbSERPIN-2 (506 aa) and AbSERPIN-3 (532 aa) were relatively larger. Bioinformatic analysis revealed the characteristic SERPIN domain in each AbSERPIN. In addition, the N-terminal region of both AbSERPIN-2 and AbSERPIN-3 contained a predicted low complexity region (LCR) and a signal peptide, suggesting that these proteins are secretory proteins and are, thus, novel peptides. Tertiary structural models of the AbSERPINs highlighted their structural and functional conservation. Ubiquitous expression of AbSERPIN transcripts was evaluated by quantitative real time PCR (qPCR) analysis in seven tissue types. AbSERPIN-1, AbSERPIN-2, and AbSERPIN-3 transcript levels were highest in mantle, hemocytes, and muscles, respectively. Temporal expression analysis revealed that AbSERPINs were significantly (P < 0.05) elevated in hemocytes during the early/middle stages following the injection of a bacterial pathogen (Vibrio parahaemolyticus or Listeria monocytogenes) or an immuno-stimulant (lipopolysaccharide). Moreover, mantle tissue injury led to significant changes in the temporal expression of AbSERPIN mRNA. Specifically, transcription of AbSERPIN-1 and AbSERPIN-3 was considerably up-regulated, while expression of AbSERPIN-2 was suppressed. These results suggest a potential role of AbSERPINs in response to pathogen invasion and tissue injury in disk abalone.

Differential expression of serine protease inhibitors 1 and 2 in Crassostrea corteziensis and C. virginica infected with Perkinsus marinus

Proliferation of Perkinsus marinus (Dermo) in vitro is inhibited by the action of 2 serine protease inhibitors belonging to the I-84 family. We compared the levels of expression of serine protease inhibitors 1 and 2 (SPI-1 and SPI-2) in 2 oyster species (Crassostrea virginica and C. corteziensis) inoculated with the parasite P. marinus. C. virginica is well known to be susceptible to this parasite, whereas C. corteziensis is apparently more tolerant. Oysters were inoculated with trophozoites (1 × 106 trophozoites oyster-1) of P. marinus while control oysters were injected with saline solution. Oysters were maintained in a closed water system for 2 wk. The oysters were then sacrificed and parasite burden, histological damage, and gene expression were evaluated. The results showed that the challenged oysters presented a significant increase in parasite burden, which generated histological damage in digestive gland and gills. Quantitative PCR detected significant differences in SPI-1 and SPI-2 expression levels in the 2 oyster species, with C. corteziensis showing higher expression levels than C. virginica as a response to P. marinus inoculation. Our results provide valuable information for the understanding of the defense response in C. corteziensis and a possible explanation for its tolerance to the parasite.

Kazal-type proteinase inhibitor from disk abalone (Haliotis discus discus): molecular characterization and transcriptional response upon immune stimulation

Proteinases and proteinase inhibitors are involved in several biological and physiological processes in all multicellular organisms. Proteinase inhibitors play a key role in regulating the activity of the respective proteinases. Among serine proteinase inhibitors, kazal-type proteinase inhibitors (KPIs) are widely found in mammals, avians, and a variety of invertebrates. In this study, we describe the identification of a kazal-type serine proteinase inhibitor (Ab-KPI) from the disk abalone, Haliotis discus discus, which is presumably involved in innate immunity. The full-length cDNA of Ab-KPI includes 600 bp nucleotides with an open reading frame (ORF) encoding a polypeptide of 143 amino acids. The deduced amino acid sequence of Ab-KPI contains a putative 17-amino acid signal peptide and two tandem kazal domains with high similarity to other kazal-type SPIs. Each kazal domain consists of reactive site (P1) residue containing a leucine (L), and a threonine (T) located in the second amino acid position after the second conserved cysteine of each domain. Temporal expression of Ab-KPI was assessed by real time quantitative PCR in hemocytes and mantle tissue following bacterial and viral hemorrhagic septicemia virus (VHSV) challenge, and tissue injury. At 6 h post-bacterial and -VHSV challenge, Ab-KPI expression in hemocytes was increased 14-fold and 4-fold, respectively, compared to control samples. The highest up-regulations upon tissue injury were shown at 9 h and 12 h in hemocytes and mantle, respectively. The transcriptional modulation of Ab-KPI following bacterial and viral challenges and tissue injury indicates that it might be involved in immune defense as well as wound healing process in abalone.

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.

Gene profiling and characterization of arginine kinase-1 (MrAK-1) from freshwater giant prawn (Macrobrachium rosenbergii)

Arginine kinase-1 (MrAK-1) was sequenced from the freshwater prawn Macrobrachium rosenbergii using Illumina Solexa Genome Analyzer Technique. MrAK-1 consisted of 1068 bp nucleotide encoded 355 polypeptide with an estimated molecular mass of 40 kDa. MrAK-1 sequence contains a potential ATP:guanido phosphotransferases active domain site. The deduced amino acid sequence of MrAK-1 was compared with other 7 homologous arginine kinase (AK) and showed the highest identity (96%) with AK-1 from cherry shrimp Neocaridina denticulate. The qRT-PCR analysis revealed a broad expression of MrAK-1 with the highest expression in the muscle and the lowest in the eyestalk. The expression of MrAK-1 after challenge with the infectious hypodermal and hematopoietic necrosis virus (IHHNV) was tested in muscle. In addition, MrAK-1 was expressed in Escherichia coli by prokaryotic expression plasmid pMAL-c2x. The optimum temperature (30 °C) and pH (8.5) was determined for the enzyme activity assay. MrAK-1 showed significant (P < 0.05) activity towards 10-50 mM ATP concentration. The enzyme activity was inhibited by α-ketoglutarate, glucose and ATP at the concentration of 10, 50 and 100 mM respectively. Conclusively, the findings of this study indicated that MrAK-1 might play an important role in the coupling of energy production and utilization and the immune response in shrimps.